24小时热门版块排行榜    

查看: 1252  |  回复: 9
当前主题已经存档。

zhaohongjie

[交流] [求助]如何读取Gaussian输出文件中轨道系数(新手)

[求助]如何从下面的Gaussian输出文件中找出轨道系数及轨道能!!(新手多谢),请帮忙标出来!!!
The electronic state is 1-A1.
Alpha  occ. eigenvalues --  -20.58265 -11.33946  -1.39265  -0.87259  -0.69715
Alpha  occ. eigenvalues --   -0.63950  -0.52294  -0.44073
Alpha virt. eigenvalues --    0.13573   0.24842   0.33338   0.37329   0.73660
Alpha virt. eigenvalues --    0.80783   0.84685   0.94689   1.10445   1.10700
Alpha virt. eigenvalues --    1.13937   1.27145   1.33529   1.62050   1.78192
Alpha virt. eigenvalues --    1.79416   1.99239   2.18347   2.23684   2.45514
Alpha virt. eigenvalues --    2.64513   2.87165   2.97616   3.27576   4.09792
Alpha virt. eigenvalues --    4.47637
     Molecular Orbital Coefficients
                           1         2         3         4         5
                        (A1)--O   (A1)--O   (A1)--O   (A1)--O   (B2)--O
EIGENVALUES --   -20.58265 -11.33946  -1.39265  -0.87259  -0.69715
   1 1   C  1S          0.00000   0.99566  -0.11060  -0.16262   0.00000
   2        2S          0.00047   0.02675   0.20981   0.33995   0.00000
   3        2PX         0.00000   0.00000   0.00000   0.00000   0.00000
   4        2PY         0.00000   0.00000   0.00000   0.00000   0.42017
   5        2PZ        -0.00007   0.00066   0.17259  -0.18451   0.00000
   6        3S         -0.00024  -0.00743   0.08051   0.31309   0.00000
   7        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
   8        3PY         0.00000   0.00000   0.00000   0.00000   0.15760
   9        3PZ        -0.00048   0.00135  -0.01160  -0.07970   0.00000
  10        4XX        -0.00002  -0.00272  -0.01628  -0.01333   0.00000
  11        4YY        -0.00006  -0.00202  -0.01365   0.03019   0.00000
  12        4ZZ        -0.00074  -0.00123   0.03302  -0.00166   0.00000
  13        4XY         0.00000   0.00000   0.00000   0.00000   0.00000
  14        4XZ         0.00000   0.00000   0.00000   0.00000   0.00000
  15        4YZ         0.00000   0.00000   0.00000   0.00000  -0.01394
  16 2   O  1S          0.99472  -0.00038  -0.19672   0.08889   0.00000
  17        2S          0.02094   0.00025   0.44184  -0.20351   0.00000
  18        2PX         0.00000   0.00000   0.00000   0.00000   0.00000
  19        2PY         0.00000   0.00000   0.00000   0.00000   0.32122
  20        2PZ        -0.00153  -0.00029  -0.13537  -0.14216   0.00000
  21        3S          0.00436  -0.00058   0.37895  -0.27048   0.00000
  22        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  23        3PY         0.00000   0.00000   0.00000   0.00000   0.17976
  24        3PZ         0.00006   0.00108  -0.04718  -0.06799   0.00000
  25        4XX        -0.00418   0.00015  -0.00022  -0.00041   0.00000
  26        4YY        -0.00383  -0.00011  -0.00073  -0.00413   0.00000
  27        4ZZ        -0.00356  -0.00019   0.01969   0.00906   0.00000
  28        4XY         0.00000   0.00000   0.00000   0.00000   0.00000
  29        4XZ         0.00000   0.00000   0.00000   0.00000   0.00000
  30        4YZ         0.00000   0.00000   0.00000   0.00000  -0.02339
  31 3   H  1S         -0.00002  -0.00020   0.03017   0.17902   0.19082
  32        2S         -0.00013   0.00210  -0.00537   0.06479   0.12026
  33 4   H  1S         -0.00002  -0.00020   0.03017   0.17902  -0.19082
  34        2S         -0.00013   0.00210  -0.00537   0.06479  -0.12026
                           6         7         8         9        10
                        (A1)--O   (B1)--O   (B2)--O   (B1)--V   (A1)--V EIGENVALUES --    -0.63950  -0.52294  -0.44073   0.13573   0.24842
   1 1   C  1S          0.01942   0.00000   0.00000   0.00000  -0.12212
   2        2S         -0.06075   0.00000   0.00000   0.00000   0.14896
   3        2PX         0.00000   0.32517   0.00000   0.40259   0.00000
   4        2PY         0.00000   0.00000  -0.19811   0.00000   0.00000
   5        2PZ        -0.37597   0.00000   0.00000   0.00000  -0.21086
   6        3S          0.03971   0.00000   0.00000   0.00000   1.98096
   7        3PX         0.00000   0.21231   0.00000   0.71124   0.00000
   8        3PY         0.00000   0.00000  -0.04477   0.00000   0.00000
   9        3PZ        -0.08856   0.00000   0.00000   0.00000  -0.74977
  10        4XX         0.00549   0.00000   0.00000   0.00000  -0.00273
  11        4YY         0.02734   0.00000   0.00000   0.00000  -0.01265
  12        4ZZ        -0.01933   0.00000   0.00000   0.00000  -0.00459
  13        4XY         0.00000   0.00000   0.00000   0.00000   0.00000
  14        4XZ         0.00000   0.03558   0.00000  -0.03288   0.00000
  15        4YZ         0.00000   0.00000   0.06035   0.00000   0.00000Sample Text
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaona158

