In this section keywords for the coupled cluster program are defined. In particular the coupled cluster model(s) and other parameters common to all submodules are specified.
.IMSKIP
.
CCSTST
disables all terms which depend on the
double excitation amplitudes or multipliers. This flag must be
set for CCS finite field calculations.
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
READ (LUCMD,*) EFIELD
READ (LUCMD,'(1X,A8)') LFIELD
Add an external finite field (operator label LFIELD
)
of strength EFIELD
to the Hamiltonian
(same input format as in Sec. 24.2.8).
These fields are only included in the CC calculation and
not in the calculation of the SCF reference state
(i.e. the orbitals).
Using this option in the calculation of numerical derivatives
thus gives so-called orbital-unrelaxed energy derivatives,
which is standard in coupled cluster response function theory.
Note that this way of adding an external field does not work for
models including triples excitations (most notably CCSD(T) and CC3).
For finite field calculations of orbital-relaxed energy
derivatives the field must be included in the SCF calculation.
For such calculations, use the input section *HAMILTONIAN (Sec. 24.2.8).
READ (LUCMD,*) NFC, NFV
Specify the number of frozen core orbitals and number of frozen virtuals (0 is recommended for the latter). The program will automatically freeze the NFC(NFV) orbitals of lowest (highest) orbital energy among the canonical Hartree-Fock orbitals. For benzene a relevant choice is for example 6 0.
READ (LUCMD,*) (NRHFFR(I),I=1,MSYM)
READ (LUCMD,*) (NVIRFR(I),I=1,MSYM)
Specify for each irreducible representation how
many orbitals should be frozen (deleted) for the coupled
cluster calculation. In calculations, the first NRHFFR(I)
orbitals will be kept frozen in symmetry class I
and
the last NVIRFR(I)
orbitals will be deleted from the
orbital list.
READ (LUCMD,*) (NRHFFR(I),I=1,MSYM)
DO ISYM = 1, MSYM
IF (NRHFFR(ISYM.NE.0) THEN
READ (LUCMD,*) (KFRRHF(J,ISYM),J=1,NRHFFR(ISYM))
END IF
END DO
READ (LUCMD,*) (NVIRFR(I),I=1,MSYM)
DO ISYM = 1, MSYM
IF (NVIRFR(ISYM.NE.0) THEN
READ (LUCMD,*) (KFRVIR(J,ISYM),J=1,NVIRFR(ISYM))
END IF
END DO
Specify explicitly for each irreducible representation the orbitals that should be frozen (deleted) in the coupled cluster calculation.
.IMSKIP
.
DIRECT
keyword was set in the general
input section of DALTON. IF HERDIR
IS NOT SPECIFIED THE ERI PROGRAM IS USED AS INTEGRAL GENERATOR.)
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
READ (LUCMD,'(I5)') MAXITE
MAXIMUM NUMBER OF ITERATIONS FOR WAVE FUNCTION OPTIMIZATION
(DEFAULT IS MAXITE = 40
).
READ (LUCMD,*) MAXRED
MAXIMUM DIMENSION OF THE REDUCED SPACE FOR THE
SOLUTION OF LINEAR EQUATIONS (DEFAULT IS MAXRED = 200
).
READ (LUCMD,*) MXDIIS
MAXIMUM NUMBER OF ITERATIONS FOR DIIS ALGORITHM
BEFORE WAVE FUNCTION INFORMATION IS DISCARDED AND A NEW DIIS
SEQUENCE IS STARTED
(DEFAULT IS MXDIIS = 8
).
READ (LUCMD, *) MXLRV
MAXIMUM NUMBER OF TRIAL VECTORS IN THE SOLUTION OF LINEAR EQUATIONS. IF THE NUMBER OF TRIAL VECTORS REACHES THIS VALUE, ALL TRIAL VECTORS AND THEIR TRANSFORMATIONS WITH THE JACOBIAN ARE SKIPPED AND THE ITERATIVE PROCEDURE FOR THE SOLUTION OF THE LINEAR (I.E. THE RESPONSE) EQUATIONS IS RESTARTED FROM THE CURRENT OPTIMAL SOLUTION.
READ (LUCMD, *) NSIMLE
SET THE MAXIMUM NUMBER OF RESPONSE EQUATIONS THAT SHOULD BE SOLVED SIMULTANEOUSLY. DEFAULT IS 0 WHICH MEANS THAT ALL COMPATIBLE EQUATIONS (SAME EQUATION TYPE AND SYMMETRY CLASS) ARE SOLVED SIMULTANEOUSLY.
READ (LUCMD,*) MSYM2
SET NUMBER OF IRREDUCIBEL REPRESENTATIONS.
.IMSKIP
.
READ (LUCMD,'(I5)') IPRINT
SET PRINT PARAMETER FOR COUPLED CLUSTER PROGRAM
(DEFAULT IS TO TAKE THE VALUE IPRUSR
, SET IN THE GENERAL
INPUT SECTION OF DALTON).
.IMSKIP
.
.IMSKIP
.
.IMSKIP
.
READ (LUCMD,*) THRENR
SET THRESHOLD FOR CONVERGENCE OF THE GROUND STATE ENERGY.
READ (LUCMD,*) THRENR
SET THRESHOLD FOR CONVERGENCE OF THE RESPONSE EQUATIONS.
.IMSKIP
.