计算化学公社

标题: 求助,Multiwfn显示的HOMO-LUMO gap 值和CP2K输出文件里的gap值不一样 [打印本页]

作者
Author: 不会扣篮的后卫 时间: 2022-3-30 11:09
标题: 求助,Multiwfn显示的HOMO-LUMO gap 值和CP2K输出文件里的gap值不一样
本帖最后由不会扣篮的后卫于 2022-3-30 11:29 编辑

请教老师，同学。
我在用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 =1.2，在这个（p 轨道alpha电子数设为3，d 轨道beta电子数设为2）的O上设置 U = 1.

输出文件显示最后结构成功完成收敛,得到的两个HOMO-LUMO gap值如图所示。

(, 下载次数 Times of downloads: 14)
(, 下载次数 Times of downloads: 13)
我想通过Multiwfn来读取CP2K产生的molden.file来考察电子结构，
但是发现载入文件后，输入0，然后回车
得到的两个HOMO-LUMO gap值与先前的输出文件的HOMO-LUMO gap值不一样。
(, 下载次数 Times of downloads: 10)
请问各位老师，同学。出现这种情况是什么原因呢？
（附上输入，输出和molden.file文件）

作者
Author: sobereva 时间: 2022-3-30 19:34
你的molden文件不是最后一帧的。在你的输入文件里搜0.29241（Multiwfn显示的alpha-HOMO能量），就会发现这只是中途一个结构的HOMO能量

你的molden文件里没记录空轨道，Multiwfn直接把LUMO能量当成了0，所以就算你载入的molden是最后一帧的，Multiwfn给出的HOMO-LUMO gap和CP2K直接显示的也不同。想和输出文件里一致，必须让CP2K把空轨道写进去。

作者
Author: 不会扣篮的后卫 时间: 2022-3-31 06:07

sobereva 发表于 2022-3-30 19:34
你的molden文件不是最后一帧的。在你的输入文件里搜0.29241（Multiwfn显示的alpha-HOMO能量），就会发现这 ...

谢谢卢老师的建议。

卢老师，我根据您的建议和您之前在这个贴子http://bbs.keinsci.com/thread-27793-1-1.html里讲的如何在OT方法下加入空轨道。
在输入文件中，加了30个空轨道的参数设置
(, 下载次数 Times of downloads: 16)
输出文件显示结构成功收敛，得到两个HOMO-LUMO值，
(, 下载次数 Times of downloads: 10)
(, 下载次数 Times of downloads: 14)
但是当把最后一帧的molden.file(1_9)
(, 下载次数 Times of downloads: 13)
载入到Multiwfn，输入0，然后回车后
(, 下载次数 Times of downloads: 11)
得到的两个HOMO-LUMO gap值与先前的输出文件的HOMO-LUMO gap值不一样。

我发现Multiwfn读取出来的的的信息是一样的（不管加不加30个空轨道参数设置），
请问卢老师，是我在输入文件中，空轨道的参数设置不对
还是什么其他原因导致Multiwfn的HOMO-LUMO gap和CP2K跑出来的HOMO-LUMO gap不同？

（附上输入，输出和molden.file文件）

作者
Author: sobereva 时间: 2022-3-31 15:52

不会扣篮的后卫发表于 2022-3-31 06:07
谢谢卢老师的建议。

卢老师，我根据您的建议和您之前在这个贴子http://bbs.keinsci.com/thread-27793- ...

molden文件里把所有含有Occup=的行列出来，就明显知道molden文件里没空轨道
建议用对角化

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

sobereva 发表于 2022-3-31 15:52
molden文件里把所有含有Occup=的行列出来，就明显知道molden文件里没空轨道
建议用对角化

谢谢卢老师的建议，我去用对角化跑一跑。

请问卢老师，要是不追求以下两种方法的HOMO-LUMO gap值结果一致。（HOMO-LUMO gap 值来自于CP2K OT方法下的输出文件和HOMO-LUMO gap值来自于Multiwfn读取CP2K 产生的molden.file）。

直接用没有空轨道的molden.file 画出CuFeO2（S,P,D轨道）PDOS图。
请问这样画出来的结果合理吗？

(, 下载次数 Times of downloads: 14)

作者
Author: sobereva 时间: 2022-4-1 00:49

不会扣篮的后卫发表于 2022-3-31 22:15
谢谢卢老师的建议，我去用对角化跑一跑。

请问卢老师，要是不追求以下两种方法的HOMO-LUMO gap值结果 ...

如果你不关心空轨道，当前是合理的

作者
Author: 不会扣篮的后卫 时间: 2022-4-1 00:55

sobereva 发表于 2022-4-1 00:49
如果你不关心空轨道，当前是合理的

谢谢卢老师！

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