---
title: q-AQUA-pol Potential Energy Surface
url: https://www.emergentmind.com/topics/q-aqua-pol-potential-energy-surface
type: topic
---

# q-AQUA-pol Potential Energy Surface

Searching arXiv for the specified papers and closely related q-AQUA/q-AQUA-pol work.
arXiv search query: "q-AQUA-pol water potential energy surface"
The q-AQUA-pol potential energy surface (PES) is a many-body water PES built from CCSD(T) reference data up to the 4-body level and used, in recent quantum vibrational work, in conjunction with an n-mode representation to enable grid-based quadrature, compressed integral storage, and selected-configuration-interaction treatments of water-cluster spectra [2204.01804; 2509.13517]. In the formulation used for cluster spectroscopy, q-AQUA-pol is expressed as a $\Delta$-correction to the classical polarizable TTM-3F model, while the underlying many-body representation resolves the total energy into 1-, 2-, 3-, and 4-body contributions, with higher-order terms treated as negligible at the CCSD(T) level [2204.01804].

## 1. Formal many-body structure

For an $N$-water system, the total potential energy in the q-AQUA framework is written as
\[
V(1,\dots,N)
=\sum_{i=1}^N V_{1\mathrm b}(i)
+\sum_{i<j} V_{2\mathrm b}(i,j)
+\sum_{i<j<k} V_{3\mathrm b}(i,j,k)
+\sum_{i<j<k<l} V_{4\mathrm b}(i,j,k,l)\,.
\]

In the q-AQUA-pol form used for the exact mid-IR vibrational calculations, the PES is built as a $\Delta$-correction to TTM-3F:
\[
V_{\mathrm{qAQUA}}(\mathbf R)
=
V_{\mathrm{TTM\text{--}3F}}(\mathbf R)
+\sum_{i<j}\Delta V_{2\mathrm B}(\mathbf r_i,\mathbf r_j)
+\sum_{i<j<k}\Delta V_{3\mathrm B}(\mathbf r_i,\mathbf r_j,\mathbf r_k)
+\sum_{i<j<k<\ell}\Delta V_{4\mathrm B}(\mathbf r_i,\mathbf r_j,\mathbf r_k,\mathbf r_\ell)\,,
\]
where $\mathbf R=(\mathbf r_1,\dots,\mathbf r_N)$ are all atomic coordinates [2509.13517].

The TTM-3F reference contains intramolecular one-body monomer potentials, permanent electrostatics with point charges placed on each H and O, induced dipole interactions with Thole damping via site polarizabilities $\alpha_O$ and $\alpha_H$, and damped dispersion and exchange-repulsion at long range. The $\Delta n\mathrm B$ terms correct the TTM-3F $n$-body energies and are constructed so that each correction goes to zero when any subset of molecules is driven apart [2509.13517].

This decomposition is significant because it combines a physically motivated polarizable baseline with explicit many-body corrections fit to high-level electronic-structure data. A plausible implication is that the model inherits long-range physical structure from TTM-3F while using explicit $\Delta$ terms to recover the CCSD(T) many-body energetics required for cluster spectroscopy and condensed-phase simulation.

## 2. Parameterization and analytic components

The one-body monomer term is taken from the Partridge–Schwenke spectroscopic monomer PES. In the q-AQUA construction it is described as analytic, full-dimensional, and permutationally invariant with respect to the two H atoms [2204.01804].

The higher-body terms are represented by permutationally invariant polynomials (PIPs) and switching constructions that separate short- and long-range physics. The essential components are as follows.

| Term | Representation | Reference data / fit quality |
|---|---|---|
| $V_{1\mathrm b}$ | Partridge–Schwenke monomer PES | Spectroscopic monomer PES |
| $V_{2\mathrm b}$ | 7th-order PIP of 42 symmetry functions + $V_{\rm dd}$ with switching | 71 892 dimer geometries; typical RMS error 25 cm$^{-1}$ |
| $V_{3\mathrm b}$ | 4th-order short-range PIP + 3rd-order longer-range PIP with switching | 45 332 trimer configurations; RMS $\approx 9$ cm$^{-1}$ and $\approx 11$ cm$^{-1}$ |
| $V_{4\mathrm b}$ | Improved CCSD(T) PIP with 200 grouped PIPs plus selected 4th-order terms | 3 692 tetramer configurations; RMS $\approx 7.2$ cm$^{-1}$ |

For the 2-body term, the short-range part is fit by a 7th-order PIP in the six inter- and intra-monomer distances, while the long-range part is described by a high-level dipole-dipole interaction $V_{\rm dd}(i,j)$ built on an accurate CCSD(T) dipole-moment surface. A smooth switching function $s_2(R_{ij})$ blends the two in a window $R_{ij}\in[R_{\rm on},R_{\rm off}]$:
\[
V_{2\mathrm b}(i,j)
=
s_2(R_{ij})\,V^{\rm PIP}_{2\mathrm b}(i,j)
+[1-s_2(R_{ij})]\,V_{\rm dd}(i,j)\,.
\]

