| 查看: 446 | 回复: 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 |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有19人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复
-
280求调剂
已经有4人回复
-
面上可以超过30页吧？
已经有12人回复
-
什么是人一生最重要的？
已经有6人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有11人回复
-
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有9人回复
-
版面费该交吗
已经有15人回复