Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2606.06584

Bootstrap bounds for Quantum Spin Systems using String Operators

Nisarg Chadha, Michael G. Scheer, Eslam Khalaf

Published Jun 4, 2026• 0 citations
arXivHTMLPDF

Abstract

Bootstrap is a numerical many-body method that provides rigorous bounds on ground-state observables by imposing a set of necessary constraints on the expectation values of operators. The quality of the resulting bounds is sensitive to the choice of operators entering the constraints. In particular, bounds on ground-state correlations are often loose in spontaneous symmetry-breaking (SSB) phases, since local operator sets cannot exclude domain-wall excitations. In this work, we introduce non-local, string-like operators into the bootstrap and show that the program can be formulated directly in thermodynamic limit. We then apply our construction to several 1D spin models. First, we obtain a significant tightening of the bounds in the SSB phase of the 1D transverse-field Ising model. Using the 1D axial next-nearest-neighbor Ising model, we further show that this tightening allows for a quantitative estimate of the locations of phase boundaries. Finally, we generalize the string operators to the 1D $\mathbb{Z}_3$ chiral clock model and track the behavior of the bounds across the phase diagram. Our results broaden the class of constraints available to the bootstrap and open a route toward bootstrapping more general symmetry-broken and topological phases, where the relevant constraints may involve non-local or extended operators.

DOI
10.48550/arxiv.2606.06584

Related papers

Defect Bootstrap: Tight Ground State Bounds in Spontaneous Symmetry Breaking Phases

Michael G. Scheer, Nisarg Chadha, Da-Chuan Lu +1 more

The recent development of bootstrap methods based on semidefinite relaxations of positivity constraints has enabled rigorous two-sided bounds on local observables directly in the thermodynamic limit. However, these bounds inevitably become loose in symmetry broken phases, where local constraints are insufficient to capture long-range order. In this work, we identify the origin of this looseness as order parameter defects which are difficult to remove using local operators. We introduce a $\textit{defect bootstrap}$ framework that resolves this limitation by embedding the system into an auxiliary $\textit{defect model}$ equipped with ancilla degrees of freedom. This construction effectively enables local operators to remove order parameter defects, yielding tighter bounds in phases with spontaneous symmetry breaking. This approach can be applied broadly to pairwise-interacting local lattice models with discrete or continuous internal symmetries that satisfy a property we call $\textit{defect diamagnetism}$, which requires that the ground state energy does not decrease upon adding any finite number of symmetry defects. Applying the method to the transverse field Ising models in 1D and 2D, we obtain significantly improved bounds on energy densities and spin correlation functions throughout the symmetry broken phase in 1D and deep within the phase in 2D. Our results demonstrate that physically motivated constraint sets can dramatically enhance the power of bootstrap methods for quantum many-body systems.

cond-mat.str-elCross-listedarXiv:2511.20860•Nov 27, 2025•0 citations

Bootstrap certification of string order in quantum spin chains

Sagnik Banerjee, Haoyu Guo, Debanjan Chowdhury

The many-body bootstrap certifies ground-state properties by minimizing energy over correlators constrained only by positivity, without any variational wavefunction. Its cost, however, grows exponentially with the size of the operators involved, placing the long-range string correlators that diagnose topological phases out of reach. We overcome this by treating string operators as primary objects: bare and endpoint-dressed strings satisfy a closed operator algebra with one another and with local words, yielding a semidefinite program whose cost grows only polynomially with string length. A single computation then yields estimates of string correlators of all lengths. Rigorous two-sided bounds are obtained for each target string once the ground-state energy is pinned within a window. When benchmarked against density-matrix renormalization group computations on the cluster Ising and spin$-1$ Heisenberg chains, the method certifies the nonlocal string order of the cluster and Haldane phases directly from the Hamiltonian. Our rigorous certification of string order parameters reveal that the sharpness of the bounds is set by which operators enter the calculation, not by their length alone. By exploiting algebraic closure of strings and local words, our framework lifts spatially extended observables into polynomially tractable bootstrap variables, opening up a route to wavefunction-free certification of nonlocal order in quantum many-body systems.

cond-mat.str-elCross-listedarXiv:2609.22420•Sep 22, 2026•0 citations

Bootstrapping Disordered Quantum Systems

