Gamess中如何设置CISD方法？

978142355 · 发表于 Post on 2017-4-21 10:16:55

如题。利用Gamess在CISD/6-31G *级别下计算甲烷的单点能。在Gamess中设关键词为：
$CONTRL RUNTYP=ENERGY CITYP=GUGA $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$CIDRT IEXCIT=2 NDOC=5 NVAL=14 $END
$DATA
cisd/6-31g* energy
C1
C 6.0             0.00000000 0.00000000 0.00000000
H 1.0             0.61776500 0.61776500 0.61776500
H 1.0             -0.61776500 -0.61776500 0.61776500
H 1.0             -0.61776500 0.61776500 -0.61776500
H 1.0             0.61776500 -0.61776500 -0.61776500
$END
但是通过高斯和gamess算出的能量差值太大（0.02hartree），是不是我gamess哪个地方弄错了呢？谢谢各位老师。

zjxitcc · 发表于 Post on 2019-3-18 22:47:53

echa 发表于 2019-3-18 21:39
请问有做过结构优化的吗？我是个GAMESS 的新手，我的结构老是出现gradient out of range. 希望有人帮帮我 ...

信息太少，最好贴出输入输出文件，便于大家分析。

sobereva · 发表于 Post on 2017-4-22 08:20:50

yflchx 发表于 2017-4-21 11:26
不用花钱、窥探黑箱、Fragment Molecular Orbital (FMO)、Effective Fragment Potential (EFP)、VSCF、Nu ...

GAMESS-US 16版还一个突出优点是LUT-IOTC有解析梯度。另外还有spin-flip、MOM等有用功能

978142355 · 发表于 Post on 2017-4-21 19:07:15

yflchx 发表于 2017-4-21 18:01
------------
EIGENVECTORS
------------

明白了，谢谢您的耐心指导和帮助

yflchx · 发表于 Post on 2017-4-21 18:01:09

978142355 发表于 2017-4-21 15:53
谢谢。原谅我基础比较薄弱………………1和4是没问题，23又是怎么回事？

      ------------
      EIGENVECTORS
      ------------
看看GAMESS输出文件的这个部分，看总数一共几个? 代表什么意思?

这样的话就明白了吧。

978142355 · 发表于 Post on 2017-4-21 15:53:22

yflchx 发表于 2017-4-21 15:20
1+4+18=23

这样就清楚了吧。

谢谢。原谅我基础比较薄弱………………1和4是没问题，23又是怎么回事？

yflchx · 发表于 Post on 2017-4-21 15:20:17

978142355 发表于 2017-4-21 15:09
谢谢您这么详细的讲解。请问这里NVAL为什么要设定18呢？

1+4+18=23

这样就清楚了吧。

978142355 · 发表于 Post on 2017-4-21 15:09:44

yflchx 发表于 2017-4-21 14:40
Sorry，是我说的不清楚。NFZC、NDOC、NVAL都很重要，要重复程序的结果需要保持严格一致。后面附的是Firefly ...

谢谢您这么详细的讲解。请问这里NVAL为什么要设定18呢？

yflchx · 发表于 Post on 2017-4-21 14:40:43

本帖最后由 yflchx 于 2017-4-21 14:45 编辑

978142355 发表于 2017-4-21 14:33
设置了还是不行……

Sorry，是我说的不清楚。NFZC、NDOC、NVAL都很重要，要重复程序的结果需要保持严格一致。后面附的是Firefly的运行结果。

$CONTRL RUNTYP=ENERGY CITYP=GUGA $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$CIDRT IEXCIT=2 NFZC=1 NDOC=4 NVAL=18 $END
$DATA
cisd/6-31g* energy
C1
C 6.0             0.00000000 0.00000000 0.00000000
H 1.0             0.61776500 0.61776500 0.61776500
H 1.0             -0.61776500 -0.61776500 0.61776500
H 1.0             -0.61776500 0.61776500 -0.61776500
H 1.0             0.61776500 -0.61776500 -0.61776500
$END

      -----------------
      ENERGY COMPONENTS
      -----------------
      WAVEFUNCTION NORMALIZATION =    1.0000000000
                     ONE ELECTRON ENERGY =    -80.1731758032
                  TWO ELECTRON ENERGY =    26.1410503269
            NUCLEAR REPULSION ENERGY =    13.6865143207
                                    ------------------
                     TOTAL ENERGY =    -40.3456111556

978142355 · 发表于 Post on 2017-4-21 14:33:30

yflchx 发表于 2017-4-21 12:07
关键是：NFZC = number of CI frozen core MOs (in $CIDRT only).

你试试是不是这样的?

设置了还是不行……

yflchx · 发表于 Post on 2017-4-21 12:07:22

978142355 发表于 2017-4-21 12:00
感觉有可能，请问这个虚轨道设置多少算是合适的？

关键是：NFZC = number of CI frozen core MOs (in $CIDRT only).

你试试是不是这样的?

978142355 · 发表于 Post on 2017-4-21 12:00:59

yflchx 发表于 2017-4-21 11:40
Gaussian: Range of M.O.s used for correlation: 2 23
NBasis= 23 NAE= 5 NBE= 5 NF ...

感觉有可能，请问这个虚轨道设置多少算是合适的？

978142355 · 发表于 Post on 2017-4-21 12:00:12

我本是个娃娃发表于 2017-4-21 10:20
怎么又开始鼓捣GAMESS了？ORCA多牛叉！

哈哈，导师想要在这个上面整点大新闻，0.0

yflchx · 发表于 Post on 2017-4-21 11:40:14

Gaussian:  Range of M.O.s used for correlation:    2 23
NBasis= 23 NAE=    5 NBE=    5 NFC=    1 NFV=    0
NROrb=    22 NOA=    4 NOB=    4 NVA= 18 NVB= 18

GAMESS: $CIDRT IEXCIT=2 NDOC=5 NVAL=14 $END

问题是不是出在这里？

yflchx · 发表于 Post on 2017-4-21 11:26:53

我本是个娃娃发表于 2017-4-21 10:40
除了能量分解，GAMESS的优点几乎就没有了

不用花钱、窥探黑箱、Fragment Molecular Orbital (FMO)、Effective Fragment Potential (EFP)、VSCF、Nuclear-Electronic Orbital (NEO)、GVB、Spin-Orbit Coupling。

用户少也就体现不出优点了。当然，有些应该算特殊点。

我本是个娃娃 · 发表于 Post on 2017-4-21 10:40:11

冰释之川发表于 2017-4-21 10:22
听说戈梅斯也出了最新版了

除了能量分解，GAMESS的优点几乎就没有了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] Gamess中如何设置CISD方法？

评分 Rate

评分 Rate

浏览过的版块