vasp优化结构报错ZBRENT: fatal error in bracketing
各位大神，刚刚使用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 !
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
金币不多，大家见谅，附上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
欢迎大家发言
這問題偶爾會出現
我想，無論是誰跳出?斫忉專琕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 是有機會可以解決這問題，
感谢你的回复，让我对 ZBRENT 这类错误有了更加清晰的认识
今天大致测试了一下，仅以本体系跟大家说一下测试的结果：
1. 仅改变了IBRION=1，但计算达到了NSW的上限（200）仍然未收敛.
2. 直接加上ADDGRID=TRUE，仍然一样的报错，计算失败
3. 改变IBRION=1 + ADDGRID=TRUE , 这个方法似乎比较有效，计算正常，几十步就收敛了
对于不同的体系可能有不同的结果，谨以此参考。
大家有什么建议或者经验，欢迎留言！
我沒有在用QQ 耶~ 呵
我沒有在用QQ 耶~ 呵