24小时热门版块排行榜    

查看: 2323  |  回复: 12
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yuhuobuku

木虫 (正式写手)

[交流] 【交流】回zhanghh1984关于HOMO、LUMO问题

发现我回不了贴了,重新开个贴吧。高斯程序计算出来的单位都是原子单位a.u.
然后自己换算吧。关于如何阅读HOMO LUMO能量的问题我以前写了个小的总结给你发过来看一下:(仅供参考)
【一】

一般选择用pop=reg来计算HOMO和LUMO轨道,分别打印出5条。如果用非限制性方法计算开窍层体系,体系有一个单电子。那么分别两套轨道,计算结果alpha,beta电子各给出了一套值。部分输出文件为:

Orbital energies and kinetic energies (alpha):
1 2
61 (A)--O -0.41907 1.15066
62 (A)--O -0.41445 1.17615
63 (A)--O -0.40228 1.46886
64 (A)--O -0.36519 1.82895
65 (A)--O -0.36450 1.83459
66 (A)--O -0.27919 1.58649    HOMO
67 (A)--V -0.16904 1.68438
68 (A)--V -0.15822 1.33172
69 (A)--V -0.13931 1.45057
70 (A)--V -0.09413 1.52481
71 (A)--V -0.09042 1.46055
Orbital energies and kinetic energies (beta):
1 2
61 (A)--O -0.40825 1.14433
62 (A)--O -0.40405 1.17127
63 (A)--O -0.38829 1.43736
64 (A)--O -0.36438 1.82856
65 (A)--O -0.36298 1.82272
66 (A)--V -0.22356 1.56890    LUMO
67 (A)--V -0.15432 1.31964
68 (A)--V -0.14677 1.69606
69 (A)--V -0.13491 1.44071
70 (A)--V -0.09264 1.49395
71 (A)--V -0.08881 1.42911

若使用限制性方法,对闭壳层体系进行计算,得出一套轨道,只有alpha电子的输出。电子全部配对。部分输出文件为:

Orbital energies and kinetic energies (alpha):
1 2
61 (A)--O -0.41191 1.14750
62 (A)--O -0.40757 1.17374
63 (A)--O -0.39302 1.45048
64 (A)--O -0.36450 1.82922
65 (A)--O -0.36345 1.82981
66 (A)--O -0.16767 1.58047     HOMO
67 (A)--V -0.16164 1.67120     LUMO
68 (A)--V -0.15696 1.34890
69 (A)--V -0.13813 1.44738
70 (A)--V -0.09351 1.51515
71 (A)--V -0.08980 1.44987


HOMO和LUMO只看占据电子的最高能量轨道和未占据电子的最低量轨道, alpha还是beta 没什么关系。

前一个:HOMO 66 (A)--O -0.27919 1.58649 (alpha)
LUMO 66 (A)--V -0.22356 1.56890 (Beta)

后一个:HOMO 66 (A)--O -0.16767 1.58047 (Alpha, beta相同,占据一个电子)
         LUMO 67 (A)--V -0.16164 1.67120  

【二】

Orbital Symmetries:
Occupied (A') (A') (A') (A') (A" (A') (A') (A"
Virtual (A') (A') (A" (A')
The electronic state is 1-A'.
Alpha occ. eigenvalues -- -20.29970 -11.10195 -1.28226 -0.81445 -0.61834
Alpha occ. eigenvalues -- -0.51763 -0.39445 -0.36126
Alpha virt. eigenvalues -- 0.23585 0.65804 0.75009 0.77751           

如上面这个输出例子(G98测试148,分子为H2CO)怎样找出HOMO和LUMO轨道?

这个计算结果存在一个问题,对称性没有判断正确,H2CO应是C2V对称性的分子。HOMO就是最高占据轨道,第一行最后一个,LUMO是最低未占据轨道,第二行第一个。
Occupied (A') (A') (A') (A') (A" (A') (A') (A" Virtual (A') (A') (A" (A') HOMO就是最高占据轨道, occupied就是占据, 当然应该是最后一个occupied就是。
LUMO是最低非占据轨道,virtual就是虚轨道,所以就是第一个virtual就是。

下面例子是个C2对称性的体系:

Alpha Molecular Orbital Coefficients:
-0.38948(773轨道--Occ.) -0.38948(774轨道--Occ.) -0.38510(775轨道--Occ.)

-0.38509(776轨道--Occ.) -0.23837(777轨道--Virt.)

Beta Molecular Orbital Coefficients:
-0.38418(773轨道--Occ.) -0.38418(774轨道--Occ.) -0.34763(775轨道--Virt.)
HOMO : -0.38418(774轨道--Occ.)
LUMO : -0.34763(775轨道--Virt.)

刚开始看这个例子的时候不是太懂,容易跟上面的弄混。但如果把这2套轨道竖立排列起来就容易懂了。

Alpha Molecular Orbital Coefficients:      Beta Molecular Orbital Coefficients:
-0.38948(773轨道--Occ.)                      -0.38418(773轨道--Occ.)

-0.38948(774轨道--Occ.)                     -0.38418(774轨道--Occ.)

-0.38510(775轨道--Occ.)                      -0.34763(775轨道--Virt.)

-0.38509(776轨道--Occ.)

-0.23837(777轨道--Virt.)

用POP=reg可以得到HOMO能量值的,你找到Orbital energies and kinetic energies (alpha)这一行,下面列出了五个最高占据轨道和五个最低未占据轨道的能量值。
至于电量,Gaussian默认输出的是Mulliken charges,你搜索它肯定能看到每个原子带的电量的!
5个轨道中,前面带o的为HOMO,v的为LUMO,O和V靠在一起.

[ Last edited by yuhuobuku on 2008-12-10 at 15:43 ]
回复此楼
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...

你是分享吗?还以为是问问题
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
4楼2008-12-10 16:16:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

erylingjet

铁杆木虫 (著名写手)

不知道最后一个什么意思?楼主再解释一下吧?

Alpha Molecular Orbital Coefficients:      Beta Molecular Orbital Coefficients:
-0.38948(773轨道--Occ.)                      -0.38418(773轨道--Occ.)

-0.38948(774轨道--Occ.)                     -0.38418(774轨道--Occ.)

-0.38510(775轨道--Occ.)                      -0.34763(775轨道--Virt.)

-0.38509(776轨道--Occ.)

-0.23837(777轨道--Virt.)

为什么不是776和777?
2楼2008-12-10 15:21:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuobuku

木虫 (正式写手)

★ ★
lei0736(金币+2,VIP+0):鼓励原创
因为占据轨道里774比776大啊 要找的就是最高占据轨道呀
同理 775比777低
我个人就是这么理解的。
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
3楼2008-12-10 15:53:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuobuku

木虫 (正式写手)

引用回帖:
Originally posted by yyx19840628 at 2008-12-10 16:16:
你是分享吗?还以为是问问题

刚刚回不了别的楼主的帖子了,我只好重新开了个贴
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
5楼2008-12-10 16:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见