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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 结构优化完结果老出现ZBRENT: fatal error in bracketing
171/2返回列表12下一页
查看: 6121  |  回复: 30
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

rainxin2008

新虫 (小有名气)

[求助] 结构优化完结果老出现ZBRENT: fatal error in bracketing

我用VASP优化完，查看结果时老出现 ZBRENT: fatal error in bracketing。please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue。刚开始INCAR中的EDIFF=1e-4，不收敛，后来改成1e-5，结果还是出现上面那句话。谁能帮我分析下原因。
INCAR文件如下：
SYSTEM    =  Azoc MD

     PREC    = Medium
     EDIFF   = 1e-5
     EDIFFG  = -3e-2     
     IALGO   = 48      
     NPAR   =  4
     LPLANE = .TRUE.
     
     ENCUT = 400
     NELMIN  = 3         
     LREAL   = On
     ROPT    = 1E-3 1E-3  1E-3  1E-3
     ISYM    = 0      
     NWRITE  = 1

  SPIN POLARIZATION SPECIFICATION
     ISPIN   = 1

IONIC RELAXATION
     NSW     = 200     
     NBLOCK  = 1      
     IBRION  = 2
     ISIF    = 0      
     POTIM   = 0.4   
     LWAVE = .F.
     LCHARG  = .F.

  ELECTRON OCCUPATION
     ISMEAR = 2  

  DOS Stuff
     SIGMA   = 0.1       check T*S to make sure
KPOINTS：3x3x1
POTCAR和POSCAR应该没问题的
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 wly的好资源专辑

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-08-10 21:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
还有，我的POSCAR也贴下来

CIF file
1.0
        9.9898004532         0.0000000000         0.0000000000
        4.9949002266         8.6514209712         0.0000000000
        0.0000000000         0.0000000000        22.0638008118
   Au    C    H    S
   48    1    4    1
Direct
0               0               0.028689999     F F F   
0.166669995     0.166669995     0.028689999     F F F   
0.333330005     0.333330005     0.028689999     F F F   
0.5             0               0.028689999     F F F   
0.666670024     0.166669995     0.028689999     F F F   
0.833329976     0.333330005     0.028689999     F F F   
0               0.5             0.028689999     F F F   
0.166669995     0.666670024     0.028689999     F F F   
0.333330005     0.833329976     0.028689999     F F F   
0.5             0.5             0.028689999     F F F   
0.666670024     0.666670024     0.028689999     F F F   
0.833329976     0.833329976     0.028689999     F F F   
0.166669995     0               0.135409996     T T T   
0.333330005     0.166669995     0.135409996     T T T   
0               0.333330005     0.135409996     T T T   
0.666670024     0               0.135409996     T T T   
0.833329976     0.166669995     0.135409996     T T T   
0.5             0.333330005     0.135409996     T T T   
0.166669995     0.5             0.135409996     T T T   
0.333330005     0.666670024     0.135409996     T T T   
0               0.833329976     0.135409996     T T T   
0.666670024     0.5             0.135409996     T T T  
0.833329976     0.666670024     0.135409996     T T T   
0.5             0.833329976     0.135409996     T T T   
0.333330005     0               0.242119998     T T T   
0               0.166669995     0.242119998     T T T   
0.166669995     0.333330005     0.242119998     T T T   
0.833329976     0               0.242119998     T T T   
0.5             0.166669995     0.242119998     T T T   
0.666670024     0.333330005     0.242119998     T T T   
0.333330005     0.5             0.242119998     T T T   
0               0.666670024     0.242119998     T T T   
0.166669995     0.833329976     0.242119998     T T T   
0.833329976     0.5             0.242119998     T T T   
0.5             0.666670024     0.242119998     T T T   
0.666670024     0.833329976     0.242119998     T T T   
0               0               0.348839998     T T T   
0.166669995     0.166669995     0.348839998     T T T   
0.333330005     0.333330005     0.348839998     T T T   
0.5             0               0.348839998     T T T   
0.666670024     0.166669995     0.348839998     T T T   
0.833329976     0.333330005     0.348839998     T T T   
0               0.5             0.348839998     T T T   
0.166669995     0.666670024     0.348839998     T T T   
0.333330005     0.833329976     0.348839998     T T T   
0.5             0.5             0.348839998     T T T   
0.666670024     0.666670024     0.348839998     T T T   
0.833329976     0.833329976     0.348839998     T T T
0.303519994     0.59246999      0.494740009     T T T
0.220149994     0.659539998     0.456409991     T T T      
0.550610006     0.570360005     0.509270012     T T T   
0.277330011     0.500820041     0.512490034     T T T   
0.289310008     0.675300002     0.532989979     T T T
0.5             0.5             0.466679999     T T T
赞一下 回复此楼
2楼2012-08-10 21:13:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
教师同学们赶紧交流啊，很着急的啊
赞一下 回复此楼
3楼2012-08-11 08:23:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
4楼: Originally posted by babaleo at 2012-08-11 09:34:01
提示里的" copy CONTCAR to POSCAR and continue" 试过吗? 一般情况这一步都能解决问题.

恩，这一步正在试，谢谢你了！
赞一下 回复此楼
5楼2012-08-11 10:37:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
4楼: Originally posted by babaleo at 2012-08-11 09:34:01
提示里的" copy CONTCAR to POSCAR and continue" 试过吗? 一般情况这一步都能解决问题.

