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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp优化结构报错ZBRENT: fatal error in bracketing
51/1返回列表
查看: 9658  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

李玥玥

铜虫 (正式写手)

[求助] vasp优化结构报错ZBRENT: fatal error in bracketing已有1人参与

各位大神，刚刚使用vasp不久，优化结构总是遇到一些问题，根据网上查询的办法似乎都没有什么效果，所以请有懂的大神能够给予一点指导，感激不尽。

ZBRENT: fatal error in bracketing
     please rerun with smaller EDIFF, or copy CONTCAR
     to POSCAR and continue

我是采用ISIF=2,ISIF=7,ISIF=2,ISIF=7的方式，来优化结构的，对于unit cell来说，优化的很好。然后超胞后，改变锂离子的数目，优化一些结构。再优化这些结构时，先采用ISIF=2优化成功，然后ISIF=7优化晶格常数，但这一步一直无法成功，总是出现上面的错误。
试过的方法：
1. 根据提示CONTCAR变POSCAR继续优化，或者使用更小的EDIFF，都仍然一样的报错。
2. IBRION改为1计算，现在已经算了120个离子步了还未收敛，优化结构后面实时更新。
3. IBRION=1，并且ADDGRID=TRUE进行优化，优化结构后面更新。（ADDGRID简单测试了一下，似乎对能量没有什么影响）
4. IBRION不变，仍然=2，加入ADDGRID=TRUE进行优化，结果后面更新。

现在通过查资料找到的就这些方法。

请问大家：
1. 这个报错的原因是什么？
2. 如何有效解决这个报错呢？
3. 看到有人说这个跟vasp的编译有关？请问这个该如何理解，具体指什么相关呢？  （本人不太相信这种原因）
4. ADDGRID的具体含义是什么，有什么影响？vasp计算好像很少看到这个参数。

本人的输入文件均已上传。

请大家能够花点时间一起讨论一下，Thanks a lot !
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : POSCAR
    • 2017-07-17 10:07:07, 12.01 K
  • 附件 2 : POTCAR
    • 2017-07-17 10:28:46, 613.53 K

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

Spintronics

» 猜你喜欢
1楼2017-07-17 10:34:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
李玥玥: 金币+15, ★★★很有帮助 2017-07-17 19:45:03
漫天飘雪: 金币+2, 谢谢交流 2017-07-17 20:45:33
這問題偶爾會出現
我想，無論是誰跳出?斫忉專琕ASP作者群的公信力仍是比較高的

http://cms.mpi.univie.ac.at/vasp ... ?f=3&t=7415
裡面的 admin (也就是 VASP 官方論壇的管理員，身份雖不確定，但~ 一定跟 VASP 作者群有直接關連)
please have a look at the OUTCAR file, at the forces of the last ionic step before Brent's algorithm fails
-- if the calculation is very well-converged already, please proceed using IBRION=1; ADDGRID=.TRUE. and a larger ENCUT
這一點是與你提出的解釋方法 3 是一致的
不過，前提是如果你的計算已經 "快要" 收斂了
-- if this error appears at an early stage of the relaxation, pelase check if your input geometry was reasonable, eg by having a look at the interatomic distances (written in OUTCAR) or the forces of the first geometries. It may then help to increase the accuracy (PREC= Accurate and a higher ENCUT) and to decrease POTIM to avoid too large stepwidhts in the first relaxation step.
若是你的計算離收斂還很遠，請檢查你的結構是否合理
檢查 OUTCAR 內第一個 ionic loop 的 FORCE，看看是否有原子出現不尋常大的受力
透過 PREC=Accurate & ENCUT 加大 & 降低 POTCAR，也許會有幫助

http://cms.mpi.univie.ac.at/vasp ... p?f=3&t=869
this message says that during geometry optimisation no reasonable next step could be found with Brent's algorithm (linear bisectioning)
please check the following
1) is the system converged already? (look at the forces in OUTCAR, especially if you use the total energy change as congvergence criterium for ionic optimisation)
2) if not: is each ionic step converged electronically ( this is needed to obtain reasonable forces)
這篇的回應剛上篇是大同小�?，此外，admin 在下面又再提到說
sometimes, if the calculation is already highly converged, vasp cannot interpolate the next step to within the numerical accuracy because it simply would be so small. (this behaviour is machine-dependent). If you have used IBRION=2, it may help to switch to IBRION=1 and set ADDGRID=.true. in addition
此點是與你說的解決方法 4 是一致

