---
title: Quantum Mechanics Bootstrap
url: https://www.emergentmind.com/topics/quantum-mechanics-bootstrap
type: topic
---

# Quantum Mechanics Bootstrap

Quantum-mechanics bootstrap is a constraint-based program for quantum systems in which spectral data, moments, and, in some extensions, dynamical observables are inferred from positivity, operator identities, and self-consistency conditions rather than from direct solution of the Schrödinger equation. The bootstrap philosophy dates back to the 1960s, but its recent formulation in quantum mechanics was shaped by analogies with conformal bootstrap and relativistic scattering-amplitude bootstrap, and it has since been applied to confining one-dimensional Hamiltonians, periodic systems, scattering on the half-line, low-energy nuclear Hamiltonians, non-Hermitian models, Calabi–Yau quantum systems, time evolution, and three-dimensional central potentials [2201.00551] [2209.14332] [2307.11724] [2409.06784] [2412.08721] [2512.09041].

## 1. Foundational framework

The basic quantum-mechanical input is the positivity statement that for any state \(|\Psi\rangle\) and any operator \(\mathcal O\),
\[
\langle \Psi|\mathcal O^\dagger \mathcal O|\Psi\rangle \ge 0.
\]
Choosing a finite operator set \(\{\mathcal O_i\}\) and writing \(\mathcal O=\sum_i \alpha_i \mathcal O_i\) turns this into positivity of a bootstrap, Gram, or moment matrix,
\[
\mathcal M_{ij}=\langle \mathcal O_i^\dagger \mathcal O_j\rangle,
\qquad
\mathcal M\succeq 0.
\]
In one-dimensional polynomial problems this often reduces to a Hankel matrix \(M^{(K)}_{ij}=x_{i+j}\), where \(x_n=\langle x^n\rangle\) and \(x_0=1\); in finite-dimensional Hamiltonian truncations, such as the deuteron in a harmonic-oscillator basis, it becomes a matrix of expectation values of projectors \(|i\rangle\langle j|\) [2209.14332] [2201.00551].

The second ingredient is the eigenstate condition. For an energy eigenstate \(|\psi\rangle\) with \(H|\psi\rangle=E|\psi\rangle\), one imposes relations such as
\[
\langle [H,O]\rangle=0,
\qquad
\langle O H\rangle = E\langle O\rangle,
\]
or their appropriate non-Hermitian generalizations. These identities generate recursions among moments and mixed correlators. Depending on the problem, the operator basis may be monomials in \(x\) and \(p\), exponentials \(e^{m\hat x+n\hat p}\), circle operators \(e^{imx}p^n\), or finite-basis projectors \(|i\rangle\langle j|\) [2209.14332] [2111.08442] [2109.02701].

In this formulation, the bootstrap can also be regarded as a generalization of uncertainty relations. For fixed-energy states, the same positivity machinery that yields the Heisenberg inequality yields state-independent bounds on observables such as \(\langle x\rangle\), \(\langle x^2\rangle\), and \(\langle p\rangle\). For the harmonic oscillator, for example, the universal bound
\[
-x_*(E)\le \langle x\rangle \le x_*(E),
\qquad
x_*(E)=\sqrt{2(E-\hbar/2)},
\]
follows from this broader positivity structure [2208.09370].

## 2. Linearization and semidefinite optimization

A central technical step is to fix the trial energy \(E\). In polynomial one-dimensional systems this linearizes the recursion relations among the moments. For a potential
\[
V(x)=\sum_{n=1}^d a_n x^n,
\]
the bootstrap reduces higher moments to a finite set of independent moments \(x_1,\dots,x_{d-1}\), and the truncated Hankel matrix can be written as
\[
M^{(K)}(E)=\sum_{n=0}^{d-1} x_n F_n(E).
\]
The spectral question then becomes a semidefinite feasibility problem: does there exist a set of low moments such that \(M^{(K)}(E)\succeq 0\)? A useful equivalent formulation introduces a slack variable \(t\) and maximizes the minimal eigenvalue through
\[
M^{(K)}(E,x)-tI \succeq 0.
\]
Positive optimal \(t\) means that the chosen \(E\) is allowed at truncation depth \(K\); negative optimal \(t\) quantifies failure of positivity [2209.14332].

