| 查看: 458 | 回复: 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 |
回复此楼» 猜你喜欢
机械工程264学硕求调剂
已经有3人回复
-
264求调剂
已经有9人回复
-
求调剂
已经有6人回复
-
22408 266求调剂
已经有7人回复
-
求调剂
已经有14人回复
-
273求调剂
已经有43人回复
-
307求调剂
已经有13人回复
-
327求调剂
已经有7人回复
-
338求调剂
已经有6人回复
-
326分，一志愿沪9，求生物学调剂
已经有3人回复