此外
http://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?t=1856
ZBRENT is an algorithm search for a root of a function by Brent's method: Numerical recipies, Section 9.3
The problem in your case might be that the Conjugate Gradient algorithm (IBRION=2) is not suitable for very small corrections of the atomic positions if your system has almost reached equilibrium (please have a look at XDATCAR to check the size of the relaxation steps done before the ZBRENT warnings show up). Usually it is sufficient to converge a system up to maximum remaining forces of about 0.01eV/A (EDIFFG=-0.01).
please try one of the following:
1) choose a different algorithm for ionic optimization (IBRION=1)
2) set ADDGRID=.True. in INCAR (only for vasp releases 4.4.5 and newer)
Admin 針對 ZBRENT 也有提出說明
ZBRENT 是數值方法內的 Brent's method
在你計算已經很接近收斂時，VASP 所做的 displacement 會變的很小
就有可能會使得這個 Brent's method 產生問題


個人經驗

這問題大多出現在 structure optimization 時
通常就是VASP在找 optimized position 產生了震盪

就好像是 (以下是想像，數字都是我隨意給的)
某個原子最佳位置是在 0.5
而現在位置是 0.4，每次VASP會移動 0.2
原子就會?� 0.4 跳到 0.6 然後又跳回 0.4 又再跳回 0.6
所以計算一直無法正常停止
結構與能量也會是在某兩個值間一直震盪
這也是為什麼改變 IBRION 是有機會可以解決這問題
赞一下(4人) 回复此楼
3楼2017-07-17 18:02:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

李玥玥

铜虫 (正式写手)
金币不多，大家见谅，附上INCAR和KPOINTS，上传出问题了。。。。
INCAR:
SYSTEM = LiFePO4_1

######Start Parameter for This Run
ISTART = 0
ICHARG = 2
INIWAV = 1

######Electronic Relaxation
ENCUT = 520
NELM = 60
EDIFF = 1E-5
GGA = PE
ISMEAR = 0
SIGMA = 0.05
PREC = Accurate
ISPIN = 2
#LORBIT = 10

######Ion Relaxation
IBRION = 2
ISIF = 2
NSW = 200
EDIFFG = -0.01
POTIM = 0.5
ALGO = Fast
LREAL = A
NPAR=4

######LDA+U
LMAXMIX = 4
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 -1 -1 -1
LDAUU = 5.3 0 0 0
LDAUJ = 1.0 0 0 0

#for out
LCHARG = .FALSE.
LWAVE = .FALSE.
LVTOT = .FALSE.
LELF = .FALSE.


KPOINTS:
k-points
0
Monkhorst Pack
4   3   4
0   0   0

欢迎大家发言
赞一下 回复此楼
2楼2017-07-17 10:39:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)
引用回帖:
3楼: Originally posted by jpchou at 2017-07-17 18:02:42
這問題偶爾會出現
我想，無論是誰跳出?斫忉專琕ASP作者群的公信力仍是比較高的

http://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?f=3&t=7415
裡面的 admin (也就是 VASP 官方論壇的管理員，身份雖不 ...

感谢你的回复，让我对 ZBRENT 这类错误有了更加清晰的认识
赞一下 回复此楼
4楼2017-07-17 19:50:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)
今天大致测试了一下，仅以本体系跟大家说一下测试的结果：
1. 仅改变了IBRION=1，但计算达到了NSW的上限（200）仍然未收敛.
2. 直接加上ADDGRID=TRUE，仍然一样的报错，计算失败
3. 改变IBRION=1 + ADDGRID=TRUE , 这个方法似乎比较有效，计算正常，几十步就收敛了

对于不同的体系可能有不同的结果，谨以此参考。

大家有什么建议或者经验，欢迎留言！
赞一下(2人) 回复此楼
5楼2017-07-17 19:55:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