From 4D to 2D: Efficient computation of the electron correlation energies of atoms in molecular systems
- 1Department of Chemistry
- 2The University of Manchester
- 3Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O.Box 91, Budapest, Hungary, H-1521
Abstract
The accurate modeling of (bio)molecular systems remains a major challenge in computational chemistry. In recent years, machine learning potentials (MLPs) have emerged as a promising approach, combining the efficiency of classical force fields with near quantum-mechanical accuracy. Among these, the physics-informed MLP FFLUX employs the Interacting Quantum Atoms (IQA) energy partitioning scheme to decompose molecular energies into atomic contributions derived from quantum mechanics for machine learning. However, incorporating two-electron correlation energies into IQA, and consequently into FFLUX, is currently unfeasible for larger systems en masse, due to the steep computational scaling associated with the necessary storage and evaluation of the two-particle density matrix (2PDM). In this work, we introduce a proof-of-concept workflow that employs density fitting (DF) to compress the 2PDM into a two-index representation, thereby accelerating the IQA partitioning of two-electron correlation energies. By expressing traditionally four-centre electron repulsion integrals in terms of auxiliary basis functions, the DF approximation enables a size reduction of the 2PDM from O(N^4) to O(NDF^2), where N is the number of basis functions and NDF is the size of the auxiliary basis set. The method proposed here is demonstrated using coupled cluster singles and doubles densities. The resulting reduction in computational cost and memory requirements enables correlated IQA analyses for significantly larger systems than previously accessible. Thus, this development provides an important step towards incorporating explicit electron correlation and hence dispersion directly into FFLUX atomic ML models, thereby improving both their physical robustness and transferability. This advance will eliminate the need to artificially add empirical potentials, or a posteriori dispersion correction.
Related
- Comparison of pair-source and correlation-function based Lévy analyses in Ar+Sc system at SPS energies
- Equality of the Jellium and the Uniform Electron Gas Energies
- Ionization Energies, Electron Affinities, Bandgaps, Exciton Binding Energies, and Polarization Energies of Orientation-Controlled Picene, [6]-Phenacene, and [7]-Phenacene Thin Films
- Accurate Atomic Correlation Energies and Virial Trends for Correlation-Consistent Effective Core Potentials III: Quantitative Separation of Fixed-Node Biases and Localization Errors in Cs-Rn
- Exchange-Correlation Potentials and Energies from Inverse Generalized Kohn-Sham Calculations
BibTeX
@misc{bandyopadhyay2026from,
title = {{From 4D to 2D: Efficient computation of the electron correlation energies of atoms in molecular systems}},
author = {Prasanta Bandyopadhyay and Gyula Samu and Mark A. Vincent and Olivia P. Aten and Peter R. Nagy and Paul L. A. Popelier},
howpublished = {ChemRxiv},
year = {2026},
month = sep,
doi = {10.26434/chemrxiv.15008932/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15008932/v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - From 4D to 2D: Efficient computation of the electron correlation energies of atoms in molecular systems
AU - Prasanta Bandyopadhyay
AU - Gyula Samu
AU - Mark A. Vincent
AU - Olivia P. Aten
AU - Peter R. Nagy
AU - Paul L. A. Popelier
PY - 2026
DA - 2026/09/16
PB - ChemRxiv
DO - 10.26434/chemrxiv.15008932/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15008932/v1
AB - The accurate modeling of (bio)molecular systems remains a major challenge in computational chemistry. In recent years, machine learning potentials (MLPs) have emerged as a promising approach, combining the efficiency of classical force fields with near quantum-mechanical accuracy. Among these, the physics-informed MLP FFLUX employs the Interacting Quantum Atoms (IQA) energy partitioning scheme to decompose molecular energies into atomic contributions derived from quantum mechanics for machine learning. However, incorporating two-electron correlation energies into IQA, and consequently into FFLUX, is currently unfeasible for larger systems en masse, due to the steep computational scaling associated with the necessary storage and evaluation of the two-particle density matrix (2PDM). In this work, we introduce a proof-of-concept workflow that employs density fitting (DF) to compress the 2PDM into a two-index representation, thereby accelerating the IQA partitioning of two-electron correlation energies. By expressing traditionally four-centre electron repulsion integrals in terms of auxiliary basis functions, the DF approximation enables a size reduction of the 2PDM from O(N^4) to O(NDF^2), where N is the number of basis functions and NDF is the size of the auxiliary basis set. The method proposed here is demonstrated using coupled cluster singles and doubles densities. The resulting reduction in computational cost and memory requirements enables correlated IQA analyses for significantly larger systems than previously accessible. Thus, this development provides an important step towards incorporating explicit electron correlation and hence dispersion directly into FFLUX atomic ML models, thereby improving both their physical robustness and transferability. This advance will eliminate the need to artificially add empirical potentials, or a posteriori dispersion correction.
ER -