24小时热门版块排行榜    

查看: 910  |  回复: 8
当前主题已经存档。

xiaoboy185

金虫 (正式写手)

[交流] 【求助】磁性计算

请教各位:我优化了一个构型   输入文件如下
   ISTART = 0  
   ICHARG = 2  
   INIWAV = 1  
   NELM   =     80;
   EDIFF  =   1E-05
   ISYM = 0        
   IBRION = 2      
   NSW = 300
   ISMEAR = 0   
   SIGMA= 0.02   
   IALGO= 48
   ISPIN = 2
   MAGMOM = 49*2
   LORBIT=11

得到的ourcar文件里面的总磁矩是4,而OSZICAR的最后写着是8,不知道怎么会差
这么多呢?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
TIGERYZZ(金币+1,VIP+0):谢谢交流 12-21 19:07
你的问题描述不是很详细,但是我想你INCAR设置的参数太说了,特别是你的LORBIT一般是计算DOS才使用。你的IBRION显示你要做relax。还有你说不一致,最好贴出来!
2楼2009-12-18 09:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)

楼上你好   我算的是一个纯碳的体系
我是想看看每个原子的磁矩,所以加上了lorbit这个关键词   
我只是做了优化,要是想算磁性的话,还需要在算个静态的计算吗?
OSZICAR最后一行
25 F= -.29309798E+03 E0= -.29309798E+03  d E =-.294688E-03  mag=     8.0000
OUTCAR磁矩部分
tot        0.716   3.813   0.000   4.529(总磁矩)
3楼2009-12-18 10:10:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice_rain

铁杆木虫 (正式写手)


TIGERYZZ(金币+1,VIP+0):谢谢交流 12-21 19:07
C本身应该很好算。为什么不试试呢。做个静态计算一下就是了,算了,可以发过总结,我们一起讨论一下。
4楼2009-12-18 13:08:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
TIGERYZZ(金币+1,VIP+0):谢谢回帖交流 12-21 19:07
你只是计算一个C吗?osicar中的因该是总磁矩,outcar中给出的应该按照各个原子给出的吧,但是磁矩相差太大,肯定是有问题的!
5楼2009-12-18 13:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chgzhong

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
TIGERYZZ(金币+1,VIP+0):谢谢回帖交流 12-21 19:08
这要看你算的什么体系,多少磁性离子,通常OSZICAR中只有总磁矩,而OUTCAR中可以看到各个原子的磁矩。
6楼2009-12-18 14:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)


zdhlover(金币+1,VIP+0):谢谢,欢迎常来交流讨论 12-18 00:29
先优化再静态计算
7楼2009-12-18 14:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)

多谢大家,能不能帮我看看这个输入文件有没有啥问题呢?
感觉算的太慢了!
8楼2009-12-18 16:04:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
TIGERYZZ(金币+1,VIP+0):谢谢回帖交流 12-21 19:08
一般加了磁性是算的比较慢!!!呵呵!
9楼2009-12-21 17:57:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoboy185 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见