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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】为什么静态计算和原先的磁矩差异如此之大！
121/2返回列表12下一页
查看: 1792  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

jugengfans

金虫 (正式写手)

[交流] 【求助】为什么静态计算和原先的磁矩差异如此之大！

RT，在结构优化后我把NSW等注释掉，进行静态计算，
但是，为什么得出的磁矩根刚才计算时的出入非常大！
结构优化时我的磁矩是4.1，静态计算后变成了6.7,这是怎么回事？
恳请各位帮忙看一下，谢谢了！

我作静态计算的的INCAR文件：

SYSTEM = layer
ISTART = 0
ICHARG = 2

ENCUT = 520
ISPIN = 2


#MAGMOM = 8*0.0 8*1.0  24*0.0
EDIFF    = 1E-4
EDFIIG = -1E-2
ISMEAR   = 2

NPAR = 4

#NELM=200
#NELMDL=8
#WEIMIN=0

AMIX     = 0.2
BMIX     = 0.00001
AMIX_MAG = 0.8
BMIX_MAG = 0.00001
LREAL = .TRUE.
Dynamic:
#IBRION = 2
#NSW = 500
#POTIM = 0.2

#LORBIT   = 11

#LDAU      = .TRUE.
#LDAUTYPE  = 2
#LDAUL     = 2 -1 -1
#LDAUU     = 3.00 0.00 0.00
#LDAUJ     = 0.87 0.00 0.00
#LDAUPRINT = 2

LMAXMIX   = 4          ! Important: mix paw occupancies up to L=4
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 第一性原理和电化学

» 猜你喜欢
1楼 2009-08-07 01:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cihinustc

★
jugengfans(金币+1):谢谢参与
期待高手
回复此楼
2楼2009-08-07 08:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六198

金虫 (正式写手)
★ ★
jugengfans(金币+1):谢谢参与
yjr(金币+1,VIP+0):谢谢！！ 8-7 09:01
结构优化时的给出的体系性质是不可信的，应该看结构优化后进行自洽计算的结果
赞一下(7人) 回复此楼
穷理求真真贯纵横皆入咏致知格物物含妙理总堪寻
3楼2009-08-07 08:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★
jugengfans(金币+1):谢谢参与
yjr(金币+1,VIP+0):谢谢！！ 8-7 09:02
我在ESPRESSO计算graphene ribbon时也遇到类似的情况。
做结构优化时，出来的是个顺磁的
做静态计算出来是反铁磁的：（

关键是文献报道的是饭铁磁的：/
赞一下(7人) 回复此楼
4楼2009-08-07 08:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
★ ★ ★
jugengfans(金币+1):谢谢参与
yjr(金币+2,VIP+0):谢谢您的热心！ 8-7 12:34
引用回帖:
Originally posted by jugengfans at 2009-8-7 01:36:
RT，在结构优化后我把NSW等注释掉，进行静态计算，
但是，为什么得出的磁矩根刚才计算时的出入非常大！
结构优化时我的磁矩是4.1，静态计算后变成了6.7,这是怎么回事？
恳请各位帮忙看一下，谢谢了！

我作静 ...

你在做静态计算的时候没有输入MAGMOM值，所以vasp默认按照铁磁序优化，有可能和你原来的磁结构不相同，而恰好体系在铁磁结构也有个稳定态，就得到了不同的结果。vasp说明里面提到可以在INCAR里设置ICHARG=1让程序从前面的CHGCAR文件读取电荷密度，同时也读取磁矩的信息，但是几乎很难以实用，因为改变了k点或者NBANDS等等参数后，可能造成读取不到合适的电荷密度，所以我一直都在静态计算的时候也设置MAGMOM值，这样和优化时得到的磁矩基本上相同。
赞一下(7人) 回复此楼
5楼2009-08-07 09:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)
引用回帖:
Originally posted by nkleof at 2009-8-7 09:49:

你在做静态计算的时候没有输入MAGMOM值，所以vasp默认按照铁磁序优化，有可能和你原来的磁结构不相同，而恰好体系在铁磁结构也有个稳定态，就得到了不同的结果。vasp说明里面提到可以在INCAR里设置ICHARG=1让程 ...

那么如果这个静态优化的能量比结构优化时给出的能量低0,可以认为这个磁矩是体系最稳态时的磁矩么？
赞一下 回复此楼
6楼2009-08-07 12:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

去日留痕

至尊木虫 (文坛精英)

小小木虫虫

优秀版主
★
jugengfans(金币+1):谢谢参与
做静态计算时没有输入MAGMOM值
赞一下(1人) 回复此楼
世间最珍贵的不是“得不到”和“已失去”，而是现在可把握的幸福....
7楼2009-08-07 12:37:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)
★
freshgirl(金币+1,VIP+0):正解，谢谢参与~ 8-7 20:51
引用回帖:
Originally posted by 去日留痕 at 2009-8-7 12:37:
做静态计算时没有输入MAGMOM值

默认的MAGMON是所有原子均为1的铁磁性排列吧。
赞一下(6人) 回复此楼
8楼2009-08-07 14:48:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六198

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by jugengfans at 2009/8/7 14:48:

默认的MAGMON是所有原子均为1的铁磁性排列吧。

是的
回复此楼
穷理求真真贯纵横皆入咏致知格物物含妙理总堪寻
9楼2009-08-07 17:06:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
fegg7502(金币+2,VIP+0):鼓励交流 8-16 17:32
楼主要是想计算FM的话  不用在设置MAGMOM了吧  默认的就是FM???
赞一下(11人) 回复此楼
10楼2009-08-08 20:58:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
关闭
请填处理意见
关闭