当前位置: 首页 > 量子化学 >(求助)如何用GAMESS通过CI方法计算He(氦)的能量

(求助)如何用GAMESS通过CI方法计算He(氦)的能量

作者 qizhuixia
来源: 小木虫 200 4 举报帖子
+关注

如何用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 返回小木虫查看更多

今日热帖
  • 精华评论
  • paramecium86

    用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,

  • beefly

    不推荐用GAMESS做多参考计算,除非你是做方法的,否则里面的坑太多。建议用molpro、molcas、orca这些“半黑箱”的程序;等你熟悉这些程序以后,可以回头再用GAMESS,就容易理解了

  • qizhuixia

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

    没办法,导师只让用这个……

  • 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

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