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.
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.
A bstract We present a new computational framework combining coarse-graining techniques with bootstrap methods to study quantum many-body systems. The method efficiently computes rigorous upper and lower bounds on both zero- and finite-temperature expectation values of any local observables of infinite quantum spin chains. This is achieved by using tensor networks to coarse-grain bootstrap constraints, including positivity, translation invariance, equations of motion, and energy-entropy balance inequalities. Coarse-graining allows access to constraints from significantly larger subsystems than previously possible, yielding tighter bounds compared to those obtained without coarse-graining.
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 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.
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.
We introduce a nullspace-guided adaptive (NGA) bootstrap method that improves the energy lower bounds of quantum many-body ground states by refining the bootstrap basis in a dynamic and incremental way. At each iteration, the optimized moment matrix reveals a nullspace of saturated positivity directions, which is intuitively interpreted as annihilators of the approximate ground-state subspace. The NGA bootstrap then prunes operators with small nullspace leverage and grows the basis along descendants of these null directions. By applying the NGA bootstrap to the transverse-field Ising chain, we obtain nearly exact energy lower bounds because the algorithm automatically discovers the eigenoperator structure in terms of Jordan-Wigner fermions from a minimal local bootstrap basis. For the Hubbard chain, it improves upon state-of-the-art energy lower bounds by up to two orders of magnitude, reaching errors ranging from 10−3 down to 10−5 in the strongly correlated regimes. We further show that the NGA framework can be used to improve the certified two-sided bounds on general observables. In addition, the bootstrap error decreases approximately as a power law with increasing computational resources. These results suggest that our method provides a practical and scalable route toward accurate bootstrap of general quantum many-body systems.
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 framework that resolves this limitation by embedding the system into an auxiliary 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, 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.