The calculation of core-polarization matrix elements is invoked by the CPP card, which can be called at an arbitrary position in the MOLPRO input, provided the usual AO integrals have been calculated before. The CPP card can have the following formats:

- CPP,INIT,
*ncentres*; - CPP,ADD[,
*factor*]; - CPP,AE[,
*record.file*]; - CPP,SET[,
*fcpp*];

CPP,INIT,;

abs() further cards will be read in the following format:

;

specifies an atomic centre with polarizable core,

fixes the form of the cutoff-function (choose 1 for Fuentealba/Stoll and 2 for Mueller/Meyer);

is the static dipole polarizability,

is the static quadrupole polarizability,

is the first non-adiabatic correction to the dipole polarizability,

is the exponential parameter of the cutoff-function, and

is the effective charge of the core with which it polarizes its surroundings
(only needed if it contributes orbitals to the CPP core (see below)).

When is lower than zero, only CPP integrals are calculated and saved in record 1490.1. Otherwise, the matrix (records 1200.1 and 1210.1) and the two-electron-integrals (record 1300.1) will be modified.

CPP,ADD,;

With this variant, previously calculated matrix elements of the polarization matrix can be added with the variable factor (default: = 1) to the -matrix as well as to the two-electron integrals. In particular, CPP,ADD,-1.; can be used to retrieve the integrals without the polarization contribution.

CPP,AE,;

This is equivalent to the CPP,ADD command (with factor 1), but
additionally sets a flag that CPP integrals projected with respect to
core orbitals should be used. The projection is necessary whenever
explicitly treated polarizable cores are present.
The core orbitals are taken from *record.file*,
and the CPP,AE command should be followed by a *core* card
specifying the numbers of these orbitals in the various irreducible representations.
Currently, the code assumes that the polarizable cores are
spherically symmetric.

CPP,SET,;

normally not necessary but may be used to tell MOLPRO after a restart, with what factor the polarization integrals are effective at the moment. Currently the CPP integrals are restricted to basis functions up to and including angular momentum 4, i.e. g functions.

molpro@molpro.net 2019-03-19