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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 结构优化完结果老出现ZBRENT: fatal error in bracketing
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
313/4返回列表上一页1234下一页
查看: 6091  |  回复: 30
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

rainxin2008

新虫 (小有名气)
引用回帖:
18楼: Originally posted by mousekingadv at 2012-08-12 19:31:31
表面的晶格大小是固定的，弛豫的是离子位置...

你说的我不太懂，我的是最下面一层固定F，其余的全是放松的T，师姐们的都是ISIF=0。不管我改ISIF=0还是2，结果都是不收敛，你能帮帮我吗？谢谢你
赞一下 回复此楼
21楼2012-08-12 19:41:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
WDD880227: 金币+2, 多谢指教 2012-08-13 19:04:45
rainxin2008: 金币+10, ★有帮助 2012-08-14 09:21:10
引用回帖:
21楼: Originally posted by rainxin2008 at 2012-08-12 19:41:47
你说的我不太懂，我的是最下面一层固定F，其余的全是放松的T，师姐们的都是ISIF=0。不管我改ISIF=0还是2，结果都是不收敛，你能帮帮我吗？谢谢你...

IBRION=2在搜寻能量最低点的时候是以下过程：
计算原子受力，乘以POTIM后让离子移动一定的距离，
{再次计算受力、能量；
  if 当前位置是这条路径上的最低点->计算原子受力，乘以POTIM后让离子移动一定的距离
  else 根据当前受力，估计可能的最低点（仅限在此一维路径上继续跑）
}
你所遇到的错误是因为计算误差或最低点附近能量曲线过于平坦之类的原因无法找到此路径上的最低点。如果确认输入文件没有问题的话，没有关系，照他说的从当前CONTCAR继续就行了（即停止寻找最低点，直接跳出循环外找下一个路径）
如果反复出现这样的问题，不妨换IBRION=1试试。
赞一下 回复此楼
22楼2012-08-13 01:20:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
22楼: Originally posted by mousekingadv at 2012-08-13 01:20:32
IBRION=2在搜寻能量最低点的时候是以下过程：
计算原子受力，乘以POTIM后让离子移动一定的距离，
{再次计算受力、能量；
  if 当前位置是这条路径上的最低点->计算原子受力，乘以POTIM后让离子移动一定的距离 ...

恩，谢谢你，我把IBRION改为1了，ISIF=0.现在在重算，看结果能不能收敛了。。。
赞一下 回复此楼
23楼2012-08-13 09:29:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
22楼: Originally posted by mousekingadv at 2012-08-13 01:20:32
IBRION=2在搜寻能量最低点的时候是以下过程：
计算原子受力，乘以POTIM后让离子移动一定的距离，
{再次计算受力、能量；
  if 当前位置是这条路径上的最低点->计算原子受力，乘以POTIM后让离子移动一定的距离 ...

你好，我把IBRION改为1算的结果如下： 195 F= -.17800925E+03 E0= -.17800569E+03  d E =-.497371E-02
BRION: g(F)=  0.255E-01 g(S)=  0.000E+00 retain N=  4 mean eig= 1.98
eig:   7.417  0.028  0.240  0.240
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.178021946329E+03   -0.12684E-01   -0.37661E+00  3273   0.301E+00    0.898E-02
RMM:   2    -0.178025353881E+03   -0.34076E-02   -0.38320E-02  3261   0.298E-01    0.686E-02
RMM:   3    -0.178025396033E+03   -0.42152E-04   -0.13120E-03  3267   0.684E-02    0.635E-02
RMM:   4    -0.178025320041E+03    0.75992E-04   -0.10529E-04  1793   0.188E-02
196 F= -.17802532E+03 E0= -.17802177E+03  d E =-.160742E-01
BRION: g(F)=  0.644E+00 g(S)=  0.000E+00 retain N=  5 mean eig= 1.59
eig:   7.444  0.235  0.235  0.022  0.022
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.178025456026E+03   -0.59993E-04   -0.41268E-03  3225   0.988E-02    0.842E-03
RMM:   2    -0.178025458590E+03   -0.25644E-05   -0.44000E-05  1737   0.115E-02    0.621E-03
RMM:   3    -0.178025457360E+03    0.12308E-05   -0.97267E-06  1648   0.719E-03
197 F= -.17802546E+03 E0= -.17802193E+03  d E =-.137319E-03
BRION: g(F)=  0.143E-01 g(S)=  0.000E+00 retain N=  1 mean eig= 0.00
eig:   0.005
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.178034835671E+03   -0.93771E-02   -0.36337E+00  3274   0.296E+00    0.937E-02
RMM:   2    -0.178038051320E+03   -0.32156E-02   -0.37021E-02  3262   0.293E-01    0.694E-02
RMM:   3    -0.178038089780E+03   -0.38460E-04   -0.12707E-03  3270   0.666E-02    0.629E-02
RMM:   4    -0.178038012742E+03    0.77038E-04   -0.10441E-04  1793   0.188E-02
198 F= -.17803801E+03 E0= -.17803453E+03  d E =-.125554E-01
BRION: g(F)=  0.175E+00 g(S)=  0.000E+00 retain N=  2 mean eig= 0.01
eig:   0.005  0.019
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.178038061183E+03    0.28597E-04   -0.51879E-02  3264   0.343E-01    0.137E-02
RMM:   2    -0.178038107135E+03   -0.45953E-04   -0.49736E-04  2332   0.336E-02    0.854E-03
RMM:   3    -0.178038108580E+03   -0.14447E-05   -0.34037E-05  1706   0.123E-02
199 F= -.17803811E+03 E0= -.17803467E+03  d E =-.958382E-04
BRION: g(F)=  0.746E-02 g(S)=  0.000E+00 retain N=  1 mean eig= 0.03
eig:   0.032
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.178043166425E+03   -0.50593E-02   -0.17555E+00  3274   0.208E+00    0.772E-02
RMM:   2    -0.178044778060E+03   -0.16116E-02   -0.17710E-02  3262   0.204E-01    0.560E-02
RMM:   3    -0.178044782047E+03   -0.39871E-05   -0.59850E-04  2607   0.463E-02
200 F= -.17804478E+03 E0= -.17804144E+03  d E =-.667347E-02
BRION: g(F)=  0.937E+00 g(S)=  0.000E+00 retain N=  2 mean eig= 0.02
eig:   0.032  0.004
start time  Mon Aug 13 14:54:12 CST 2012
end   time  Tue Aug 14 01:51:52 CST 2012
。。。。。。。。。。。。。。。。
是不是说还没有算完？算到200步就结束了，我的NSW设的是200.
还是就是我的是结构优化，IBRION设为1和2算的结果对我的结构有影响没？是不是为
赞一下 回复此楼
24楼2012-08-14 09:15:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)
引用回帖:
24楼: Originally posted by rainxin2008 at 2012-08-14 09:15:21
你好，我把IBRION改为1算的结果如下： 195 F= -.17800925E+03 E0= -.17800569E+03  d E =-.497371E-02
BRION: g(F)=  0.255E-01 g(S)=  0.000E+00 retain N=  4 mean eig= 1.98
eig:   7.417  0.028  0.240  0 ...

确实没算完。只是因为NSW到了就停了
IBRION和ISMEAR以及各种MIX，只要最后收敛了，收敛到的都是局域的能量最低点，不会因为方法不同而使基态变化。IBRION=2（CG方法）常用于大范围地寻找能量低点，IBRION=1（quasi-Newton）在搜寻范围较小但能得到精确的结果（CG基本EDIFFG=-0.01就是极限了，当然0.01已经很精确了）
为什么能量比你之前贴出来的下降了那么多（2eV）？五个原子吸附上去不可能有这么大的能量差吧。你的衬底原子间距是怎么得到的？
另外建议过一段时间就用VESTA看看CONTCAR，确认在跑的东西是你想要的。
赞一下 回复此楼
25楼2012-08-14 10:30:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainxin2008

新虫 (小有名气)
引用回帖:
25楼: Originally posted by mousekingadv at 2012-08-14 10:30:14
确实没算完。只是因为NSW到了就停了
IBRION和ISMEAR以及各种MIX，只要最后收敛了，收敛到的都是局域的能量最低点，不会因为方法不同而使基态变化。IBRION=2（CG方法）常用于大范围地寻找能量低点，IBRION=1（quas ...

我的衬底原子间距就是直接在MS里导入的Au,然后建的模型，间距也没有手动修改，不知道为什么能量下降了那么多。怎么这么多问题啊，哎，太多不懂了
赞一下 回复此楼
26楼2012-08-14 10:37:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2012-08-14 18:59:27
引用回帖:
26楼: Originally posted by rainxin2008 at 2012-08-14 10:37:15
我的衬底原子间距就是直接在MS里导入的Au,然后建的模型，间距也没有手动修改，不知道为什么能量下降了那么多。怎么这么多问题啊，哎，太多不懂了...

晶格常数不能使用实验数据，必须用相同的方法和赝势优化一遍块材Au。不需要真空层。ENCUT和K点不同也会有影响，不过差距不是很大。
赞一下 回复此楼
27楼2012-08-14 12:39:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)
改IBRION=1，就可以了，算到200步还不收敛，就先把POTIM改成0.5，收敛后再改成0.3左右继续算，试试吧
赞一下 回复此楼
28楼2012-11-15 16:09:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李加伟

荣誉版主 (文坛精英)

不过凡人

优秀版主优秀版主
遇到相同问题，请问楼主后来怎么处理的
赞一下 回复此楼
人生海海，在于内心
29楼2014-03-25 08:26:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微6凉

铁虫 (小有名气)
引用回帖:
6楼: Originally posted by rainxin2008 at 2012-08-11 11:03:37
我的用CONTCAR替换POSCAR又算了一遍，还是不收敛。怎么办啊啊？纳了闷了...

请问楼主问题解决了吗？我也遇到了同样的问题
赞一下 回复此楼
30楼2015-08-11 16:47:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rainxin2008 的主题更新
313/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