24小时热门版块排行榜    

查看: 3588  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

hn19870519

木虫 (正式写手)

[求助] vasp计算磁矩的困惑

最近计算了一个体系,掺杂
弛豫的KPOINTS用4*4*4,包含Gamma点,计算出总磁矩是2。静态计算增大K点,8*8*8,计算出来的总磁矩为5.1。
感觉这个差别太大了也。
弛豫的INCAR
SYSTEM=common
ISTART=0;ICHARG=2
ENCUT=300
ISMEAR=0
SIGMA=0.01
IBRION=2
NSW=200
ALGO=Fast
EDIFF=1E-05;EDIFFG=-0.02
ISPIN=2
MAGMOM=0.0 62*0 5.0
LWAVE=.FALSE.;LCHARG=.FALSE.
LPLANE=.TRUE.
NPAR=4;NSIM=4
LREAL=Auto
静态计算的INCAR
SYSTEM=common
ISTART=0;ICHARG=2
ENCUT=300
ISMEAR=-5
IBRION=-1
ALGO=Fast
EDIFF=1E-05;EDIFFG=-0.02
ISPIN=2
MAGMOM=0.0 62*0 5.0
LWAVE=.FALSE.;LCHARG=.FALSE.
LPLANE=.TRUE.
NPAR=4;NSIM=4
LREAL=Auto
希望高手能解惑。。。
还有就是问一下大家在计算每个原子磁矩的时候是怎么算的?用不用在弛豫的时候指定RWIGS,这样是否可以得到一个比较准确的值?指定RWIGS的时候具体标准是什么?
问题比较多,先谢谢大家了。。。
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 vasp

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


hn19870519(金币+2): 换了算法会对磁矩有影响么?这种情况倒没有遇到过。有什么根据么? 2011-06-04 17:39:40
youzhizhe(金币+1): 谢谢交流。 2011-06-04 19:13:05
hn19870519(金币+6): 谢谢您热心的指导。 2011-06-05 10:30:01
引用回帖:
Originally posted by uuv2010 at 2011-06-04 15:03:13:
为什么要“对静态计算去掉ALGO=Fast” ?

静态计算的计算量小,加了这么用的算法不一样,变快了得出的结果可能不太好吧,加了也节省不了什么时间
嘟啦啦嘟啦啦嘟
10楼2011-06-04 16:00:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-06-02 22:17:52
hn19870519(金币+5): 好的。我先去试试。 2011-06-03 08:25:56
所有的计算ENCUT=300加大点,对静态计算去掉ALGO=Fast,你先试试看吧
嘟啦啦嘟啦啦嘟
2楼2011-06-02 18:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

自己顶一下。希望大家再讨论一下。
3楼2011-06-03 08:31:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

★ ★
gzqdyouxia(金币+2): 谢谢指教 2011-06-03 16:18:53
上次ENCUT去掉就是300ev的截断能太小了。另外一方面,要看你的体系的大小,如果是个比较小的体系,或许对于4x4x4你取的k点就不够,你可以试着加大k点的网格密度试试,如果还有变化那就是你的K点取少
嘟啦啦嘟啦啦嘟
4楼2011-06-03 15:51:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见