请问，使用CP2K计算轨道能量时，LUMO的能量小于HOMO的能量，这种现象正常吗？

junyuzhang

请问，使用CP2K计算轨道能量时，LUMO的能量小于HOMO的能量，这种现象正常吗？

我将 Gd Al Ga O 组成的体系中的两个Gd替换成了Ce，多了两个价电子，因此结果多了两个约5.3eV的轨道，但能量较高的一个是占据轨道，这是什么原因呢？


Orb:  2238 Ene(au/eV):     0.116622       3.1734 Occ: 2.000000 Type:A+B (    )
Orb:  2239 Ene(au/eV):     0.119056       3.2397 Occ: 2.000000 Type:A+B (    )
Orb:  2240 Ene(au/eV):     0.121170       3.2972 Occ: 2.000000 Type:A+B (    )
Orb:  2241 Ene(au/eV):     0.195917       5.3312 Occ: 2.000000 Type:A+B (    )
Orb:  2242 Ene(au/eV):     0.195831       5.3288 Occ: 0.000000 Type:A+B (    )
Orb:  2243 Ene(au/eV):     0.310951       8.4614 Occ: 0.000000 Type:A+B (    )
Orb:  2244 Ene(au/eV):     0.327857       8.9214 Occ: 0.000000 Type:A+B (    )
Orb:  2245 Ene(au/eV):     0.000000       0.0000 Occ: 0.000000 Type:A+B (    )



另外，我参考《详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件》里面关于轨道的记录里的方法2读取了空轨道的能量，那么如何才能通过命令行的方式将这些信息保存下来呢？

sobereva

junyuzhang 发表于 2024-10-26 21:38
老师，我计算了3个空轨道的能量，分别是2242，2243，2244号，但是2242空轨道的能量小于2241占据轨道，这 ...

大抵是OT在post-SCF阶段的求解精度问题，可以到邮件列表里问开发者
建议用对角化结合ADDED_MOS得到空轨道能量，这样肯定准确、没这个问题

junyuzhang

sobereva 发表于 2024-10-25 23:45
你没要求计算空轨道能量，所以空轨道能量自动都为0

详谈Multiwfn的命令行方式运行和批量运行的方法

老师，我计算了3个空轨道的能量，分别是2242，2243，2244号，但是2242空轨道的能量小于2241占据轨道，这是正常的还是哪里计算可能有问题？

sobereva

你没要求计算空轨道能量，所以空轨道能量自动都为0

详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612（http://bbs.keinsci.com/thread-24929-1-1.html）