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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算磁矩的困惑
111/2返回列表12下一页
查看: 3723  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +7 1624447980 2026-04-08 8/400 2026-04-11 00:45 by dc19971224
[考研] 284求调剂 +11 archer.. 2026-04-10 12/600 2026-04-10 18:57 by HPUCZ
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 299求调剂 +6 ZVVZ13 2026-04-08 6/300 2026-04-09 15:48 by wp06
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 专硕310求调剂 +7 捞捞我…. 2026-04-04 8/400 2026-04-08 20:24 by yutian743
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 求助 +3 卡卡东88 2026-04-06 4/200 2026-04-06 15:28 by going home
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
信息提示
关闭
请填处理意见
关闭