24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1866  |  回复: 8
【悬赏金币】回答本帖问题,作者秃头娃将赠送您 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.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注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的回帖
查看全部 9 个回答

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕283求调剂 +12 试试看呗 2026-04-04 13/650 2026-04-05 15:45 by wwytracy
[考研] 315求调剂 +10 欣喜777 2026-04-04 11/550 2026-04-05 12:55 by Hdyxbekcb
[考研] 081700,311,求调剂 +12 冬十三 2026-04-04 13/650 2026-04-05 08:43 by qlm5820
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 070300一志愿211,312分求调剂院校 +16 小黄鸭宝 2026-03-30 16/800 2026-04-03 19:53 by lijunpoly
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
[考研] 江苏苏北高校诚邀调剂同学 +3 zzll406 2026-03-31 3/150 2026-03-31 16:54 by 及时行乐fan
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
信息提示
请填处理意见