The `POLY` program allows for the transformation of the potential energy surface
and property surfaces from a grid representation to a polynomial representation. Once
a polynomial representation has been chosen, the corresponding VSCF, VCI or VMP2
programs need to be selected (see below). The `POLY` program is fully parallelized
in terms of the `MPPX` scheme. An example for the use of the `POLY` program
can be found in chapter 55.1.3.

The following *options* are available:

`NDIM`=*n*- The keyword
`NDIM=`terminates the transformation after the D terms within the -mode expansion of the surfaces. The default is set to 3. The transformation of the 4D terms can be very time consuming. `NDIMDIP`=*n*- Term after which the -body expansions of the dipole surfaces are truncated. The default is set to 3.
Note that
`NDIMDIP`has to be lower or equal than`NDIM`. `NDIMPOL`=*n*- Term after which the -body expansions of the polarizability tensor surfaces are truncated. The default is 0.
`NDIMPOL`has to be lower or equal than`NDIM`and must be smaller than 4. `NGRID`=*n*- Once the value of the
`NGRID=`keyword for controlling the number of grid points has been changed in the`SURF`program, this information will be passed to the`POLY`program automatically. `AUTOPOL`=*n*`AUTOPOL=1`(default) determines the optimal parameters for`PMAX`and`PSUM`in an iterative procedure. In combination with`AUTOPOL=1`the`PMAX`and`PSUM`keywords specify the maximum values within the automated determination of these parameters.`PMAX`=*n*- This keyword controls the order of the polynomial to be determined in the fitting
procedure. The default is set to 8. The definition of
this keyword changes in combination with the
`AUTOPOL`keyword (see above). `PSUM`=*n*- In multidimensional fitting the maximum order of the polynomial is not given
by
`PMAX`times the dimension, but rather by`PSUM`. The default is set to 24. This reduces the computational cost for fitting tremendously. However, for accurate calculations`PMAX`=12 and`PSUM`=30 are recommended. The definition of this keyword changes in combination with`AUTOPOL`keyword (see above). `TYPE`=*variable*- Once polynomials have been generated, the expansion of the potential can be restricted
to the harmonic approximation or a semi-quartic force field, i.e.
`TYPE`=*HARM*or`TYPE=`=*QFF*. `DIPOLE`=*n*- In case that dipole surfaces have been computed, they need to be transformed
to a polynomial representation as well. This can be accomplished by
`DIPOLE=1`and will be set by default, if dipole surfaces are available.`DIPOLE=0`excludes existing dipole surfaces from the transformation. `POLAR`=*n*- Once set to one this keyword switches the transformation of the polarizability tensor surfaces on.
The default is
`POLAR`=0. `VAM`=*n*- The Watson correction term, i.e.
`VAM=1`, is absorbed in the potential by default. Once the polynomials to be generated should exclude this correction, this needs to be specified by`VAM=0`. Any other values, i.e. 2 or 3, will be ignored. `DELAUTO`=*value*- In order to delete troublesome 2D, 3D or 4D surfaces from the multi-mode expansion of
the potential, the
`DELAUTO`keyword sets all coefficients of the polynomial expansion of these surfaces to zero. The threshold passed to the`DELAUTO`keyword corresponds to the value of the least squares fitting procedure. It acts on both, the energy and the dipole surfaces. The default is set to 9.d99, i.e. all surfaces with values below this threshold will not be affected. `SOLVER`=*n*- By default,
`SOLVER=1`, an LU decomposition will be used within the fitting procedure. Alternatively, a Cholesky decomposition is available, which sometimes leads to more stable solutions (`SOLVER=2`). `NVARC`=*n*- For the calculation of vibrationally averaged rotational constants
-tensor surfaces will be generated and transformed to polynomials up to 2nd order (default).
The order can be changed by the integer passed to the
`NVARC`keyword. This transformation can also be extended within the`VSCF`program. `VRC`=*n*- Once vibrational-rotational coupling surfaces have been computed in the
`SURF`program, these couplings can be considered (`VRC=1`) or excluded (`VRC=0`, default) in the`POLY`program. `SHOW`=*n*- The coefficients of the polynomial representation can be printed. In order to
identify quartic potentials, it is recommended to use
`SHOW=1`. Higher values will lead to a very long output file. `START`=*record*- This keyword controls the record, from which the potential shall be read. The default is 5600.2. This keyword usually is only needed within the calculation of vibronic spectra.
`SAVE`=*record*- Specification of the record, where to dump the polynomials. The default is 5750.2. This keyword usually is only needed within the calculation of vibronic spectra.
`INFO`=*n*`INFO=1`provides a list of the values of all relevant program parameters (options).