IKKT/IIB Matrix Model

• The IKKT/IIB matrix model is a maximally supersymmetric, zero-dimensional formulation derived from 10D super Yang–Mills theory that defines type IIB superstring theory nonperturbatively.
• It demonstrates spontaneous symmetry breaking where matrix eigenvalues dynamically generate emergent spacetime dimensions and a holographic relation to type IIB supergravity.
• Numerical and analytic methods reveal transitions between fuzzy and commuting phases, underpinning gravity emergence and chiral gauge sector formation.

The IKKT or IIB matrix model is a maximally supersymmetric, zero-dimensional matrix model resulting from the dimensional reduction of ten-dimensional $\mathcal{N}=1$ super Yang–Mills theory to a point. Introduced by Ishibashi, Kawai, Kitazawa, and Tsuchiya in 1996 as a nonperturbative definition of type IIB superstring theory, it is widely regarded as a central framework for studying emergent spacetime, holography, and the quantum structure of string theory. The defining action couples ten Hermitian $N\times N$ matrices and their Majorana–Weyl fermionic partners. In the large-$N$ limit, the model is conjectured to yield a dynamically emergent spacetime, with gravity and string interactions arising from the collective dynamics of matrices. Both Euclidean and Lorentzian versions are actively studied, encompassing analytic, numerical, and holographic analyses.

1. Model Definition, Symmetries, and Action

The IKKT matrix model is formulated with fundamental dynamical variables

• ten Hermitian matrices $A_M$ ($M=0,\dots,9$) transforming as SO(9,1) or SO(10) vectors,
• a Majorana–Weyl spinor $\Psi$ (also realized as $N\times N$ traceless Hermitian matrices), transforming as a chiral spinor of SO(10).

The action is

$S_{\rm IKKT} = -\frac{1}{4g^2} \, \Tr\bigl([A^M, A^N][A_M, A_N]\bigr) - \frac{1}{2g^2} \, \Tr(\bar\Psi\, \Gamma^M [A_M, \Psi])$

where $g$ is the coupling, $\Gamma^M$ denote 10D gamma matrices, and the trace is over $N\times N$ matrices (Ydri, 2017, Anagnostopoulos et al., 2022).

The crucial symmetries are:

• Global SO(9,1) Lorentz (Euclidean models use SO(10)),
• Gauge invariance under $A_M \to U A_M U^{-1}$, $\Psi \to U \Psi U^{-1}$,
• $\mathcal{N}=2$ supersymmetry in 10D (32 real supercharges).

Interpretation of the model posits that the eigenvalues of $A_M$ parametrize spacetime points—i.e., spacetime is not pre-specified, but arises dynamically as an emergent phenomenon (Ydri, 2017, Anagnostopoulos et al., 2022).

The IKKT path integral,

$Z = \int dA\, d\Psi\; e^{-S_{\rm IKKT}}$

is argued to define the nonperturbative completion of type IIB superstrings (Asano, 2024, Anagnostopoulos et al., 2022).

2. Emergent Spacetime and Spontaneous Symmetry Breaking

The central dynamical question concerns emergent dimensions and the spontaneous symmetry breaking (SSB) of the original rotational invariance.

• In Euclidean models, order parameters such as

$T_{\mu\nu} = \frac{1}{N} \Tr(A_\mu A_\nu)$

are studied; SSB occurs if the eigenvalues $\lambda_1 \gg \lambda_{2,\dots,10}$ in the large-$N$ limit (Anagnostopoulos et al., 2020, Anagnostopoulos et al., 2012). Numerical and analytic studies (Gaussian Expansion Method—GEM, Monte Carlo, and Complex Langevin Method—CLM) demonstrate that the SO(10) symmetry is spontaneously broken to SO(3) in the Euclidean model: three extended directions emerge dynamically, with seven compactified (Anagnostopoulos et al., 2020, Kumar et al., 2022, Anagnostopoulos et al., 2022).

