今天遇到 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 Optimization 的 Max.iterations是500没错，但是Elastic Constants 的 Max.iterations还是默认的100。

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

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

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

没办法，只能试试看了。

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