对于二维材料的结构驰豫，在真空层方向一直不收敛

静心68

我设置的力的收敛精度为1e-3 eV/A ， 对于块体材料很快就可以达到收敛精度，但是对于含有真空层的二维材料，在运行了800步结构驰豫以后，各原子在真空层方向的力一直也保持在1e-1到1e-2 eV/A 的水平（从一开始差不多就是这个水平）。
我的输入文件与log文件在如下链接中，论坛暂时有文件上传大小限制。
https://github.com/abacusmodeling/abacus-develop/issues/246

zjdai

我测试了你计算的二维MoS构型，使用的ABACUS版本为abacus.3.0.4。测试发现使用：
1. cg算法进行结构优化能量收敛已经受力收敛优化很慢，计算700个离子步才达到收敛精度(<1e-3eV/A);
2. 使用md的fire算法能量优化和受力收敛会大大加快，100个离子步可以达到收敛精度；
3. 看了你的INPUT文件，你的自洽收敛精度为1e-5，对于自洽计算和结构优化过低，导致你每一个离子步电荷密度已经能量结果不够好，受力误差就会很大，这也许是你计算800个离子步驰豫不能收敛原因；
4. ABACUS计算速度比较快，我是用56个进程计算，使用md中fire算法，56min就可以计算完100个离子步计算；
5. 至于cg方法计算收敛比较慢点原因，我不能给你很好的解释，也许和你计算的体系相关，或者设置的参数不够合理，这个需要进一步测试和讨论；
6. 附上我计算的INPUT文件
INPUT_PARAMETERS
#Parameters     (General)
ntype                   2
#Parameters (Accuracy)
ecutwfc                 100
scf_nmax                100
scf_thr                 1e-8
basis_type              lcao
smearing_method         gauss
smearing_sigma          0.01
mixing_type             pulay
mixing_beta             0.2


calculation             md
md_type                 -1
md_nstep                800
force_thr_ev            1.0e-3
md_tfirst               0
cal_force               1
cal_stress              0
vdw_method              d3_0

symmetry                0
6

hermit

建议直接使用 ABACUS 的默认电荷密度收敛参数。

静心68

zjdai 发表于 2022-12-16 22:55
我测试了你计算的二维MoS构型，使用的ABACUS版本为abacus.3.0.4。测试发现使用：
1. cg算法进行结构优化能 ...

好的，十分感谢，祖建
我按照你的提示再尝试一下驰豫

静心68

hermit 发表于 2022-12-17 09:34
建议直接使用 ABACUS 的默认电荷密度收敛参数。

好的，我现在设置到推荐的精度