• In the Lorentzian model, solutions show dynamical emergence of a single time and three spatial dimensions, with the remaining six "frozen"—indicative of a (3+1) dimensional expanding universe (Nishimura, 2022, Anagnostopoulos et al., 2022). Notably, a precise Lorentz-invariant mass deformation regulates the model and selects real spacetime signature at late times; for sufficiently large mass parameter $\gamma$, the model exhibits real time slices and exponential late-time spatial expansion [a la Nagata–Nishimura, Aarts et al.].
• The SSB is critically dependent on the fermion Pfaffian phase. The purely bosonic model (phase-quenched) does not exhibit SSB, but restoration of the phase yields spontaneous breaking consistent with full supersymmetry (Anagnostopoulos et al., 2020, Kumar et al., 2022, Anagnostopoulos et al., 2012).
• Numerical methods: CLM with carefully chosen fermionic mass deformation avoids singular-drift issues and confirms SO(10)$\to$SO(3) breaking (Anagnostopoulos et al., 2020, Kumar et al., 2022); factorization Monte Carlo provides complementary results (Anagnostopoulos et al., 2012). Earlier studies using "approximate Lorentzian dynamics" yielded partially unphysical vacua, now superseded by CLM analyses (Anagnostopoulos et al., 2022).

3. Gravity and Effective Field Theory from Matrix Dynamics

A major achievement is the demonstration that (3+1)-dimensional gravity—specifically, the Einstein–Hilbert action—emerges at one loop in appropriate backgrounds:

• For backgrounds $$\mathcal{M}^{3,1}\times \mathcal{K} \subset \mathbb{R}^{9,1}$$, where $\mathcal{K}$ is a compact fuzzy internal space, the one-loop effective action contains a term

$$\int d^4x\; \sqrt{|G|} \frac{1}{16\pi G_N} \left( \mathcal{R} + \tfrac34\, \tilde{T}_\mu\tilde{T}^\mu \right)$$

with $G_N$ set by the Kaluza–Klein scale of $\mathcal{K}$ (Steinacker, 2023, Steinacker, 2022). The torsion $\tilde{T}_\mu$ is associated with matrix background structure.

• The effective Newton constant is

$$\frac{1}{G_N} = 2\pi^3 c_\mathcal{K}^2 \rho^{-2} m_\mathcal{K}^2$$

where $c_\mathcal{K}^2$ and $m_\mathcal{K}^2$ parameterize fuzzy extra dimension structure (Steinacker, 2023).

• The vacuum energy does not produce an effective four-dimensional cosmological constant due to its dependence on the symplectic volume form, not the Riemannian measure. Instead, it stabilizes the vacuum modulus of the brane plus compact-factor configuration (Steinacker, 2023).
• Modifications of gravity at large distances are also predicted: additional geometric tensor contributions ("mirage matter") appear in the effective gravitational equations, with non-Ricci-flat infrared modes and dark-matter-like phenomenology (Steinacker, 12 Jan 2026).

4. Compactification, Mass Deformations, and Matrix Model Phases

Explicit compactification and mass/gauge deformations provide controlled mechanisms for breaking rotational symmetry and stabilizing moduli:

• Toroidal compactification with anti-periodic fermion boundary conditions for six spatial directions induces an effective mass matrix at one loop: $S_{\rm eff}^0 = -\frac{1}{4g_{\rm eff}^2} \Tr([A^M, A^N]^2) + \frac12 M_{MN}^2 \Tr(A^M A^N) + \cdots$ with block-diagonal $M_{MN}^2$ breaking SO(1,9) to SO(1,3)$\times$SO(6) (Laliberte, 2024). The mechanism single out four large directions and six compactified ones, providing a dynamical symmetry-breaking origin for observed spacetime structure.
• Supersymmetry-preserving mass (or Myers-type) deformations lead to the so-called polarized IKKT model. In the large deformation (mass) limit, the vacuum is dominated by a fuzzy-sphere solution, corresponding to polarized D-instantons (D1-brane emergent) in a background flux (Hartnoll et al., 2024, Hartnoll et al., 8 Apr 2025). At intermediate or low mass, the dominant matrix configurations can transition in a first-order manner between fuzzy-sphere and commuting-matrix phases; this phase structure is sharply visible in the associated partition function and Monte Carlo studies (Hartnoll et al., 8 Apr 2025, Chou et al., 24 Jul 2025).
• Numerical simulations for small $N$ via Monte Carlo and at larger $N$ via localization and analytic methods confirm that the original IKKT vacuum is smoothly connected to the fuzzy sphere solution as deformation increases, with noncommutative matrix structure giving way to commuting (clustered instanton) phases at small mass (Chou et al., 24 Jul 2025).

