24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2432  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程 085701,267求调剂 +7 minht 2026-03-29 7/350 2026-03-30 20:46 by dophin1985
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +9 XX小邓 2026-03-29 9/450 2026-03-30 17:18 by limeifeng
[考研] 307求调剂 +16 超级伊昂大王 2026-03-24 17/850 2026-03-30 17:15 by 544594351
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 070300一志愿211,312分求调剂院校 +4 小黄鸭宝 2026-03-30 4/200 2026-03-30 11:49 by 棒棒球手
[考研] 化学工程085602 305分求调剂 +25 RichLi_ 2026-03-25 25/1250 2026-03-30 11:02 by klasasda
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +6 @taotao 2026-03-26 7/350 2026-03-30 10:43 by 我是小康
[考研] 0856求调剂 +8 楒桉 2026-03-28 8/400 2026-03-30 10:00 by wzy-lxz
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 压国家一区线,求导师收留,有恩必谢! +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 求调剂推荐 材料 304 +15 荷包蛋hyj 2026-03-26 15/750 2026-03-28 04:13 by fmesaito
[考研] 321求调剂 +6 Ymlll 2026-03-24 6/300 2026-03-26 20:50 by 不吃魚的貓
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
信息提示
请填处理意见