24小时热门版块排行榜    

查看: 1548  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见