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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算磁矩的困惑
111/2返回列表12下一页
查看: 3664  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 292求调剂 +5 yhk_819 2026-02-28 5/250 2026-03-01 21:38 by 公瑾逍遥
[考研] 0856化工专硕求调剂 +12 董boxing 2026-03-01 12/600 2026-03-01 19:45 by 材子momo
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 化工专硕348，一志愿985求调剂 +5 弗格个 2026-02-28 8/400 2026-03-01 17:25 by sunny81
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[基金申请] 刚录用，没有期刊号，但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
关闭
请填处理意见
关闭