We present an analytic proof of the existence of phase transition in the large N limit of certain random noncommutative geometries. These geometries can be expressed as ensembles of Dirac operators. When they reduce to single matrix ensembles, one can apply the Coulomb gas method to find the empirical spectral distribution. We elaborate on the nature of the large N spectral distribution of the Dirac operator itself. Furthermore, we show that these models exhibit both a single and double cut region for certain values of the order parameter and find the exact value where the transition occurs.

We apply the bootstrap technique to find the moments of certain multi-trace and multi-matrix random matrix models suggested by noncommutative geometry. Using bootstrapping we are able to find the relationships between the coupling constant of these models and their second moments. Using the Schwinger-Dyson equations, all other moments can be expressed in terms of the coupling constant and the second moment. Explicit relations for higher mixed moments are obtained.

In this paper, we find spectral properties in the large N limit of Dirac operators that come from random finite noncommutative geometries. In particular, for a Gaussian potential, the limiting eigenvalue spectrum is shown to be universal, regardless of the geometry, and is given by the convolution of the semicircle law with itself. For simple non-Gaussian models, this convolution property is also evident. In order to prove these results, we show that a wide class of multi-trace multimatrix models have a genus expansion.

We review recent progress in the analytic study of random matrix models suggested by noncommutative geometry. One considers fuzzy spectral triples where the space of possible Dirac operators is assigned a probability distribution. These ensembles of Dirac operators are constructed as toy models of Euclidean quantum gravity on finite noncommutative spaces and display many interesting properties. The ensembles exhibit spectral phase transitions, and near these phase transitions they show manifold-like behavior. In certain cases one can recover Liouville quantum gravity in the double scaling limit. We highlight examples where bootstrap techniques, Coulomb gas methods, and Topological Recursion are applicable.

In this paper we study ensembles of finite real spectral triples equipped with a path integral over the space of possible Dirac operators. In the noncommutative geometric setting of spectral triples, Dirac operators take the center stage as a replacement for a metric on a manifold. Thus, this path integral serves as a noncommutative analogue of integration over metrics, a key feature of a theory of quantum gravity. From these integrals in the so-called double scaling limit we derive critical exponents of minimal models from Liouville conformal field theory coupled with gravity. Additionally, the asymptotics of the partition function of these models satisfy differential equations such as Painlevé I, as a reduction of the KDV hierarchy, which is predicted by conformal field theory. This is all proven using well-established and rigorous techniques from random matrix theory.