用Multiwfn分析CP2K产生的molden文件，HOMO和LUMO能量显示都为0的问题求助

ljc050512 · 发表于 Post on 2022-7-17 18:30:27

诸位老师，我计算的是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，这个能量跨度这么大正常吗?还是我哪里设置的不对导致算错了？

sobereva · 发表于 Post on 2022-7-17 18:57:56

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

ljc050512 · 发表于 Post on 2022-7-17 20:41:35

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

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

ljc050512 · 发表于 Post on 2022-7-17 21:45:35

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

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

sobereva · 发表于 Post on 2022-7-18 10:15:59

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

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

sobereva · 发表于 Post on 2022-7-18 10:24:16

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

别用ADMM和OT了

ljc050512 · 发表于 Post on 2022-7-18 11:04:04

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

再次感谢老师的解答！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Multiwfn使用咨询] 用Multiwfn分析CP2K产生的molden文件，HOMO和LUMO能量显示都为0的问题求助

浏览过的版块