
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
twosphere and its noncommutative gauge theory. This provides a novel and a
much simpler alternative to Connes noncommutative 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
YangMills 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 foursphere ${\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 YangMills 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 foursphere
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 nonperturbative 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 phifour 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 disorderedtononuniformordered transition. The WilsonFisher
fixed point at $\theta=0$ which describes the disorderedtouniformordered
transition, and a noncommutative WilsonFisher fixed point at a maximum value
of $\theta$ which is associated with the transition between nonuniformorder
and uniformorder 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 MoyalWeyl plane ${\bf R}^{2}_{\theta, \Omega}$,
with a nonzero harmonic oscillator term added. The doubletrace matrix model
symmetric under $M\longrightarrow M$ is solved in closed form. An analytical
prediction for the disorderedtononuniformordered 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
YangMills matrix models is briefly reviewed. In particular, the eigenvalues
distributions of YangMills matrix models in lower dimensions in the commuting
(matrix or YangMills) phase of these models are discussed.

We propose a new algorithm for simulating noncommutative phifour 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 phifour 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 YangMills matrix model with a Myers term. As an application,
the phase diagram, including the triple point, of noncommutative phifour
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
nonzero 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 phifour 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 noncommutative vector sigma model with the most
general \phi^4 interaction on MoyalWeyl spaces. We compute the 2 and 4point
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 noncommutative WilsonFisher fixed point
interpolates between the commutative WilsonFisher 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 noncommutativity beyond which the twosheeted
structure of a_* as a function of the dilation parameter disappears.

We present a study of D=4 supersymmetric YangMills 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 nonpolynomial 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 nonpolynomial law.

We present a study of D=4 supersymmetric YangMills 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 nonpolynomial 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 nonpolynomial law. We also
discuss the D=3 models and by using cohomological deformation, localization
techniques and the saddlepoint 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 UVIR 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 twosphere 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
twosphere, 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 nonperturbative 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 UVIR 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 0dimensional 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) YangMills 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 weakcoupling and strongcoupling
phases of the gauge theory. In the ``matrix'' phase we have a U(N) gauge theory
on a single point.