应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 计算orbital composition
111/2返回列表12下一页
查看: 1885  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

meteoric30

木虫 (正式写手)

[求助] 计算orbital composition

我简单计算了甲醇的单点能 基组6-31g** pop=reg得到
Molecular Orbital Coefficients
                           5         6         7         8         9
                        (A)--O    (A)--O    (A)--O    (A)--O    (A)--O
     EIGENVALUES --    -0.50123  -0.43327  -0.41890  -0.32941  -0.26492
   1 1   C  1S          0.02599   0.00000  -0.00543   0.01363   0.00000
   2        2S         -0.04962   0.00001   0.01119  -0.03358   0.00000
   3        2PX        -0.23161   0.00004  -0.27326   0.11697   0.00000
   4        2PY         0.21363   0.00003  -0.30666  -0.26989   0.00001
   5        2PZ         0.00001   0.42015   0.00004   0.00000  -0.17876
   6        3S         -0.05786   0.00002   0.02227  -0.03258   0.00000
   7        3PX        -0.07912   0.00002  -0.11207   0.03513   0.00000
   8        3PY         0.09332   0.00000  -0.15399  -0.07930   0.00000
   9        3PZ        -0.00001   0.17326   0.00003   0.00000  -0.02244
  10        4XX         0.02178   0.00000   0.01845  -0.00832   0.00000
  11        4YY        -0.00188  -0.00002  -0.02006  -0.00488   0.00002
  12        4ZZ        -0.01434   0.00002   0.00052   0.01518  -0.00002
  13        4XY        -0.00986   0.00000  -0.00836  -0.02462   0.00000
  14        4XZ         0.00000   0.00249   0.00000   0.00000  -0.03162
  15        4YZ        -0.00001  -0.01587   0.00002   0.00002   0.01328
  16 2   H  1S          0.02750  -0.00013  -0.23672  -0.17711   0.00011
  17        2S          0.01577  -0.00009  -0.18684  -0.18524   0.00015
  18        3PX        -0.00412   0.00000   0.00113   0.00336   0.00000
  19        3PY         0.00185   0.00001   0.00617   0.00333  -0.00001
  20        3PZ         0.00000   0.00779  -0.00001  -0.00001  -0.00486
  21 3   H  1S         -0.11948  -0.20300   0.03706   0.11400   0.14852
  22        2S         -0.07696  -0.17200   0.02421   0.13849   0.19822
  23        3PX        -0.00066   0.00461  -0.00480   0.00065   0.00034
  24        3PY         0.00006  -0.00539  -0.00425  -0.00360   0.00311
  25        3PZ        -0.00584  -0.00226   0.00174   0.00534   0.00115
  26 4   H  1S         -0.11936   0.20314   0.03689   0.11378  -0.14863
  27        2S         -0.07687   0.17210   0.02406   0.13825  -0.19837
  28        3PX        -0.00066  -0.00461  -0.00480   0.00065  -0.00034
  29        3PY         0.00007   0.00539  -0.00425  -0.00361  -0.00311
  30        3PZ         0.00584  -0.00227  -0.00174  -0.00533   0.00116
  31 5   O  1S         -0.03947   0.00000   0.03131  -0.06665   0.00000
  32        2S          0.09216   0.00000  -0.05801   0.12995   0.00000
  33        2PX         0.28346  -0.00003   0.34108  -0.19977   0.00000
  34        2PY         0.31776   0.00001  -0.12932   0.42084   0.00000
  35        2PZ         0.00000   0.25642   0.00003   0.00001   0.59486
  36        3S          0.14668   0.00001  -0.14062   0.32712  -0.00001
  37        3PX         0.12762  -0.00002   0.19338  -0.14321   0.00000
  38        3PY         0.16321   0.00001  -0.08316   0.28146   0.00000
  39        3PZ         0.00000   0.16208   0.00002   0.00000   0.45628
  40        4XX         0.01207   0.00000   0.00965  -0.00899   0.00000
  41        4YY        -0.02116   0.00000   0.00558  -0.03088   0.00000
  42        4ZZ         0.00150   0.00000   0.00106  -0.00012   0.00000
  43        4XY         0.00010   0.00000  -0.02405   0.00764   0.00000
  44        4XZ         0.00000   0.01420   0.00000   0.00000   0.00830
  45        4YZ         0.00000  -0.01107   0.00000   0.00000  -0.02046
  46 6   H  1S         -0.21191   0.00001  -0.04498  -0.15220   0.00000
  47        2S         -0.13754   0.00000  -0.04682  -0.12187   0.00001
  48        3PX        -0.00194   0.00000   0.01073  -0.01014   0.00000
  49        3PY        -0.01180   0.00000  -0.00761   0.00213   0.00000
  50        3PZ         0.00000   0.00855   0.00000   0.00000   0.01974
可是如何才能得到原子或者基团在HOMO LUMO中所占百分比？
如看文献中
Ru   HOMO 25.3%   NSC    32.8%
        LOMO 6.96%              0.727%
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-24 17:29:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)

【答案】应助回帖

meteoric30(金币+10): 能详细些么 2011-05-26 11:30:16
pop=nbo

111111111111111111111111111111111111111111111111111111111111111111111111
赞一下 回复此楼
更痛苦的是梦没醒路已经不能走了
2楼2011-05-26 05:13:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

