| 查看: 780 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 两相流模拟，残差图一段时间后一直跌宕着，没有收敛的倾向。。。
已经有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