pith. sign in

arxiv: 2601.14141 · v2 · submitted 2026-01-20 · 🧮 math-ph · gr-qc· math.MP

Symmetry Breaking and Phase Transitions in Random Non-Commutative Geometries and Related Random-Matrix Ensembles

Mauro D'Arcangelo , Sven Gnutzmann This is my paper

Pith reviewed 2026-05-16 12:22 UTC · model grok-4.3

classification 🧮 math-ph gr-qcmath.MP
keywords random matrix ensemblesnon-commutative geometryphase transitionssymmetry breakingBarrett-Glaser ensemblesquartic potentialDirac operatorslarge N limit
0
0 comments X

The pith

Quartic one-matrix ensembles of non-commutative geometries exhibit fully analytic phase transitions and symmetry breaking as N tends to infinity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper supplies a complete analytic account of crossovers, phase transitions, and symmetry breaking for one-parameter quartic families of Barrett-Glaser random-matrix ensembles in the two one-matrix cases. It works in the N to infinity limit using saddle-point methods on the underlying Dirac operators and verifies the resulting phase diagram against Monte Carlo data. A sympathetic reader would care because these ensembles function as solvable toy models for quantum gravity, so exact control over their large-scale behavior clarifies when and how random geometries change character. The agreement with numerics indicates that the saddle-point treatment captures the essential physics without missing corrections.

Core claim

For the quartic Barrett-Glaser ensembles restricted to the (1,0) and (0,1) one-matrix cases, the large-N limit is governed by saddle-point equations whose solutions yield an explicit phase diagram in the single coupling constant g, including the locations of all crossovers and symmetry-breaking transitions.

What carries the argument

The large-N saddle-point analysis of the eigenvalue distributions arising from the quartic potential on the finite Dirac operators.

If this is right

  • The phase boundaries and symmetry-breaking patterns are given by explicit analytic expressions in g for both one-matrix cases.
  • Different ranges of g correspond to distinct phases separated by crossovers or sharp transitions.
  • The saddle-point results become exact in the N to infinity limit and already match finite-N Monte Carlo data.
  • The same methods furnish a benchmark for treating higher-matrix or multi-parameter extensions of the ensembles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same saddle-point technique may locate phase structure in multi-matrix or higher-dimensional versions of these geometry models.
  • The identified symmetry-breaking patterns could translate into distinct continuum geometric or topological phases once a suitable scaling limit is taken.
  • Analogous random-matrix constructions with other potentials may admit equally explicit large-N phase diagrams.

Load-bearing premise

The probability measure on the finite Dirac operators remains well-defined and the large-N saddle-point analysis captures every relevant phase transition without non-perturbative effects that would invalidate the analytic expressions.

What would settle it

Monte Carlo simulations at successively larger N that locate a transition or symmetry-breaking point at a value of g different from the saddle-point prediction would falsify the analytic picture.

Figures

Figures reproduced from arXiv: 2601.14141 by Mauro D'Arcangelo, Sven Gnutzmann.

