A review of M-(atrix) theory (the BFFS matrix quantum mechanics), type IIB matrix model (the IKKT matrix model) and Matrix String Theory (the DVV matrix gauge theory) is presented.

This article is divided into three parts. First, a systematic derivation of the Hawking radiation is given in three different ways. The information loss problem is then discussed in great detail. The last part contains a concise discussion of black hole thermodynamics. This article was published as chapter $6$ of the IOP book "Lectures on General Relativity, Cosmology and Quantum Black Holes" (July $2017$).

We present a new model of quantum gravity as a theory of random geometries given explicitly in terms of a multitrace matrix model. This is a generalization of the usual discretized random surfaces of 2D quantum gravity which works away from two dimensions and captures a large class of spaces admiting a finite spectral triple. These multitrace matrix models sustain emergent geometry as well as growing dimensions and topology change.

We present a new multitrace matrix model, which is a generalization of the real quartic one matrix model, exhibiting dynamical emergence of a fuzzy two-sphere and its non-commutative gauge theory. This provides a novel and a much simpler alternative to Connes non-commutative geometry and to the IKKT matrix model for emergent geometry in two dimensions.

A detailed Monte Carlo calculation of the phase diagram of bosonic IKKT Yang-Mills matrix models in three and six dimensions with quartic mass deformations is given. Background emergent fuzzy geometries in two and four dimensions are observed with a fluctuation given by a noncommutative $U(1)$ gauge theory very weakly coupled to normal scalar fields. The geometry, which is determined dynamically, is given by the fuzzy spheres ${\bf S}^2_N$ and ${\bf S}^2_N\times{\bf S}^2_N$ respectively. The three and six matrix models are in the same universality class with some differences. For example, in two dimensions the geometry is completely stable, whereas in four dimensions the geometry is stable only in the limit $M\longrightarrow \infty$, where $M$ is the mass of the normal fluctuations. The behavior of the eigenvalue distribution in the two theories is also different. We also sketch how we can obtain a stable fuzzy four-sphere ${\bf S}^2_N\times{\bf S}^2_N$ in the large $N$ limit for all values of $M$ as well as models of topology change in which the transition between spheres of different dimensions is observed. The stable fuzzy spheres in two and four dimensions act precisely as regulators which is the original goal of fuzzy geometry and fuzzy physics. Fuzzy physics and fuzzy field theory on these spaces are briefly discussed.

We study a six matrix model with global $SO(3)\times SO(3)$ symmetry containing at most quartic powers of the matrices. This theory exhibits a phase transition from a geometrical phase at low temperature to a Yang-Mills matrix phase with no background geometrical structure at high temperature. This is an exotic phase transition in the same universality class as the three matrix model but with important differences. The geometrical phase is determined dynamically, as the system cools, and is given by a fuzzy four-sphere background ${\bf S}^2_N\times{\bf S}^2_N$, with an Abelian gauge field which is very weakly coupled to two normal scalar fields playing the role of dark energy.

This book is divided into two parts. In the first part we give an elementary introduction to computational physics consisting of 21 simulations which originated from a formal course of lectures and laboratory simulations delivered since 2010 to physics students at Annaba University. The second part is much more advanced and deals with the problem of how to set up working Monte Carlo simulations of matrix field theories which involve finite dimensional matrix regularizations of noncommutative and fuzzy field theories, fuzzy spaces and matrix geometry. The study of matrix field theory in its own right has also become very important to the proper understanding of all noncommutative, fuzzy and matrix phenomena. The second part, which consists of 9 simulations, was delivered informally to doctoral students who are working on various problems in matrix field theory. Sample codes as well as sample key solutions are also provided for convenience and completness. An appendix containing an executive arabic summary of the first part is added at the end of the book.

This book includes my lectures, together with their problem sets and solutions, on 1) classical mechanics (one semester), 2) thermodynamics and statistical mechanics (one semester), and 3) quantum mechanics (one semester), which I have been giving to graduate students of theoretical physics at Annaba University since 2010 .

The subject of matrix field theory involves matrix models, noncommutative geometry, fuzzy physics and noncommutative field theory and their interplay. In these lectures, a lot of emphasis is placed on the matrix formulation of noncommutative and fuzzy spaces, and on the non-perturbative treatment of the corresponding field theories. In particular, the phase structure of noncommutative $\phi^4$ theory is treated in great detail, and an introduction to noncommutative gauge theory is given.

We present a study of phi-four theory on noncommutative spaces using a combination of the Wilson renormalization group recursion formula and the solution to the zero dimensional vector/matrix models at large $N$. Three fixed points are identified. The matrix model $\theta=\infty$ fixed point which describes the disordered-to-non-uniform-ordered transition. The Wilson-Fisher fixed point at $\theta=0$ which describes the disordered-to-uniform-ordered transition, and a noncommutative Wilson-Fisher fixed point at a maximum value of $\theta$ which is associated with the transition between non-uniform-order and uniform-order phases.

