求助，FERWE 和FERDO的设置

guyizhang

我现在想计算dimer的电子转移重组能，具体的反应是 A(4+)---(anion)---(A5+)     ->      A(5+)---(anion)---(A4+)
我在用vasp计算，我把A(4+)---(anion)---(A5+) 结构优化后，想通过设置FERWE和FERDO计算A(4+)---(anion)---(A5+) 这个几何结构下电子态为A(5+)---(anion)---(A4+)的单点能。


   87       -6.141256     -6.142310   1.000000   1.000000
   88       -5.605999(delocalized O)     -5.605999(delocalized O)   1.000000   1.000000
   89       -5.094215(34号原子,4价)     -4.049428 (35,5价)  1.000000   0.000000
   90       -4.049797(35号原子,5价)     -3.184751(34)   0.000000   0.000000
   91       -3.377757     -3.054222   0.000000   0.000000
   92       -3.365011     -2.915274   0.000000   0.00000
上面我vasp 计算结构优化后的eigenval 文件和其中原子的占据情况，现在我想单点计算后让34号5价，35号4价，我通过设置FERWE交换了spin up 的band 89 和90 的占据情况， 但是并没有得到我想要的结果，最后的电子态还是A(4+)---(anion)---(A5+) ， 请教这种情况是什么原因，是因为spin up的band 90和spin down 的band 89 能量接近吗，这种情况应该怎么做？谢谢！

guyizhang

LDIAG = .FALSE
ALGO=all

在INCAR里添加上面的关键词，就可以解决了


For ALGO = Damped and All, the final subspace diagonalization is simply skipped if LDIAG = .FALSE. is set. Generally using LDIAG = .FALSE. is only advised, if one wants to maintain a certain orbital order, for instance, when reading the orbitals from an existing WAVECAR file.