northern_mar 发表于 2026-4-20 21:34
个人感觉还有几个小地方也许可以一起讨论一下。
1. x方向和y方向没设置k点，但是晶胞又不够大；
2. 不知 ...

谢谢答复
不加smearing，SCF容易收敛，结构也可以收敛，只是会散的比例厉害
我会尝试加K 点

如果加221 K点，全pop dzvp rev2  ，收敛不太好
是tep     Update method      Time    Convergence         Total energy    Change
  ------------------------------------------------------------------------------
     1 NoMix/Diag. 0.40E+00   12.3     4.15358412    -94736.7063089387 -9.47E+04
     2 Broy./Diag. 0.40E+00   13.9     1.32952109    -94619.6904807796  1.17E+02
     3 Broy./Diag. 0.40E+00   13.7    17.51187137    -94727.3492622163 -1.08E+02
     4 Broy./Diag. 0.40E+00   13.5   120.26863828    -94467.8404698225  2.60E+02
     5 Broy./Diag. 0.40E+00   13.5   162.11694664    -93633.4100162796  8.34E+02
     6 Broy./Diag. 0.40E+00   13.7   231.03799664    -92704.3994921492  9.29E+02
     7 Broy./Diag. 0.40E+00   13.6   231.06615429    -92242.9112017088  4.61E+02
     8 Broy./Diag. 0.40E+00   13.7   163.87664050    -90074.2296181857  2.17E+03
     9 Broy./Diag. 0.40E+00   13.7   122.28343960    -89881.3651945907  1.93E+02
    10 Broy./Diag. 0.40E+00   13.4   175.46421865    -85390.7740296135  4.49E+03
    11 Broy./Diag. 0.40E+00   13.6   227.75519169    -70859.0561701974  1.45E+04
    12 Broy./Diag. 0.40E+00   13.7   227.75536503    -65671.5757085649  5.19E+03
    13 Broy./Diag. 0.40E+00   13.8   217.70651379    -65908.6503469109 -2.37E+02

个人感觉还有几个小地方也许可以一起讨论一下。
1. x方向和y方向没设置k点，但是晶胞又不够大；
2. 不知道这个例子的收敛情况如何，自旋多重度不一定设置合理；
3. 如同上一位老师说的，含Ni的体系CUTOFF设置太低了。如果耗时太长可以考虑OT+全电子基组；
4. 除非特别大的体系，不然一开始一般不用LBFGS，而是优先使用BFGS和CG；

Uus/pMeC6H4-/キ 发表于 2026-4-20 14:35
先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从40 ...

谢谢
不开smearing SCF收敛没问题
我试试开smearing ，看看结构是否正常

Uus/pMeC6H4-/キ 发表于 2026-4-20 14:35
先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从40 ...

谢谢！我先试试开smearing

先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从400 Ry提高），确定把IGNORE_CONVERGENCE_FAILURE给去掉也能正常收敛了再说几何优化的问题。