| 查看: 444 | 回复: 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 |
回复此楼» 猜你喜欢
面上项目申报
已经有3人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有5人回复
-
酰胺脱乙酰基
已经有9人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复