| 查看: 499 | 回复: 0 | ||
| 当前主题已经存档。 | ||
[资源]
【分享】Trust radius
|
||
|
Today, I try to opitimize carbon atom chain with 4 atoms per cell. When I got the following atom position, it failed. ATOMIC_POSITIONS (crystal) C 0.001769728 0.000000000 0.000000000 C 0.251731600 0.000000000 0.000000000 C 0.501768886 0.000000000 0.000000000 C 0.751729786 0.000000000 0.000000000 From the above data, it is easy to see that these four atoms almost equally divide the primitive vector of x axis. Besides, I find more information in .rx.out file: trust_radius < trust_radius_min Total force = 0.000080 Total SCF correction = 0.000222 SCF correction compared to forces is too large, reduce conv_thr history already reset at previous step Furthermore, I find useful guide at http://www.democritos.it/piperma ... ctober/005249.html: "My advise is that you should not reduce trust_radius_min. Instead, increase the accuracy of your forces by reducing conv_thr even more (or by increasing upscale). You can monitor the error on the forces due to lack of self-consistency by checking the size of total scf correction (it must be much smaller than the total force). However this is not the only source of numerical errors. Since some terms are computed in real space on a equispaced grid, the energy of the system is not invariant for arbitrary translations. This can easily confuse BFGS algorithm when the convergence criteria are too tight.To reduce this error you should increase ecutrho. However I don't see the reason to waste cpu-time by using such a tight criteria. Also, consider that the self-consistency itself is going to fail if you ask conv_thr < 1.D-12 because the iterative diagonalization schemes we use can suffer numerical instabilities when the required accuracy for eigenvalues/eigenvectors is to high." From the expert's explanation, I know the accuracy of forces is close related with the conv_thr of scf. By the way,if atom positions are almost equal,trust_radius_ini should be reduced. 以上内容来自我的博客：blog.sina.com.cn/nkasir |
回复此楼» 猜你喜欢
面上提前没消息，有中的吗
已经有6人回复
-
2026年国自然面上资助率
已经有21人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有10人回复
-
2027年申博
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复