木虫 (著名写手)

自由人

EIGENVALUES 后面的数字就是这个轨道对应的能量,但是它的单位是HF,一般使的时候需要换成电子福特,用这个系数乘27.2116就可以了。
2楼2008-04-03 08:13:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer

★ ★
zzgyb(金币+2,VIP+0):谢谢你的帮助,欢迎继续关注计算模拟版!
那个应该是基函数对应的能量,现在量子化学中的基组已经超出了原子轨道的范围。一般来说只有HOMO和LUMO是感兴趣的
聪明人知道自己是傻子,傻子以为自己是聪明人
3楼2008-04-03 08:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer

另外,注意#后面的关键字,如果选简略输出,会少很多信息
聪明人知道自己是傻子,傻子以为自己是聪明人
4楼2008-04-03 08:42:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliracial

银虫 (小有名气)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的帮助,欢迎继续关注计算模拟版!
以 1  (A1)--O分子轨道为例,O表示已占分子轨道,其下对应的EIGENVALUES --   -20.58265,就是该轨道的本征值,即轨道能量。在该值下面对应一长串的各原子各轨道对应的数值,就是各原子轨道组成该分子轨道的系数。原子轨道系数绝对值越大,说明该原子轨道对该分子轨道的贡献越大。对1  (A1)--O分子轨道而言,很明显主要是由O  1S 原子轨道组成,其轨道系数达到0.99472。
修炼成一块愚钝的大石头
5楼2008-04-03 08:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaohongjie

非常感谢大家的帮助!!!!
6楼2008-04-03 09:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaohongjie

问题

若是我还想算一下原子轨道贡献率呢?比如说算 (A1)--O分子轨道中碳C原子的,我要把它后面对应的1S、2S、。。。。。。直到4YZ对应的系数都加在一块吗?谢谢!
7楼2008-04-03 09:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaochao85

铜虫 (正式写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的帮助,欢迎继续关注计算模拟版!
那你就做一个nbo 基本上你想要的轨道信息 一般都会有
8楼2008-04-03 10:17:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaohongjie

nbo请问是什么呢?

不好意思,我是新手
9楼2008-04-03 10:22:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliracial

银虫 (小有名气)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的帮助,欢迎继续关注计算模拟版!
引用回帖:
Originally posted by zhaohongjie at 2008-4-3 09:37:
若是我还想算一下原子轨道贡献率呢?比如说算 (A1)--O分子轨道中碳C原子的,我要把它后面对应的1S、2S、。。。。。。直到4YZ对应的系数都加在一块吗?谢谢!

根据原子轨道线形组合成分子轨道,各个系数的平方之和应该等于1。要算贡献也应该是计算各系数的平方,而不是简单的加和。你可以看看结构化学或者量子化学这方面的资料。
nbo,是自然键轨道,Guassian03手册里面有这个关键字,可以去查查手册,看怎么用。
修炼成一块愚钝的大石头
10楼2008-04-03 12:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhaohongjie 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见