Yaprak Önder, Michael G. Scheer, Minjae Cho +1 more

We present a bootstrap framework that yields rigorous two-sided bounds on ground-state observables of quantum systems with quenched disorder. The key idea is to extend the operator algebra to include classical disorder variables and to encode the disorder distribution through a set of moments. The method directly bounds disorder-averaged expectation values for finite or infinite lattice systems and for arbitrary disorder distributions. In contrast to conventional approaches, we do not need to sample disorder realizations, which enables us to make thorough use of symmetries that hold only on average rather than for each disorder realization. Additionally, we show that multi-replica operators can be used to access disorder fluctuations such as the variance of an observable. We apply the method to the one-dimensional random transverse-field Ising model (RTFIM) with discrete, uniform, and Gaussian disorder, obtaining tight bounds on the disorder-averaged ground-state energy density and short-range spin-spin correlators, as well as bounds on the disorder variance of the energy density. Additionally, we show that the difference between the upper and lower bounds on the disorder-averaged ground-state energy density, which we term the bootstrap gap, can be used to map out the phase diagram of the RTFIM. The bootstrap gap is small and flat within the paramagnetic and ferromagnetic phases at weak disorder, grows upon entering the Griffiths regions, and peaks along the critical line.

arXiv:2609.20916•Sep 17, 2026•0 citations

Quantum Many-body Bootstrap

Xizhi Han

A numerical bootstrap method is proposed to provide rigorous and nontrivial bounds in general quantum many-body systems with locality. In particular, lower bounds on ground state energies of local lattice systems are obtained by imposing positivity constraints on certain operator expectation values. Complemented with variational upper bounds, ground state observables are constrained to be within a narrow range. The method is demonstrated with the Hubbard model in one and two dimensions, and bounds on ground state double occupancy and magnetization are discussed.

arXiv:2006.06002•Jun 10, 2020•6 citations

Bootstrapping the Quantum Hall Problem

Qiang Gao, Ryan A. Lanzetta, Patrick Ledwith +2 more

The bootstrap method aims to solve problems by imposing constraints on the space of physical observables, which often follow from physical assumptions such as positivity and symmetry. Here, we employ a bootstrap approach to study interacting electrons in the lowest Landau level by minimizing the energy as a function of the static structure factor subject to a set of constraints, bypassing the need to construct the full many-body wave function. This approach rigorously lower bounds the ground state energy, making it complementary to conventional variational upper bounds. We show that the lower bound we obtain is relatively tight, within at most a few percent from the ground state energy computed with exact diagonalization at small system sizes, and generally gets tighter as we include more constraints. We also show that, by combining the bootstrap lower bounds with variational Monte Carlo calculations for various quantum Hall states, we can obtain two-sided bounds on the ground state energy that rigorously bounds the error to below a few percentage for large system sizes where exact diagonalization is not accessible. Beyond the energetics, our results reproduce the correct power law dependence of the pair correlation function at short distances and the existence of a large entanglement gap in the two-particle entanglement spectra for the Laughlin states at ν = 1 / 3 . We further identify signatures of the composite Fermi liquid state close to half filling. This shows that the bootstrap approach is capable, in principle, of describing nontrivial gapped topologically ordered, as well as gapless, phases. We further study the evolution of the ground state energy and correlation functions as we interpolate between the lowest and first Landau levels at half filling and find indications of a phase transition from a composite Fermi liquid to a non-Abelian state. We also discuss more generally how the geometry of the allowed set of structure factors can be used to diagnose phase transitions. At the end, we discuss possible extensions and current limitations of this approach and how they can be overcome in future studies. Our work establishes numerical bootstrap as a promising method to study many-body phases in topological bands, paving the way to its application in moiré platforms where the energetic competition between fractional quantum anomalous Hall, symmetry broken, and gapless states remains poorly understood.

arXiv:2409.10619•Physical Review X•Jul 8, 2025•11 citations

Bootstrapping the gap in quantum spin systems

Colin Oscar Nancarrow, Yuan Xin

