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.
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.
We introduce a Basis Adaptive (BA) algorithm for hybrid quantum-classical simulation of correlated quantum many-body systems. Starting from a small set of physically motivated bitstrings, the algorithm iteratively applies a single-step first-order Trotterized circuit on a quantum processor, filters the sampled configurations by enforcing U(1) spin conservation and lattice reflection symmetry, and classically diagonalizes the Hamiltonian in the resulting reduced Hilbert space. This design avoids the variational optimization overhead of VQE, the deep coherent circuits required by QPE, and the symmetry-violating subspaces that arise in SKQD. The ground-state energy error is bounded analytically by 8∥H∥(1−αDT)1/2, where αDT is the probability weight captured by the DT sampled basis states. This bound connects algorithm performance directly to ground-state sparsity and explains the observed accuracy hierarchy across different phases. Benchmarked on the spin-1/2 Heisenberg XXZ chain (up to N=62 qubits on the IBM Heron processor), the algorithm achieves a 3.5% energy error in the gapped Neel phase (Δ=2.0) and below 0.5% at the ferromagnetic boundary (Δ=−1.0). The accuracy degrades to 28.7% in the strongly quasi-long-range-ordered regime (Δ=0.5). Spin-spin correlation functions are reproduced across all regimes, confirming that symmetry-filtered real-time sampling provides a practical and noise-resilient pathway to ground-state properties on near-term quantum hardware.
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 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.
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.
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.
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.