The same logic extends beyond the simplest Hankel setup. In half-line problems one must impose positivity of both \(M_{ij}=x_{i+j}\) and \(M'_{ij}=x_{1+i+j}\), conveniently encoded by a block-diagonal semidefinite constraint. In finite harmonic-oscillator spaces, as in the deuteron problem, one solves linear expectation-value constraints from normalization and the discrete eigenvalue equation and then checks positivity of the resulting finite bootstrap matrix [2307.11724] [2201.00551].

A compact summary of recurrent formulations is given below.

| Setting | Bootstrap data | Representative papers |
|---|---|---|
| 1D polynomial Schrödinger operators | \(x_n=\langle x^n\rangle\), Hankel matrix \(M^{(K)}_{ij}=x_{i+j}\) | [2209.14332] |
| Finite-basis nuclear Hamiltonians | \(X_{ij}=\langle |i\rangle\langle j| \rangle\), projector bootstrap matrix | [2201.00551] |
| Circle and periodic systems | \(\langle e^{imx}p^n\rangle\) or mixed position-momentum moments | [2109.02701], [2209.09919] |
| Half-line scattering | \(M_{ij}=x_{i+j}\), \(M'_{ij}=x_{1+i+j}\) with boundary anomalies | [2307.11724] |

This optimization viewpoint also clarifies why the method often produces bounds rather than exact values at finite truncation. The admissible set is convex, and increasing the operator basis or matrix size monotonically tightens the feasible region.

## 3. Benchmark systems and exact solvability

The earliest systematic tests focused on analytically tractable systems and simple confining Hamiltonians. For the hydrogen atom and the harmonic oscillator, the bootstrap resolves many energy levels and the spectra converge exponentially fast as the matrix size increases [2108.08757]. For the quartic anharmonic oscillator and the double-well potential, the method captures non-perturbative structure; in particular, the double-well study found that the bootstrap correctly captures non-perturbative aspects, while supersymmetric partner potentials reproduce the expected paired spectra, and the singlet sector of the \(O(N)\) vector model agrees with large-\(N\) saddle-point analysis [2108.11416].

A distinct line of work concerns exactly solvable systems. For shape-invariant one-dimensional Hamiltonians, the bootstrap can derive exact energy eigenvalues analytically, and the information of the annihilation operators is also obtained naturally. This was demonstrated numerically for harmonic oscillators, Morse potentials, Rosen–Morse potentials, and hyperbolic Scarf potentials, leading to the claim that the numerical bootstrap can determine the solvability of a given unknown system if it satisfies shape invariance [2504.08586].

Exactness can also hinge on operator ordering. In the interval model with Hamiltonian
\[
H=SZ(1-Z)S,
\]
defined on \(L^2([0,1],dz)\), one ordering of the bootstrap constraints restricts the energy only into bands, whereas an alternative ordering makes a finite number of constraints sufficient to fix the low-lying energy levels exactly. Since the exact spectrum is
\[
E_n=n(n+1),
\qquad n=0,1,2,\dots,
\]
this model provides a concrete example in which finite bootstrap data reproduces exact results rather than asymptotic bounds [2402.03434].

For generic polynomial confining potentials, exactness is replaced by rapid numerical convergence. In the eighth-degree test potential
\[
V(x)=\frac{1}{2}x^2-x^4+\frac{1}{8}x^8,
\]
the interval widths for the first five excited states shrink exponentially with truncation depth,
\[
\bar w(K)\propto e^{-0.83 K},
\]
and the bootstrap values at \(K=30\) match finite-element benchmarks to the quoted precision [2209.14332].

## 4. Periodicity, boundary conditions, and scattering

Periodic and compact-configuration systems require operator bases adapted to translation or winding structure. In the Kronig–Penney problem, the bootstrap was applied to the Schrödinger equation with periodic potentials, with an operator basis involving position and momenta. In this setting the method efficiently computes band gaps of the energy spectrum but has trouble effectively constraining the minimum energy; more complex constraints involving higher powers of momenta were proposed as a possible remedy. The same study also proposed an approach for analytically constructing the dispersion relation associated with the Bloch momentum of the system [2209.09919].

A related compact system is quantum mechanics on a circle with a \(\theta\)-term, implemented as a charged particle on \(S^1\) in a constant gauge potential. Here the bootstrap correctly reproduces correlations among observables for energy eigenstates for any \(\theta\), but it is hard to determine physical quantities as functions of \(\theta\), such as \(E(\theta)\), except at \(\theta=0\) and \(\pi\). This identifies a structural difficulty associated with gauge choice and periodicity rather than a failure of positivity itself [2109.02701].

Scattering on the half-line requires a further modification because boundary conditions become part of the bootstrap data. For
\[
H=-\partial_x^2+V(x)
\]
on \(x\ge 0\), the basic operator identity becomes
\[
[H,O] + (H^\dagger-H)O = 0,
\]
where the second term is an anomaly reflecting the fact that \(O\) may not preserve the domain of the self-adjoint Hamiltonian. For Robin boundary conditions
\[
\psi(0)+a\,\psi'(0)=0,
\qquad
a=-\frac{\psi(0)}{\psi'(0)},
\]
the bootstrap imposes positivity of both \(M_{ij}=x_{i+j}\) and \(M'_{ij}=x_{1+i+j}\). Varying the Robin parameter \(a\) traces the discrete half-line spectrum \(E_n(a)\), and for a purely reflecting one-dimensional problem the reflection coefficient is recovered from
\[
R(k)=\frac{ika+1}{ika-1}=e^{i\delta(k)},
\qquad
\delta(k)=2\tan^{-1}(ka).
\]
This framework was applied to half-harmonic walls, a metastable quartic potential, and the exponential potential of Liouville theory, with WKB used to connect the bootstrap data to the asymptotic scattering phase [2307.11724].

## 5. Generalizations to realistic, geometric, non-Hermitian, and dynamical settings

One prominent extension is low-energy nuclear physics. “Bootstrapping the deuteron” formulates the deuteron in the \({}^3S_1\) channel within pionless effective field theory in a harmonic-oscillator basis and reports the first bootstrap results in low-energy nuclear physics. The Hamiltonian is a finite real symmetric matrix, the operator basis is built from projectors \(\mathcal O_{ij}=|i\rangle\langle j|\), and positivity of the corresponding bootstrap matrix, together with linear self-consistency conditions, reproduces the deuteron ground-state energy and excited levels in excellent agreement with exact diagonalization [2201.00551].

Another extension concerns mirror-curve quantum mechanics from local toric Calabi–Yau geometries. For the local \(\mathbb P^1\times\mathbb P^1\) and local \(\mathbb P^2\) Hamiltonians, the bootstrap is built from exponentials \(e^{m\hat x+n\hat p}\), and an important improvement is to use a larger set of two-dimensional operators instead of one-dimensional ones. The same improved strategy also gives better numerical accuracies for the two-body non-relativistic Toda system and the quartic anharmonic oscillator [2111.08442].

The bootstrap has also been generalized beyond Hermitian stationary spectra. For generic complex polynomial potentials,
\[
H=p^2+V_1(x)+iV_2(x),
\]
one replaces the Hermitian constraints by
\[
\langle R_n|O^\dagger O|R_n\rangle \ge 0,
\]
\[
\langle R_n|(OH-H^\dagger O)|R_n\rangle = 2i E_{\mathcal I}^n \langle R_n|O|R_n\rangle,
\]
\[
\langle R_n|OH|R_n\rangle = (E_{\mathcal R}^n+iE_{\mathcal I}^n)\langle R_n|O|R_n\rangle,
\]
and then bootstraps complex spectra directly. This formulation handles generic non-Hermitian systems, includes PT-symmetric Hamiltonians as a special case, and captures PT-symmetric phase transitions such as the onset of complex eigenvalues in \(x^4+iax\) [2409.06784].

Time evolution admits a different but closely related bootstrap. Instead of a static moment matrix, one considers a time-dependent matrix
\[
M_{ij}(t)=\langle \mathcal O_i^\dagger(t)\mathcal O_j(t)\rangle
\]
subject to positivity, operator identities, known initial data, and Heisenberg equations. This yields a hierarchy of rigorous bounds on observables at later times, systematically generalizing Mandelstam–Tamm-like relations. For any fixed hierarchy level the bounds are tightest at short times and loosen over time, while for fixed time the evidence indicates that increasing the hierarchy can make them arbitrarily tight; the computational effort scales polynomially with the number of degrees of freedom at fixed level [2412.08721].

Three-dimensional central potentials provide a further major generalization. The bootstrap has now been applied to Coulomb, Yukawa, Gaussian, Cornell, and conformal quantum mechanics, including non-algebraic potentials such as Yukawa and Gaussian. In this setting the method is particularly effective for ground-state lower bounds: the Cornell critical coupling is determined to better than one part in \(10^7\), and lower bounds on energies are occasionally accurate to one part in greater than \(10^8\). The same study also analyzes when meaningful upper bounds can and cannot be obtained [2512.09041].

| Domain | Characteristic feature | Representative papers |
|---|---|---|
| Low-energy nuclear physics | Finite-basis projector bootstrap for pionless EFT deuteron | [2201.00551] |
| Mirror-curve / Calabi–Yau QM | Exponential operator basis \(e^{m\hat x+n\hat p}\) and 2D bootstrap matrices | [2111.08442] |
| Non-Hermitian QM | Complex spectra and PT phase transitions | [2409.06784] |
| Time evolution | Hierarchy of rigorous bounds on \(M_{ij}(t)\) | [2412.08721] |
| 3D central potentials | High-precision ground-state bounds for Coulomb, Yukawa, Gaussian, Cornell | [2512.09041] |

## 6. Limitations, misconceptions, and open directions

A recurring limitation is that the strongest results are usually lower bounds on low-lying states, especially the ground state. In three-dimensional central potentials, meaningful upper bounds are obtained only under specific circumstances, whereas lower bounds are routinely much sharper [2512.09041]. In periodic problems the method efficiently constrains band gaps but may leave the minimum energy weakly constrained unless higher-momentum constraints are added [2209.09919].

Another limitation concerns global data that are not naturally encoded by moment positivity. For systems with a \(\theta\)-term, the bootstrap reproduces correlations among observables in energy eigenstates but has difficulty reconstructing \(E(\theta)\) except at \(\theta=0\) and \(\pi\) [2109.02701]. For identical particles, the method faces a more basic obstruction: it has difficulty distinguishing bosons from fermions at the level of expectation-value constraints, and its predictive power in multi-particle systems is therefore limited in the derivation of observables including energy eigenstates. In thermal applications, temperature and entropy cannot be handled, although some quantities in micro-canonical ensembles of integrable systems governed by generalized Gibbs ensembles can still be derived [2208.09370].

Boundary conditions and operator orderings are also structural, not merely technical, inputs. Half-line and interval problems require anomaly terms or domain-sensitive identities, and in the interval model \(H=SZ(1-Z)S\) one operator ordering gives only bands while another produces exact low-lying eigenvalues from finite constraints [2307.11724] [2402.03434]. This suggests that bootstrap performance depends strongly on how the operator algebra is represented in the positivity problem.

The current literature nevertheless delineates a coherent research program. Bootstrap methods in quantum mechanics now cover discrete spectra, band structure, reflection phases, realistic finite-basis nuclear systems, generic complex spectra, exact solvability via shape invariance, and rigorous bounds on real-time evolution [2504.08586]. What remains open is not the existence of useful bootstrap constraints, but the extent to which basis choice, boundary data, and truncation hierarchies can be systematized so that exactness, when present, can be detected reliably and high-precision bounds can be made routine across broader classes of quantum systems.

Source: https://www.emergentmind.com/topics/quantum-mechanics-bootstrap