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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » 分子轨道数 和 本征值数目 的关系
41/1返回列表
查看: 856  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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朵）

fwang2011

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-03-04 00:34:37
已阅   回复此楼   关注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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