24小时热门版块排行榜    

查看: 1346  |  回复: 11
当前主题已经存档。

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 第一性原理和电化学

» 猜你喜欢

已阅   回复此楼   关注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
结构优化时的给出的体系性质是不可信的,应该看结构优化后进行自洽计算的结果
穷理求真真贯纵横皆入咏致知格物物含妙理总堪寻
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时也遇到类似的情况。
做结构优化时,出来的是个顺磁的
做静态计算出来是反铁磁的:(

关键是文献报道的是饭铁磁的:/
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值,这样和优化时得到的磁矩基本上相同。
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值
世间最珍贵的不是“得不到”和“已失去”,而是现在可把握的幸福....
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的铁磁性排列吧。
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???
10楼2009-08-08 20:58:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见