应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算磁矩的困惑
111/2返回列表12下一页
查看: 3671  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-02 17:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-06-02 22:17:52
hn19870519(金币+5): 好的。我先去试试。 2011-06-03 08:25:56
所有的计算ENCUT=300加大点，对静态计算去掉ALGO=Fast，你先试试看吧
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
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点取少
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
4楼2011-06-03 15:51:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)
★ ★
gzqdyouxia(金币+2): 谢谢指教 2011-06-03 21:32:24
引用回帖:
Originally posted by 贺仪 at 2011-06-03 15:51:50:
上次ENCUT去掉就是300ev的截断能太小了。另外一方面，要看你的体系的大小，如果是个比较小的体系，或许对于4x4x4你取的k点就不够，你可以试着加大k点的网格密度试试，如果还有变化那就是你的K点取少

可能是K点取小的原因，我试试。
有个新问题就是，ENCUT增大到400，离子步跑到60还没收敛。ENCUT=300时跑到20步左右就收敛了。
赞一下(1人) 回复此楼
5楼2011-06-03 17:29:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★
hn19870519(金币+5): 文献上的ENCUT也是300左右的。我觉得不会是ENCUT太小吧。 2011-06-04 13:07:34
zzy870720z(金币+1): 谢谢提示，这个看测试吧 2011-06-04 14:32:51
这样看的话ENCUT是小了，我一般都是设置成500的
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
6楼2011-06-04 09:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)
引用回帖:
Originally posted by 贺仪 at 2011-06-04 09:14:10:
这样看的话ENCUT是小了，我一般都是设置成500的

SMEAR和SIGMA的设置是否会影响磁矩的计算呢？
赞一下 回复此楼
7楼2011-06-04 13:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
引用回帖:
Originally posted by 贺仪 at 2011-06-02 18:17:17:
所有的计算ENCUT=300加大点，对静态计算去掉ALGO=Fast，你先试试看吧

为什么要“对静态计算去掉ALGO=Fast” ？
赞一下 回复此楼
8楼2011-06-04 15:03:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-06-04 19:12:59
hn19870519(金币+2): 您说的更物理能不能解释一下呢？ 2011-06-05 10:29:30
引用回帖:
Originally posted by hn19870519 at 2011-06-04 13:08:41:
SMEAR和SIGMA的设置是否会影响磁矩的计算呢？

没做过这方面的研究，对于此行我一般用ismear=-1，因为这样更物理，其他的情况没研究过
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
9楼2011-06-04 15:59:25
已阅   回复此楼   关注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” ？

静态计算的计算量小，加了这么用的算法不一样，变快了得出的结果可能不太好吧，加了也节省不了什么时间
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
10楼2011-06-04 16:00:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hn19870519 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
关闭
请填处理意见
关闭