24小时热门版块排行榜    

查看: 1301  |  回复: 4

qizhuixia

新虫 (初入文坛)

[求助] (求助)如何用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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

不推荐用GAMESS做多参考计算,除非你是做方法的,否则里面的坑太多。建议用molpro、molcas、orca这些“半黑箱”的程序;等你熟悉这些程序以后,可以回头再用GAMESS,就容易理解了
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2019-03-14 18:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by beefly at 2019-03-14 18:18:21
不推荐用GAMESS做多参考计算,除非你是做方法的,否则里面的坑太多。建议用molpro、molcas、orca这些“半黑箱”的程序;等你熟悉这些程序以后,可以回头再用GAMESS,就容易理解了

没办法,导师只让用这个……
4楼2019-03-14 20:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by paramecium86 at 2019-03-14 04:47:37
用GAMESS自带的例子修改是可以的。只是要把基组改一下。不要用 STO-3G 这个最小基。这个算完 RHF 是没有virtual orbitals信息的。如果没有virtual orbital信息 就没法设置CI了。因为He本身就是一个1s轨道
想算CI就 ...

厉害老哥,STO-3G无法用CI计算,我说我搞算了半天都是error
5楼2019-03-14 20:18:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qizhuixia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 320求调剂 +5 魏zy 2026-03-08 5/250 2026-03-14 03:41 by JourneyLucky
[考研] 一志愿中国科学技术大学,086000,生物与医药,求调剂 +3 12341234789 2026-03-08 4/200 2026-03-14 03:14 by JourneyLucky
[考研] 一志愿天大化工专调剂总分338 +3 蔡大美女 2026-03-08 3/150 2026-03-14 03:12 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见