Machine learning the two-electron reduced density matrix in molecules and condensed phases

TL;DR

This paper demonstrates that machine learning models can accurately predict the two-electron reduced density matrix (2-RDM) in molecules and condensed phases, enabling efficient electronic structure calculations beyond traditional methods.

Contribution

The authors develop ML surrogates for correlated wavefunction methods to predict 2-RDMs with high fidelity, facilitating large-scale, energy-conserving molecular dynamics and condensed-phase simulations.

Findings

ML models accurately predict 2-RDMs for correlated methods

Enabled large solvated system calculations at Hartree-Fock cost

Achieved coupled-cluster-quality electronic structure for glucose in water

Abstract

Machine learning is rapidly accelerating materials and chemical discovery, but most current models target energies, forces, or selected molecular properties rather than the underlying many-body electronic structure. Learning electronic-structure proxies, such as reduced density matrices, offers a path to surrogates that can predict a broad range of observables from a single ML model. Short of learning the full wavefunction, the two-electron reduced density matrix (2-RDM) is among the most information-rich, minimally lossy targets, providing direct access to expectation values of arbitrary one- and two-electron operators regardless of the strength of the underlying electron correlation. Here we show that learning the 2-RDM is a feasible goal, yielding exceptionally accurate models. We develop surrogates for correlated wavefunction methods (including configuration interaction and coupled cluster) that yield 2-RDMs with sufficient fidelity to provide direct, training-free access to energies and forces for driving energy-conserving molecular dynamics. To tackle realistic molecular condensed phases, we leverage a many-body expansion of the 2-RDM, using our ML models to supply the expansion terms and enabling ML-powered, coupled-cluster-quality electronic structure and energetics for large solvated systems. As a demonstration, we showcase a coupled-cluster-level electronic-structure calculation of glucose solvated by 500 water molecules achieved at Hartree-Fock cost. This work establishes a general framework for learning correlated electronic structure with high fidelity and deploying it to systems beyond the reach of conventional ab initio methods.