aiqbench analysis
This research introduces local relaxation hierarchies and parallel message passing algorithms that offer polynomial-time classical methods for estimating quantum ground state energies with rigorous convergence guarantees. By achieving exponential convergence for specific Hamiltonian perturbations, the work addresses critical scalability challenges in solving many-body quantum systems without requiring full quantum hardware. The authors provide a software library to benchmark these approaches against standard convex optimizers, potentially advancing classical simulation capabilities for complex quantum lattices.
Local Relaxation Hierarchies for Quantum Ground State Energies: Convergence Guarantees and Message Passing Algorithms
arXiv:2609.40336v1 Announce Type: new Abstract: Convex relaxation hierarchies provide lower bounds to the ground state energy of quantum many-body systems that can be computed in polynomial time on a classical computer, at any fixed hierarchy level. However, scaling these methods to large systems and accurate approximations remains challenging due to the computational cost of traditional solvers and the scarcity of efficiency guarantees. In this work, we develop local relaxation hierarchies and efficient, highly parallelisable message passing algorithms for estimating the relaxed ground state energies. We show that the first level of the hierarchy---based on local consistency of pairwise reduced density matrices---is exact for commuting Hamiltonians on trees. We further establish that another hierarchy, based on consistent intervals, converges exponentially fast in the interval size to the ground state energy for weak perturbations of separable Hamiltonians on a chain, thereby providing an efficient classical algorithm for these systems. Then, we introduce two variants of message passing algorithms that run in $\mathcal{O}(n/\epsilon^2)$ and $\mathcal{O}(n/\epsilon)$ time for any fixed level of the local hierarchy on bounded-degree graphs, where $\epsilon$ is the precision for the relaxed ground state energy per site. This assumes that the optimal messages have $\mathcal{O}(1)$ norm---a condition we observe in practical settings in our experiments. These algorithms are based on the subgradient method and the Nesterov-type accelerated gradient descent method applied to an entropy-smoothed objective. Finally, we benchmark the message passing algorithms across different quantum Hamiltonians, lattice geometries, and relaxation levels, validating the theoretical predictions and their potential to surpass standard convex optimisation solvers for this problem. We release the resulting library at github.com/rick1924/gse-message-passing.
Abstract:Convex relaxation hierarchies provide lower bounds to the ground state energy of quantum many-body systems that can be computed in polynomial time on a classical computer, at any fixed hierarchy level. However, scaling these methods to large systems and accurate approximations remains challenging due to the computational cost of traditional solvers and the scarcity of efficiency guarantees. In this work, we develop local relaxation hierarchies and efficient, highly parallelisable message passing algorithms for estimating the relaxed ground state energies. We show that the first level of the hierarchy---based on local consistency of pairwise reduced density matrices---is exact for commuting Hamiltonians on trees. We further establish that another hierarchy, based on consistent intervals, converges exponentially fast in the interval size to the ground state energy for weak perturbations of separable Hamiltonians on a chain, thereby providing an efficient classical algorithm for these systems. Then, we introduce two variants of message passing algorithms that run in $\mathcal{O}(n/\epsilon^2)$ and $\mathcal{O}(n/\epsilon)$ time for any fixed level of the local hierarchy on bounded-degree graphs, where $\epsilon$ is the precision for the relaxed ground state energy per site. This assumes that the optimal messages have $\mathcal{O}(1)$ norm---a condition we observe in practical settings in our experiments. These algorithms are based on the subgradient method and the Nesterov-type accelerated gradient descent method applied to an entropy-smoothed objective. Finally, we benchmark the message passing algorithms across different quantum Hamiltonians, lattice geometries, and relaxation levels, validating the theoretical predictions and their potential to surpass standard convex optimisation solvers for this problem. We release the resulting library at this http URL.
Quantum Physics (quant-ph); Distributed, Parallel, and Cluster Computing (cs.DC)
arXiv:2609.40336 [quant-ph]
(or arXiv:2609.40336v1 [quant-ph] for this version)
https://doi.org/10.48550/arXiv.2609.40336 Focus to learn more arXiv-issued DOI via DataCite (pending registration)
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.
Full text under the publisher's open licence · 361 words
Related coverage

Quantum Max d-Cut via qudit swap operators
Quantum Journal · Sep 3, 2026
Dynamical BEC formation and vanishing entropy for trapped bosons in non-relativistic QED via nonlinear resonance cascades
arXiv · Oct 14, 2026
A convergent hierarchy of spectral gap certificates for qubit Hamiltonians
arXiv · Oct 5, 2026
Convergence rates of Sum-of-Hermitian-Squares Hierarchies for the Pauli algebra
arXiv · Jun 4, 2026