今天遇到 Materials Studio cij 计算失败,提示基态电子最小化失败。检查输出文件,发现蛛丝马迹。

错误信息

1
2
3
4
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
Current trace stack:
 check_elec_ground_state
 castep

检查castep 文件发现下面的片段

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
.
.
.
     99  -2.09594982E+004  6.31070540E-001   3.18129596E-006    4013.01  <-- SCF
    100  -2.09594980E+004  6.31017470E-001  -2.41328571E-006    4053.95  <-- SCF
------------------------------------------------------------------------ <-- SCF
 *Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy, E     =  -20959.33865776     eV
Current free energy (E-TS)  =  -20959.49797103     eV
(energies not corrected for finite basis set)
 .
 .
 .

意思大概是我设置的 SCF 迭代次数不够,100次后仍然未能找到自恰解。

回想使用步骤皆是以前的所用的,马上翻看之前成功的castep文件,

1
2
3
4
5
6
7
8
     13  -2.78030685E+004 -1.99686373E+000  -7.49554981E-007     898.41  <-- SCF
     14  -2.78030685E+004 -1.99941165E+000   4.13854537E-008     943.16  <-- SCF
     15  -2.78030685E+004 -2.00267847E+000  -2.65786251E-007     988.60  <-- SCF
------------------------------------------------------------------------ <-- SCF
 
Final energy, E             =  -27803.06847596     eV
Final free energy (E-TS)    =  -27803.06847596     eV
(energies not corrected for finite basis set)

几乎都是十几次得到自恰解。

找到问题所在。

迭代次数我不是设置的500 么

回到 windows 版的MS,检查设置。我勒个去,Geometry OptimizationMax.iterations是500没错,但是Elastic ConstantsMax.iterations还是默认的100。

还有 SCF 迭代次数可以设置。

问题有可能就在这里,提高迭代次数上线,有可能能够收敛,也有可能不能,白白浪费时间。

截断能、k点、自旋、对称性等都有可能影响收敛。

没办法,只能试试看了。

后续如果解决,会详细记录。