In this article we provide a multitrace analysis of the theory of noncommutative $\Phi^4$ in two dimensions on the fuzzy sphere ${\bf S}^2_{N,\Omega}$, and on the Moyal-Weyl plane ${\bf R}^{2}_{\theta, \Omega}$, with a non-zero harmonic oscillator term added. The doubletrace matrix model symmetric under $M\longrightarrow -M$ is solved in closed form. An analytical prediction for the disordered-to-non-uniform-ordered phase transition and an estimation of the triple point, from the termination point of the critical boundary, are derived and compared with previous Monte Carlo measurement.

The phenomena of emergent fuzzy geometry and noncommutative gauge theory from Yang-Mills matrix models is briefly reviewed. In particular, the eigenvalues distributions of Yang-Mills matrix models in lower dimensions in the commuting (matrix or Yang-Mills) phase of these models are discussed.

We propose a new algorithm for simulating noncommutative phi-four theory on the fuzzy sphere based on, i) coupling the scalar field to a U(1) gauge field, in such a way that in the commutative limit N\longrightarrow \infty, the two modes decouple and we are left with pure scalar phi-four on the sphere, and ii) diagonalizing the scalar field by means of a U(N) unitary matrix, and then integrating out the unitary group from the partition function. The number of degrees of freedom in the scalar sector reduces, therefore, from N^2 to the N eigenvalues of the scalar field, whereas the dynamics of the U(1) gauge field, is given by D=3 Yang-Mills matrix model with a Myers term. As an application, the phase diagram, including the triple point, of noncommutative phi-four theory on the fuzzy sphere, is reconstructed with small values of N up to N=10, and large numbers of statistics.

We present a (hopefully) novel calculation of the vacuum energy in expanding FLRW spacetimes based on the renormalization of quantum field theory in non-zero backgrounds. We compute the renormalized effective action up to the $2-$point function and then apply the formalism to the cosmological backgrounds of interest. As an example we calculate for quasi de Sitter spacetimes the leading correction to the vacuum energy given by the tadpole diagram and show that it behaves as $\sim H_0^2 \Lambda_{\rm pl}$ where $H_0$ is the Hubble constant and $\Lambda_{\rm pl}$ is the Planck constant. This is of the same order of magnitude as the observed dark energy density in the universe.

In this article we exhibit explicitly the matrix model ($\theta=\infty$) fixed point of phi-four theory on noncommutative spacetime with only two noncommuting directions using the Wilson renormalization group recursion formula and the 1/N expansion of the zero dimensional reduction and then calculate the mass critical exponent $\nu$ and the anomalous dimension $\eta$ in various dimensions .

In this article we study non-commutative vector sigma model with the most general \phi^4 interaction on Moyal-Weyl spaces. We compute the 2- and 4-point functions to all orders in the large N limit and then apply the approximate Wilson renormalization group recursion formula to study the renormalized coupling constants of the theory. The non-commutative Wilson-Fisher fixed point interpolates between the commutative Wilson-Fisher fixed point of the Ising universality class which is found to lie at zero value of the critical coupling constant a_* of the zero dimensional reduction of the theory, and a novel strongly interacting fixed point which lies at infinite value of a_* corresponding to maximal non-commutativity beyond which the two-sheeted structure of a_* as a function of the dilation parameter disappears.

We present a study of D=4 supersymmetric Yang-Mills matrix models with SO(3) mass terms based on the Monte Carlo method. In the bosonic models we show the existence of an exotic first/second order transition from a phase with a well defined background geometry (the fuzzy sphere) to a phase with commuting matrices with no geometry in the sense of Connes. At the transition point the sphere expands abruptly to infinite size then it evaporates as we increase the temperature (the gauge coupling constant). The transition looks first order due to the discontinuity in the action whereas it looks second order due to the divergent peak in the specific heat. The fuzzy sphere is stable for the supersymmetric models in the sense that the bosonic phase transition is turned into a very slow crossover transition. The transition point is found to scale to zero with N. We conjecture that the transition from the background sphere to the phase of commuting matrices is associated with spontaneous supersymmetry breaking. The eigenvalues distribution of any of the bosonic matrices in the matrix phase is found to be given by a non-polynomial law obtained from the fact that the joint probability distribution of the four matrices is uniform inside a solid ball with radius R. The eigenvalues of the gauge field on the background geometry are also found to be distributed according to this non-polynomial law.