Figure 1
Figure 1. Figure 1: FIG. 1. Equilibrium density (equivalent to expected density of states) for the [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Left: Dependence of the second moment [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Free energy of various candidates for the equilibrium density of states in [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Equilibrium densities in the [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Left: Support of the equilibrium density in the [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Comparison of Monte-Carlo simulations for the equilibrium density at matrix size [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
read the original abstract

Ensembles of random fuzzy non-commutative geometries may be described in terms of finite (\(N^2\)-dimensional) Dirac operators and a probability measure. Dirac operators of type \((p,q)\) are defined in terms of commutators and anti-commutators of \(2^{p+q-1}\) hermitian matrices \(H_k\) and tensor products with a representation of a Clifford algebra. Ensembles based on this idea have recently been used as a toy model for quantum gravity, and they are interesting random-matrix ensembles in their own right. We provide a complete theoretical picture of crossovers, phase transitions, and symmetry breaking in the \(N \to \infty \) limit of 1-parameter families of quartic Barrett-Glaser ensembles in the one-matrix cases \((1,0)\) and \((0,1)\) that depend on one coupling constant \(g\). Our theoretical results are in full agreement with previous and new Monte-Carlo simulations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper analyzes ensembles of random fuzzy non-commutative geometries via finite Dirac operators and associated probability measures. It focuses on one-parameter quartic Barrett-Glaser ensembles in the one-matrix cases (1,0) and (0,1) depending on a single coupling g. The authors derive a complete theoretical description of crossovers, phase transitions, and symmetry breaking in the N to infinity limit using large-N saddle-point analysis, claiming exact agreement with both prior and new Monte-Carlo simulations.

Significance. If the large-N saddle-point results are rigorous and capture the full set of transitions without missing non-perturbative contributions, the work would provide a valuable analytic framework for these random-matrix ensembles and their use as toy models in quantum gravity. Explicit agreement with simulations for the one-parameter families strengthens the claims and offers falsifiable predictions for the order parameters and critical values of g.

major comments (1)
  1. [Large-N saddle-point analysis] The central large-N analysis (detailed in the sections deriving the saddle-point equations for the effective potential) assumes that the saddle-point equations determine all phase transitions and symmetry breaking without additional non-perturbative effects such as eigenvalue tunneling or instantons altering the global minimum for any g. While Monte-Carlo agreement is asserted, the manuscript should include an explicit argument or bound showing that such effects are absent or negligible in the N to infinity limit for the (1,0) and (0,1) cases.
minor comments (2)
  1. [Introduction and setup] Clarify the precise definition of the probability measure on the finite Dirac operators and confirm its normalizability for all g in the one-parameter families, as this underpins the validity of the ensemble.
  2. [Numerical results] In the simulation comparison sections, specify the range of N and g values used in the Monte-Carlo runs and any criteria for identifying phase transitions from the data.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comment on strengthening the large-N analysis. We address the point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Large-N saddle-point analysis] The central large-N analysis (detailed in the sections deriving the saddle-point equations for the effective potential) assumes that the saddle-point equations determine all phase transitions and symmetry breaking without additional non-perturbative effects such as eigenvalue tunneling or instantons altering the global minimum for any g. While Monte-Carlo agreement is asserted, the manuscript should include an explicit argument or bound showing that such effects are absent or negligible in the N to infinity limit for the (1,0) and (0,1) cases.

    Authors: We thank the referee for highlighting this point. In the one-matrix quartic Barrett-Glaser ensembles studied here, the effective potential obtained after integrating out the angular degrees of freedom is a quartic polynomial in a single order parameter (the trace of the matrix or the eigenvalue moment). For all values of the coupling g, the relevant minima are separated by barriers whose height scales as N². Standard large-deviation estimates then imply that the tunneling (instanton) contribution is of order exp(−c N²) and is therefore absent in the strict N → ∞ limit. This is the same mechanism that guarantees the validity of the saddle-point approximation in classical Hermitian matrix models with polynomial potentials. We will add a short paragraph (or subsection) immediately after the derivation of the saddle-point equations, spelling out this scaling argument and citing the relevant random-matrix literature on the suppression of non-perturbative corrections in the planar limit. The existing Monte-Carlo data for successively larger N already show that the measured critical values of g converge to the saddle-point predictions without detectable shifts, consistent with the exponential suppression. revision: yes

Circularity Check

0 steps flagged

Large-N saddle-point derivation of phase transitions is independent of inputs

full rationale

The paper starts from the explicit probability measure on finite Dirac operators of type (p,q) defined via commutators and anticommutators of Hermitian matrices, then applies standard large-N saddle-point analysis to the quartic Barrett-Glaser ensembles in the (1,0) and (0,1) cases. The resulting effective potential and order-parameter equations are derived directly from the action without fitting parameters to the target phase diagram or invoking self-citations as load-bearing uniqueness theorems. Monte-Carlo simulations are performed independently and serve as external validation rather than input. No step reduces the claimed crossovers, transitions, or symmetry breaking to a tautological renaming or fitted quantity by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard large-N saddle-point analysis of random-matrix ensembles with a quartic potential and on the assumption that the Dirac-operator construction yields a well-defined probability measure. No new free parameters beyond the coupling g are introduced; no invented entities are postulated.

free parameters (1)
  • g
    Single coupling constant of the quartic potential; its value controls the location of the phase transitions.
axioms (2)
  • domain assumption Large-N limit of the matrix ensembles admits a saddle-point description that captures all phase transitions.
    Invoked when taking N → ∞ to obtain closed-form expressions for the eigenvalue density and order parameters.
  • domain assumption The finite-dimensional Dirac operators define a valid probability measure on non-commutative geometries.
    Stated in the opening paragraph as the starting point for the ensembles.

pith-pipeline@v0.9.0 · 5508 in / 1491 out tokens · 38884 ms · 2026-05-16T12:22:18.113978+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    The effective potential in this case is always symmetricW −,g(λ) =W −,g(−λ)

    The case of(0,1)geometries We start with the simpler case of(0,1)geometries. The effective potential in this case is always symmetricW −,g(λ) =W −,g(−λ). The variational conditions in the(0,1)case are equivalent to random-matrix models where uniqueness hold. One just has to add as a constraint that the parameterm2 that appears inW −,g(λ)is equal to the se...

    work page

  2. [2]

    Forg→ ∞the density approaches a rescaled semi-circle law

    This extends the existence of these solutions to some negative values of the coupling constant. Forg→ ∞the density approaches a rescaled semi-circle law. Next we consider the 2-cut solutions (see Appendix A for details). Again we can assume a symmetric supportΩ = [−b, a]∪[a, b]with0< a < b. These simplify to G2-cut −,g (z) =− 2(6m2 +g)z+ 8z 3 2πi + 8z 2πi...

    work page

  3. [3]

    In Ap- pendix A we give the general equilibrium densities and their Borel transforms for a 4th-order polynomial potential

    The case of(1,0)geometries In the(1,0)geometries, the effective potential (36) is generally not symmetric. In Ap- pendix A we give the general equilibrium densities and their Borel transforms for a 4th-order polynomial potential. For each of these solutions the boundary of the support are solutions to a set of simultaneous non-linear equations that follow...

    work page

  4. [4]

    In the(1,0)case the additional broken-symmetry solution minimizes the free energy for g < g +,cr. The inset graph shows where the free energies of the symmetric 1-cut solution and the broken symmetry 2-cut solutions cross (the broken symmetry solution exists beyond the crossing which cannot be seen in the large graph due to the width of the curves which h...

    work page

  5. [5]

    5 −2 −1 0 1 2 λ ρ +, −3(λ ) 0

    work page

  6. [6]

    5 −2 −1 0 1 2 λ ρ −, −7(λ ) 0

    work page

  7. [7]

    5 −2 −1 0 1 2 λ ρ +, −4(λ ) FIG. 6. Comparison of Monte-Carlo simulations for the equilibrium density at matrix sizeN= 1024 (blue and red curves) with the theoretical equilibrium density (black) that we obtained using the Riemann-Hilbert approach. All plots are free of any free fitting parameters. Left: 1-cut spectral density in the(1,0)case at coupling c...

    work page 2048

  8. [8]

    J. W. Barrett and L. Glaser, J. Phys. A49, 245001 (2016)

    work page 2016

  9. [9]

    J. W. Barrett, J. Math. Phys.56, 082301 (2015)

    work page 2015

  10. [10]

    M. L. Mehta,Random Matrices(Academic Press, 2004)

    work page 2004

  11. [11]

    P. J. Forrester,Log-Gases and Random Matrices(Princeton University Press, 2010)

    work page 2010

  12. [12]

    Akemann, J

    G. Akemann, J. Baik, and P. Di Francesco, eds.,The Oxford Handbook of Random Matrix Theory(Oxford University Press, 2011)

    work page 2011

  13. [13]

    Haake, S

    F. Haake, S. Gnutzmann, and M. Kuś,Quantum Signatures of Chaos, 4th ed. (Springer, 2018)

    work page 2018

  14. [14]

    Glaser, J

    L. Glaser, J. Phys. A50, 275201 (2017)

    work page 2017

  15. [15]

    J. W. Barrett, P. Druce, and L. Glaser, J. Phys. A52, 275203 (2019)

    work page 2019

  16. [16]

    Khalkhali and N

    M. Khalkhali and N. Pagliaroli, J. Phys. A54, 035202 (2020)

    work page 2020

  17. [17]

    Khalkhali and N

    M. Khalkhali and N. Pagliaroli, J. Math. Phys.63, 053504 (2022)

    work page 2022

  18. [18]

    Hessam, M

    H. Hessam, M. Khalkhali, and N. Pagliaroli, J. Phys. A55, 335204 (2022)

    work page 2022

  19. [19]

    Hessam, M

    H. Hessam, M. Khalkhali, N. Pagliaroli, and L. S. Verhoeven, J. Phys. A55, 413002 (2022)

    work page 2022

  20. [20]

    Azarfar and M

    S. Azarfar and M. Khalkhali, Annales de l’Institut H. Poincaré D11, 409 (2024)

    work page 2024

  21. [21]

    J. W. Barrett, J. Phys. A57, 455201 (2024)

    work page 2024

  22. [22]

    Khalkhali, N

    M. Khalkhali, N. Pagliaroli, and L. S. Verhoeven, J. Math. Phys.66, 053502 (2025)

    work page 2025

  23. [23]

    Jurkiewicz, Phys

    J. Jurkiewicz, Phys. Lett. B261, 260 (1991). 28

    work page 1991

  24. [24]

    Ambjørn, L

    J. Ambjørn, L. Chekhov, C. Kristjansen, and Y. Makeenko, Nucl. Phys. B404, 127 (1993)

    work page 1993

  25. [25]

    Akemann and J

    G. Akemann and J. Ambjørn, J. Phys. A29, L555 (1996)

    work page 1996

  26. [26]

    Akemann, Nucl

    G. Akemann, Nucl. Phys. B482, 403 (1996)

    work page 1996

  27. [27]

    Bonnet, F

    G. Bonnet, F. David, and B. Eynard, J. Phys. A33, 6739 (2000)

    work page 2000

  28. [28]

    Deift,Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach, (Amer- ican Mathematical Soc., 2000)

    P. Deift,Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach, (Amer- ican Mathematical Soc., 2000)

    work page 2000

  29. [29]

    Iso and A

    S. Iso and A. Kavalov, Nucl. Phys. B501, 670 (1997)

    work page 1997

  30. [30]

    Boutet de Monvel, L

    A. Boutet de Monvel, L. Pastur, and M. Shcherbina, J. Stat. Phys.79, 585 (1995)

    work page 1995

  31. [31]

    Johansson, Duke Math

    K. Johansson, Duke Math. J.91, 151 (1996)

    work page 1996

  32. [32]

    Wigner, Ann

    E. Wigner, Ann. Math.62, 548 (1955)

    work page 1955

  33. [33]

    D’Arcangelo,Numerical simulation of random Dirac operators, Phd thesis, University of Nottingham (2022)

    M. D’Arcangelo,Numerical simulation of random Dirac operators, Phd thesis, University of Nottingham (2022)

    work page 2022

  34. [34]

    D’Arcangelo,Random Fuzzy Library(2020), available athttps://github.com/ darcangelomauro/RFL

    M. D’Arcangelo,Random Fuzzy Library(2020), available athttps://github.com/ darcangelomauro/RFL

    work page 2020

  35. [35]

    Duane, A

    S. Duane, A. Kennedy, B. J. Pendleton, and D. Roweth, Physics Letters B195, 216 (1987)

    work page 1987

  36. [36]

    Pastur and M

    L. Pastur and M. Shcherbina, J. Stat. Phys.86, 109 (1997)

    work page 1997