计算化学公社

标题: 用Multiwfn分析CP2K产生的molden文件，HOMO和LUMO能量显示都为0的问题求助 [打印本页]

作者
Author: ljc050512 时间: 2022-7-17 18:30
标题: 用Multiwfn分析CP2K产生的molden文件，HOMO和LUMO能量显示都为0的问题求助
诸位老师，我计算的是K参杂的C3N4，想在输出文件中读取体系 HOMO-LUMO的gap。设置自旋多重度为3的UKS，得到输出文件显示如下结果： Eigenvalues of the occupied subspace spin          1
---------------------------------------------
   -1.24062624    -1.23961847    -1.02236521    -1.01987027
...
   -0.13488807    -0.13206934    -0.12865901    -0.12627031
   -0.04171729    -0.03571244
Fermi Energy [eV] : -0.971785

Eigenvalues of the occupied subspace spin          2
---------------------------------------------
   -1.24057191    -1.23956607    -1.02144788    -1.01897067
...
   -0.13543226    -0.13362275    -0.12984849    -0.12783345
Fermi Energy [eV] : -3.478525

  Lowest Eigenvalues of the unoccupied subspace spin          1
-----------------------------------------------------
OT| Eigensolver reached convergence in 243 iterations
   -0.02901051

  Lowest Eigenvalues of the unoccupied subspace spin          2
-----------------------------------------------------
OT| Eigensolver reached convergence in 159 iterations
   -0.02508385

HOMO - LUMO gap [eV] : 0.182369
HOMO - LUMO gap [eV] : 2.795959

我在用Multiwfn分析.molden时候，却发现给出的HOMO和LUMO能量都是0（如下图）

为什么会出现这个情况呢？
应该选用哪个HOMO - LUMO gap和费米能级数值呢？
还有
在  Eigenvalues of the occupied subspace spin          1
---------------------------------------------
   -1.24062624    -1.23961847    -1.02236521    -1.01987027
...
   -0.13488807    -0.13206934    -0.12865901    -0.12627031
   -0.04171729    -0.03571244
这里，-0.12627031下一个能量直接跳到了-0.04171729，这个能量跨度这么大正常吗?还是我哪里设置的不对导致算错了？

作者
Author: sobereva 时间: 2022-7-17 18:57
用OT时不计算轨道能量，用对角化
用对角化时默认也不计算空轨道。如果想输出HOMO、LUMO能级以及gap，必须在&DFT里加上
&PRINT
   &MO_CUBES
      STRIDE 1 #Stride of exported cube file
      NHOMO    1
      NLUMO    1
   &END MO_CUBES
&END PRINT

作者
Author: ljc050512 时间: 2022-7-17 20:41

sobereva 发表于 2022-7-17 18:57
用OT时不计算轨道能量，用对角化
用对角化时默认也不计算空轨道。如果想输出HOMO、LUMO能级以及gap，必须 ...

非常感谢老师的解答！那么就是说用HSE06_ADMM算HOMO - LUMO是不合适的，对吗？

作者
Author: ljc050512 时间: 2022-7-17 21:45

sobereva 发表于 2022-7-17 18:57
用OT时不计算轨道能量，用对角化
用对角化时默认也不计算空轨道。如果想输出HOMO、LUMO能级以及gap，必须 ...

还想请教老师一个问题，在Multiwfn里面，CP2K部分列出的泛函里面，哪一个纯泛函比较合适计算轨道能量或者HOMO - LUMO？杂化泛函太贵了，先请您推荐个相对好一些的纯泛函。

作者
Author: sobereva 时间: 2022-7-18 10:15

ljc050512 发表于 2022-7-17 21:45
还想请教老师一个问题，在Multiwfn里面，CP2K部分列出的泛函里面，哪一个纯泛函比较合适计算轨道能量或者 ...

哪个都不理想，都半斤八两，普遍低估gap

作者
Author: sobereva 时间: 2022-7-18 10:24

ljc050512 发表于 2022-7-17 20:41
非常感谢老师的解答！那么就是说用HSE06_ADMM算HOMO - LUMO是不合适的，对吗？

别用ADMM和OT了

作者
Author: ljc050512 时间: 2022-7-18 11:04

sobereva 发表于 2022-7-18 10:24
别用ADMM和OT了

再次感谢老师的解答！

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