The paper establishes a microscopic realization of an interacting 3D CFT with emergent N=1 supersymmetry using a Gross–Neveu–Yukawa model coupled with a gauged Majorana fermion.
It employs fuzzy sphere regularization to discretize spatial S² while preserving SO(3) symmetry, enabling precise extraction of operator scaling dimensions and a clear operator-state correspondence.
The findings reveal distinct supermultiplet structures and continuous RG flow of the operator spectrum, overcoming traditional challenges in lattice and perturbative methods.
The three-dimensional (3D) gauged Majorana conformal field theory (CFT) constitutes a nontrivial example of an interacting, strongly coupled CFT with emergent N=1 supersymmetry. Its microscopic realization involves coupling a 3D Ising CFT to a gauged Majorana fermion via a Yukawa interaction, capturing the Gross–Neveu–Yukawa universality class at criticality. The model exhibits the full structure of superconformal multiplets, including the signatures of spacetime supersymmetry such as precise relations between the scaling dimensions of bosonic and fermionic operators. The fuzzy sphere regularization provides a non-perturbative and numerically exact route to study these phenomena, enabling the extraction of operator scaling dimensions and the renormalization-group (RG) flow at the operator level (Tang et al., 31 Dec 2025).
1. Microscopic Formulation and Hamiltonian Structure
The continuum action of the 3D gauged Majorana CFT is given by the Gross–Neveu–Yukawa (GNY) Lagrangian: LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ4
where ψ is a two-component Majorana fermion, ϕ is a real Z2-odd scalar field, g is the Yukawa coupling, and u is the self-interaction of ϕ.
On the fuzzy sphere, the full model Hamiltonian combines separate bosonic (Ising) and fermionic (gauged Majorana) sectors with a short-range Yukawa coupling: H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]
Here, HIsing and LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ40 are fuzzy-sphere Hamiltonians for the concentric spheres with monopole charge LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ41; LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ42 and LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ43 denote bosonic and fermionic density operators; couplings LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ44 and LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ45 are expanded in Haldane pseudopotentials. The parameter LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ46 tunes the relative "speed of light" for the noninteracting CFTs, ensuring Lorentz invariance at criticality. Setting LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ47 yields the decoupled Ising LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ48 Majorana theory, while tuning to LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ49 brings the system to the GNY (ψ0 supersymmetric) critical point (Tang et al., 31 Dec 2025).
2. Fuzzy-Sphere Regularization and Operator-State Correspondence
The fuzzy sphere regularization discretizes the spatial ψ1 while preserving exact ψ2 symmetry. In this approach:
The presence of a monopole with charge ψ3 results in a lowest Landau level (LLL) of dimension ψ4, corresponding to the spin-ψ5 representation of ψ6. The coordinates ψ7 are replaced by rescaled angular-momentum operators satisfying ψ8, maintaining the sphere's geometry through ψ9.
Scalar and fermionic fields become ϕ0 matrices acting on the LLL.
Through radial quantization, energy eigenstates of the Hamiltonian map directly onto local scaling operators of the CFT. The scaling dimension is extracted as
ϕ1
with ϕ2 the speed of light and ϕ3 the ground-state energy. The normalization is set such that the stress tensor has ϕ4, or equivalently through the integer spacing of conformal towers predicted by the algebra (Tang et al., 31 Dec 2025).
3. Scaling Dimensions and Operator Spectrum at Fixed Points
The operator spectrum is resolved both at the decoupled and the interacting fixed points:
Decoupled theory (ϕ5):
Ising primaries: ϕ6, ϕ7
Free Majorana: ϕ8, ϕ9
Interacting SCFT (Z20), using fuzzy-sphere exact diagonalization:
Leading supermultiplet: Z21, Z22, Z23
Subleading supermultiplet: Z24, Z25, Z26
Supercurrent: Z27 at Z28
Stress tensor: Z29 (by normalization)
Additional spin-2 operator: g0
The hallmark of supersymmetry is manifest in the spacing g1, g2, and in the integer spacing of the corresponding conformal towers (Tang et al., 31 Dec 2025).
4. Multiplet Structure and Supersymmetry Relations
Three main supermultiplets are established, connected by the emergent supercharge g3 with scaling dimension g4:
g5, with g6 and g7.
Higher multiplet g8 at larger scaling dimensions.
Superconformal current multiplet with stress tensor g9 (u0) and supercurrent u1 (u2).
Integer-u3 sectors reproduce conformal towers predicted by the 3D conformal algebra, with each primary operator succeeded by descendants at u4 with the expected degeneracies (Tang et al., 31 Dec 2025).
5. Operator-Level Renormalization Group Flow
Tracking the operator spectrum as the microscopic pseudopotentials u5 are varied from u6 to u7 reveals a continuous RG trajectory:
u8: u9 flows from ϕ0 to ϕ1
ϕ2: ϕ3 from ϕ4 to ϕ5
ϕ6: ϕ7 from ϕ8 to ϕ9
The Yukawa operator H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]0 is initially relevant but becomes irrelevant at the interacting fixed point, reflecting SCFT stability.
Composite operators H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]1 and H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]2 hybridize along the flow, obeying the critical equation of motion H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]3.
Individual eigenstates evolve smoothly between the decoupled and interacting spectra, making the RG flow at the operator level directly visible in the spectrum (Tang et al., 31 Dec 2025).
6. Evidence for Emergent H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]4 Supersymmetry and Physical Consequences
Several features confirm the emergence of supersymmetry at the critical point:
Protected gaps: The low-lying spectrum at H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]5 features exact integer spacing characteristic of conformal multiplets and uniform SUSY spacing of H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]6 and H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]7.
Supercharge algebra: H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]8, with numerical agreement between scaling dimensions at the H=HIsing⊗I+I⊗rHMajorana+∫dΩfdΩbnfz(Ωf)[λz0(Ωfb)nb0(Ωb)+λzx(Ωfb)nbx(Ωb)]9–HIsing0 level.
Supercurrent and stress tensor: Observation of HIsing1 at the unitarity bound (HIsing2) and HIsing3 at HIsing4.
Central charge HIsing5: Not extracted directly, but fixable from the HIsing6 gap, with HIsing7, using the constraint HIsing8.
Non-perturbative control: The fuzzy sphere approach, in contrast to lattice or perturbative HIsing9-expansion, exactly preserves LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ400 symmetry and permits direct measurement of scaling dimensions (Tang et al., 31 Dec 2025).
7. Significance and Methodological Advancements
The fuzzy-sphere regularized Gross–Neveu–Yukawa model provides the first controlled, microscopic realization of a 3D LGNY=21(∂μϕ)2+iψˉ∂ψ+gϕψˉψ+2mϕ2+4uϕ401 superconformal Ising CFT. It allows for non-perturbative extraction of operator dimensions, transparent demonstration of operator-level RG flow, and identification of the full supermultiplet and conformal multiplet structure in the spectrum. This approach overcomes the challenges faced by lattice Monte Carlo simulations (notably the fermion sign problem) and perturbative expansions, establishing a pathway for future investigations of strong-coupling phenomena and emergent spacetime supersymmetry in higher-dimensional CFTs (Tang et al., 31 Dec 2025).