24小时热门版块排行榜    

Znn3bq.jpeg
查看: 977  |  回复: 3

fwang2011

金虫 (小有名气)

[求助] 分子轨道数 和 本征值数目 的关系

可能这并不仅仅是个GAMESS的问题,更像一个基本的量化计算问题。

以水为例,if we run the following:
$CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE MAXIT=30 MULT=1 MOLPLT=.TRUE. $END
$SYSTEM TIMLIM=525600 MEMORY=1000000 $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$SCF DIRSCF=.TRUE. $END
$STATPT OPTTOL=0.0001 NSTEP=200 HSSEND=.t. $END
$FORCE METHOD=ANALYTIC VIBANL=.TRUE. $END
$DATA
Test run of H2O
CNV 2

O     8.0     0.00000     0.00000     0.00000
H     1.0    -0.75800     0.00000     0.54500
$END

在输出文件的最后,we will find the following output in the final MO:
                     11         12         13         14         15
                    1.1785     1.3851     1.4312     2.0206     2.0306
                     A1         B1         A1         A1         A2
    1  O  1  S   -0.051055   0.000000  -0.085666  -0.007482   0.000000
    2  O  1  S   -0.052893   0.000000  -1.443495  -0.078523   0.000000
    3  O  1  X    0.000000  -1.040316   0.000000   0.000000   0.000000
    4  O  1  Y    0.000000   0.000000   0.000000   0.000000   0.000000
    5  O  1  Z    0.731788   0.000000  -0.501204  -0.002830   0.000000
    6  O  1  S    0.411818   0.000000   3.637788   0.234068   0.000000
    7  O  1  X    0.000000   1.540982   0.000000   0.000000   0.000000
    8  O  1  Y    0.000000   0.000000   0.000000   0.000000   0.000000
    9  O  1  Z   -0.315514   0.000000   1.158296   0.130567   0.000000
   10  O  1 XX   -0.249651   0.000000  -0.643896   0.563979   0.000000
   11  O  1 YY    0.080110   0.000000  -0.293083   0.385254   0.000000
   12  O  1 ZZ   -0.063729   0.000000  -0.403630  -1.011299   0.000000
   13  O  1 XY    0.000000   0.000000   0.000000   0.000000   1.000000
   14  O  1 XZ    0.000000   0.026841   0.000000   0.000000   0.000000
   15  O  1 YZ    0.000000   0.000000   0.000000   0.000000   0.000000
   16  H  2  S   -0.694183   0.108435  -0.301211  -0.080583   0.000000
   17  H  2  S    0.368089  -0.934728  -0.812287  -0.037357   0.000000
   18  H  3  S   -0.694183  -0.108435  -0.301211  -0.080583   0.000000
   19  H  3  S    0.368089   0.934728  -0.812287  -0.037357   0.000000

也就是说这里只有15个MO。

但是,检查一下我们知道:
TOTAL NUMBER OF BASIS SET SHELLS             =    8
NUMBER OF CARTESIAN GAUSSIAN BASIS FUNCTIONS =   19
这里”NUMBER OF CARTESIAN GAUSSIAN BASIS FUNCTIONS“ 等于19,也就是将有19个本征值。这也可以从输出文件中验证。

我的问题是:
如何解释这个不同?因为不同的basis选择就会有不同的NUMBER OF CARTESIAN GAUSSIAN BASIS FUNCTIONS,也就会有不同的数量的”eigenvalues“(最终解释为分子轨道)。

相关书籍推荐?

谢谢。
回复此楼

» 收录本帖的淘帖专辑推荐

我的收藏

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990(金币+2): 2012-03-04 21:09:38
分子轨道数根据使用的是球函数还是直角函数而有所不同。一般而言,纵向的原子轨道数也应该使用同样的球函数/直角函数设置。因此两个值应该是相同的。

但是GAMESS不同,它永远打印直角函数原子轨道,所以一旦出现了D,F,G函数,两个轨道数就不一样了。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-03-04 10:35:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fwang2011

金虫 (小有名气)

送鲜花一朵
多谢。

---
不知金币有没有送出。
3楼2012-03-04 13:18:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fwang2011

金虫 (小有名气)

For more information in this topic, please go to (you have to a member of it):
https://groups.google.com/d/topic/gamess/OBhW8maYlyE/discussion

For those who cannot/don't access the page, I copied the one which I think as answer for your reference.

And the reference is
in statp.scr
C*MODULE STATPT  *DECK SIGEND
       SUBROUTINE SIGEND
line 5359
NMOPRT = MIN(L0,NOCC+10)

in human language min(size of basis, number of occ+10)
min(19,10+5)=15
Robert Balawender
btw Sean

>>In 6-31G(d), 2*H has 2*4 primitives and O has 11 primitives
>> (including
>> the extra d function), giving a total of 19.

19 is a number of the contracted gaussian
and all your answer is wrong

- show quoted text -

--
dr Robert Balawender
Instytut Chemii Fizycznej PAN
ul. Kasprzaka 44/52
01-224 Warszawa
4楼2012-05-18 00:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fwang2011 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见