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:2603.17364

Finite-NNN Bootstrap Constraints in Matrix and Tensor Models

Samuel Laliberte, Reiko Toriumi

Published Mar 18, 2026• 0 citations
arXivHTMLPDF

Abstract

We explore how matrix bootstrap techniques can be used to constrain matrix and tensor models at finite $N$, where $N$ is the dimension of the matrix/tensor, taking a Gaussian model with a quartic interaction as example. For matrix models, we find further evidence that bounds do not depend explicitly on $N$, but rather on properties of multi-trace expectation values. For tensor models, the structure of the Schwinger-Dyson equations allow for bounds that vary as a function of $N$, admitting a broader scan of the parameter space of the theory. In the latter case, we find novel bounds on the two-point function as a function of the quartic coupling of the theory.

DOI
10.48550/arxiv.2603.17364

Related papers

Bootstrapping Tensor Integrals

Nathan Pagliaroli, Carlos I. Pérez-Sánchez, Brayden Smith

This work proposes a bootstrapping with positivity methodology to study random $U(N)^{D}$ invariant tensors in the large $N$ limit. As has been done for $U(N)$ invariant random matrices, we combine the Dyson-Schwinger equations and positivity constraints of moments to approximate the moments of such tensor models. As examples, we bootstrap the quartic and two hexic rank three tensor models. All models studied converge quickly, and for those which have known analytic formulae, they converge to such solutions. We conjecture new explicit formulae for all moments of the rank three quartic model and support this conjecture using bootstrapped results and explicit double-series computations with 'feyntensor'.

arXiv:2604.19714•Apr 21, 2026•0 citations

Toward bootstrapping tensor-network contractions

Seishiro Ono, Yanbai Zhang, Hoi Chun Po

Accurate contraction of tensor networks beyond one dimension is essential in various fields including quantum many-body physics. Existing approaches typically rely on approximate contraction schemes and do not provide certified error bars. We introduce a numerical bootstrap framework which casts the problem of tensor-network contractions into a convex optimization problem, thereby yielding certified lower and upper bounds on expectation values of physical observables. As a proof-of-principle, we construct such constraints explicitly for translationally invariant matrix product states and demonstrate that, assuming a canonical form, second-order-cone relaxation can provide tight bounds on the contraction result. We further demonstrate that when the requirement on canonical form is lifted, a more general semidefinite-programming approach could yield similar tight bounds at higher but still polynomial computational cost. Our work suggests numerical bootstrap could be a possible way forward for the rigorous contractions of tensor networks.

cond-mat.str-elCross-listedarXiv:2603.17856•Mar 19, 2026•0 citations

Redundant moments of non-melonic random tensor models: simplifying the tensor bootstrap

Carlos I. Pérez-Sánchez

We prove the redundancy of a family (the so-called melonic family) of large-$N$ moments in arbitrary tensor models. This reduces the number of independent entries in the positive semidefinite matrices that build the core of the tensor bootstrap (for tensor models there are already several generalizations of the positive semidefinite Toeplitz and Hankel matrices of lattice gauge theory and random matrix bootstrap, respectively). More concretely, we prove that any melonic operator $C$ at large-$N$ has an expectation value that depends on $C$ exclusively through its degree. A similar statement is shown here to hold for the Schwinger-Dyson equations of melonic moments. The strength of our result is also its scope, which is not limited to melonic tensor models.

arXiv:2607.28518•Jul 30, 2026•0 citations

The Discrete Bootstrap for Quantum Thermodynamics

Sandor, Joshua

A novel bootstrap method is exemplified by lattice theories and quantum mechanics. Positivity of expectation values combined with sets of non-trivial constraint equations and semidefinite optimization rigorously and precisely bound correlation functions. Starting with the example of a one-dimensional integral, the technique is built upon to extract limit cycles in a two-dimensional classical dynamical system and spin correlation functions in the Statistical Ising Model. Further positivity constraints and simplifications arising in a large rank limit are implemented to bound the energy expectation value as a function of temperature in Matrix Quantum Mechanics. Along the way, the S-matrix of Chern-Simons-matter theory is computed as an example of computational simplification at large gauge group rank.

•Jan 1, 2025•0 citations

Bosonic tensor models at large N and small ε

Simone Giombi, Igor R. Klebanov, Grigory Tarnopolsky

