| 查看: 1229 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
（求助）如何用GAMESS通过CI方法计算He（氦）的能量已有1人参与
|
|||
|
如何用GAMESS通过CI方法计算He（氦）的能量 GAMESS小白，导师直接让我算氦的能量，用CI方法，不知道如何输入…… 下面是CH2的例子，改了半天也不对，求救！ $CONTRL SCFTYP=RHF CITYP=GUGA RUNTYP=ENERGY $END $SYSTEM TIMLIM=1 $END $BASIS GBASIS=STO NGAUSS=3 $END $GUESS GUESS=HUCKEL $END ! look at all state symmetries, by using C1 symmetry $CIDRT GROUP=C1 IEXCIT=2 NFZC=1 NDOC=3 NVAL=3 $END ! lowest singlet is 1-A-1, 1st excited singlet is 1-B-1 $GUGDIA NSTATE=2 $END $DATA Methylene...CI...STO-3G basis Cnv 2 Carbon 6.0 Hydrogen 1.0 0.0 0.82884 0.7079 $END |
回复此楼» 收录本帖的淘帖专辑推荐
 Gaussian |
» 猜你喜欢
北京师范大学胡德平课题组以及袁萌伟课题组考研招生
已经有0人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有218人回复
-
祈福论文小修
已经有21人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12738.4
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2019-03-14 18:18:21
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 586 (博士)
- 贵宾: 0.01
- 金币: 26316.1
- 红花: 348
- 帖子: 2021
- 在线: 1440.5小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2019-03-14 08:46:02
qizhuixia: 金币+30, ★★★★★最佳答案, 厉害，我用STO-3G主要是要一步一步计算He的能量，从最小basis set开始，然后对比计算出的能量的差异 2019-03-14 20:11:04
感谢参与，应助指数 +1
zhou2009: 金币+5 2019-03-14 08:46:02
qizhuixia: 金币+30, ★★★★★最佳答案, 厉害，我用STO-3G主要是要一步一步计算He的能量，从最小basis set开始，然后对比计算出的能量的差异 2019-03-14 20:11:04
|
用GAMESS自带的例子修改是可以的。只是要把基组改一下。不要用 STO-3G 这个最小基。这个算完 RHF 是没有virtual orbitals信息的。如果没有virtual orbital信息 就没法设置CI了。因为He本身就是一个1s轨道 想算CI就得包括虚轨道信息。大概改成这样就行了 $CONTRL SCFTYP=RHF CITYP=GUGA RUNTYP=ENERGY $END $basis gbasis=n31 ngauss=6 ndfunc=0 npfunc=1 diffsp=.f. $end 这个地方改成了 6-31g(p) $GUESS GUESS=HUCKEL $END $CIDRT GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=1 $END 这行设置的意思是 C1 对称， 计算singles and doubles CI，没有Frozen core，有1个全满轨道 ，有1个空轨道 $GUGDIA NSTATE=2 $END 计算两个CI states $GUGDM NFLGDM(1)=1,1 IROOT=1 $END 储存输出第一个state（也就是基态）信息 $GUGDM2 WSTATE(1)=1.0,0.0 $END $DATA He...CI...6-31g(p) basis C1 He 2.0 $END |
2楼2019-03-14 04:47:37
4楼2019-03-14 20:11:30
5楼2019-03-14 20:18:28