【答案】应助回帖

meteoric30(金币+5): 3q a lot 2011-05-26 11:20:40
meteoric30(金币+5): 2011-05-26 16:13:52
引用回帖:
Originally posted by meteoric30 at 2011-05-24 17:29:15:
我简单计算了甲醇的单点能 基组6-31g** pop=reg得到
Molecular Orbital Coefficients
                           5         6         7         8         9
                        (A)--O    (A)--O    ( ...

http://web.rfeda.net/dp-bbsthread-793103.html
以供参考
赞一下 回复此楼
3楼2011-05-26 10:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
引用回帖:
Originally posted by mchen10 at 2011-05-26 05:13:33:
pop=nbo

111111111111111111111111111111111111111111111111111111111111111111111111

我用甲烷分子简略的做了一下HOMO分析：它的HOMO是第５个轨道 Molecular Orbital Coefficients 1 2 3 4 5 (A1)--O (A1)--O (T2)--O (T2)--O (T2)--O EIGENVALUES -- -10.14326-0.70258-0.39844-0.39844-0.39844 1 1 C1S 0.56289-0.10958 0.00000 0.00000 0.00000 2 2S 0.46308-0.16939 0.00000 0.00000 0.00000 3 2PX 0.00000 0.00000 0.20725 0.00000 0.00000 4 2PY 0.00000 0.00000 0.00000 0.00000 0.20725 5 2PZ 0.00000 0.00000 0.00000 0.20725 0.00000 6 3S 0.02002 0.43848 0.00000 0.00000 0.00000 7 3PX 0.00000 0.00000 0.29256 0.00000 0.00000 8 3PY 0.00000 0.00000 0.00000 0.00000 0.29256 9 3PZ 0.00000 0.00000 0.00000 0.29256 0.0000010 4S -0.00309 0.25481 0.00000 0.00000 0.00000（后面的就不贴出来了）　　我把第HOMO（第５列）的所有数平方后求和，然后分辨把1C，2H，3H，4H，5H的轨道的平方和除以HOMO（第5列）的所有系数的平方和，从而得到各原子在甲烷分子的HOMO中所占比例： 1C占了0.527013397 2H占了0.118246651 3H占了0.118246651 4H占了0.118246651 5H占了0.118246651请教大侠们，具体是不是这样算的？

大侠， 这样对吗？
我做了 pop=nbo，但是输出文件中如何才能得到homo中每个原子的所占的比例呢？ 请大侠指明。
赞一下 回复此楼
4楼2011-05-26 11:29:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
引用回帖:
Originally posted by kaegi at 2011-05-26 10:54:31:
http://web.rfeda.net/dp-bbsthread-793103.html
以供参考

我觉得 他的平方和不等于一 是因为有些轨道所占比例太小.因此被省略 导致 但是我计算了甲醇的 才0.729，  如果用部分除以不等于1的那一项（0.729） 总感觉不是太对   更倾向于直接平方和所得的数。
赞一下 回复此楼
5楼2011-05-26 11:33:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)
引用回帖:
Originally posted by meteoric30 at 2011-05-26 11:33:51:
我觉得 他的平方和不等于一 是因为有些轨道所占比例太小.因此被省略 导致 但是我计算了甲醇的 才0.729，  如果用部分除以不等于1的那一项（0.729） 总感觉不是太对   更倾向于直接平方和所得的数。

我猜是因为所列出的轨道并不是所有轨道，楼主如果用nbo=full看看会怎么样。不过这个其实我也不会，期待有人可以研究出来。
赞一下 回复此楼
6楼2011-05-26 14:14:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
引用回帖:
Originally posted by kaegi at 2011-05-26 14:14:47:
我猜是因为所列出的轨道并不是所有轨道，楼主如果用nbo=full看看会怎么样。不过这个其实我也不会，期待有人可以研究出来。

full的话 输出的轨道比reg多  nbo 我实在是不懂如何寻找前线轨道组成。
赞一下 回复此楼
7楼2011-05-26 15:15:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)
http://muchong.com/html/200811/1078422.html
这个不知道能否符合你的要求
赞一下 回复此楼
8楼2011-05-26 15:28:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
引用回帖:
Originally posted by kaegi at 2011-05-26 15:28:02:
http://muchong.com/html/200811/1078422.html
这个不知道能否符合你的要求

nice 原理懂了就ok 其实没必要用软件  excel 足矣  多谢你啦
赞一下 回复此楼
9楼2011-05-26 16:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
引用回帖:
Originally posted by kaegi at 2011-05-26 15:28:02:
http://muchong.com/html/200811/1078422.html
这个不知道能否符合你的要求

nice 原理懂了就ok 其实没必要用软件  excel 足矣  多谢你啦
赞一下 回复此楼
10楼2011-05-26 16:10:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meteoric30 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291 求调剂 +4 化工2026届毕业� 2026-03-21 4/200 2026-03-23 09:59 by bingchuan
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123！
[考研] 277材料科学与工程080500求调剂 +7 自由煎饼果子 2026-03-16 7/350 2026-03-22 22:40 by ACS Nano——
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 289材料与化工（085600）B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
关闭
请填处理意见
关闭