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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP优化结构时，E0发生突变？
181/1返回列表
查看: 1594  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