(求助)如何用GAMESS通过CI方法计算He(氦)的能量
如何用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 返回小木虫查看更多
今日热帖
用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,
不推荐用GAMESS做多参考计算,除非你是做方法的,否则里面的坑太多。建议用molpro、molcas、orca这些“半黑箱”的程序;等你熟悉这些程序以后,可以回头再用GAMESS,就容易理解了
没办法,导师只让用这个……
厉害老哥,STO-3G无法用CI计算,我说我搞算了半天都是error