| 查看: 2679 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 raman |
» 猜你喜欢
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有183人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
yyx19840628
木虫 (著名写手)
鹰的眼睛，狼的耳朵，豹的速 ...
- 应助: 4 (幼儿园)
- 贵宾: 0.41
- 金币: 3521.9
- 散金: 50
- 红花: 8
- 帖子: 2321
- 在线: 129小时
- 虫号: 508615
- 注册: 2008-02-21
- 性别: MM
- 专业: 环境污染化学
4楼2008-12-10 16:16:23
erylingjet
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 5.372
- 金币: 9933.4
- 散金: 752
- 红花: 2
- 帖子: 2248
- 在线: 241.5小时
- 虫号: 415302
- 注册: 2007-06-28
- 性别: GG
- 专业: 计算化学
2楼2008-12-10 15:21:44
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
3楼2008-12-10 15:53:54
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
5楼2008-12-10 16:43:04
