13.11.1 Theory

XPol is an approximate, fragment-based molecular orbital method that was developed as a “next-generation” force field.^{285, 1026, 1028, 1027} The basic idea of the method is to treat a molecular liquid, solid, or cluster as a collection of fragments, where each fragment is a molecule. Intra-molecular interactions are treated with a self-consistent field method (Hartree-Fock or DFT), but each fragment is embedded in a field of point charges that represent electrostatic interactions with the other fragments. These charges are updated self-consistently by collapsing each fragment’s electron density onto a set of atom-centered point charges, using charge analysis procedures (Mulliken, Löwdin, or ChElPG, for example; see Section 11.2.1). This approach incorporates many-body polarization, at a cost that scales linearly with the number of fragments, but neglects the anti-symmetry requirement of the total electronic wave function. As a result, intermolecular exchange-repulsion is neglected, as is dispersion since the latter is an electron correlation effect. As such, the XPol treatment of polarization must be augmented with empirical, Lennard–Jones-type intermolecular potentials in order to obtain meaningful optimized geometries, vibrational frequencies or dynamics.

The XPol method is based upon an ansatz in which the super-system wave function is written as a direct product of fragment wave functions,

$$|\mathrm{\Psi }⟩=\prod _A^{N_{\mathrm{frag}}}|{\mathrm{\Psi }}_A⟩,$$

(13.29)

where $N_{\mathrm{frag}}$ is the number of fragments. We assume here that the fragments are molecules and that covalent bonds remain intact. The fragment wave functions are anti-symmetric with respect to exchange of electrons within a fragment, but not to exchange between fragments. For closed-shell fragments described by Hartree-Fock theory, the XPol total energy is^{1028, 417}

$$E_{\mathrm{XPol}}=\sum _A\left[2\sum _a{𝐜}_a^{†}\left({𝐡}^A+{𝐉}^A-\frac{1}{2}{𝐊}^A\right){𝐜}_a+E_{\mathrm{nuc}}^A\right]+E_{\mathrm{embed}}.$$

(13.30)

The term in square brackets is the ordinary Hartree-Fock energy expression for fragment $A$. Thus, ${𝐜}_a$ is a vector of occupied MO expansion coefficients (in the AO basis) for the occupied MO $a\in A$; ${𝐡}^A$ consists of the one-electron integrals; and ${𝐉}^A$ and ${𝐊}^A$ are the Coulomb and exchange matrices, respectively, constructed from the density matrix for fragment $A$. The additional terms in Eq. (13.30),

$$E_{\mathrm{embed}}=\frac{1}{2}\sum _A\sum _{B\ne A}\sum _{J\in B}\left(-2\sum _a{𝐜}_a^{†}{𝐈}_J{𝐜}_a+\sum _{I\in A}{L}_{IJ}\right)q_J,$$

(13.31)

arise from the electrostatic embedding. The matrix ${𝐈}_J$ is defined by its AO matrix elements,

$${\left({𝐈}_J\right)}_{\mu \nu }=⟨\mu \left|\frac{1}{\left|\overrightarrow{r}-{\overrightarrow{R}}_J\right|}\right|\nu ⟩,$$

(13.32)

and $L_{IJ}$ is given by

$$L_{IJ}=\frac{Z_I}{\left|{\overrightarrow{R}}_I-{\overrightarrow{R}}_J\right|}.$$

(13.33)

According to Eqs. (13.30) and (13.31), each fragment is embedded in the electrostatic potential arising from a set of point charges, $\{q_J\}$, on all of the other fragments; the factor of $1/2$ in Eq. (13.31) avoids double-counting. Exchange interactions between fragments are ignored, and the electrostatic interactions between fragments are approximated by interactions between the charge density of one fragment and point charges on the other fragments.

Crucially, the vectors ${𝐜}_a$ are constructed within the ALMO ansatz,⁴⁶⁷ so that MOs for each fragment are represented in terms of only those AOs that are centered on atoms in the same fragment. This choice affords a method whose cost grows linearly with respect to $N_{\mathrm{frag}}$, and where basis set superposition error is excluded by construction. In compact basis sets, the ALMO ansatz excludes inter- fragment charge transfer as well.

The original XPol method of Xie et al.^{1026, 1028, 1027} uses Mulliken charges for the embedding charges $q_J$ in Eq. (13.31), though other charge schemes could be envisaged. In non-minimal basis sets, the use of Mulliken charges is beset by severe convergence problems,⁴¹⁷ and Q-Chem’s implementation of XPol offers the alternative of using either Löwdin charges or ChElPG charges,¹²⁰ the latter being derived from the electrostatic potential as discussed in Section 11.2.1. The ChElPG charges are found to be stable and robust, albeit with a somewhat larger computational cost as compared to Mulliken or Löwdin charges.^{417, 372} An algorithm to compute ChElPG charges using atom-centered Lebedev grids rather than traditional Cartesian grids is available (see Section 11.2.1),³⁹⁰ which uses far fewer grid points and thus can significantly improve the performance for the XPol/ChElPG method, where these charges must be iteratively updated.

Researchers who use Q-Chem’s XPol code are asked to cite Refs. 417, 372.