我的用CONTCAR替换POSCAR又算了一遍，还是不收敛。怎么办啊啊？纳了闷了
赞一下 回复此楼
6楼2012-08-11 11:03:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
4楼: Originally posted by babaleo at 2012-08-11 09:34:01
提示里的" copy CONTCAR to POSCAR and continue" 试过吗? 一般情况这一步都能解决问题.

第一次生成的CONTCAR用软件看了后结构变的微小。而我把用CONTCAR改为POSCAR后又算了一遍，生成的CONTCAR再用软件看了一下图发现结构变的比较厉害了。这是怎么回事？
赞一下 回复此楼
7楼2012-08-11 11:20:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
8楼: Originally posted by mousekingadv at 2012-08-11 16:09:32
为什么会是ISIF=0？

我的INCAR 是直接从师姐那里拷的，用的原子都一样，所以没改。我这是第一步的结构优化，之前从没有算过这些，也是刚接触VASP第一次提交计算，谁知道就老不收敛，着急的很。IBRION=2,ISIF=0不合适吗？你帮我看看应该改为多少？
赞一下 回复此楼
9楼2012-08-11 21:47:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
11楼: Originally posted by forestwolf9291 at 2012-08-11 23:17:50
网上搜到官方解答说，查看outcar里的力常数看是否已经收敛，如不是，给了一组参数设置让尝试一下。楼主自己搜索看看。
另外，我也几次遇到，印象中将SYMPREC = 0.00001改为0。0001就可以了。对结果有无大的影响也没 ...

这是我OUTCAR里最后几行： 155 F= -.17784496E+03 E0= -.17783980E+03  d E =-.429362E-03
curvature:   0.00 expect dE=-0.932E-07 dE for cont linesearch -0.314E-09
ZBRENT:  can not reach accuracy
opt :   0.0000  next Energy=     0.000000 (dE= 0.178E+03)
bond charge predicted
       N       E                     dE             d eps       ncg     rms     
     rms(c)
RMM:   1    -0.177844506136E+03    0.45596E-03   -0.23214E-02  3276   0.237E-01
   0.778E-03
RMM:   2    -0.177844525942E+03   -0.19806E-04   -0.19793E-04  1670   0.236E-02
   0.464E-03
RMM:   3    -0.177844529056E+03   -0.31147E-05   -0.32983E-05  1661   0.121E-02
156 F= -.17784453E+03 E0= -.17783937E+03  d E =0.368263E-05
curvature:   0.00 expect dE= 0.122-306 dE for cont linesearch  0.143-307
ZBRENT: fatal error in bracketing
     please rerun with smaller EDIFF, or copy CONTCAR
     to POSCAR and continue
start time  Sat Aug 11 21:59:02 CST 2012
end   time  Sun Aug 12 00:39:49 CST 2012
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
你说的SYMPREC是哪个输入文件里的参数，我怎么没找到？谢谢
赞一下 回复此楼
13楼2012-08-12 08:42:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
12楼: Originally posted by mousekingadv at 2012-08-12 00:26:18
我感觉这个POSCAR是在Au表面吸附CH4S
那么应该保持原胞大小不变，优化原子位置，使用ISIF=2
IBRION=2没有问题...

恩，昨晚我尝试着改了ISIF=2，然后刚刚查看结果还是一样，OUTCAR
赞一下 回复此楼
14楼2012-08-12 08:42:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
13楼: Originally posted by rainxin2008 at 2012-08-12 08:42:11
这是我OUTCAR里最后几行： 155 F= -.17784496E+03 E0= -.17783980E+03  d E =-.429362E-03
curvature:   0.00 expect dE=-0.932E-07 dE for cont linesearch -0.314E-09
ZBRENT:  can not reach accuracy
op ...

恩，昨晚我改了ISIF=2重算了一下，早上起来查看结果还是一样的。OUTCAR最后几行如下： 155 F= -.17784496E+03 E0= -.17783980E+03  d E =-.429362E-03
curvature:   0.00 expect dE=-0.932E-07 dE for cont linesearch -0.314E-09
ZBRENT:  can not reach accuracy
opt :   0.0000  next Energy=     0.000000 (dE= 0.178E+03)
bond charge predicted
       N       E                     dE             d eps       ncg     rms     
     rms(c)
RMM:   1    -0.177844506136E+03    0.45596E-03   -0.23214E-02  3276   0.237E-01
   0.778E-03
RMM:   2    -0.177844525942E+03   -0.19806E-04   -0.19793E-04  1670   0.236E-02
   0.464E-03
RMM:   3    -0.177844529056E+03   -0.31147E-05   -0.32983E-05  1661   0.121E-02
156 F= -.17784453E+03 E0= -.17783937E+03  d E =0.368263E-05
curvature:   0.00 expect dE= 0.122-306 dE for cont linesearch  0.143-307
ZBRENT: fatal error in bracketing
     please rerun with smaller EDIFF, or copy CONTCAR
     to POSCAR and continue
start time  Sat Aug 11 21:59:02 CST 2012
end   time  Sun Aug 12 00:39:49 CST 2012
赞一下 回复此楼
15楼2012-08-12 08:44:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rainxin2008 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