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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP优化结构时，E0发生突变？
51/1返回列表
查看: 1607  |  回复: 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的回帖

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的回帖
查看全部 18 个回答

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的回帖
简单回复
lijuanshu13楼
2013-01-09 20:40   回复  
引用回帖:
12楼: Originally posted by fanarsenal at 2013-01-09 11:01:18 手册上有写这个问题，你看一下去吧……

3Q
查看全部 18 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