24小时热门版块排行榜    

查看: 869  |  回复: 10
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

powerful77

铁虫 (初入文坛)

[交流] 【求助】不用gv等软件画CO2的HOMO

现已经采用3-21G基组计算得到CO2的HOMO, 由于我想得到分子轨道波函数的具体值所以不能用gv等软件。那么由计算果HOMO应该是:
0.36823*第一个氧的2PY轨道+0.44979*第一个氧的3PY轨道-0.36823第二个氧的2PY轨道-0.44979*第二个氧的3PY轨道。
那么现在就需要知道氧在3-21G基组下的2PY和3PY轨道,那么这些具体值可以到哪里去找呢?

                                  11         12          13           14            15
                              (PIG)--O  (PIU)--V  (PIU)--V  (SGG)--V  (SGU)--V
EIGENVALUES --  -0.51585   0.13583   0.13583   0.22990   0.68522
   1 1   C  1S         0.00000   0.00000   0.00000  -0.13819   0.00000
   2        2S           0.00000   0.00000   0.00000  -0.00031   0.00000
   3        2PX         0.00000   0.00000   0.46934   0.00000   0.00000
   4        2PY         0.00000   0.46934   0.00000   0.00000   0.00000
   5        2PZ         0.00000   0.00000   0.00000   0.00000  -0.02926
   6        3S           0.00000   0.00000   0.00000   1.59057   0.00000
   7        3PX         0.00000   0.00000   0.66656   0.00000   0.00000
   8        3PY         0.00000   0.66656   0.00000   0.00000   0.00000
   9        3PZ         0.00000   0.00000   0.00000   0.00000   2.88905
  10 2   O  1S       0.00000   0.00000   0.00000   0.07537   0.11154
  11        2S          0.00000   0.00000   0.00000  -0.05164  -0.03222
  12        2PX        0.00000   0.00000  -0.26811   0.00000   0.00000
  13        2PY        0.36823  -0.26811   0.00000   0.00000   0.00000
  14        2PZ        0.00000   0.00000   0.00000   0.28155   0.02754
  15        3S          0.00000   0.00000   0.00000  -0.61727  -1.74799
  16        3PX        0.00000   0.00000  -0.40972   0.00000   0.00000
  17        3PY        0.44979  -0.40972   0.00000   0.00000   0.00000
  18        3PZ        0.00000   0.00000   0.00000   0.56399   0.65528
  19 3   O  1S       0.00000   0.00000   0.00000   0.07537  -0.11154
  20        2S          0.00000   0.00000   0.00000  -0.05164   0.03222
  21        2PX        0.00000   0.00000  -0.26811   0.00000   0.00000
  22        2PY       -0.36823  -0.26811   0.00000   0.00000   0.00000
  23        2PZ        0.00000   0.00000   0.00000  -0.28155   0.02754
  24        3S          0.00000   0.00000   0.00000  -0.61727   1.74799
  25        3PX        0.00000   0.00000  -0.40972   0.00000   0.00000
  26        3PY       -0.44979  -0.40972   0.00000   0.00000   0.00000
  27        3PZ        0.00000   0.00000   0.00000  -0.56399   0.65528

[ Last edited by yjcmwgk on 2009-11-13 at 15:28 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★
powerful77(金币+2,VIP+0):谢谢! 9-16 22:02
erylingjet(金币+2,VIP+0):感谢交流! 9-18 10:56
用cubegen获取HOMO的空间格点文件,获取后应能把数据导入fortran,matlab之类.

这个HOMO空间各点的密度可以自己取,比如,沿x轴1埃取10个点,y轴1埃取5个点,z轴1埃取3个点。
5楼2009-09-16 17:54:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

powerful77

铁虫 (初入文坛)

没有人知道吗?
2楼2009-09-16 14:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★
powerful77(金币+2,VIP+0):很感谢你参加讨论,最近我一直头疼着个问题,还希望能继续向你请教 9-16 17:24
erylingjet(金币+3,VIP+0):感谢交流 9-18 10:56
你的问题不够明确,所以匆匆而过的人不能回答:

1、        你标题说“不用gv等软件画CO2的HOMO”,是个画HOMO图的问题。
这可以用cubegen作这HOMO的空间格点文件,用gsgrid程序作一个HOMO的截面,在sigmaplot中就可以作出HOMO的等值线图来。
为了理解你的提问,我也算了一下CO2,你的这个HOMO,应是在安排两个O的孤对电子(部分),且二O的孤对电子是相位相反的,体现二O的孤对电子是排斥的。但不是反键轨道,因为中间有C隔开了,相距甚远。

2、你说:我想得到分子轨道波函数的具体值。什么具体值?分子轨道的能级值、波函数系数值不是都有了吗?

3、说需要知道氧在3-21G基组下的2PY和3PY轨道这些具体值,是指这些原子的轨道的值吗?一个原子的轨道,它占有一个电子时的能级值是可算的,但在分子中,只是给它分配了系数,且这系数还不是归一的,是一些分数,这时的2PY和3PY是没有具体值可输出的。
3楼2009-09-16 15:32:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

powerful77

铁虫 (初入文坛)

非常感谢你的回答,你的回答使我受益匪浅。我所说的波函数具体值是指gauss计算后,假设HOMO的波函数是Ψ,我怎么知道Ψ(x,y,z)的值,x,y,z是笛卡尔空间坐标,而且x,y,z可以在空间内任取。

因此我原来的猜想是从高斯的out文件中,对于CO2我只能看出HOMO主要是由,2氧的2Py,3Py轨道贡献,因此我猜测如果知道氧的2py设为Ψ1,3py设为Ψ2,

[ Last edited by powerful77 on 2009-9-17 at 10:49 ]
4楼2009-09-16 17:22:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见