8.8.1 Variables set by the program

A number of variables are predefined by the program. The following variables can be used to convert between atomic units and other units:

EV=1.d0/27.2113839d0 HARTREE KELVIN=1.d0/3.157747d5 HARTREE KJOULE=1.d0/2625.500d0 HARTREE KCAL=1.d0/627.5096d0 HARTREE CM=1.d0/219474.63067d0 HARTREE CM-1=1.d0/219474.63067d0 HARTREE HZ=1.d0/6.5796839207d15 HARTREE HERTZ=1.d0/6.5796839207d15 HARTREE ANG=1.d0/0.529177209d0 BOHR ANGSTROM=1.d0/0.529177209d0 BOHR TOEV=27.2113839d0 EV TOK=3.157747d5 K TOKELVIN=3.157747d5 K TOCM=219474.63067d0 CM-1 TOHERTZ=6.5796839207d15 HZ TOHZ=6.5796839207d15 HZ TOKJ=2625.500d0 KJ/MOL TOKJOULE=2625.500d0 KJ/MOL TOKCAL=627.5096d0 KCAL/MOL TOA=0.529177209d0 ANGSTROM TOANG=0.529177209d0 ANGSTROM TODEBYE=2.54158d0 DEBYE

Further variables which are set during execution of the program:

`INTYP`- defines integral program to be used. Either
`INTS`(Seward) or`INTP`(Argos). `INTDONE`- has the value
`.true.`if the integrals are done for the current geometry. `CARTESIAN`- Set to one if Cartesian basis functions are used.
`SCFDONE`- has the value
`.true.`if an SCF calculation has been done for the current geometry. `NUMVAR`- number of variables presently defined
`STATUS`- status of last step (1=no error, -1=error or no convergence)
`CHARGE`- Total charge of the molecule
`NELEC`- number of electrons in last wavefunction
`SPIN`- spin multiplicity minus one of last wavefunction
`ORBITAL`- record of last optimized orbitals (set but never used in the program)
`LASTORB`- Type of last optimized orbitals (
`RHF`,`UHF`,`UHFNAT`, or`MCSCF`. `LASTSYM`- Symmetry of wavefunction for last optimized orbitals.
`LASTSPIN`- for wavefunctions for last optimized orbitals.
`LASTNELEC`- Number of electrons in wavefunction for last optimized orbitals.
`ENERGR(istate)`- Reference energy for state
*istate*in MRCI and CCSD. `ENERGY(istate)`- last computed total energy for state
*istate*for the method specified in the input (e.g.,`HF`,`MULTI`,`CCSD(T)`, or`CCSD[T]`. `ENERGY_METHOD`- String variable holding name of the method
used for calculating
`ENERGY` `ENERGY_BASIS`- String variable holding name of the orbital
basis-set used for calculating
`ENERGY` `GEOMETRY_METHOD`- Equal to the value of
`ENERGY_METHOD`for the most recent geometry optimisation `GEOMETRY_BASIS`- Equal to the value of
`ENERGY_BASIS`for the most recent geometry optimisation `ENERGD(istate)`- Total energy for state
*istate*including Davidson correction (set only in`CI`). `ENERGP(istate)`- Total energy for state
*istate*including Pople correction (set only in`CI`). `ENERGT(1)`- Total energy including perturbative triples
`(T)`correction (set only in`CCSD(T), QCI(T)`). `ENERGT(2)`- Total energy including perturbative triples
`[T]`correction (set only in`CCSD(T), QCI(T)`). `ENERGT(3)`- Total energy including perturbative triples
`-t`correction (set only in`CCSD(T), QCI(T)`). `EMP2`- holds MP2 energy in MPn, CCSD, BCCD, or QCISD calculations, and RS2 energy in MRPT2 (CASPT2) calculations.
`EMP3`- holds MP3 energy in MP3 and MP4 calculations, and RS3 energy in MRPR3 (CASPT3) calculations.
`EMP4`- holds MP4(SDQ) energy in MP4 calculations. The MP4(SDTQ) energy is stored in variable
`ENERGY`. `METHODC`- String variable holding name of the methods used for
`ENERGC`, e.g.,`CCSD, BCCD, QCI`. `METHODT(1)`- String variable holding name of the methods used for
`ENERGT(1)`, e.g.,`CCSD(T), BCCD(T), QCI(T)`. `METHODT(2)`- String variable holding name of the methods used for
`ENERGT(2)`, e.g.,`CCSD[T], BCCD[T], QCI[T]`. `METHODT(3)`- String variable holding name of the methods used for
`ENERGT(3)`, e.g.,`CCSD-T, BCCD-T, QCI-T`. `ENERGC`- Total energy excluding perturbative triples correction (set only in
`QCI`or`CCSD`with triples correction enabled). `DFTFUN`- total value of density functional in
`DFT`or`KS`. `DFTFUNS(ifun)`- value of
`ifun`'th component of density functional in`DFT`or`KS`. `DFTNAME(ifun)`- name of
`ifun`'th component of density functional in`DFT`or`KS`. `DFTFAC(ifun)`- factor multiplying
`ifun`'th component of density functional in`DFT`or`KS`. `DFTEXFAC`- factor multiplying exact exchange in
`KS`. `PROP(istate)`- computed property for state
*istate*. See below for the names`PROP`of various properties. `PROGRAM`- last program called, as specified in the input (e.g.,
`HF`,`CCSD(T)`, etc.) `ITERATIONS`- Number of iterations used. Set negative if no convergence or max number of iterations reached.
`CPUSTEP`- User-CPU time in seconds for last program called.
`SYSSTEP`- System-CPU time in seconds for last program called.
`WALLSTEP`- Elapsed time in seconds for last program called.

The variable names for properties are the same as used on the `EXPEC` input cards.

`OV`- Overlap
`EKIN`- Kinetic energy
`POT`- Potential
`DELTA`- Delta function
`DEL4``DARWIN`- Darwin term of relativistic correction
`MASSV`- Mass-velocity term of relativistic correction
`EREL`- Total relativistic correction
`DMX, DMY, DMZ`- Dipole moments
`XX, YY, ZZ, XY, XZ, XY`- Second moments
`XXX, XXY, XXZ, XYY, XYZ, XZZ, YYY, YYZ, YZZ, ZZZ`- Third moments
`QMXX, QMYY, QMZZ, QMXY, QMXZ, QMXY`- Quadrupole moments
`EFX, EFY, EFZ`- Electric field
`FGXX, FGYY, FGZZ, FGXY, FGXZ, FGXY`- Electric field gradients
`D/DX, D/DY, D/DZ`- Velocity
`LSX, LSY, LSZ`- One-electron spin-orbit
`LL`- Total angular momentum squared
`LX, LY, LZ`- Electronic angular momentum
`LXLX, LYLY, LZLZ, LXLY, LXLZ, LYLZ`- Two-electron angular momentum

By default, only the dipole moments are computed and defined. The values of other
properties are only
stored in variables if they are requested by `EXPEC` cards. If more than
one state is computed (e.g., in state-averaged MCSCF, corresponding arrays
`PROP(istate)` are returned. If properties are computed for more than
one center, the center number is appended to the name, e.g. `EFX1`, `EFX2` etc.

If transition properties are computed, their values are stored in corresponding
variables with prefix `TR`, e.g., `TRDMX, TRDMY, TRDMZ` for transition
dipole moments. If more than two states are computed, the index is
, where are state numbers. In a state-averaged
calculation, states are counted sequentially for all state symmetries.

For instance, in the following state-averaged MCSCF

`MULTI;WF,14,1,0;STATE,3;WF,14,2,0;STATE,2;WF,14,3,0`

the states are counted as

1 | 2 | 3 | 4 | 5 | 6 | ||

Symmetry | 1 | 1 | 1 | 2 | 2 | 3 | |

Root in Sym. | 1 | 2 | 3 | 1 | 2 | 1 |

molpro@molpro.net 2019-08-22