求助，在用CP2K跑结构优化任务时，为什么最后会得到两个HOMO-LUMO gap值？

不会扣篮的后卫

请教老师，同学。
我在用CP2K 跑12个Chemical Unit：CuFeO2的结构。 设置了以下输入文件参数。

• run_type:Oprimizingstructure（cell is fixed）
• 自旋多重度设为62的情况下
  

（在&BS下面把12个 +3 价Fe离子的d 轨道alpha 电子数设为5，d 轨道beta电子数设为0.
                   把1个 -2 价O离子的p 轨道alpha电子数设为3，d 轨道beta电子数设为2.
                   自旋多重度=12*5+1+1=62）

• charge = +1 (hole polaron---这个（p 轨道alpha电子数设为3，d 轨道beta电子数设为2）的O离子上移除一个电子)
• 使用DFT+U 方法，在12个Fe上设置U =2，在这个（p 轨道alpha电子数设为3，d 轨道beta电子数设为2）的O上设置 U = 2.5.
  

输出文件显示最后结构成功完成收敛。



但是发现出现了输出文件在收敛过程中，以及最后完成收敛的结果中，一直出现两个HOMO-LUMO gap 的值。

请问各位老师，同学，为什么会出现这种情况？（附上输入和输出文件）

highlight

你这问题和 http://bbs.keinsci.com/thread-28439-1-1.html 简直一毛一样
显然是开壳层下两种自旋各自对应的 HOMO-LUMO gap
输出文件前面写的就是不同自旋的占据非占据本征值，自己减一下就能算出来

不会扣篮的后卫

highlight 发表于 2022-3-23 11:39
你这问题和 http://bbs.keinsci.com/thread-28439-1-1.html 简直一毛一样
显然是开壳层下两种自旋各自对应 ...

老师您好，我发现我的输出文件里只有最低非占据轨道本征值。

请问老师，在输入文件里，输入什么关键词，可以在输出文件里得到最高占据轨道本征值？

highlight

不会扣篮的后卫 发表于 2022-3-24 00:07
老师您好，我发现我的输出文件里只有最低非占据轨道本征值。

请问老师，在输入文件里，输入什么关键词 ...

最高占据轨道不就是占据里本征值最大的那个么，你截图部分的前面就有了
或者你直接用 Fermi Energy，但要注意单位换算

不会扣篮的后卫

highlight 发表于 2022-3-24 08:38
最高占据轨道不就是占据里本征值最大的那个么，你截图部分的前面就有了
或者你直接用 Fermi Energy，但 ...

谢谢老师！

再请教老师一个问题，因为实验值CuFeO2的bandgap是1.5eV, 我想通过调整不同的DFT+U的U值来使结构成功收敛后，输出文件最后的HOMO-LUMO gap的值接近1.5eV.



请问老师，由于我最后有两个HOMO-LUMO gap 的值，
所以我该选择哪一个HOMO-LUMO gap呢？

highlight

不会扣篮的后卫 发表于 2022-3-24 11:51
谢谢老师！

再请教老师一个问题，因为实验值CuFeO2的bandgap是1.5eV, 我想通过调整不同的DFT+U的U值来 ...

http://bbs.keinsci.com/thread-28439-1-1.html

如果你是指非限制性计算给了两套轨道能量，应当取两种自旋中最高的占据的那个

不会扣篮的后卫

highlight 发表于 2022-3-24 13:42
http://bbs.keinsci.com/thread-28439-1-1.html

谢谢老师！