用NWChem做SODFT在DFT计算中考虑旋轨耦合效应

beefly

冰释之川 发表于 2017-4-18 13:01
问题解决了，是Linear dependencies的锅：

Unfortunately, the spin-orbit code does not correctly h ...

nwchem判断基组线性依赖的默认阈值是1.0e-5。由于sodft不支持消除冗余积函数的算法，这么高的阈值显然不适用。可以用手册里没有讲的一个命令降低阈值。根据出错信息里的“Smallest S eigenvalue”，假如此值在1.0e-8以上，可以用

set lindep:tol 1.0e-8

但是如果出现更低的特征值，降低阈值的做法就不保险了，可能会带来其他问题，比如scf不收敛。