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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算反铁磁构型能量变化不大是什么原因？
91/1返回列表
查看: 1806  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者秃头娃将赠送您 15 个金币

秃头娃

新虫 (初入文坛)

[求助] vasp计算反铁磁构型能量变化不大是什么原因？ 已有1人参与

请教一下小木虫上的各位老师同学，最近在计算一种材料的反铁磁构型，之前没有计算经验，经过尝试后，发现计算后不同的磁构型计算出来的能量有变化但是不明显，这种是正常现象还是我参数设置上有问题？附上我的INCAR文件。麻烦有相关计算经验的朋友帮忙解答一下啦！谢谢

#Startparameter for this run:
  4   PREC = Accurate
  5   ISTART = 0
  6   ICHARG = 2
  7   LREAL = .FALSE.
  8 #Ionic Relaxation
  9   IBRION = -1
10   ISIF = 2
11   NSW = 0
12   EDIFFG = 1.0E-2
13   POTIM = 0.25
14 #Electronic Relaxation   
15   EDIFF = 1.0E-05
16   ENCUT = 800 eV
17   ALGO = Fast
18 #magnetism and SOC
19   ISPIN = 2
20   MAGMOM=24*0 16*0 2 -2
21   LORBIT=11
22 #Force constants
23  # ADDGRID = .TRUE.
24 #DOS related values
25   ISMEAR = 0        
26   SIGMA = 0.05
27 #Writing items
28   LWAVE = .TRUE.
29   LCHARG = .TRUE.
30   LVTOT = .FALSE.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2021-12-31 11:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

1. 你需要检查OUTCAR中的磁矩结果，是否符合你的预期
2. 不太清楚你是什么元素？初始磁矩对元素是不同的
3. 能量具体相差多少？有无实验值对照？有的体系差别是比较小的
赞一下 回复此楼
2楼2022-01-03 22:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秃头娃

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 043114076 at 2022-01-03 22:52:45
1. 你需要检查OUTCAR中的磁矩结果，是否符合你的预期
2. 不太清楚你是什么元素？初始磁矩对元素是不同的
3. 能量具体相差多少？有无实验值对照？有的体系差别是比较小的

1.初始结构是否符合预期是指顺磁结构状态下计算后OUTCAR里面的磁矩结果与计算时在INCAR文件里面磁矩设置做对比吗？我设置的是1，后续OUTCAR里面的magnetization (x)只有0.5左右，这是否存在问题？
2.我的体系之前是不含磁性的，后期掺杂的氢元素之后就有磁性了，在计算结果OSZICAR里面发现氢元素的磁矩大概就是1了，所以我设置的是1
3.我计算的这个体系没有参数对比，能量差的话原始在169.55eV.我计算之后的能量差距也就0.03eV，我感觉这个差距实在是太小了。
赞一下 回复此楼
3楼2022-01-04 09:54:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 秃头娃 at 2022-01-04 09:54:46
1.初始结构是否符合预期是指顺磁结构状态下计算后OUTCAR里面的磁矩结果与计算时在INCAR文件里面磁矩设置做对比吗？我设置的是1，后续OUTCAR里面的magnetization (x)只有0.5左右，这是否存在问题？
2.我的体系之前 ...

如果材料本身不含磁性，我的经验是加了H以后产生的磁性变化是正常的，能量差距很小也算正常。说明这个体系没有明显的反铁磁性质。仅供参考。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

秃头娃
4楼2022-01-04 23:04:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秃头娃

新虫 (初入文坛)
送红花一朵
引用回帖:
4楼: Originally posted by 043114076 at 2022-01-04 23:04:18
如果材料本身不含磁性，我的经验是加了H以后产生的磁性变化是正常的，能量差距很小也算正常。说明这个体系没有明显的反铁磁性质。仅供参考。...

好的，谢谢
回复此楼
5楼2022-01-05 14:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
3楼: Originally posted by 秃头娃 at 2022-01-04 09:54:46
1.初始结构是否符合预期是指顺磁结构状态下计算后OUTCAR里面的磁矩结果与计算时在INCAR文件里面磁矩设置做对比吗？我设置的是1，后续OUTCAR里面的magnetization (x)只有0.5左右，这是否存在问题？
2.我的体系之前 ...

1、magnetization (x)小于1也正常，因为这个是通过WS半径划分得来的，你看看OSZICAR总的磁矩显示多少？
2、没磁性的体系加上一个H，因为H自身带一个电子，如果这个电子能够局域在H上，整个体系应该是有磁矩的，MAGMOM可以按照H带一个磁矩的情形进行猜测；反之有些体系其实也会出现这个电子不局域在H上的情况，比如有些金属吸附H可能还是会出现没磁矩的能量更低，这点建议经过计算之后确认。
3、磁性引起的能量变化确实不大，这种情况建议使用小一些的展宽看一下，以及尽可能多copy CONTCAR优化几次，vasp的结果对初始结构和电荷猜测有些敏感。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

