GAMESS轨道读取 (以及用于ORMAS计算提取分子轨道下积分)
(2012-01-12 10:18:43)
标签:
gamess
轨道读取
casci
杂谈
|
分类:
学习科研
|
忙活了好几天,一直还纳闷为什么用GAMESS算的CAS-CI的方法算不对。因为用自己编写的full-CI程序来做CAS-CI的话可以得到合理的结果,但是自己的程序还是很慢的,full-CI部分离实际应用还是很远很远的。(不过QC_orb_opt部分还是可以接近实际应用的)今天过来又翻GAMESS的INPUT.DOC,终于在里面发现错误的原因了。就是如果$GAUSS里面采用MOREAD读入轨道的话,默认是只读取占据轨道的系数的,空轨道按照C'SC=1随便给。如果想同时读入空轨道的话,需要在$GAUSS里面同时指定读入的轨道数(NORB=?)。记录一下,毕竟GAMESS不是自己写的,防止忘记。
+++++++++++++++++++++++++++++++++++++++++++++++++++++++
GAMESS里面INPUT.DOC关于此部分的说明
+++++++++++++++++++++++++++++++++++++++++++++++++++++++
NORB = The number of orbitals
to be read in the $VEC
group. This applies only to GUESS=MOREAD.
For -RHF-, -UHF-, -ROHF-, and -GVB-, NORB defaults to the
number of occupied
orbitals. NORB must be given for
-CI-
and -MCSCF-. For -UHF-, if NORB is not given,
only the
occupied alpha and beta orbitals should be given, back to
back. Otherwise, both alpha and beta orbitals
must
consist of NORB vectors.
NORB may be larger than the number of occupied MOs, if you
wish to read in the virtual orbitals. If NORB is
less
than the number of atomic orbitals, the remaining orbitals
are generated as the orthogonal complement to those read.
++++++++++++++++++++ sample input file
+++++++++++++++++++++++
$CONTRL SCFTYP=NONE CITYP=ORMAS RUNTYP=ENERGY
$END
$BASIS GBASIS=STO NGAUSS=6
$END
$GUESS GUESS=MOREAD
NORB=6 $END
$SYSTEM MWORDS=1080 $END
$CIDET
ncore=0 nact=6 nels=6 NSTATE=4 SZ=0.0
$END
$ORMAS
NSPACE=6 MSTART(1)=
1 MINE(1)=2 MAXE(1)=2
MSTART(2)= 2 MINE(2)=2
MAXE(2)=2
MSTART(3)= 3 MINE(3)=0
MAXE(3)=2
MSTART(4)= 4 MINE(4)=0
MAXE(4)=0
MSTART(5)= 5 MINE(5)=0
MAXE(5)=0
MSTART(6)= 6 MINE(6)=0
MAXE(6)=2
$END
$DATA
H_chainsX4
C1
H
1
0.000000
0.000000
0.000000
H
1
0.000000
0.741400
0.000000
H
1
0.000000
3.707000
0.000000
H
1
0.000000
4.448400
0.000000
H
1
0.000000
7.414000
0.000000
H
1
0.000000
8.155400
0.000000
$END
$VEC
1 1 0.54889004E+00
0.54927915E+00-0.11744229E-01 0.57920941E-02-0.74667599E-04
1 2 0.35148626E-04
2 1
0.57893939E-02-0.11741529E-01 0.54915193E+00
0.54915193E+00-0.11741529E-01
2 2 0.57893939E-02
3 1
0.35148626E-04-0.74667599E-04 0.57920941E-02-0.11744229E-01
0.54927915E+00
3 2 0.54889004E+00
4 1
0.12107377E+01-0.12109142E+01 0.15331675E-01-0.73810645E-02
0.94148530E-04
4 2-0.44344763E-04
5