Correction exponents in the chiral Heisenberg model at $1/N^2$: singular contributions and operator mixing

This paper calculates the $1/N^2$ correction exponents for the chiral Heisenberg model, identifies a divergence in three dimensions caused by operator mixing among four-fermion operators, and proposes a resummation procedure that yields results consistent with direct three-dimensional calculations and the $\epsilon$ expansion.

The Gist

Imagine you are trying to predict how a crowd of people behaves when they are about to panic. In physics, we call this a "critical point." Whether it's a magnet losing its magnetism, a fluid turning into a gas, or electrons in a sheet of graphene (a super-thin material), the rules of the game change drastically right at that tipping point.

Physicists use a set of "rules" called critical exponents to describe exactly how things change near this tipping point. Think of these exponents as the "speed limits" or "acceleration rates" of the system as it approaches chaos.

This paper by Alexander Manashov and Leonid Shumilov is like a high-stakes engineering report on a specific type of crowd: the Chiral Heisenberg Model. This model describes a complex dance between particles called fermions (like electrons) and fields (like magnetic waves).

Here is the story of their discovery, broken down into simple concepts:

1. The "Big N" Strategy

Calculating these rules for a real-world system is incredibly hard because there are too many particles interacting. To make it manageable, the authors use a trick called the $1/N$ expansion.