秃头娃
不一定挂在论坛，计算问题问题欢迎留言。
6楼2022-01-05 14:50:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秃头娃

新虫 (初入文坛)
送红花一朵
引用回帖:
6楼: Originally posted by 卡开发发 at 2022-01-05 14:50:49
1、magnetization (x)小于1也正常，因为这个是通过WS半径划分得来的，你看看OSZICAR总的磁矩显示多少？
2、没磁性的体系加上一个H，因为H自身带一个电子，如果这个电子能够局域在H上，整个体系应该是有磁矩的，MA ...

1.掺杂H元素之后，OSZICAR显示的磁矩是1.但是我后续掺杂的F进去，OSZICAR 显示只有0.9164，我计算的时候还是选择MAGMOM=1.这个应该不影响吧？
2.我后期做了一下自旋电荷密度图，显示的整个物质仅仅只在H周围有一圈，应该就是您所说的该电子局域在H上面吧。
3.我接下来尝试一下多计算几次CONTCAR,谢谢！我也书上看到说一般而言不同磁结构的能量差都比较小，磁性结构还依赖与经验参数U,但是我好像又看到的资料显示结构+U的话可能更针对哪些原子半径更大的物质有效果，我计算的物质主要集中与第二周期，感觉有点不太需要。是不是意义不太大啊？
赞一下 回复此楼
7楼2022-01-07 16:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by 秃头娃 at 2022-01-07 16:23:10
1.掺杂H元素之后，OSZICAR显示的磁矩是1.但是我后续掺杂的F进去，OSZICAR 显示只有0.9164，我计算的时候还是选择MAGMOM=1.这个应该不影响吧？
2.我后期做了一下自旋电荷密度图，显示的整个物质仅仅只在H周围有一圈 ...

1、不影响，不过你可以也检查下ISPIN=1的情况两者谁能量更低。
2、对，按道理是。
3、U主要是针对d或f这些体系的自作用问题导致的低估局域化的修正。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

秃头娃
不一定挂在论坛，计算问题问题欢迎留言。
8楼2022-01-07 23:11:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秃头娃

新虫 (初入文坛)
送红花一朵
引用回帖:
8楼: Originally posted by 卡开发发 at 2022-01-07 23:11:20
1、不影响，不过你可以也检查下ISPIN=1的情况两者谁能量更低。
2、对，按道理是。
3、U主要是针对d或f这些体系的自作用问题导致的低估局域化的修正。...

好的，谢谢指导！

发自小木虫Android客户端
回复此楼
9楼2022-01-09 11:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 秃头娃 的主题更新
91/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 267一志愿南京工业大学0817化工求调剂 +4 SUICHILD 2026-03-12 4/200 2026-03-13 14:14 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 复试调剂 +7 Copy267 2026-03-10 7/350 2026-03-13 10:37 by houyaoxu
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 2026博士申请 +7 起泡酒 2026-03-08 7/350 2026-03-13 09:56 by mumu真真
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 290 材料与化工求调剂 +8 Nebulala 2026-03-08 8/400 2026-03-12 15:19 by wyb0328
[考研] 材料工程调剂 +3 咪咪空空 2026-03-11 3/150 2026-03-12 15:02 by 每天只摆一小会
[考研] 2026考研求调剂-材料类-本科211一志愿985-初试301分 +10 虫友233 2026-03-07 11/550 2026-03-12 14:52 by ruiyingmiao
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子！
[考研] 材料工程专硕调剂自荐信，初试323 +6 李白26 2026-03-07 6/300 2026-03-10 16:16 by 无际的草原
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
[考研] 307求调剂 +3 辛仁豆腐 2026-03-08 5/250 2026-03-09 14:09 by macy2011
[考研] 083000环境科学与工程调剂 +5 加油呀fxy 2026-03-07 6/300 2026-03-08 21:07 by sdgdsdgd
[考研] 288求调剂（0703）一志愿东北大学 +5 好好- 2026-03-07 5/250 2026-03-08 18:22 by zzpnuli111
[考研] 求调剂，一志愿华中科大0702，数一英一，293 +4 小罗露一二 2026-03-07 4/200 2026-03-08 16:36 by 星空星月
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 求调剂 +3 snw石 2026-03-07 3/150 2026-03-07 20:01 by houyaoxu
信息提示
关闭
请填处理意见
关闭