Abstract In this work we report on a new bootstrap method for quantum mechanical problems that closely mirrors the setup from conformal field theory (CFT). We use the equations of motion to develop an analogue of the conformal block expansion for matrix elements and impose crossing symmetry in order to place bounds on their values. The method can be applied to any quantum mechanical system with a local Hamiltonian, and we test it on an anharmonic oscillator model as well as the (1 + 1)-dimensional transverse field Ising model (TFIM). For the anharmonic oscillator model we show that a small number of crossing equations provides an accurate solution to the spectrum and matrix elements. For the TFIM we show that the Hamiltonian equations of motion, translational invariance and global symmetry selection rules imposes a rigorous bound on the gap and the matrix elements of TFIM in the thermodynamic limit. The bound improves as we consider larger systems of crossing equations, ruling out more finite-volume solutions. Our method provides a way to probe the low energy spectrum of an infinite lattice from the Hamiltonian rigorously and without approximation.

arXiv:2211.03819•Journal of High Energy Physics•Aug 11, 2023•20 citations
Bootstrap is a numerical many-body method that provides rigorous bounds on ground-state observables by imposing a set of necessary constraints on the expectation values of operators. The quality of the resulting bounds is sensitive to the choice of operators entering the constraints. In particular, bounds on ground-state correlations are often loose in spontaneous symmetry-breaking (SSB) phases, since local operator sets cannot exclude domain-wall excitations. In this work, we introduce non-local, string-like operators into the bootstrap and show that the program can be formulated directly in thermodynamic limit. We then apply our construction to several 1D spin models. First, we obtain a significant tightening of the bounds in the SSB phase of the 1D transverse-field Ising model. Using the 1D axial next-nearest-neighbor Ising model, we further show that this tightening allows for a quantitative estimate of the locations of phase boundaries. Finally, we generalize the string operators to the 1D Z3\mathbb{Z}_3Z3​ chiral clock model and track the behavior of the bounds across the phase diagram. Our results broaden the class of constraints available to the bootstrap and open a route toward bootstrapping more general symmetry-broken and topological phases, where the relevant constraints may involve non-local or extended operators.
The recent development of bootstrap methods based on semidefinite relaxations of positivity constraints has enabled rigorous two-sided bounds on local observables directly in the thermodynamic limit. However, these bounds inevitably become loose in symmetry broken phases, where local constraints are insufficient to capture long-range order. In this work, we identify the origin of this looseness as order parameter defects which are difficult to remove using local operators. We introduce a defect bootstrap\textit{defect bootstrap}defect bootstrap framework that resolves this limitation by embedding the system into an auxiliary defect model\textit{defect model}defect model equipped with ancilla degrees of freedom. This construction effectively enables local operators to remove order parameter defects, yielding tighter bounds in phases with spontaneous symmetry breaking. This approach can be applied broadly to pairwise-interacting local lattice models with discrete or continuous internal symmetries that satisfy a property we call defect diamagnetism\textit{defect diamagnetism}defect diamagnetism, which requires that the ground state energy does not decrease upon adding any finite number of symmetry defects. Applying the method to the transverse field Ising models in 1D and 2D, we obtain significantly improved bounds on energy densities and spin correlation functions throughout the symmetry broken phase in 1D and deep within the phase in 2D. Our results demonstrate that physically motivated constraint sets can dramatically enhance the power of bootstrap methods for quantum many-body systems.
The many-body bootstrap certifies ground-state properties by minimizing energy over correlators constrained only by positivity, without any variational wavefunction. Its cost, however, grows exponentially with the size of the operators involved, placing the long-range string correlators that diagnose topological phases out of reach. We overcome this by treating string operators as primary objects: bare and endpoint-dressed strings satisfy a closed operator algebra with one another and with local words, yielding a semidefinite program whose cost grows only polynomially with string length. A single computation then yields estimates of string correlators of all lengths. Rigorous two-sided bounds are obtained for each target string once the ground-state energy is pinned within a window. When benchmarked against density-matrix renormalization group computations on the cluster Ising and spin−1-1−1 Heisenberg chains, the method certifies the nonlocal string order of the cluster and Haldane phases directly from the Hamiltonian. Our rigorous certification of string order parameters reveal that the sharpness of the bounds is set by which operators enter the calculation, not by their length alone. By exploiting algebraic closure of strings and local words, our framework lifts spatially extended observables into polynomially tractable bootstrap variables, opening up a route to wavefunction-free certification of nonlocal order in quantum many-body systems.
We present a bootstrap framework that yields rigorous two-sided bounds on ground-state observables of quantum systems with quenched disorder. The key idea is to extend the operator algebra to include classical disorder variables and to encode the disorder distribution through a set of moments. The method directly bounds disorder-averaged expectation values for finite or infinite lattice systems and for arbitrary disorder distributions. In contrast to conventional approaches, we do not need to sample disorder realizations, which enables us to make thorough use of symmetries that hold only on average rather than for each disorder realization. Additionally, we show that multi-replica operators can be used to access disorder fluctuations such as the variance of an observable. We apply the method to the one-dimensional random transverse-field Ising model (RTFIM) with discrete, uniform, and Gaussian disorder, obtaining tight bounds on the disorder-averaged ground-state energy density and short-range spin-spin correlators, as well as bounds on the disorder variance of the energy density. Additionally, we show that the difference between the upper and lower bounds on the disorder-averaged ground-state energy density, which we term the bootstrap gap, can be used to map out the phase diagram of the RTFIM. The bootstrap gap is small and flat within the paramagnetic and ferromagnetic phases at weak disorder, grows upon entering the Griffiths regions, and peaks along the critical line.
A numerical bootstrap method is proposed to provide rigorous and nontrivial bounds in general quantum many-body systems with locality. In particular, lower bounds on ground state energies of local lattice systems are obtained by imposing positivity constraints on certain operator expectation values. Complemented with variational upper bounds, ground state observables are constrained to be within a narrow range. The method is demonstrated with the Hubbard model in one and two dimensions, and bounds on ground state double occupancy and magnetization are discussed.
The bootstrap method aims to solve problems by imposing constraints on the space of physical observables, which often follow from physical assumptions such as positivity and symmetry. Here, we employ a bootstrap approach to study interacting electrons in the lowest Landau level by minimizing the energy as a function of the static structure factor subject to a set of constraints, bypassing the need to construct the full many-body wave function. This approach rigorously lower bounds the ground state energy, making it complementary to conventional variational upper bounds. We show that the lower bound we obtain is relatively tight, within at most a few percent from the ground state energy computed with exact diagonalization at small system sizes, and generally gets tighter as we include more constraints. We also show that, by combining the bootstrap lower bounds with variational Monte Carlo calculations for various quantum Hall states, we can obtain two-sided bounds on the ground state energy that rigorously bounds the error to below a few percentage for large system sizes where exact diagonalization is not accessible. Beyond the energetics, our results reproduce the correct power law dependence of the pair correlation function at short distances and the existence of a large entanglement gap in the two-particle entanglement spectra for the Laughlin states at ν = 1 / 3 . We further identify signatures of the composite Fermi liquid state close to half filling. This shows that the bootstrap approach is capable, in principle, of describing nontrivial gapped topologically ordered, as well as gapless, phases. We further study the evolution of the ground state energy and correlation functions as we interpolate between the lowest and first Landau levels at half filling and find indications of a phase transition from a composite Fermi liquid to a non-Abelian state. We also discuss more generally how the geometry of the allowed set of structure factors can be used to diagnose phase transitions. At the end, we discuss possible extensions and current limitations of this approach and how they can be overcome in future studies. Our work establishes numerical bootstrap as a promising method to study many-body phases in topological bands, paving the way to its application in moiré platforms where the energetic competition between fractional quantum anomalous Hall, symmetry broken, and gapless states remains poorly understood.
Abstract In this work we report on a new bootstrap method for quantum mechanical problems that closely mirrors the setup from conformal field theory (CFT). We use the equations of motion to develop an analogue of the conformal block expansion for matrix elements and impose crossing symmetry in order to place bounds on their values. The method can be applied to any quantum mechanical system with a local Hamiltonian, and we test it on an anharmonic oscillator model as well as the (1 + 1)-dimensional transverse field Ising model (TFIM). For the anharmonic oscillator model we show that a small number of crossing equations provides an accurate solution to the spectrum and matrix elements. For the TFIM we show that the Hamiltonian equations of motion, translational invariance and global symmetry selection rules imposes a rigorous bound on the gap and the matrix elements of TFIM in the thermodynamic limit. The bound improves as we consider larger systems of crossing equations, ruling out more finite-volume solutions. Our method provides a way to probe the low energy spectrum of an infinite lattice from the Hamiltonian rigorously and without approximation.