We present a study of D=4 supersymmetric Yang-Mills matrix models with SO(3) mass terms based on the cohomological approach and the Monte Carlo method. In the bosonic models we show the existence of an exotic first/second order transition from a phase with a well defined background geometry (the fuzzy sphere) to a phase with commuting matrices with no geometry in the sense of Connes. At the transition point the sphere expands abruptly to infinite size then it evaporates as we increase the temperature (the gauge coupling constant). The transition looks first order due to the discontinuity in the action whereas it looks second order due to the divergent peak in the specific heat. The fuzzy sphere is stable for the supersymmetric models in the sense that the bosonic phase transition is turned into a very slow crossover transition. The transition point is found to scale to zero with N. We conjecture that the transition from the background sphere to the phase of commuting matrices is associated with spontaneous supersymmetry breaking. The eigenvalues distribution of any of the bosonic matrices in the matrix phase is found to be given by a non-polynomial law obtained from the fact that the joint probability distribution of the four matrices is uniform inside a solid ball with radius R. The eigenvalues of the gauge field on the background geometry are also found to be distributed according to this non-polynomial law. We also discuss the D=3 models and by using cohomological deformation, localization techniques and the saddle-point method we give a derivation of the D=3 eigenvalues distribution starting from a particular D=4 model.

We give a brief review of two nonperturbative phenomena typical of noncommutative field theory which are known to lead to the perturbative instability known as the UV-IR mixing. The first phenomena concerns the emergence/evaporation of spacetime geometry in matrix models which describe perturbative noncommutative gauge theory on fuzzy backgrounds. In particular we show that the transition from a geometrical background to a matrix phase makes the description of noncommutative gauge theory in terms of fields via the Weyl map only valid below a critical value g_*. The second phenomena concerns the appearance of a nonuniform ordered phase in noncommutative scalar \phi^4 field theory and the spontaneous symmetry breaking of translational/rotational invariance which happens even in two dimensions. We argue that this phenomena also originates in the underlying matrix degrees of freedom of the noncommutative field theory. Furthermore it is conjectured that in addition to the usual WF fixed point at $\theta=0$ there must exist a novel fixed point at \theta=\infty corresponding to the quartic hermitian matrix model.

We present, theoretical predictions and Monte Carlo simulations, for a simple three matrix model that exhibits an exotic phase transition. The nature of the transition is very different if approached from the high or low temperature side. The high temperature phase is described by three self interacting random matrices with no background spacetime geometry. As the system cools there is a phase transition in which a classical two-sphere condenses to form the background geometry. The transition has an entropy jump or latent heat, yet the specific heat diverges as the transition is approached from low temperatures. We find no divergence or evidence of critical fluctuations when the transition is approached from the high temperature phase. At sufficiently low temperatures the system is described by small fluctuations, on a background classical two-sphere, of a U(1) gauge field coupled to a massive scalar field. The critical temperature is pushed upwards as the scalar field mass is increased. Once the geometrical phase is well established the specific heat takes the value 1 with the gauge and scalar fields each contributing 1/2.

We study a three matrix model with global SO(3) symmetry containing at most quartic powers of the matrices. We find an exotic line of discontinuous transitions with a jump in the entropy, characteristic of a 1st order transition, yet with divergent critical fluctuations and a divergent specific heat with critical exponent $\alpha=1/2$. The low temperature phase is a geometrical one with gauge fields fluctuating on a round sphere. As the temperature increased the sphere evaporates in a transition to a pure matrix phase with no background geometrical structure. Both the geometry and gauge fields are determined dynamically. It is not difficult to invent higher dimensional models with essentially similar phenomenology. The model presents an appealing picture of a geometrical phase emerging as the system cools and suggests a scenario for the emergence of geometry in the early universe.

In these notes we review Klimcik's construction of noncommutative gauge theory on the fuzzy supersphere. This theory has an exact SUSY gauge symmetry with a finite number of degrees of freedom and thus in principle it is amenable to the methods of matrix models and Monte Carlo numerical simulations. We also write down in this article a novel fuzzy supersymmetric scalar action on the fuzzy supersphere.

In this letter we show that supersymmetry like geometry can be approximated using finite dimensional matrix models and fuzzy manifolds. In particular we propose a non-perturbative regularization of {\cal N}=2 supersymmetric U(n) gauge action in 4D. In some planar large N limits we recover exact SUSY together with the smooth geometry of R^4_{theta}.

Many gauge theory models on fuzzy complex projective spaces will contain a strong instability in the quantum field theory leading to topology change. This can be thought of as due to the interaction between spacetime via its noncommutativity and the fields (matrices) and it is related to the perturbative UV-IR mixing. We work out in detail the example of fuzzy CP^2 and discuss at the level of the phase diagram the quantum transitions between the 3 spaces (spacetimes) CP^2, S^2 and the 0-dimensional space consisting of a single point {0}.

We construct noncommutative U(1) gauge theory on the fuzzy sphere S^2_N as a unitary 2N x 2N matrix model. In the quantum theory the model is equivalent to a nonabelian U(N) Yang-Mills theory on a 2 dimensional lattice with 2 plaquettes. This equivalence holds in the " fuzzy sphere" phase where we observe a 3rd order phase transition between weak-coupling and strong-coupling phases of the gauge theory. In the ``matrix'' phase we have a U(N) gauge theory on a single point.