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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算磁矩的困惑
111/2返回列表12下一页
查看: 2981  |  回复: 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下一页
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 连续两年医学口青年项目初审体会 +4 进击的荣耀 2024-06-09 7/350 2024-06-10 13:19 by 进击的荣耀
[基金申请] 今年香江初审出结果了吗 +5 悲催科研狗 2024-06-06 14/700 2024-06-10 11:29 by Pickfoot
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话，在编者有什么补偿？ +7 河西夜郎 2024-06-09 8/400 2024-06-10 08:31 by zzq770204
[教师之家] 公办双非，学生论文升学就业都不行。一本大学的论文升学就业没问题吧？ +4 河西夜郎 2024-06-09 4/200 2024-06-10 08:11 by Quakerbird
[考博] 希望能25博士入学，可提前一年做科研助理 +3 干饭版小太阳 2024-06-09 7/350 2024-06-10 00:05 by 干饭版小太阳
[找工作] 今年的就业形势真的很差吗 +4 215747018 2024-06-07 5/250 2024-06-09 23:25 by ffffjjjj
[基金申请] C口食品 +4 好好不过分 2024-06-09 16/800 2024-06-09 22:08 by Pickfoot
[基金申请] 看到大家一直在说代表作NC，不明觉厉 +23 翰海2022 2024-06-07 23/1150 2024-06-09 18:55 by hhs666
[基金申请] G口yq发了吗 +9 nnyyxx 2024-06-03 11/550 2024-06-09 14:34 by oylb2
[基金申请] B01化学口的，近几年要什么样的评价才能上会呀？ +4 salmon95 2024-06-06 5/250 2024-06-09 09:45 by chem-cai
[论文投稿] Analytica Chimica Acta 投稿系统中 abstract 字符数限制600什么意思？ 3+3 3921 2024-06-05 5/250 2024-06-09 09:45 by bobvan
[论文投稿] 工作后评职称需要发表论文，想自己写，求帮助 50+3 上官逸夜 2024-06-04 9/450 2024-06-07 22:57 by xy66xy
[基金申请] B口yq答辩 +5 学员qYccD0 2024-06-05 6/300 2024-06-06 22:39 by 2673228047
[基金申请] 5B能上会么？ +12 Tenfu 2024-06-05 14/700 2024-06-06 22:34 by mch20120521
[论文投稿] 电化学投稿 +4 爱安沐梅梅 2024-06-06 14/700 2024-06-06 21:51 by 爱安沐梅梅
[论文投稿] 返修之后，今天with editor状态没变，但是日期更新了，有人知道是怎么回事么？ 50+3 huanpo116 2024-06-04 8/400 2024-06-05 21:44 by lucky_my2010
[论文投稿] 不要投这个杂志Journal of materials science +5 蛋蛋198956 2024-06-04 5/250 2024-06-05 18:18 by nandi2212
[基金申请] E口YQ答辩 +3 firepick 2024-06-05 3/150 2024-06-05 17:12 by kaka198888
[基金申请] E口人才项目发通知了么？ +3 linan2011 2024-06-04 3/150 2024-06-04 22:34 by hncxu
[论文投稿] 硅酸盐学报审稿费晚交有啥影响没 +4 赵杉嫚 2024-06-03 4/200 2024-06-04 21:20 by topedit
信息提示
关闭
请填处理意见
关闭