计算化学公社

标题: Gamess中如何设置CISD方法？ [打印本页]

作者
Author: 978142355 时间: 2017-4-21 10:16
标题: Gamess中如何设置CISD方法？
如题。利用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哪个地方弄错了呢？谢谢各位老师。

作者
Author: 我本是个娃娃 时间: 2017-4-21 10:20

怎么又开始鼓捣GAMESS了？ORCA多牛叉！

作者
Author: 冰释之川 时间: 2017-4-21 10:22

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

听说戈梅斯也出了最新版了

作者
Author: 我本是个娃娃 时间: 2017-4-21 10:40

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

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

作者
Author: yflchx 时间: 2017-4-21 11:26

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

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

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

作者
Author: yflchx 时间: 2017-4-21 11:40
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

问题是不是出在这里？

作者
Author: 978142355 时间: 2017-4-21 12:00

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

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

作者
Author: 978142355 时间: 2017-4-21 12:00

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

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

作者
Author: yflchx 时间: 2017-4-21 12:07

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

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

你试试是不是这样的?

作者
Author: 978142355 时间: 2017-4-21 14:33

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

你试试是不是这样的?

设置了还是不行……

作者
Author: yflchx 时间: 2017-4-21 14:40
本帖最后由 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

作者
Author: 978142355 时间: 2017-4-21 15:09

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

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

作者
Author: yflchx 时间: 2017-4-21 15:20

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

1+4+18=23

这样就清楚了吧。

作者
Author: 978142355 时间: 2017-4-21 15:53

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

这样就清楚了吧。

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

作者
Author: yflchx 时间: 2017-4-21 18:01

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

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

这样的话就明白了吧。

作者
Author: 978142355 时间: 2017-4-21 19:07

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

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

作者
Author: sobereva 时间: 2017-4-22 08:20

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

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

作者
Author: zjxitcc 时间: 2019-3-18 22:47

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)