| 查看: 728 | 回复: 2 | ||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[求助]
CNEB收敛到0.05后，提高精度不收敛
|
|
（1）q-m算法收敛到1eV/A (2) LBFGS算法熟练到0.05 （3）按照（2）的设置，修改EDIFFG=-0.03，很难收敛。 下面一段话是早前Graeme给某个帖子的回复。看他的建议，似乎是 q-m算法收敛到1,然后CG或者LBFGS熟练到0.05，这时EDIFFG=1E-4或者1E-5；如果要达到更高的精度，需要更高的力的精度，要求EDIFFG = 1E-5或者更低。这时可以退回q-m算法。 及CNEB的熟练过程的算法切换是：Q-M （<1 eV/A）->CG/LBFGS (<0.05 eV/A)->Q-M (<0.03 eV/A). 这样的话，如果（有人报道过这样的结果）力精度提高到0.001,EDIFFG岂不是要设置的很小很小？ Quickmin is a first-order optimizer which is relatively stable at any point during optimization. Near the minimum where the potential can be expanded using a harmonic approximation, it is not as efficient as second-order optimizers such as CG or LBFGS. There are a few caveats. CG and LBFGS are not particularly stable far from the minimum where forces are high. So using quick-min until the forces drop below 1 eV/Ang is a sensible way to start. You can also make sure that the path at this point makes physical sense. Then it will be more efficient to switch to a second order method. When I talk about forces being accurate enough near the minimum, I'm referring to the errors in the forces which are controlled by the electronic convergence criteria, ediff. When the forces are larger than the errors, curvatures can be determined with sufficient accuracy. Near the minimum, the errors are a higher portion of the force. At this point, you can either go back to quick-min (or fire), or decrease ediff. Your value of ediff=1e-8 is extremely tight. Using a value of 1e-5 will generally speed up your calculations without causing problems unless you are trying to get forces down to ~0.05 eV/Ang. Also, as a side note, your cutoff is very high and I have never seen a case where you need to use prec=accurate for the calculation of a barrier. [ Last edited by cenwanglai on 2011-9-27 at 20:45 ] |
回复此楼» 猜你喜欢
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有6人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有4人回复
-
考博
已经有5人回复
-
26年申博自荐-计算机视觉
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 两相流模拟，残差图一段时间后一直跌宕着，没有收敛的倾向。。。
已经有11人回复
- 优选opt后，能量一直不收敛且停在同个状态
已经有5人回复
- TD-DFT激发态优化不收敛
已经有4人回复
- CASTEP 加自旋不收敛
已经有6人回复
- NEB过渡态搜索遇到的问题
已经有33人回复
- 解决非线性分析不收敛的技巧
已经有37人回复
- DMol3计算频率不收敛
已经有4人回复
- siesta 表面不收敛
已经有3人回复
- 【求助】用Orca计算一个片簇的单点能，不收敛。
已经有7人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
04nylxb
木虫 (正式写手)
- 应助: 33 (小学生)
- 金币: 2321.9
- 散金: 46
- 红花: 4
- 帖子: 824
- 在线: 262.6小时
- 虫号: 817223
- 注册: 2009-07-28
- 性别: GG
- 专业: 工程热物理相关交叉领域
2楼2012-01-16 11:10:25
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
3楼2012-01-16 12:46:20