Reliable Structures and Electronic Energies of Small Water Clusters using Density Functional and Local Correlation Coupled Cluster Model Chemistries

Author

Benjamin T. Petty
Vance R. Fowler
Audrey Ryu
Caroline S. Glick
Carly A. Rock
Qihang Wang
Gregory S. Tschumper, Missouri University of Science and TechnologyFollow
George C. Shields

Abstract

In this paper we have assessed the ability of the domain-based local pair natural orbital (DLPNO)-CCSD(T) method to match the explicitly correlated CCSD(T) relative energies of (H₂O)n=3–7isomers along with the impact of the level of theory utilized to optimize the water cluster geometries. The benchmark structures were optimized using a 2-body:Many-body procedure in which all of the 1- and 2-body contributions are computed using CCSD(T) while all of the higher-order many-body interactions are computed using MP2 (denoted CCSD(T):MP2). Benchmark relative energies were computed for these CCSD(T):MP2 optimized geometries with explicitly correlated CCSD(T)-F12b single point energies (SPEs) using the cc-pVQZ-F12 and cc-pV5Z-F12 basis sets augmented with diffuse functions on the O atoms. The benchmark structures and energies were used to gauge the performance of less demanding computational protocols. For example, DLPNO–CCSD(T) computations on the 31 benchmark structures with the analogous family of correlation consistent basis sets (cc-pVNZ for H and aug-cc-pVNZ for O, or simply haNZ where N = D-6) were used to estimate relative energies at the complete basis set (CBS) limit via three-point extrapolations. When compared to the CCSD(T)-F12 benchmark data, the mean absolute differences (MADs) were ≤ 0.13 kcal/mol when triple-ζ and larger basis sets were employed. Using these DLPNO–CCSD(T) results, we demonstrate that 2 less-demanding geometry optimization procedures, specifically the ωB97X-D density functional theory (DFT) method paired with the 6–31++G(d,p) basis set and the density-fitted MP2 method paired with the haTZ basis set, give structures that yield nearly identical relative energies (MADs of only 0.07 and 0.02 kcal/mol, respectively, when comparing DLPNO–CCSD(T)/ha6Z data). In addition, we show how the presence or absence of diffuse functions in the basis sets used for DLPNO–CCSD(T) SPEs impact the quality of the relative energies. The protocol that combines ωB97X-D/6–31++G(d,p) optimized structures with DLPNO–CCSD(T) SPEs using triple-ζ or higher Dunning basis sets that include augmentation with diffuse functions on the oxygen atoms provides a fast and accurate method for determining the relative electronic energies of (H₂O)n=3–7water cluster isomers.

Recommended Citation

B. T. Petty et al., "Reliable Structures and Electronic Energies of Small Water Clusters using Density Functional and Local Correlation Coupled Cluster Model Chemistries," Journal of Physical Chemistry A, vol. 129, no. 40, pp. 9291 - 9302, American Chemical Society, Oct 2025.

The definitive version is available at https://doi.org/10.1021/acs.jpca.5c04923

Department(s)

Chemistry

Publication Status

Open Access

Comments

Gauss Centre for Supercomputing, Grant CHE-2452726

International Standard Serial Number (ISSN)

1520-5215; 1089-5639

Document Type

Article - Journal

Document Version

Citation

File Type

text

Language(s)

English

Rights

© 2025 American Chemical Society, All rights reserved.

Creative Commons Licensing

This work is licensed under a Creative Commons Attribution 4.0 License.

Publication Date

09 Oct 2025

PubMed ID

40994268