计算化学公社

标题: 求助，FERWE 和FERDO的设置 [打印本页]

作者
Author: guyizhang 时间: 15 hour ago
标题: 求助，FERWE 和FERDO的设置
我现在想计算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 能量接近吗，这种情况应该怎么做？谢谢！

作者
Author: guyizhang 时间: 7 hour ago
本帖最后由 guyizhang 于 2026-4-27 16:05 编辑

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.

欢迎光临计算化学公社 (http://bbs.keinsci.com/)