24小时热门版块排行榜    

查看: 1655  |  回复: 4

算来一梦浮生

新虫 (初入文坛)

[求助] vasp磁性计算优化步与自洽步得出的磁矩数值差距很大 已有1人参与

求助各位大神!
我在用vasp计算材料磁性时,优化步得到的磁矩和自洽步得到的磁矩数值相差14uB左右,正常情况下优化步和自洽步得到的磁矩不是应该相差不大吗?这个问题一直未得到解决,还请大家不吝赐教!多谢多谢!以下是我优化自洽时的INCAR和KPOINTS:
优化步INCAR
ISTART =  0
  ICHARG =  2
  PREC   =  Accurate
  ENCUT  =  500
  EDIFF  =  1E-05
  EDIFFG =  -0.01
  NSW    = 500
  IBRION =  2
  ISIF   =  2
  #ISYM   =  2
  ISMEAR =  0
  LREAL=Auto
  NPAR=4
  ISPIN = 2
  SIGMA = 0.05
  LWAVE = F
  IVDW = 11  
自洽步INCAR
ISTART =  0
  ICHARG =  2
  PREC   =  Accurate
  ENCUT  =  500
  EDIFF  =  1E-05
  EDIFFG =  -0.01
  #NSW    = 500
  #IBRION =  2
  #ISIF   =  2
  #ISYM   =  2
  ISMEAR =  0
  LREAL=Auto
  NPAR=4
  ISPIN = 2
  SIGMA = 0.05
  LAECHG = .True.
  IVDW = 11  
KPOINTS均为:
Automatic generation
0
Gamma
5    5     1
0.0  0.0   0.0
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
算来一梦浮生(蛋蛋小童鞋代发): 金币+3, 谢谢帮助 2017-09-15 18:20:33
我的疑问有
1. 你第二部自洽的POSCAR 是第一步的CONTCAR吗
2. 一般自洽电子需要用LCHARGE=TURE ISTART=1,使用之前的电荷文件
如果这俩是对的 总能和磁性差别很小的
路漫漫其修远兮
2楼2017-09-15 17:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

算来一梦浮生

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by hs0218 at 2017-09-15 17:02:14
我的疑问有
1. 你第二部自洽的POSCAR 是第一步的CONTCAR吗
2. 一般自洽电子需要用LCHARGE=TURE ISTART=1,使用之前的电荷文件
如果这俩是对的 总能和磁性差别很小的

非常感谢您!
1.我第二部自洽的POSCAR是用的第一步的CONTCAR
2.我自洽步用的ISTART=0,计算态密度和能带时才用的ISTART=1
您的意思是在自洽步就用ISTART=1吗?
3楼2017-09-15 20:05:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feng27

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 算来一梦浮生 at 2017-09-15 20:05:11
非常感谢您!
1.我第二部自洽的POSCAR是用的第一步的CONTCAR
2.我自洽步用的ISTART=0,计算态密度和能带时才用的ISTART=1
您的意思是在自洽步就用ISTART=1吗?...

我也遇到同样问题
4楼2018-08-29 08:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bujiesan

新虫 (正式写手)

楼主你好,我最近刚接触vasp,你可以把你这次计算的完整程序给我看看吗?万分感谢

发自小木虫Android客户端
好好学物理!
5楼2018-09-04 23:18:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 算来一梦浮生 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 13:54 by GkCVyODlx4mR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见