以前做过dft或hf下的势能面扫描,一般都能成功。现在在尝试casscf下势能面扫描,但没成功,不知道为什么,来这里咨询下高手。
我的体系是一个Oh构型的结构CrCl6_-3,稳定结构下Cr-Cl键长为2.4674A。我用hf计算获得了该结构下的电子结构,并将其三个d轨道B1g(dxy)、B2g(dyz)、B3g(dxz)以及对应dz2和d_x2-y2两个空轨道放入活性空间进行单点计算,即cas(3,5),成功了。
现在我想利用do循环做一下其基态能量E(4Ag)随bond(Cr-Cl)的变化曲线,计算时用restart命令读取hf产生的轨道信息作为初猜。但是在结构偏离稳定结构0.1A时又出现了如下错误,而且终止。
? Warning
? Non zero Pspace value in the residuum
我不知道这个问题该如何解决,也不知道文献中的势能面扫描是不是这样做的?还是在每个结构下单独计算获得其各个组态的能量?(如果是后者,的确有点麻烦)
所以想请教一下高手们,cas的势能面扫描应该如何做?