Basis-Set Limit Wavefunction Energies of Massive Water Clusters: Resolution of a Challenge for Many-Body Expansion
5 affiliations
- 1School of Chemistry
- 2University of New South Wales
- 3School of Computing
- 4Australian National University
- 5Monash Institute of Pharmaceutical Sciences
Abstract
The many-body expansion (MBE) offers a route to high-level wavefunction energies of systems beyond the reach of conventional calculations, but its practical convergence is entangled with basis set superposition error (BSSE), electrostatic embedding, and the scarcity of direct supersystem benchmarks at large size. Here, we show that as the complete-basis-set (CBS) limit is approached, electrostatically embedded three-body MBE, EE-MBE(3), provides accurate absolute and relative HF, RI-MP2, and CCSD(T) energies across water clusters of varying sizes. Critically, the largesystem conclusions are tested against direct same-level supersystem references rather than lowerlevel or extrapolated benchmarks. A distributed canonical RI-MP2 implementation enables a systematic HF/RI-MP2 campaign for 64 and 128 waters with cc-pVXZ and aug-cc-pVXZ basis sets (X = D, T, Q), with three configurations at each cluster size/basis set point. The largest calculation, (H2O)128/aug-cc-pVQZ, contains 22,016 spherical primary basis functions and, to our knowledge, is the largest reported MP2 calculation performed without a local-correlation approximation. At augcc-pVQZ, EE-MBE(3) maintains per-water errors within the 0.09 kcal mol-1 dynamical-accuracy
Related
- Band structure of two-dimensional graphene and of one-dimensional nanoribbons via variational quantum algorithms: from energies to wavefunction quality
- High‐Resolution Angle‐Resolved Photoelectron Spectroscopy of Cu5− Clusters
- Static correlation, a coupled-cluster diagnostic blind spot, and its quantum-computing resolution in silver–dinitrogen clusters
- Full complex wavefunction reconstruction with direct-interference ptychography
- An Algorithm for computing Three- and Four-Body Density Matrices from CASCI Wavefunction and AI-Assisted Code Optimization
BibTeX
@misc{jiang2026basisset,
title = {{Basis-Set Limit Wavefunction Energies of Massive Water Clusters: Resolution of a Challenge for Many-Body Expansion}},
author = {Yuhong Jiang and Jiarui Wang and Giuseppe M. J. Barca and Junming Ho},
howpublished = {ChemRxiv},
year = {2026},
month = oct,
doi = {10.26434/chemrxiv.15009940/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15009940/v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Basis-Set Limit Wavefunction Energies of Massive Water Clusters: Resolution of a Challenge for Many-Body Expansion
AU - Yuhong Jiang
AU - Jiarui Wang
AU - Giuseppe M. J. Barca
AU - Junming Ho
PY - 2026
DA - 2026/10/05
PB - ChemRxiv
DO - 10.26434/chemrxiv.15009940/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15009940/v1
AB - The many-body expansion (MBE) offers a route to high-level wavefunction energies of systems beyond the reach of conventional calculations, but its practical convergence is entangled with basis set superposition error (BSSE), electrostatic embedding, and the scarcity of direct supersystem benchmarks at large size. Here, we show that as the complete-basis-set (CBS) limit is approached, electrostatically embedded three-body MBE, EE-MBE(3), provides accurate absolute and relative HF, RI-MP2, and CCSD(T) energies across water clusters of varying sizes. Critically, the largesystem conclusions are tested against direct same-level supersystem references rather than lowerlevel or extrapolated benchmarks. A distributed canonical RI-MP2 implementation enables a systematic HF/RI-MP2 campaign for 64 and 128 waters with cc-pVXZ and aug-cc-pVXZ basis sets (X = D, T, Q), with three configurations at each cluster size/basis set point. The largest calculation, (H2O)128/aug-cc-pVQZ, contains 22,016 spherical primary basis functions and, to our knowledge, is the largest reported MP2 calculation performed without a local-correlation approximation. At augcc-pVQZ, EE-MBE(3) maintains per-water errors within the 0.09 kcal mol-1 dynamical-accuracy
ER -