• The Analogy: Imagine trying to predict the traffic flow of a single car (hard to generalize) versus a city with a million cars (easier to find patterns).
• The Trick: They pretend there are a huge number of particle types (let's call this number $N$). When $N$ is huge, the math becomes much simpler. They then calculate the answer for "infinite cars" and add small corrections for when $N$ is smaller (like our real world, where $N$ might be small).
• The Goal: They wanted to calculate the "correction exponents"—the second-order rules that tell us how the system adjusts as it gets closer to the tipping point.

2. The Unexpected Glitch (The Pole)

The authors did their math for a general number of dimensions (let's call it $d$). Usually, these calculations work smoothly whether you are in 2D, 3D, or 4D.

However, when they looked at the result for 3 dimensions (our physical world), they found a mathematical disaster: a pole.

• The Analogy: Imagine you are driving a car and calculating your speed. For every road condition, your speed is a normal number. But suddenly, when you hit a specific patch of road (3D), your speed calculation says "Infinity!" or "Divide by Zero!"
• What it means: In physics, a "pole" (a division by zero) usually means the math has broken down. It suggests that the way the particles interact in 3D is fundamentally different from how they interact in 4D or 2D. The standard "small correction" method they were using stopped working.

3. The Root Cause: The "Mixing" Problem

Why did the math break? The authors realized it was due to Operator Mixing.

• The Analogy: Imagine you have two different types of musical instruments: a Flute (a simple particle) and a Drum (a complex particle).
  • In 4D, the Flute and the Drum play different songs. They don't interfere with each other.
  • In 3D, the laws of physics change slightly. Suddenly, the Flute and the Drum start playing the exact same note at the same time. They "mix."
  • Because they are mixing, you can no longer treat them as separate instruments. If you try to calculate the sound of just the Flute, the Drum's sound is secretly messing up your math, causing that "Infinity" error.

In this paper, they found that in 3D, a specific particle interaction (involving four fermions) suddenly becomes identical in "weight" to another interaction. This causes a mathematical collision that breaks the standard calculation.

4. The Fix: Resummation

You can't just ignore the infinity. The authors had to invent a new way to add up the numbers, a process called Resummation.

• The Analogy: Imagine you are trying to add up a list of numbers: $1 + 0.5 + 0.25 + 0.125...$
  • If you stop early, you get a rough guess.
  • But if the numbers start behaving weirdly (like jumping to infinity) at a certain point, you need a special formula to "re-sum" the whole infinite series to get the true answer.
• The Result: By using this new "resummation" technique, they fixed the broken math. They found that the "speed limit" (the exponent) for the system in 3D isn't what the simple calculation suggested. It's actually a completely different value, determined by how the Flute and Drum mix together.

5. Why This Matters

This isn't just abstract math. The Chiral Heisenberg model is believed to describe Graphene and other exotic materials.

• The Takeaway: If you want to build better electronics using these materials, you need to know exactly how they behave at their critical points.
• The Discovery: The authors showed that if you use the old, simple math, you get the wrong answer for 3D materials. You must account for this "mixing" effect. They provided the corrected "speed limits" (exponents) that match perfectly with other methods, proving that their new way of fixing the math works.

Summary

The paper is a detective story where the authors:

1. Tried to calculate the rules of a complex particle system.
2. Found a "divide by zero" error when looking at our 3D world.
3. Realized the error was caused by two different particle interactions suddenly becoming identical (mixing).
4. Invented a new mathematical "glue" (resummation) to fix the error.
5. Confirmed that the corrected rules match reality, ensuring our understanding of materials like graphene is accurate.

It's a reminder that in the quantum world, sometimes things that look different in one dimension suddenly become twins in another, and you have to change your whole way of thinking to solve the puzzle.

Technical Summary

1. Problem Statement

The paper addresses the calculation of correction exponents (also known as critical exponents related to the slopes of $\beta$-functions) in the Chiral Heisenberg (CH) model.

• Context: The CH model describes scalar fields ($\pi_a$) interacting with $N$-component fermion fields ($q_i$). It is a candidate for describing critical phenomena in systems like graphene and belongs to the Gross-Neveu-Yukawa (GNY) universality class.
• Gap: While basic critical indices (like the anomalous dimension $\eta$) were known up to order $1/N^2$ (and $\eta$ to $1/N^3$), the correction exponents $\omega_{\pm}$ were previously only known at order $1/N$.
• Specific Challenge: The authors aim to calculate these exponents at order $1/N^2$. A critical issue arises when taking the limit $d \to 3$ (three dimensions). Previous calculations suggested a divergence (a pole) in one of the exponents at $d=3$, which contradicts the expectation that physical observables should be well-behaved. The paper seeks to explain the origin of this pole and propose a method to resolve it.

2. Methodology

The authors employ the Large $N$ expansion technique combined with conformal bootstrap methods and operator mixing analysis.

• Model Formulation: They utilize the critical equivalence between the CH model in $d=4-2\epsilon$ dimensions and a simplified effective action in $d=2\mu$ dimensions involving a non-local kernel.
• Operator Identification: The correction exponents are identified with the scaling dimensions of specific composite operators:
  • $O_- = \pi \partial^2 \pi$
  • $O_+ = \pi^4$
  • These operators mix under renormalization.
• Calculation Framework:
  • The anomalous dimensions are calculated as a series in $1/N$.
  • The authors use a specific renormalization technique where the propagator is modified by an auxiliary parameter $u$ to extract anomalous dimensions from renormalization factors $Z(u)$.
  • They compute Feynman diagrams up to order $1/N^2$, including single-box, double-box, self-energy, and vertex correction diagrams.
• Analysis of Singularities: The authors investigate the behavior of the scaling dimensions as $d \to 3$. They identify that the pole arises due to operator mixing between operators of different canonical dimensions that coincide only at specific dimensions (specifically, the mixing between the scalar operator $\pi \partial^2 \pi$ and a system of four-fermion operators).

3. Key Contributions

A. Calculation of $1/N^2$ Correction Exponents

The authors successfully derived the correction exponents $\omega_{\pm}$ at order $1/N^2$ for arbitrary dimension $d$.

• They calculated the anomalous dimensions for a general set of scalar operators $\pi^m$ and traceless symmetric operators.
• Consistency Check: When expanded in $\epsilon$ (near $d=4$), their results match perfectly with existing four-loop perturbative $\epsilon$-expansion results, validating the Large $N$ approach.

B. Identification of the $d=3$ Pole

The authors found that the correction exponent $\omega_-$ (associated with the operator $\pi \partial^2 \pi$) contains a pole at $d=3$ at order $1/N^2$:
$$\gamma_-(d) \sim \frac{1}{d-3}$$
They argue this is not a mathematical error but a physical phenomenon caused by operator mixing.

• Mechanism: In $d \neq 3$, the operator $\pi \partial^2 \pi$ (canonical dimension 4) does not mix with four-fermion operators (canonical dimension $2d-2$). However, at $d=3$, their canonical dimensions coincide ($4 = 2(3)-2$).
• Consequence: Diagrams that are finite for $d \neq 3$ become divergent at $d=3$ because the subgraphs start to diverge when the dimensions align. This leads to a breakdown of the standard perturbative expansion at the intersection point.

C. Resummation Procedure

To resolve the divergence and obtain physical results at $d=3$, the authors propose a resummation procedure:

1. Instead of treating the mixing perturbatively, they construct the full characteristic equation (determinant of the anomalous dimension matrix) for the mixing system involving $\pi \partial^2 \pi$ and the relevant four-fermion operators.
2. They demonstrate that while individual matrix elements ($\gamma_{ij}$) have poles at $d=3$, the coefficients of the characteristic polynomial ($A, B, C$) remain finite (regular) at $d=3$ up to high orders in $1/N$.
3. By solving the characteristic equation at $d=3$ using these regular coefficients, they derive finite, resummed scaling dimensions.

4. Results

• Analytical Expressions: The paper provides explicit, complex analytical formulas for $\omega_{\pm}$ at order $1/N^2$ in terms of $d$ (or $\mu$) and the digamma function $\psi$.
• Agreement with $\epsilon$-expansion: The expansion of these results near $d=4$ confirms agreement with known 4-loop results.
• Behavior at $d=3$:
  • The "naive" $1/N$ expansion (ignoring mixing) predicts a divergence.
  • The resummed exponents are finite and continuous at $d=3$.
  • The resummed value for the scaling dimension of the operator $\pi \partial^2 \pi$ at $d=3$ shifts significantly from the naive leading-order prediction. Specifically, the correction changes from $0$ to $-8\frac{\eta_1}{N}\sqrt{\frac{2}{3}}$.
• Numerical Impact: For $N=2$ (relevant for graphene), the deviation between the naive and resummed exponents is substantial, indicating that the mixing effect is crucial for accurate physical predictions in 3D.

5. Significance

• Theoretical Insight: The paper clarifies a subtle but critical issue in Large $N$ expansions: the breakdown of perturbation theory when operators of different canonical dimensions cross at specific spacetime dimensions. It generalizes the understanding of operator mixing beyond simple diagonalization.
• Methodological Advance: The proposed resummation technique provides a robust framework for calculating critical exponents in $d=3$ without needing to perform separate, complex 3D calculations that lose connection with the $4-2\epsilon$ expansion.
• Physical Relevance: Since the Chiral Heisenberg model is a primary candidate for describing the critical behavior of graphene and other Dirac materials, accurate determination of correction exponents is vital for predicting experimental observables. The finding that $1/N^2$ corrections significantly alter the $d=3$ results suggests that previous estimates might have been inaccurate due to unaccounted operator mixing.
• Universality: The authors argue that this pole-and-resummation phenomenon is a general feature of field theories with large $N$ expansions, applicable to other models like the Gross-Neveu model.

In summary, this work bridges the gap between $4-2\epsilon$ and $1/N$ expansions by resolving a singularity at $d=3$ through a sophisticated analysis of operator mixing, providing corrected and physically consistent critical exponents for the Chiral Heisenberg model.