|||
结构优化过程错误提示: ZBRENT: fatal error in bracketing please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue ############################################################################################### 根据输出文件给出优化方法仍不收敛,采用以下几法。单独更改一种参数即可。 该错误一般是由于电荷重叠引起的或对称性不匹配引起的 INCAR内容如下: SYSTEM = Various- local optimisation PREC = Accurate ENCUT = 400.0 【增加ENCUT,如 460 eV】 (不建议,收敛性测试做好!) EDIFF = 1e-7 【改变EDIFF,如EDIFF= 1e-7 】 #测试通过,能量略高 EDIFFG = -1e-3 #SYMPREC=1e-6 IBRION = 2 【改变IBRION方式,如IBRION=1】 #测试通过,能量最低 POTIM = 0.05 【增加POTIM如,POTIM=0.1】【为加快收敛,可以减小KPOINTS】 #测试通过,不建议改KPOINTS ISIF = 3 【对于包含20几个原子的体系,8个核优化,一个结构耗时~20h正常】 NSW = 100 【或改动/调整POSCAR中atomi position,例如改动小数点后的第5位及其以后的数值】 PSTRESS = $i ISMEAR = 1 SIGMA = 0.2 ISYM = 0 LCHARG = FALSE LWAVE = FALSE #LREAL= FALSE IBRION=2在搜寻能量最低点的时候是以下过程: {再次计算受力、能量; if 当前位置是这条路径上的最低点->计算原子受力,乘以POTIM后让离子移动一定的距离 else 根据当前受力,估计可能的最低点(仅限在此一维路径上继续跑) } 你所遇到的错误是因为计算误差或最低点附近能量曲线过于平坦之类的原因无法找到此路径上的最低点。如果确认输入文件没有问题的话,没有关系,照他说的从当前CONTCAR继续就行了(即停止寻找最低点,直接跳出循环外找下一个路径) 如果反复出现这样的问题,不妨换IBRION=1试试。 tips: IBRION和ISMEAR以及各种MIX,只要最后收敛了,收敛到的都是局域的能量最低点,不会因为方法不同而使基态变化。IBRION=2(CG方法)常用于大范围地寻找能量低点,IBRION=1(quasi-Newton)在搜寻范围较小但能得到精确的结果(CG基本EDIFFG=-0.01就是极限了,当然0.01已经很精确了) http://muchong.com/t-4816372-3 |
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-19 21:24
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社