Automated generation and execution of large-scale semidefinite programs for variational two-electron reduced density matrix theory
Abstract
The variational two-electron reduced density matrix (v2RDM) method provides a wave-function-free approach to the description of strong electron correlation whose computational cost scales polynomially with the number of correlated orbitals. Here, we present an automated pipeline for the generation and execution of semidefinite programs (SDPs) associated with the v2RDM problem. Algebraic constraints encoding ensemble N-representability conditions are generated via the p†q second-quantized algebra package, and the p†q-generated expressions are automatically translated into a sparse-matrix representation that can be digested by an appropriate SDP solver. We have implemented a boundary-point SDP (BPSDP) solver in the libsdp package that can execute the v2RDM SDPs on either CPU or GPU hardware. We apply the code in several contexts. First, two-particle conditions are used to describe the 76-orbital model of the iron-molybdenum cofactor (FeMoco) of nitrogenase, the simulation of which requires roughly 25 hours of compute time on a single Nvidia H100 GPU. Second, the performance of high-order (four-particle positivity) conditions is evaluated on a challenging two-dimensional H10 cluster problem. Lastly, new classes of positivity constraints for mixed fermion-boson systems are generated and applied within the context of ab initio cavity quantum electrodynamics simulations.
Related
- Automated transition state generation for mechanistic exploration in organic synthesis
- Graph Theory-Based De Novo Strategy for Automated Generation of Organic Molecular Cage
- From Execution to Discovery : Physical Verification as the Key Gap in AI Agents for Science
- ETALON: Verified Execution and Persistent Research Memory for Agentic Computer-Aided Drug Design
- vibe-queue: provenance-aware job execution for computational chemistry across workstations and clusters
BibTeX
@misc{deprince2026automated,
title = {{Automated generation and execution of large-scale semidefinite programs for variational two-electron reduced density matrix theory}},
author = {A. Eugene DePrince III and Nicholas C. Rubin},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.30278},
eprint = {2610.30278},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.30278},
}
RIS
TY - UNPB
TI - Automated generation and execution of large-scale semidefinite programs for variational two-electron reduced density matrix theory
AU - A. Eugene DePrince III
AU - Nicholas C. Rubin
PY - 2026
DA - 2026/10/21
PB - arXiv
DO - 10.48550/arxiv.2610.30278
UR - https://arxiv.org/abs/2610.30278
AB - The variational two-electron reduced density matrix (v2RDM) method provides a wave-function-free approach to the description of strong electron correlation whose computational cost scales polynomially with the number of correlated orbitals. Here, we present an automated pipeline for the generation and execution of semidefinite programs (SDPs) associated with the v2RDM problem. Algebraic constraints encoding ensemble N-representability conditions are generated via the p^†q second-quantized algebra package, and the p^†q-generated expressions are automatically translated into a sparse-matrix representation that can be digested by an appropriate SDP solver. We have implemented a boundary-point SDP (BPSDP) solver in the libsdp package that can execute the v2RDM SDPs on either CPU or GPU hardware. We apply the code in several contexts. First, two-particle conditions are used to describe the 76-orbital model of the iron-molybdenum cofactor (FeMoco) of nitrogenase, the simulation of which requires roughly 25 hours of compute time on a single Nvidia H100 GPU. Second, the performance of high-order (four-particle positivity) conditions is evaluated on a challenging two-dimensional H_10 cluster problem. Lastly, new classes of positivity constraints for mixed fermion-boson systems are generated and applied within the context of ab initio cavity quantum electrodynamics simulations.
ER -