In order to obtain physical results, one must either supplement the XPol energy expression with either empirical intermolecular potentials or else with an ab initio treatment of intermolecular interactions. The latter approach is described in Section 13.13. Here, we describe how to add Lennard-Jones or Buckingham potentials to the XPol energy, using the $xpol_mm and $xpol_params sections described below.
The Lennard-Jones potential is
where represents the distance between atoms and . This potential is characterized by two parameters, a well depth and a length scale . Although quite common, the repulsion is unrealistically steep. The Buckingham potential replaces this with an exponential function,
Here, , , and are additional (dimensionless) constants, independent of atom type. In both Eq. (13.34) and Eq. (13.35), the parameters and are determined using the geometric mean of atomic well-depth and length-scale parameters. For example,
The atomic parameters and must be specified using a $xpol_mm section in the Q-Chem input file. The format is a molecular mechanics-like specification of atom types and connectivities. All atoms specified in the $molecule section must also be specified in the $xpol_mm section. Each line must contain an atom number, atomic symbol, Cartesian coordinates, integer atom type, and any connectivity data. The $xpol_params section specifies, for each atom type, a value for in kcal/mol and a value for in Ångstroms. A Lennard-Jones potential is used by default; if a Buckingham potential is desired, then the first line of the $xpol_params section should contain the string BUCKINGHAM followed by values for the , , and parameters.