🔬 mesoscale physics

Interacting Chern insulator transition on the sphere: revealing the Gross-Neveu-Yukawa criticality

Using exact diagonalization of Dirac fermions on a sphere to leverage full $SO(3)$ symmetry and circumvent the parity anomaly, this study demonstrates that the interaction-driven multicritical point of the two-dimensional Chern insulator transition belongs to the $N=2$ Gross-Neveu-Yukawa universality class, with extracted operator scaling dimensions showing strong agreement with conformal bootstrap predictions.

Original authors: Zhi-Qiang Gao, Taige Wang, Dung-Hai Lee

Published 2026-08-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zhi-Qiang Gao, Taige Wang, Dung-Hai Lee

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world of solid materials, electrons usually settle into predictable patterns. Sometimes they flow freely like a river, creating electricity; other times they get stuck in place, turning the material into an insulator. For decades, physicists understood these changes through the lens of symmetry breaking, a concept where a material loses a specific symmetry as it cools down, much like water freezing into ice crystals. However, a newer class of materials defies this old rule. These are topological insulators, where the electrons arrange themselves in a way that is fundamentally different from ordinary insulators, not because they broke a symmetry, but because their internal structure twisted into a new shape. This twist creates a special state where electricity can flow perfectly along the edges while the interior remains blocked. The moment a material switches from this twisted, topological state to a normal, non-twisted state is a mystery that has puzzled scientists. It is a transition that happens without any obvious symmetry breaking, making it a clean example of a "beyond-Landau" phase change, a frontier where our standard theories of matter struggle to keep up.

A team of researchers at the University of California, Berkeley, has now taken a significant step toward solving this puzzle by simulating the transition in a controlled, mathematical environment. They focused on a specific type of topological material known as a Chern insulator, which can be tuned from a topological state to an ordinary insulator by adjusting a single control knob. The challenge has always been that when electrons interact with each other strongly, the math becomes incredibly difficult, and standard computer simulations often fail or produce misleading results due to inherent flaws in how they are built. To bypass these obstacles, the researchers placed their model of interacting electrons onto the surface of a sphere. This curved geometry acts as a natural shield against certain mathematical errors that plague flat-grid simulations, allowing them to study the system with high precision using a method called exact diagonalization, which solves the quantum equations directly for small groups of electrons.

By running these calculations on spheres with different sizes, the team was able to observe the exact moment the material undergoes its phase transition. They found that when the electrons interact strongly, the transition does not follow the simple, free-particle rules one might expect. Instead, it settles into a complex, critical state known as the Gross–Neveu–Yukawa fixed point. This state is a benchmark problem in physics, appearing in diverse fields from the study of subatomic particles in high-energy physics to the behavior of electrons in graphene. The researchers were able to map out the full spectrum of excitations, or the specific ways the system can vibrate and fluctuate, right at this critical point. Their results matched the most advanced theoretical predictions available, confirming that their spherical setup accurately captures the physics of this elusive transition.

Perhaps the most exciting part of their discovery was finding new types of collective behaviors that had never been seen in previous studies. In the language of quantum physics, these are called "operators," which describe how the system responds to changes. The team identified several new modes, including a nearly conserved current that behaves almost like a perfect flow of charge, and a specific pairing of electrons that acts like a high-spin wave. These findings are significant because they provide a detailed, numerical confirmation of theories that were previously only guessed at through approximations. The researchers also uncovered a new type of electron pairing that sits very close to a theoretical limit, suggesting it is a nearly conserved quantity, a feature that had been inferred from other mathematical techniques but never directly observed in a simulation of this kind.

The study demonstrates that placing a model of interacting electrons on a curved surface is a surprisingly powerful tool. It offers a compact and transparent way to look at strongly interacting topological transitions, cutting through the noise that usually obscures these phenomena. The team's work confirms that the transition is governed by a specific set of rules that govern how the electrons and their collective fields behave together. While the simulations were limited to systems with fewer than ten electrons, the patterns they revealed are robust and align perfectly with the best theoretical bounds known to science. This suggests that the spherical approach is not just a mathematical trick, but a genuine window into the deep structure of matter.

Looking forward, this work opens the door to exploring even more complex scenarios. The researchers suggest that with larger systems and different types of particles, this method could be used to extract even more detailed data about how these materials work. It could help scientists investigate whether certain critical points in nature possess hidden symmetries, such as supersymmetry, which would link particles of different types in a profound way. For now, the study stands as a clear, concrete demonstration that by changing the geometry of the problem, physicists can finally see the hidden machinery behind some of the most mysterious phase transitions in the universe. The transition from a topological insulator to a normal one is no longer just a theoretical concept; it has been mapped, measured, and understood in a way that brings us closer to a complete theory of how matter changes its fundamental nature.

Technical Summary: Interacting Chern Insulator Transition on the Sphere

Problem Statement
The transition from a trivial band insulator to a Chern insulator in two dimensions represents a prototypical "beyond-Landau" criticality, where topology changes without symmetry breaking. The low-energy physics of this transition is governed by a single Dirac cone. While the non-interacting limit ( $u_0=0$ ) is described by a free-fermion conformal field theory (CFT), electron interactions ( $u_0 > 0$ ) drive the system toward a Gross–Neveu–Yukawa (GNY) fixed point. This multicritical point, where fermions and the mass operator simultaneously become gapless, is a benchmark problem of broad interest in high-energy and condensed-matter physics.

However, obtaining a faithful numerical treatment of an interacting single Dirac cone is notoriously difficult due to four primary obstacles:

