| 查看: 877 | 回复: 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 |
回复此楼
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 564 (博士)
- 贵宾: 0.01
- 金币: 24607.4
- 红花: 338
- 帖子: 1969
- 在线: 1335.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
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10769.2
- 散金: 27307
- 红花: 295
- 帖子: 4586
- 在线: 916.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2019-03-14 18:18:21
4楼2019-03-14 20:11:30
5楼2019-03-14 20:18:28
15