计算化学公社

标题: 求助，在用CP2K跑结构优化任务时，为什么最后会得到两个HOMO-LUMO gap值？ [打印本页]

作者
Author: 不会扣篮的后卫 时间: 2022-3-23 11:25
标题: 求助，在用CP2K跑结构优化任务时，为什么最后会得到两个HOMO-LUMO gap值？
本帖最后由不会扣篮的后卫于 2022-3-23 11:26 编辑

请教老师，同学。
我在用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.

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

(, 下载次数 Times of downloads: 16)
但是发现出现了输出文件在收敛过程中，以及最后完成收敛的结果中，一直出现两个HOMO-LUMO gap 的值。
(, 下载次数 Times of downloads: 19)
请问各位老师，同学，为什么会出现这种情况？（附上输入和输出文件）

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

作者
Author: 不会扣篮的后卫 时间: 2022-3-24 00:07
本帖最后由不会扣篮的后卫于 2022-3-24 05:28 编辑

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

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

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

作者
Author: highlight 时间: 2022-3-24 08:38

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

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

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

作者
Author: 不会扣篮的后卫 时间: 2022-3-24 11:51

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呢？

作者
Author: highlight 时间: 2022-3-24 13:42

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

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

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

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

作者
Author: 不会扣篮的后卫 时间: 2022-3-24 22:23

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

谢谢老师！

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