Bootstrapping the critical behavior of multi-matrix models
Masoud Khalkhali, Nathan Pagliaroli, Andrei Parfeni, Brayden Smith
TL;DR
The authors introduce bootstrapping with positivity to bound large-$N$ moments in multi-matrix models by coupling Schwinger-Dyson relations with positivity constraints from the Hamburger moment problem. They validate the approach on the quartic single-matrix model and apply it to several unsolved 2-matrix models, obtaining evidence for a string susceptibility exponent $\gamma=1/2$ in some cases (suggesting the Continuum Random Tree) and identifying other exponents that may signal new continuum limits. They also derive explicit first terms of the large-$N$ free energy as a power series in the coupling constants and study a 3-matrix model, demonstrating the method’s utility for extracting critical behavior and generating series expansions. The work broadens analytical access to multi-matrix models, offering a practical route to characterize phase structure and potential continuum limits, with potential extensions to tensor models and map enumerations."
Abstract
Given a matrix model, by combining the Schwinger-Dyson equations with positivity constraints on its solutions, in the large $N$ limit one is able to obtain explicit and numerical bounds on its moments. This technique is known as bootstrapping with positivity. In this paper we use this technique to estimate the critical points and exponents of several multi-matrix models. As a proof of concept, we first show it can be used to find the well-studied quartic single matrix model's critical phenomena. We then apply the method to several similar ``unsolved" 2-matrix models with various quartic interactions. We conjecture and present strong evidence for the string susceptibility exponent for some of these models to be $γ= 1/2$, which heuristically indicates that the continuum limit will likely be the Continuum Random Tree. For the other 2-matrix models, we find estimates of new string susceptibility exponents that may indicate a new continuum limit. We then study an unsolved 3-matrix model that generalizes the 3-colour model with cubic interactions. Additionally, for all of these models, we are able to derive explicitly the first several terms of the free energy in the large $N$ limit as a power series expansion in the coupling constants at zero by exploiting the structure of the Schwinger-Dyson equations.
