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:
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.
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.
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.
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.