24小时热门版块排行榜    

查看: 2359  |  回复: 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的回帖

zhanghh1984

木虫 (正式写手)

我也找到了你说的那部分,我就是不明白在这个文件里这些数据又是什么意思?
谢谢啊!
12楼2008-12-13 09:22: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的回帖

yyx19840628

木虫 (著名写手)

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

你是分享吗?还以为是问问题
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
4楼2008-12-10 16:16:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +4 瞬息宇宙 2026-02-15 5/250 2026-02-15 22:28 by baiboxie
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见