Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Bauman, Nicholas P. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000339626265); Pathak, Himadri [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000259198002); Liebenthal, Marcus D. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)]; Panyala, Ajay R. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000208463347); Mejia-Rodriguez, Daniel [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000203502941); Govind, Niranjan [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000033625366X); Kowalski, Karol [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000016357785X)

doi:10.1021/acs.jctc.5c01599

DOE OSTI · 3015982

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Bauman, Nicholas P. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000339626265)·Pathak, Himadri [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000259198002)·Liebenthal, Marcus D. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)]·Panyala, Ajay R. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000208463347)·Mejia-Rodriguez, Daniel [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000203502941)·Govind, Niranjan [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000033625366X)·Kowalski, Karol [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000016357785X)

View original source ↗Read source document ↗Export metadata ↓

Abstract

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Explore related subjects

Basis sets Cavities Chemical calculations Chemical reactions Transfer reactions

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Bauman, Nicholas P. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000339626265), Pathak, Himadri [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000259198002), Liebenthal, Marcus D. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)], Panyala, Ajay R. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000208463347), Mejia-Rodriguez, Daniel [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000203502941), Govind, Niranjan [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000033625366X), Kowalski, Karol [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Univ. of Washington, Seattle, WA (United States)] (ORCID:000000016357785X). 2025-12-15. Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems. https://doi.org/10.1021/acs.jctc.5c01599

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗

DOE OSTI2025

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗

DOE OSTI2025

Implementation of McMurchie–Davidson Algorithm for Gaussian AO Integrals Suited for SIMD Processors

We report an implementation of the McMurchie− Davidson evaluation scheme for 1- and 2-particle Gaussian AO integrals designed for processors with Single Instruction Multiple Data (SIMD) instruction sets. Like in our recent MD implementation for graphical processing units (GPUs) [Asadchev, A.; Valeev, E. F.. J. Chem. Phys. 2024, 160, 244109.], variable-sized batches of shellsets of integrals are evaluated at a time. By optimizing for the floating point instruction throughput rather than minimizing the number of operations, this approach achieves up to 50% of the theoretical hardware peak FP64 performance for many common SIMD-equipped platforms (AVX2, AVX512, NEON), which translates to speedups of up to 30 over the state-of-the-art one-shellset-at-a-time implementation of Obara−Saika-type schemes in Libint for a variety of primitive and contracted integrals. As with our previous work, we rely on the standard C++ programming language such as the std::simd standard library feature to be included in the 2026 ISO C++ standard without any explicit code generation to keep the code base small and portable. The implementation is part of the open source LibintX library freely available at https://github.com/ValeevGroup/libintx.

Basis sets↗