| 查看: 1250 | 回复: 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人回复
-
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有192人回复
-
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有1人回复
-
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘
已经有0人回复
-
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复
-
福州大学新能源材料与工程研究院招收2026年入学博士
已经有0人回复
-
海南师范大学2026年博士研究生招收 (在职想提升学历人员可报考) 申请考核制
已经有0人回复
-
锂离子电池循环寿命衰减过快，求机理分析与改进思路
已经有0人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 587 (博士)
- 贵宾: 0.01
- 金币: 26352.6
- 红花: 348
- 帖子: 2024
- 在线: 1441.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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12877.4
- 散金: 27307
- 红花: 300
- 帖子: 4668
- 在线: 925.8小时
- 虫号: 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