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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP优化结构时，E0发生突变？
北京石油化工学院2026年研究生招生接收调剂公告
181/1返回列表
查看: 1704  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lijuanshu

铜虫 (正式写手)

[交流] VASP优化结构时，E0发生突变？

我在用VASP优化结构时，前4个离子步的E0都在175 eV左右，但是从5F后能量就发生了突变，变成了73 eV左右，并且显示错误信息。如果是正确的结果能量应该在176eV左右的。
  1 F= -.17431755E+03 E0= -.17432851E+03  d E =-.174318E+03  mag=    33.8873
   2 F= -.17441190E+03 E0= -.17442320E+03  d E =-.943533E-01  mag=    33.8831
   3 F= -.17468180E+03 E0= -.17469457E+03  d E =-.364249E+00  mag=    33.8813
   4 F= -.17520759E+03 E0= -.17522201E+03  d E =-.890035E+00  mag=    33.8579
   5 F= -.70424051E+02 E0= -.70427134E+02  d E =0.103893E+03  mag=    32.8561
   6 F= -.72780584E+02 E0= -.72785187E+02  d E =0.101537E+03  mag=    32.5818
   7 F= -.72810234E+02 E0= -.72815120E+02  d E =0.101507E+03  mag=    32.5715
   8 F= -.72829505E+02 E0= -.72834576E+02  d E =0.101488E+03  mag=    32.5578
   9 F= -.72838108E+02 E0= -.72843141E+02  d E =0.101479E+03  mag=    32.5506
  10 F= -.72842228E+02 E0= -.72847536E+02  d E =0.101475E+03  mag=    32.5461
  11 F= -.72843606E+02 E0= -.72848723E+02  d E =0.101474E+03  mag=    32.5419
  12 F= -.72844565E+02 E0= -.72849802E+02  d E =0.101473E+03  mag=    32.5403
  13 F= -.72844971E+02 E0= -.72850165E+02  d E =0.101473E+03  mag=    32.5386
  14 F= -.72845060E+02 E0= -.72850245E+02  d E =0.101472E+03  mag=    32.5367
  15 F= -.72845127E+02 E0= -.72850308E+02  d E =0.101472E+03  mag=    32.5358
  16 F= -.72845215E+02 E0= -.72850408E+02  d E =0.101472E+03  mag=    32.5353
  17 F= -.72845162E+02 E0= -.72850352E+02  d E =0.101472E+03  mag=    32.5346
  18 F= -.72845144E+02 E0= -.72850316E+02  d E =0.101472E+03  mag=    32.5341
  19 F= -.72845109E+02 E0= -.72850304E+02  d E =0.101472E+03  mag=    32.5336
  20 F= -.72699958E+02 E0= -.72708257E+02  d E =0.101618E+03  mag=    32.2532   
curvature: ****** expect dE=-0.275E+02 dE for cont linesearch -0.695E+01
ZBRENT: fatal internal in brackting
         system-shutdown; contact gK immediately
我的INCAR如下：
SYSTEM = Ni2
ENCUT = 280
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIGMA = 0.2
NSW = 200
IBRION = 2
ISPIN = 2
NELM = 2000
POTIM = 0.2
EDIFF = 1E-4
EDIFFG = -0.04
IALGO = 48
LREAL = Auto
GGA = 91
VOSKOWN = 1
LWAVE = .F.
LCHARG = .F.
还希望虫友们给些修改的意见
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2013-01-05 09:53:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwangan

