$ORMAS group (required by MCSCF if CISTEP=ORMAS)
(required for CITYP=ORMAS or SFORM)
This group partitions an active space, defined in $DET
or $CIDET, into Occupation Restricted Multiple Active
Spaces (ORMAS). All possible determinants satisfying the
occupation restrictions (and of course the space symmetry
restriction given in $DET/$CIDET) will be generated. This
group's usefulness lies in reducing the large number of
determinants present in full CI calculations with large
active spaces.
There are no sensible defaults for these inputs, but if
the group is entirely omitted, a full CI calculation will
be performed. That is, the defaults are
NSPACE=1, MSTART(1)=NCORE+1, MINE(1)=NELS, MAXE(1)=NELS
meaning all active orbitals are in one partition.
NSPACE = number of orbital groups you wish to partition
the active space (NACT in $DET/$CIDET) into.
MSTART = an array of NSPACE integers. These specify where
each orbital group starts in the full list. You
must not overlook the NCORE core orbitals in
computing MSTART values. Space I runs from
orbital MSTART(I) up to orbital MSTART(I+1)-1,
or NACT+NCORE if I is the last space, I=NSPACE.
IMPORTANT !!!! Remember to make sure your orbitals have
been reordered to suit MSTART, using NORDER in $GUESS.
MINE = an array of NSPACE integers. These specify the
minimum numbers of electrons that must always
occupy the orbital groups. In other words,
MINE(I) is the minimum number of electrons that
can occupy space I in any of the determinants.
MAXE = an array of NSPACE integers. These specify the
maximum numbers of electrons that must always
occupy the orbital groups. In other words,
MAXE(I) is the maximum number of electrons that
can occupy space I in any of the determinants.
The number of active electrons is NELS in $DET or $CIDET,
and the program will check that MINE/MAXE values are
consistent with this total number.
Input for Coupled Electron Pair Approximations
These two keywords provide size-extensivity corrections for
singly and doubly excited CI calculations. At present,
these approximate corrections may only be applied to a
single state. The CISD may follow RHF (single reference,
SR) or complete active space (multireference, MR)
calculations. Usually, the final ORMAS orbital space will
be the entire external space of the SR-CISD or MR-CISD
calculations. The excitation level (MAXE) into this space
should only be two. The number of electrons being
correlated (NELS in $CIDET) is denoted Ne. The value of a
parameter G determines the weight of the approximate size-
extensivity correction, with the literature containing
several possible choices: GVAL here sometimes is written
as -(1-G) in the literature. The first order density
matrix for the chosen CEPA-style correction will be
generated, and used for property calculation.
CEPA = NONE ordinary SR-CISD or MR-CISD (default)
= CEPA0 Coupled Electron Pair Approximation,
whose GVAL= 1.0
= ACPF Averaged Coupled Pair Functional,
whose GVAL= (Ne-2)/Ne
= AQCC Averaged Quadratic Coupled-Cluster,
whose GVAL= (Ne-3)(Ne-2)/[Ne(Ne-1)]
= READ user will supply the desired GVAL.
GVAL = is given only when CEPA=READ
An example of the SR-AQCC size-extensivity correction for a
molecule with 10 chemical cores, 8 occupied valence
orbitals, and 163 total MOs is
$CONTRL SCFTYP=RHF CITYP=ORMAS RUNTYP=ENERGY
$CIDET NCORE=10 NACT=153 NELS=16
$ORMAS CEPA=AQCC NSPACE=2
MSPACE(1)=11,19 MINE(1)=14,0 MAXE(1)=16,2
A MR-CEPA case with N active electrons and M electrons in
filled valence orbitals should use three spaces, with NELS
M+N; MINE M-2,N-2,0; and MAXE M,N+2,2.
Review: P.G.Szalay, in Modern Ideas in Coupled-Cluster
methods, R.J.Bartlett (ed), World Scientific, Singapore
(1997), pp 81-123.
* * * * *
BLOCK = a flag to request that natural orbital generation
for CISTEP=ORMAS MCSCF or CITYP=ORMAS CI runs
should prevent mixing between the NSPACE orbital
subspaces. This means the NOs only diagonalize
the diagonal blocks of the density, and are thus
not the genuine NOs. However, these approximate
NOs can be used with MOREAD to exactly reproduce
the ORMAS energy, which is invariang to rotations
within the orbital subspaces. (Default=.FALSE.)
QCORR = a flag to request Davidson-style +Q corrections.
If this is not sensible for your CI choice, the
program will not print this correction, anyway.
The default is .TRUE.
FDIRCT = a flag to choose storage in memory of some
intermediates. This is very large, and slower in
the case of many occupied orbitals, but helpful
with a smaller number of orbitals. Therefore the
default for this is .TRUE. for MCSCF runs, but
.FALSE. during CI computations.
*** See REFS.DOC for more information on using ORMAS ***
121 lines are written.
Edited by Shiro KOSEKI on Tue May 17 15:19:38 2022.