For the 3-body term, two complementary PIP fits are used: a short-range fit for $\max R_{OO}<7$ Å in Morse-type variables of OO, OH, and HH distances with 222111 symmetry, and a longer-range fit for $\max R_{OO}\in[5,9.5]$ Å in inverse distances with the same symmetry; the two are joined by a second switching function $s_3(\max R_{OO})$ that ensures analytic continuity. The 4-body term is based on an improved version of a CCSD(T) PIP, fitted with a compact basis of 200 grouped PIPs plus a selected set of 4th-order terms, and multiplied by a final switch $s_4(R_{\rm max})$ at large tetramer size [2204.01804].

The fitting protocol uses PIP bases that enforce permutation invariance of monomers and H atoms. “Purification” removes basis functions that are exactly redundant under the molecular symmetry group, and “pruning” discards high-order monomials that contribute negligibly to the fit or lead to ill-conditioning. Coefficients are obtained by linear least squares, analytic gradients are generated efficiently via reverse-mode differentiation, and no empirical regularization is used; fit quality is judged by RMS errors on large test sets [2204.01804].

## 3. n-mode representation and quadrature realization

In the vibrational calculations, the q-AQUA-pol PES is not used directly on a full-dimensional grid. Instead, it is expanded in a normal-mode coordinate basis $\{q_i\}$ by an n-mode representation truncated at $n=3$:
\[
V(q_1,\dots,q_N)
\simeq
V_0
+\sum_i V^{(i)}(q_i)
+\sum_{i<j} V^{(ij)}(q_i,q_j)
+\sum_{i<j<k} V^{(ijk)}(q_i,q_j,q_k)\,.
\]

More formally, if $U\subseteq\{1,\dots,N\}$ with $|U|\le 3$, then
\[
V(q)\simeq \sum_{U:|U|\le 3} V^{(U)}(\{q_i\}_{i\in U})\,,
\]
with each higher-order contribution defined by subtraction of all lower-order terms. For example,
\[
V^{(ij)}(q_i,q_j)
=
V(q_i,q_j;\text{ all other }q=0)
-
V^{(i)}(q_i)
-
V^{(j)}(q_j)
-
V_0\,.
\]

This truncation avoids the full-dimensional grid while still capturing anharmonic coupling up to three modes simultaneously. The matrix elements required for the vibrational Hamiltonian,
\[
\langle \phi_{n_i}(q_i)\phi_{n_j}(q_j)\dots|V^{(U)}(q_U)|\phi_{m_i}(q_i)\phi_{m_j}(q_j)\dots\rangle,
\]
factorize into products of one-, two-, or three-dimensional integrals because of the n-mode form [2509.13517].

One-dimensional integrals are evaluated by Gauss–Hermite quadrature,
\[
\int_{-\infty}^{\infty}\phi_n(q)\phi_m(q)V^{(i)}(q)\,dq
\approx
\sum_{\alpha=1}^{N_q} w_\alpha \phi_n(q_\alpha)\phi_m(q_\alpha)V^{(i)}(q_\alpha),
\]
and two-mode and three-mode terms are treated by corresponding product quadratures. Within selected CI, the two-mode integrals
\[
V^{(ij)}_{m_i m_j\,n_i n_j}
=
\langle \phi_{m_i}(q_i)\phi_{m_j}(q_j)|V^{(ij)}|\phi_{n_i}(q_i)\phi_{n_j}(q_j)\rangle
\]
are precomputed once, sorted in descending magnitude, and stored in compressed form. This “heat-bath” ordering allows the HCI selection step to test excitations in order of largest coupling rather than scanning all possible combinations; one-mode integrals are stored analogously [2509.13517].

## 4. Benchmark performance from clusters to liquid water

The original q-AQUA validation emphasizes both cluster benchmarks and condensed-phase observables. For water hexamer isomers, the dissociation-energy mean absolute error versus CCSD(T)/CBS benchmarks is under $0.3$ kcal/mol, and harmonic-frequency MAEs for Prism, Cage, Book, and Ring are $\sim 10$–$15$ cm$^{-1}$, with q-AQUA generally matching or slightly improving over MB-pol [2204.01804].

Unrestricted diffusion Monte Carlo (DMC) zero-point energies relative to the Prism equilibrium were reported for Prism, Cage, and Book1 as $32\,647\pm 9$ cm$^{-1}$, $32\,553\pm 19$ cm$^{-1}$, and $32\,652\pm 12$ cm$^{-1}$, respectively. Removing the 4-body term shifts these by up to $\sim 75$ cm$^{-1}$, which was used to confirm the physical importance of the 4-body interaction. For $(\mathrm H_2\mathrm O)_{20}$ isomers A3, A2d, and 9, q-AQUA gives binding energies of $-199.8$, $-201.7$, and $-200.5$ kcal/mol, within $1$ kcal of the stated MP2/CBS benchmarks, whereas MB-pol underbinds by $\approx 6$ kcal. Preliminary DMC on the 20-mer shows no “holes” and stable propagation [2204.01804].

