| 查看: 1083 | 回复: 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朵）
fwang2011» 猜你喜欢
申请2027年材料/化学类博士
已经有6人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有105人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
投稿咨询
已经有6人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用Multiwfn做拓扑分析以及计算孤对电子角度
已经有18人回复
- 怎么得到旋轨耦合中不同J值的本征值和波函数？
已经有6人回复
- 【求助】用Gaussview画分子轨道图和数据分析的结果不一致的问题
已经有8人回复
- 【求助】VASP怎样精确计算出费米能级处的态密度值？
已经有22人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- 【讨论】深问：原子轨道对分子轨道的贡献，函数如何如何归属到原子轨道里？
已经有3人回复
- 【重点讨论】Castep中的Empty Band的用处和设置？（参与讨论就有机会赢取大礼包！）
已经有54人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13331
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2012-03-04 10:35:40
fwang2011
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1682.9
- 红花: 1
- 帖子: 57
- 在线: 21.7小时
- 虫号: 1356013
- 注册: 2011-07-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-03-04 13:18:50
fwang2011
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1682.9
- 红花: 1
- 帖子: 57
- 在线: 21.7小时
- 虫号: 1356013
- 注册: 2011-07-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
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
40