Critical Ripples and Dirac Fermions in Crystalline Membranes

Technical Summary: Critical Ripples and Dirac Fermions in Crystalline Membranes

Problem Statement
This work investigates the low-energy field theory of crystalline membranes hosting massless Dirac fermions, with graphene serving as the paradigmatic example. The system presents a fundamental theoretical challenge: the coupling of two sectors with sharply different dynamical scalings. The flexural (out-of-plane) phonon modes of the membrane are nonrelativistic, governed by a quadratic dispersion ($\omega \sim k^2$) and a dynamical exponent $z_h = 2$. In contrast, the Dirac quasiparticles are relativistic-like, with a linear dispersion ($\omega \sim k$) and $z_\psi = 1$. The central question is how this dynamical mismatch controls the coupling between the membrane's elastic degrees of freedom and the electronic density, specifically determining whether electronic feedback destabilizes the flat phase or leads to new critical phenomena at finite momentum.

Methodology
The authors develop a controlled low-energy effective field theory, analyzing the system across three distinct scaling regimes:

1. Classical Static Membrane: Used as a benchmark for tensor algebra, with upper critical dimension $D=4$.
2. Zero-Temperature Quantum Membrane: Analyzed with $z=2$ scaling, lowering the upper critical dimension to $D=2$.
3. Finite-Momentum Order Parameter: Treated as a relativistic theory with $z=1$, relevant for rippling instabilities at wavevectors $\pm Q$, with upper critical dimension $D=3$.

The analysis relies on Renormalization Group (RG) flow equations derived via one-loop perturbation theory. Key theoretical tools include:

• Ward Identities: The rotational Ward identity is utilized to fix the scaling dimension of the isotropic tension operator (conjugate to the scalar strain $H_T$) without independent vertex renormalization.
• Weighted Power Counting: To handle the $z_h=2$ vs. $z_\psi=1$ mismatch, the authors expand the fermionic determinant in the parameter $s \sim \omega_h / (v_F k)$, which vanishes in the infrared.
• Composite Operator Analysis: The scaling of mixed vertices (e.g., strain-density coupling) is determined by projecting onto tree-level tensor structures and analyzing the irrelevance of fermionic feedback loops.
• Universality Classes: The finite-momentum transitions are classified based on symmetry matching between the structural ripple field and Dirac bilinears, mapping them to either the Wilson-Fisher fixed point or the Gross-Neveu-Yukawa (GNY) universality class.

Key Results

1. Stability of the Flat Phase (Charge Neutrality)
In the long-wavelength flat phase, the authors demonstrate an asymptotic decoupling between flexural modes and Dirac fermions.

• Renormalization of Coupling: The dimensionless strength of the leading local scalar strain–density coupling ($\bar{g}$) decreases logarithmically at long wavelengths. This flow is dictated entirely by the pure-membrane renormalization of the bending rigidity ($\kappa_b$) via the Ward identity.
• Irrelevance of Feedback: Fermionic corrections to the mixed vertex are suppressed by the running coupling $\alpha_e \sim \bar{g}^2 s$. Since the flexural modes become parametrically slower than the Dirac fermions ($s \to 0$) as the scale increases, the fermionic feedback vanishes as a power law.
• Dynamics: The Dirac polarization does not generate overdamped flexural dynamics, nor does it modify the leading infrared scaling of the Dirac kinetic term. Consequently, the flat phase remains stable against this local perturbation.

2. Finite-Momentum Rippling Instabilities
When microscopic physics (elasticity, lattice anisotropy, or electronic softening) selects isolated minima at finite wavevectors $\pm Q$, the system undergoes a structural transition. The nature of this transition depends on the symmetry of the ripple field relative to the Dirac bilinears:

• Mechanical (Wilson-Fisher) Route: If no low-energy Dirac bilinear shares the momentum and quantum numbers (including horizontal-reflection parity) of the ripple field, the leading fermionic interaction is a coupling between the conserved fermion density and the bosonic tuning operator ($|\Phi_Q|^2 \rho_\psi$).

  • In $D=2$, the scaling dimension of this coupling is $y_\lambda = 1/\nu - D < 0$.
  • The perturbation is irrelevant. The transition is governed by the purely bosonic Wilson-Fisher fixed point, with Dirac fermions acting as "spectators" at criticality.
  • In the ordered phase, the coupling shifts the Dirac-point energy relative to the chemical potential but does not alter the critical exponents.

• Symmetry-Locked (Gross-Neveu-Yukawa) Route: If the structural field transforms as a mass-type Dirac bilinear under the full microscopic symmetry group (requiring broken horizontal reflection or a symmetry-compatible bond component), a linear Yukawa coupling is allowed.

  • The critical mode becomes a hybrid electronic–structural order parameter.
  • The transition falls into the chiral-XY Gross-Neveu-Yukawa (GNY) universality class.
  • Using known one-loop critical exponents, the authors characterize the transition (e.g., $\nu = 1/2 + 3\epsilon/10$ for the chiral-XY case with $N_f=N_b=2$).
  • Velocity Locking: A significant result is that the fermionic velocity ($v_F$) and bosonic velocity ($c$) flow toward a common terminal value within the isotropic continuum limit, restoring emergent Lorentz symmetry at the critical point.

Significance and Claims
The paper claims to establish the universal infrared behavior of crystalline membranes coupled to Dirac fermions, distinguishing between asymptotic stability and finite-momentum instabilities.

• It resolves the question of electronic softening in the flat phase, proving that the "dynamical mismatch" renders electronic feedback irrelevant, thereby preserving the stability of the flat phase against local scalar perturbations.
• It provides a rigorous classification of finite-Q rippling transitions, identifying two distinct universality classes: a mechanically driven transition with spectator fermions and a hybrid GNY transition where electronic and structural degrees of freedom are inextricably linked.
• The work explicitly constructs the graphene realization of the chiral-XY GNY transition via Kekulé mass matrices, clarifying the specific symmetry conditions (broken horizontal reflection) required for this hybrid criticality to occur.

The authors note that their analysis is restricted to the local scalar projection and the isotropic continuum limit; extensions to the full projected deformation potential, shear channels, and finite-density regimes are left for future work.