金虫 (初入文坛)
★
lijuanshu(金币+1): 谢谢参与
把POTIM取小点试试
赞一下 回复此楼
2楼2013-01-05 13:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
★ ★ ★
lijuanshu(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流，， 2013-01-07 15:34:06
数值前面有一个负号还有一个小数点儿，能量是在降低的，只是能量在之后就没有收敛迹象了。
楼主可以先检查一下初始磁矩设置。
赞一下 回复此楼
3楼2013-01-05 19:32:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinhanqing

金虫 (著名写手)
★ ★ ★
lijuanshu(金币+1): 谢谢参与
franch: 金币+2, 鼓励交流，， 2013-01-07 15:34:48
引用回帖:
4楼: Originally posted by yingying1588 at 2013-01-06 18:15:39

我都是把CONTCAR改成POSCAR重新提交一遍就解决的。。这个错误经常遇见。。。
赞一下 回复此楼

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

lijuanshu
5楼2013-01-06 21:45:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
★ ★ ★
lijuanshu(金币+1): 谢谢参与
franch: 金币+2, 鼓励交流，， 2013-01-07 15:34:55
本帖仅楼主可见
赞一下
6楼2013-01-06 23:53:18
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

lijuanshu

铜虫 (正式写手)
引用回帖:
6楼: Originally posted by 旭旭宝宝 at 2013-01-06 23:53:18
看看你的运行脚本 先运行CONTCAR  还是POSCAR

什么意思啊？这有什么区别吗？
赞一下 回复此楼
7楼2013-01-08 08:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
3楼: Originally posted by fanarsenal at 2013-01-05 19:32:54
数值前面有一个负号还有一个小数点儿，能量是在降低的，只是能量在之后就没有收敛迹象了。
楼主可以先检查一下初始磁矩设置。

初始磁矩取的是默认值，这个应该怎么设置啊？
赞一下 回复此楼
8楼2013-01-08 08:55:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
MAGMOM
回复此楼
9楼2013-01-08 12:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
9楼: Originally posted by fanarsenal at 2013-01-08 12:01:50
MAGMOM

我知道是这个参数，但是我在所有的计算中都没有设置它，取的都是默认值，这个有什么影响吗
赞一下 回复此楼
10楼2013-01-09 08:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
送鲜花一朵
引用回帖:
5楼: Originally posted by yinhanqing at 2013-01-06 21:45:34
我都是把CONTCAR改成POSCAR重新提交一遍就解决的。。这个错误经常遇见。。。...

我按照你说的方法，问题解决了，非常感谢啊
赞一下 回复此楼
11楼2013-01-09 08:40:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
手册上有写这个问题，你看一下去吧……
赞一下 回复此楼
12楼2013-01-09 11:01:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Lathander

新虫 (小有名气)
★
lijuanshu(金币+1): 谢谢参与
显然是磁的问题。你之前的magmom是33.x，然后变成32.x，说明至少有一个电子发生转移了。比如
Fe3+Mn2+ magmom=10
-->
Fe2+Mn3+ magmom=8

这当然会引起剧烈的能量变化
赞一下 回复此楼
14楼2013-01-19 14:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
14楼: Originally posted by Lathander at 2013-01-19 14:15:31
显然是磁的问题。你之前的magmom是33.x，然后变成32.x，说明至少有一个电子发生转移了。比如
Fe3+Mn2+ magmom=10
-->
Fe2+Mn3+ magmom=8

这当然会引起剧烈的能量变化

那该如何解决呢？是设置MAGMOM这个参数吗？
赞一下 回复此楼
15楼2013-01-21 08:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
12楼: Originally posted by fanarsenal at 2013-01-09 11:01:18
手册上有写这个问题，你看一下去吧……

我看了下手册，手册上说MAGMOM初始值取NIONS的值，这个值我在OUTCAR中找到了，是不是把MAGMOM取为1.5倍的NIONS就行了
赞一下 回复此楼
16楼2013-01-21 09:39:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by lijuanshu at 2013-01-21 09:39:23
我看了下手册，手册上说MAGMOM初始值取NIONS的值，这个值我在OUTCAR中找到了，是不是把MAGMOM取为1.5倍的NIONS就行了...

手册的意思是默认磁矩为1，它还有根据实验值设置的部分啊，因为初始值偏小可能导致收敛不到（VASP老家那个东西，自己搜一下）。你提到的这个1.5是说每个原子磁矩都为1.5……
赞一下 回复此楼
17楼2013-01-21 10:08:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
17楼: Originally posted by fanarsenal at 2013-01-21 10:08:34
手册的意思是默认磁矩为1，它还有根据实验值设置的部分啊，因为初始值偏小可能导致收敛不到（VASP老家那个东西，自己搜一下）。你提到的这个1.5是说每个原子磁矩都为1.5……...

是不是找到该材料磁矩的实验值，MAGMOM取实验值的1.2或1.5倍？但悲剧的是找不到实验值啊
赞一下 回复此楼
18楼2013-01-21 12:02:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
yingying15884楼
2013-01-06 18:15   回复  
lijuanshu(金币+1): 谢谢参与
lijuanshu13楼
2013-01-09 20:40   回复  
引用回帖:
12楼: Originally posted by fanarsenal at 2013-01-09 11:01:18 手册上有写这个问题，你看一下去吧……

3Q
相关版块跳转 我要订阅楼主 lijuanshu 的主题更新
181/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 食品工程专硕一志愿中海洋309求调剂 +5 小张zxy张 2026-03-26 10/500 2026-03-31 00:29 by jp9609
[考研] 311求调剂一志愿合肥工业大学 +3 秋二十二 2026-03-30 3/150 2026-03-30 23:59 by jp9609
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 286求调剂 +5 Faune 2026-03-30 5/250 2026-03-30 17:15 by 无际的草原
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 290求调剂 +3 dfffsar 2026-03-29 3/150 2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 275求调剂 +15 Micky11223 2026-03-25 20/1000 2026-03-29 20:44 by 唐沐儿
[考研] 0856材料化工调剂 总分330 +14 zhubinhao 2026-03-27 14/700 2026-03-29 10:01 by Sjndkwm
[考研] 一志愿太原理工安全工程300分，求调剂 +5 0857求调剂. 2026-03-24 6/300 2026-03-28 22:04 by zhq0425
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 086000调剂 +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
信息提示
关闭
请填处理意见
关闭