24小时热门版块排行榜    

查看: 911  |  回复: 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的回帖

xiaoboy185

金虫 (正式写手)

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

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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见