Condensed-phase validation was carried out with classical MD, PIMD, and RPMD in i-PI for 256 waters under periodic boundary conditions. At 298 K, a model truncated at 1- and 2-body gives severe over-structuring and peak misplacement in the radial distribution function $g_{OO}(r)$, inclusion up to 3-body gives correct peak positions but amplitudes that are too large, and the full 4-body model with PIMD yields nearly quantitative agreement with the x-ray data of Skinner et al. For self-diffusion, the reported values at 298 K are $0.145\pm 0.012$ Å$^2$/ps for classical MD, $0.226\pm 0.020$ Å$^2$/ps for RPMD, and $0.230$ Å$^2$/ps for experiment; nuclear quantum effects increase $D$ by roughly $50\%$ at 298 K. The orientational relaxation time $\tau_2$ at 298 K is $3.2$ ps classically and $2.4$ ps in RPMD, compared with an experimental $2.5$ ps. Turning the 4-body term off raises $D$ by $\sim 5$–$10\%$ [2204.01804].

## 5. Role in exact mid-IR quantum vibrational spectroscopy

The q-AQUA-pol PES is the central interaction model in selected-configuration-interaction calculations of the zero-temperature mid-infrared ($2800$–$3800$ cm$^{-1}$) vibrational spectra of the water monomer, dimer, trimer, and hexamer in its cage and prism geometries [2509.13517]. In that setting, the PES is combined with the n-mode representation specifically to facilitate grid-based quadrature and integral storage.

Within selected configuration interaction, a new spectral strategy is introduced that is complementary to eigenstate enumeration. Instead of enumerating excited states directly, the spectrum is calculated with the response-vector method, and the resulting system of linear equations is solved using a basis of configurations that is optimally selected at each frequency of interest. The pre-sorted one- and two-mode integrals inherited from the q-AQUA-pol/n-mode representation are reused in this frequency-resolved HCI procedure [2509.13517].

The reported computational savings are explicit: the spectral HCI approach yields sizable savings in both storage, by a factor of $\sim 4$ in the trimer, and CPU time, with minutes per frequency point on a single node for $N=9$. The same study compares the resulting spectra to previous work and highlights limitations of the local monomer approximation. It also states that, to the best of its authors’ knowledge, the hexamer spectra are the most accurate ones reported to date [2509.13517].

## 6. Accuracy bounds, limitations, and computational significance

Several distinct accuracy statements are attached to q-AQUA-pol as used in vibrational applications. The underlying four-body-corrected fit reproduces benchmark CCSD(T) harmonic frequencies of the water hexamer to within $\approx 5$ cm$^{-1}$ on average, as cited to Howard and Tschumper. In the present vibrational calculations, after freezing low-frequency modes and truncating at three-mode coupling, absolute fundamental frequencies are converged to better than $10$–$20$ cm$^{-1}$ compared to previous VCI on the HBB PES, with systematic behavior [2509.13517].

The main limitations stated for the vibrational implementation are associated with the n-mode truncation. The three-mode truncation can introduce spurious wells or unbound behavior in very floppy, low-frequency modes, especially torsions. To avoid this, all modes below $\approx 400$ cm$^{-1}$ were frozen in the reported work. As with all many-body expansions, $\Delta V_{5\mathrm B}$ and higher effects are neglected, although for water they are described as very small beyond four bodies [2509.13517].

The computational advantages are equally specific. The n-mode expansion breaks a $3N$-dimensional integral into sums of at most three-dimensional integrals, making full grid quadrature tractable up to $N\approx 18$ modes. Using Gauss–Hermite quadrature and storing only one-, two-, and three-mode integrals minimizes memory and flops, while heat-bath ordering makes variational-space construction efficient. In the original q-AQUA implementation for periodic simulations, the total energy and gradient evaluation for 256 waters up to 4-body with periodic boundary conditions is $\sim 4.5$ s per step on a single 2.4 GHz Intel Xeon core, and with 8-core OpenMP this drops by a factor of $\approx 6.5$; the 4-body work is the most expensive, at $\sim 115\,000$ tetramers evaluated per step, and truncating at the 3-body level roughly halves the cost [2204.01804].

Taken together, the q-AQUA-pol many-body fit, the n-mode/Gauss–Hermite representation, and the heat-bath integral strategy make fully quantum, near-exact mid-IR spectra of water clusters up to the hexamer feasible and reliable at the $5$–$20$ cm$^{-1}$ level, while the underlying q-AQUA model also reproduces cluster energetics, liquid structure, and diffusion with accuracy that remains anchored to CCSD(T) reference data [2509.13517; 2204.01804].

Source: https://www.emergentmind.com/topics/q-aqua-pol-potential-energy-surface