We study the spectrum of the large $N$ quantum field theory of bosonic rank-3 tensors, the quartic interactions of which are such that the perturbative expansion is dominated by the melonic diagrams. We use the Schwinger-Dyson equations to determine the scaling dimensions of the bilinear operators of arbitrary spin. Using the fact that the theory is renormalizable in $d=4$, we compare some of these results with the $4\ensuremath{-}\ensuremath{\epsilon}$ expansion, finding perfect agreement. This helps elucidate why the dimension of operator ${\ensuremath{\phi}}^{abc}{\ensuremath{\phi}}^{abc}$ is complex for $d<4$: the large $N$ fixed point in $d=4\ensuremath{-}\ensuremath{\epsilon}$ has complex values of the couplings for some of the $O(N{)}^{3}$ invariant operators. We show that a similar phenomenon holds in the $O(N{)}^{2}$ symmetric theory of a matrix field ${\ensuremath{\phi}}^{ab}$, where the double-trace operator has a complex coupling in $4\ensuremath{-}\ensuremath{\epsilon}$ dimensions. We also study the spectra of bosonic theories of rank-$q\ensuremath{-}1$ tensors with ${\ensuremath{\phi}}^{q}$ interactions. In dimensions $d>1.93$, there is a critical value of $q$, above which we have not found any complex scaling dimensions. The critical value is a decreasing function of $d$, and it becomes 6 in $d\ensuremath{\approx}2.97$. This raises a possibility that the large $N$ theory of rank-5 tensors with sextic potential has an IR fixed point which is free of perturbative instabilities for $2.97<d<3$. This theory may be studied using renormalized perturbation theory in $d=3\ensuremath{-}\ensuremath{\epsilon}$.

arXiv:1707.03866•Physical review. D/Physical review. D.•Nov 21, 2017•103 citations

Bootstrapping Euclidean Two-point Correlators

Minjae Cho, Barak Gabai, Henry W. Lin +2 more

We develop a bootstrap approach to Euclidean two-point correlators, in the thermal or ground state of quantum mechanical systems. We formulate the problem of bounding the two-point correlator as a semidefinite programming problem, subject to the constraints of reflection positivity, the Heisenberg equations of motion, and the Kubo-Martin-Schwinger condition or ground-state positivity. In the dual formulation, the Heisenberg equations of motion become "inequalities of motion" on the Lagrange multipliers that enforce the constraints. This enables us to derive rigorous bounds on continuous-time two-point correlators using a finite-dimensional semidefinite or polynomial matrix program. We illustrate this method by bootstrapping the two-point correlators of the ungauged one-matrix quantum mechanics, from which we extract the spectrum and matrix elements of the low-lying adjoint states. Along the way, we provide a new derivation of the energy-entropy balance inequality and establish a connection between the high-temperature two-point correlator bootstrap and the matrix integral bootstrap.

