| 查看: 715 | 回复: 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 ] |
回复此楼» 猜你喜欢
298求调剂
已经有3人回复
-
博士自荐
已经有7人回复
-
博士推荐
已经有4人回复
-
求环氧树脂研发1名
已经有10人回复
-
280求调剂
已经有5人回复
-
什么是人一生最重要的？
已经有10人回复
-
面上可以超过30页吧？
已经有13人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有13人回复
-
版面费该交吗
已经有17人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 两相流模拟，残差图一段时间后一直跌宕着，没有收敛的倾向。。。
已经有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
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
3楼2012-01-16 12:46:20