计算化学公社

标题: HOMO-LUMO gap 的方法与基组选择求助 [打印本页]

作者
Author: lei234 时间: 2016-4-15 02:44
标题: HOMO-LUMO gap 的方法与基组选择求助
请问要算HOMO-LUMO gap一般选什么方法和基组？文献中比较推荐B3PW91，不知道怎么样？我的体系中有Ti，C，H 和O, Ti 用lanl2dz基组，C,H,O用6-31+g**可以吗？还有，很重要的一个问题，就是，必须用B3PW91优化后再用pop=full读轨道能量，还是用别的方法优化过，用B3PW91算个pop=full之后直接读轨道能量就可以？

谢谢~

作者
Author: sobereva 时间: 2016-4-15 03:14
分子体系用B3LYP就行。6-31G**即可，不用弥散。优化完了直接读gap就完了，用不着pop=full。

作者
Author: scf 时间: 2016-4-15 05:53
反正homo-lumo gap实验上也测不到，用哪个泛函好如何评估呢？和fundamental gap/optical gap比？

作者
Author: 冰释之川 时间: 2016-4-15 08:45
本帖最后由冰释之川于 2016-4-15 08:48 编辑

scf 发表于 2016-4-15 05:53
反正homo-lumo gap实验上也测不到，用哪个泛函好如何评估呢？和fundamental gap/optical gap比？

HOMO-LUMO gap可以近似的和 fundamental gap 比较，只能说是“近似”的； optical gap 是 fundamental gap 减去 elelctron-hole binding energy

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

作者
Author: sobereva 时间: 2016-4-15 13:09

scf 发表于 2016-4-15 05:53
反正homo-lumo gap实验上也测不到，用哪个泛函好如何评估呢？和fundamental gap/optical gap比？

对于导电寡聚物（很长的链），是有有物理意义的HOMO-LUMO gap的，B3LYP此时表现尤佳。虽然小分子HOMO-LUMO gap只是fundamental gap和optical gap的很粗略的近似，物理意义不明确，但说用B3LYP算这个合适并不为过。

作者
Author: lei234 时间: 2016-4-16 00:42

sobereva 发表于 2016-4-15 03:14
分子体系用B3LYP就行。6-31G**即可，不用弥散。优化完了直接读gap就完了，用不着pop=full。

请问直接读gap的意思是读图中红色框出的最后一次的LUMO-HOMO吗？
Alpha  occ. eigenvalues -- -0.38447  -0.35254  -0.29446  -0.28547
Alpha virt. eigenvalues -- 0.04077 0.04651 0.11846 0.14770 0.17393
那HOMO-LUMO gap=0.04077-（-0.28547） = 0.32624(a.u)=8.8ev
如果是开壳层体系，有两个单电子，非限制的，那α轨道和β轨道的HOMO和LUMO值是不是不一样，这样哪个算是gap呢？(以O2为例)
Alpha  occ. eigenvalues -- -0.67577  -0.64845  -0.39975  -0.39975
Alpha virt. eigenvalues -- 0.10632 0.10856 0.14784 0.14784 0.20805
α-HOMO-LUMO gap=0.10632- （-0.39975）=0.50607 a.u =13.7 ev
Beta  occ. eigenvalues -- -0.54281  -0.54281
Beta virt. eigenvalues -- -0.05841  -0.05841 0.10805 0.10830 0.15350 （为什么O2的β轨道LUMO是负的？）
β-HOMO-LUMO gap=-0.05841- （-0.54281）=0.60122 a.u=16.4 ev

谢谢~~

作者
Author: sobereva 时间: 2016-4-16 00:48
是。
非限制性开壳层的时候alpha gap和beta gap单独计算。

作者
Author: Scarlett-ww 时间: 2016-4-18 10:09

sobereva 发表于 2016-4-15 03:14
分子体系用B3LYP就行。6-31G**即可，不用弥散。优化完了直接读gap就完了，用不着pop=full。

Sob老师，如果分子体系不是用B3LYP优化的，gap计算和优化方法可以不用保持一致？

作者
Author: sobereva 时间: 2016-4-18 14:37

Scarlett-ww 发表于 2016-4-18 10:09
Sob老师，如果分子体系不是用B3LYP优化的，gap计算和优化方法可以不用保持一致？

可以不一致

作者
Author: Scarlett-ww 时间: 2016-4-19 11:21

sobereva 发表于 2016-4-18 14:37
可以不一致

恩好的，谢谢老师~

作者
Author: sukey 时间: 2019-10-30 21:17
请问energy计算之后如何查看HOMO-LUMO gap

作者
Author: sobereva 时间: 2019-10-31 06:31

sukey 发表于 2019-10-30 21:17
请问energy计算之后如何查看HOMO-LUMO gap

把fch文件载入Multiwfn，进入主功能0，文本窗口里直接给出gap值

这里提了
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

作者
Author: sukey 时间: 2019-11-1 10:08
请问别人给的example为什么用gaussview打不开，显示如图

作者
Author: sobereva 时间: 2019-11-2 08:47

sukey 发表于 2019-11-1 10:08
请问别人给的example为什么用gaussview打不开，显示如图

路径有中文

和帖子主题无关的问题不要在回帖里问，否则扰乱论坛帖子秩序，去置顶的简单量化答疑转帖里问

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