Parity Anomaly: Local, time-reversal-symmetric lattice regularizations suffer from the parity anomaly, often requiring enlarged flavor numbers or broken symmetries that obscure the intended universality class.
Strong Correlation: The critical interacting Dirac fermion cannot be reached perturbatively from the free Dirac CFT because the four-fermion interaction is irrelevant at the free fixed point.
Symmetry Limitations: Without $SO(3)$ symmetry, most microscopic calculations can only extract critical exponents from order parameters rather than utilizing state-operator correspondence to determine the full operator content of the CFT.
Dispersion Issues: Standard fuzzy sphere constructions, which preserve $SO(3)$ , typically generate Landau levels due to the monopole at the sphere center, failing to reproduce the target linear Dirac dispersion.

Methodology
The authors propose a framework that formulates the interacting Dirac cone directly on the surface of a sphere ( $S^2$ ) to circumvent these obstacles.

Geometric Setup: By placing the Dirac operator on a sphere of radius $R=1$ , the necessary spin connection is supplied by the curvature. The Dirac operator takes the form $D = \sigma \cdot L + 1$ , where $L$ is orbital angular momentum and $\sigma$ is spin. This setup preserves full rotational $SO(3)$ symmetry and discrete symmetries (Time-reversal $T$ , Charge conjugation $C$ , and Parity $P$ ).
Hamiltonian Construction: The system is described by a Hamiltonian $H = H_0 + H_{int}$ . The interaction term preserves $SO(3)$ and discrete symmetries by utilizing a local spin operator $\sigma_\Omega = \hat{n}_\Omega \cdot \sigma$ aligned with the radial direction. The interaction potential is expanded in Haldane pseudopotentials $U_l$ , allowing for contact interactions ( $U_0$ ) and higher-order terms ( $U_1, U_2$ ).
Exact Diagonalization (ED): The authors perform ED on systems with a finite angular momentum cutoff $j_{max}$ . Due to the rapid growth of the Hilbert space with $j_{max}$ , the study focuses on $j_{max} = 3/2$ and $j_{max} = 5/2$ , corresponding to systems with fewer than ten electrons.
Spectral Analysis: Using radial quantization, the energy spectrum on $S^2 \times \mathbb{R}$ is mapped to the scaling dimensions of CFT operators via $\delta E = v \Delta$ . The velocity $v$ is calibrated by matching the scaling dimension of the lowest scalar operator $\sigma$ to the conformal bootstrap value ( $\Delta_\sigma = 0.66$ ). The pseudopotentials $(U_0, U_1, U_2)$ are fine-tuned to align the scaling dimensions of the Dirac fermion $\psi$ and its descendants with bootstrap predictions.

Key Contributions and Results

Realization of GNY Criticality: The study demonstrates that a robust GNY critical point can be realized in a modest Hilbert space on a sphere. By tuning the interaction parameters, the authors pinpoint a critical regime where the system exhibits the expected conformal symmetry.
Conformal Spectrum Extraction: The ED spectrum reveals eight distinct conformal multiplets spanning charge sectors $q=0, 1, 2$ . The leading scaling dimensions for the mass operator $\sigma$ , the Dirac spinor $\psi$ , and the singlet pairing operator $\sigma_T$ match state-of-the-art conformal bootstrap bounds within a few percent (e.g., $\sim 3\%$ error for $\sigma_T$ at $j_{max}=5/2$ ).
Discovery of New Primaries: The authors identify two previously unseen primary operators in the GNY spectrum:
- $\tilde{J}_\mu$ : A non-conserved current with the same quantum numbers as the conserved current $J_\mu$ but admitting derivative descendants.
- $D_{\mu\nu}$ : A spin-2 operator in the $q=2$ sector with a scaling dimension $\Delta \approx 3.03$ (at $j_{max}=5/2$ ). This value nearly saturates the unitarity bound $\Delta = \ell + 1$ for spin- $\ell$ operators in $(2+1)$ D CFTs.
Evidence for Weakly Broken Higher-Spin Symmetries: The near-saturation of the unitarity bound by $D_{\mu\nu}$ , combined with the absence of lower-spin descendants, provides direct numerical evidence that $D_{\mu\nu}$ acts as a nearly conserved current. The authors conjecture this corresponds to an unconventional $p$ -wave pairing operator $\psi^T R (\gamma_\mu \partial_\nu + \gamma_\nu \partial_\mu) \psi$ , representing a weakly broken higher-spin symmetry expected in the GNY CFT.
Validation of Spherical ED: The results confirm that the low-lying excitations form representations of the special-conformal group, verifying the emergence of 3D conformal symmetry. The method successfully reproduces the restricted descendant structures of conserved operators (where $\partial_\mu J^\mu = 0$ forbids lower-spin derivatives) and the characteristic structure of the spinor multiplet.

Significance
The paper claims that placing a Dirac cone on a curved geometry provides a "surprisingly simple and accurate setup" for studying strongly interacting topological transitions. By generalizing the fuzzy sphere framework to interacting Dirac fermions, this approach offers a compact and transparent numerical microscope for the GNY fixed point. It resolves long-standing difficulties in simulating single-flavor Dirac fermions by naturally avoiding the parity anomaly and preserving the full rotational symmetry required to extract the complete operator content of the CFT. The identification of new primaries and the confirmation of weakly broken higher-spin symmetries enrich the CFT data set for the GNY universality class, offering an independent numerical benchmark for theoretical analyses that previously relied on $\epsilon$ -expansions or large- $N$ techniques.

More like this