| 查看: 759 | 回复: 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 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 两相流模拟，残差图一段时间后一直跌宕着，没有收敛的倾向。。。
已经有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
