| 查看: 1299 | 回复: 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 |
» 猜你喜欢
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
-
化学化工学院 招收化工、化学、材料等相关方向研究生（学硕、专硕都有调剂名额）
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有258人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
温州大学新能源材料微观溯源课题组招收2026级博士研究生（材料、化学方向）
已经有6人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流
已经有5人回复
-
26申博
已经有2人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12986.2
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2019-03-14 18:18:21
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 587 (博士)
- 贵宾: 0.01
- 金币: 26383.1
- 红花: 349
- 帖子: 2025
- 在线: 1442小时
- 虫号: 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