5. Holography and Dual Supergravity Realizations

The IKKT model is now firmly integrated into the landscape of string holography:

• It is the $p=-1$ endpoint of the D$p$/SYM family, dual (in the large-$N$ and strong-coupling limit) to type IIB supergravity on the near-horizon background of stacked D-instantons (Ciceri et al., 11 Mar 2025, Ciceri et al., 28 Nov 2025).
• The holographic dictionary pairs matrix operators—traces of matrices and their commutators—with Kaluza–Klein fluctuations of IIB fields on $S^9$. Notably, the protected $\mathcal{B}_2$ multiplet contains trace bilinears, fermionic traces, and commutator traces matched to SO(10) harmonics (Ciceri et al., 11 Mar 2025, Ciceri et al., 28 Nov 2025).
• The one-dimensional maximal $\mathrm{SO}(10)$ gauged supergravity Lagrangian, constructed as the dual bulk theory, describes all nonlinear dynamics of these KK modes, with explicit bosonic and fermionic content corresponding to the protected multiplet (Ciceri et al., 28 Nov 2025). Half-BPS solutions of the 1D supergravity describe flows to various symmetry-breaking vacua, whose uplifts yield families of ten-dimensional backgrounds (including the D(–1) instanton, spherical brane solutions, and ultimately 12D pp-waves) (Ciceri et al., 11 Mar 2025, Ciceri et al., 28 Nov 2025).
• This framework provides a setting for "timeless" holography: since the IKKT model lacks an explicit time, geometrical and holographic flows are realized in terms of matrix eigenvalue and moduli space distributions, rather than geometric radial directions (Hartnoll et al., 2024, Chou et al., 24 Jul 2025, Ciceri et al., 11 Mar 2025).

6. Matrix Model Regularization of IIB Superstrings

From a first-principles string theory perspective, the IKKT model emerges as the matrix regularization of type IIB string theory in Schild gauge. Rigorous path-integral analyses demonstrate the precise equivalence between the Polyakov, Schild, and matrix-model (IKKT) formulations, and between Minkowskian/Euclidean signatures (Asano, 2024). In the Schild gauge, quantization and matrix discretization lead directly to the (Euclidean or Lorentzian) IKKT model, with the matrix path integral realizing causality via eigenvalue prescriptions and s-sector summations (Asano, 2024).

7. Phenomenology, Chiral Matter, and Model Building

The IKKT framework enables first-principles constructions of phenomenologically relevant backgrounds:

• Noncommutative tori, realized as infinite-dimensional matrices, support Dirac operators whose zero-modes in magnetic flux backgrounds possess chirality and generation structures analogous to commutative toroidal compactifications (Honda, 2019).
• Yukawa couplings among chiral matter fields can be computed exactly as traces of matrix products, with explicit theta-function dependence, mirroring intersecting D-brane constructions (Honda, 2019). This provides a fully background-independent origin for chiral gauge sectors from the underlying matrix dynamics.

---

In summary, the IKKT or IIB matrix model provides a mathematically precise, nonperturbative formulation of type IIB superstring theory in which both geometry and field content—spacetime, gravity, and gauge/matter spectra—emerge dynamically from matrix degrees of freedom. It unifies string field theory, gauge/gravity duality, and the theory of emergent geometry, and, through a rich network of analytic, numerical, and holographic results, has become a central platform for the study of quantum spacetime, symmetry breaking, and nonperturbative string dynamics (Ydri, 2017, Anagnostopoulos et al., 2022, Anagnostopoulos et al., 2020, Steinacker, 2023, Laliberte, 2024, Hartnoll et al., 2024, Ciceri et al., 11 Mar 2025, Chou et al., 24 Jul 2025, Steinacker, 12 Jan 2026, Ciceri et al., 28 Nov 2025, Asano, 2024, Honda, 2019).