24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2425  |  回复: 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的回帖

yuhuobuku

木虫 (正式写手)

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

yyx19840628

木虫 (著名写手)

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

你是分享吗?还以为是问问题
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
4楼2008-12-10 16:16:23
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +10 张zz111 2026-03-27 11/550 2026-03-30 09:17 by 无际的草原
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-26 6/300 2026-03-30 09:10 by 123xxb
[考研] 材料工程专硕求调剂 +3 hyl3153942 2026-03-29 3/150 2026-03-30 07:41 by melodiousnow
[考研] 求调剂 +6 家佳佳佳佳佳 2026-03-29 6/300 2026-03-30 06:48 by zhyzzh
[考研] 南京大学化学调剂 +9 景随风 2026-03-29 14/700 2026-03-29 21:22 by chixmc
[考研] 070305高分子化学与物理 304分求调剂 +12 c297914 2026-03-28 12/600 2026-03-29 16:21 by Serene1974
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 071000生物学求调剂,初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +7 Micky11223 2026-03-25 7/350 2026-03-28 18:34 by 无际的草原
[考研] 299求调剂 +7 嗯嗯嗯嗯2 2026-03-27 7/350 2026-03-28 13:09 by 唐沐儿
[考研] 070300求调剂306分 +4 26要上岸 2026-03-27 4/200 2026-03-28 13:06 by 唐沐儿
[考研] 339求调剂,想调回江苏 +6 烤麦芽 2026-03-27 8/400 2026-03-28 10:40 by 烤麦芽
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 316求调剂 +5 Pigcasso 2026-03-24 5/250 2026-03-27 12:10 by zhshch
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 324求调剂 +8 hanamiko 2026-03-26 10/500 2026-03-27 08:06 by hypershenger
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
信息提示
请填处理意见