hep-thCross-listedarXiv:2511.08560•Journal of High Energy Physics•Nov 12, 2025•2 citations
We explore how matrix bootstrap techniques can be used to constrain matrix and tensor models at finite NNN, where NNN is the dimension of the matrix/tensor, taking a Gaussian model with a quartic interaction as example. For matrix models, we find further evidence that bounds do not depend explicitly on NNN, but rather on properties of multi-trace expectation values. For tensor models, the structure of the Schwinger-Dyson equations allow for bounds that vary as a function of NNN, admitting a broader scan of the parameter space of the theory. In the latter case, we find novel bounds on the two-point function as a function of the quartic coupling of the theory.
This work proposes a bootstrapping with positivity methodology to study random U(N)DU(N)^{D}U(N)D invariant tensors in the large NNN limit. As has been done for U(N)U(N)U(N) invariant random matrices, we combine the Dyson-Schwinger equations and positivity constraints of moments to approximate the moments of such tensor models. As examples, we bootstrap the quartic and two hexic rank three tensor models. All models studied converge quickly, and for those which have known analytic formulae, they converge to such solutions. We conjecture new explicit formulae for all moments of the rank three quartic model and support this conjecture using bootstrapped results and explicit double-series computations with 'feyntensor'.
Accurate contraction of tensor networks beyond one dimension is essential in various fields including quantum many-body physics. Existing approaches typically rely on approximate contraction schemes and do not provide certified error bars. We introduce a numerical bootstrap framework which casts the problem of tensor-network contractions into a convex optimization problem, thereby yielding certified lower and upper bounds on expectation values of physical observables. As a proof-of-principle, we construct such constraints explicitly for translationally invariant matrix product states and demonstrate that, assuming a canonical form, second-order-cone relaxation can provide tight bounds on the contraction result. We further demonstrate that when the requirement on canonical form is lifted, a more general semidefinite-programming approach could yield similar tight bounds at higher but still polynomial computational cost. Our work suggests numerical bootstrap could be a possible way forward for the rigorous contractions of tensor networks.
We prove the redundancy of a family (the so-called melonic family) of large-NNN moments in arbitrary tensor models. This reduces the number of independent entries in the positive semidefinite matrices that build the core of the tensor bootstrap (for tensor models there are already several generalizations of the positive semidefinite Toeplitz and Hankel matrices of lattice gauge theory and random matrix bootstrap, respectively). More concretely, we prove that any melonic operator CCC at large-NNN has an expectation value that depends on CCC exclusively through its degree. A similar statement is shown here to hold for the Schwinger-Dyson equations of melonic moments. The strength of our result is also its scope, which is not limited to melonic tensor models.
A novel bootstrap method is exemplified by lattice theories and quantum mechanics. Positivity of expectation values combined with sets of non-trivial constraint equations and semidefinite optimization rigorously and precisely bound correlation functions. Starting with the example of a one-dimensional integral, the technique is built upon to extract limit cycles in a two-dimensional classical dynamical system and spin correlation functions in the Statistical Ising Model. Further positivity constraints and simplifications arising in a large rank limit are implemented to bound the energy expectation value as a function of temperature in Matrix Quantum Mechanics. Along the way, the S-matrix of Chern-Simons-matter theory is computed as an example of computational simplification at large gauge group rank.
We study the spectrum of the large NNN quantum field theory of bosonic rank-3 tensors, the quartic interactions of which are such that the perturbative expansion is dominated by the melonic diagrams. We use the Schwinger-Dyson equations to determine the scaling dimensions of the bilinear operators of arbitrary spin. Using the fact that the theory is renormalizable in d=4d=4d=4, we compare some of these results with the $4\ensuremath{-}\ensuremath{\epsilon}$ expansion, finding perfect agreement. This helps elucidate why the dimension of operator ${\ensuremath{\phi}}^{abc}{\ensuremath{\phi}}^{abc}$ is complex for d<4d<4d<4: the large NNN fixed point in $d=4\ensuremath{-}\ensuremath{\epsilon}$ has complex values of the couplings for some of the O(N)3O(N{)}^{3}O(N)3 invariant operators. We show that a similar phenomenon holds in the O(N)2O(N{)}^{2}O(N)2 symmetric theory of a matrix field ${\ensuremath{\phi}}^{ab}$, where the double-trace operator has a complex coupling in $4\ensuremath{-}\ensuremath{\epsilon}$ dimensions. We also study the spectra of bosonic theories of rank-$q\ensuremath{-}1$ tensors with ${\ensuremath{\phi}}^{q}$ interactions. In dimensions d>1.93d>1.93d>1.93, there is a critical value of qqq, above which we have not found any complex scaling dimensions. The critical value is a decreasing function of ddd, and it becomes 6 in $d\ensuremath{\approx}2.97$. This raises a possibility that the large NNN theory of rank-5 tensors with sextic potential has an IR fixed point which is free of perturbative instabilities for 2.97<d<32.97<d<32.97<d<3. This theory may be studied using renormalized perturbation theory in $d=3\ensuremath{-}\ensuremath{\epsilon}$.
We develop a bootstrap approach to Euclidean two-point correlators, in the thermal or ground state of quantum mechanical systems. We formulate the problem of bounding the two-point correlator as a semidefinite programming problem, subject to the constraints of reflection positivity, the Heisenberg equations of motion, and the Kubo-Martin-Schwinger condition or ground-state positivity. In the dual formulation, the Heisenberg equations of motion become "inequalities of motion" on the Lagrange multipliers that enforce the constraints. This enables us to derive rigorous bounds on continuous-time two-point correlators using a finite-dimensional semidefinite or polynomial matrix program. We illustrate this method by bootstrapping the two-point correlators of the ungauged one-matrix quantum mechanics, from which we extract the spectrum and matrix elements of the low-lying adjoint states. Along the way, we provide a new derivation of the energy-entropy balance inequality and establish a connection between the high-temperature two-point correlator bootstrap and the matrix integral bootstrap.