Constraint effective action and critical correlation functions at fixed magnetization

This paper extends the functional renormalization group framework to compute momentum-dependent critical observables at fixed magnetization for the Ising universality class, demonstrating that the second-order derivative expansion accurately reproduces universal rate functions and correlation functions in three dimensions and qualitatively agrees with simulations in two dimensions, where lower-order approximations fail.

The Gist

The Big Picture: Predicting the "Mood" of a Crowd

Imagine a giant crowd of people (like the atoms in a magnet) standing in a room. Each person can be in one of two moods: happy (spin up) or sad (spin down).

Usually, if you look at the whole crowd, the moods balance out, and you get a mix of happy and sad people. But sometimes, the crowd reaches a "critical moment" (like a phase transition). At this exact moment, everyone starts influencing everyone else. The whole room becomes a single, giant entity where a change in one corner ripples through the entire space.

Scientists want to know: What does the probability distribution of this crowd's mood look like?

• Is it a standard bell curve (most people are neutral, few are extreme)?
• Or is it something weird and non-Gaussian (lots of extreme moods)?

This paper is about a new, more powerful way to calculate that "mood distribution" and how the people in the crowd are connected to each other, specifically when we force the total mood of the room to be fixed at a certain level.

---

The Problem: The "Fixed Magnetization" Puzzle

In physics, this "total mood" is called magnetization.

• The Old Way: Scientists used a tool called the Functional Renormalization Group (FRG). Think of FRG as a high-powered microscope that zooms out to see how the system behaves at different scales.
• The Limitation: Previous versions of this microscope were a bit "blurry." They used a simple approximation (called LPA) that assumed the crowd was perfectly smooth and ignored how the "texture" of the connections between people changed. This worked okay for 3D systems (like a cube of atoms) but failed completely for 2D systems (like a flat sheet of atoms) because the 2D crowd is much more chaotic and "wiggly."

The Goal of this Paper:
The authors wanted to upgrade the microscope. They wanted to:

1. Fix the "blur" by adding more detail (calculating up to the "second order" of complexity).
2. Apply this to a specific, tricky scenario: What happens if we lock the total magnetization of the system to a specific value?
3. See if this new, sharper tool works for both 2D and 3D systems and matches real computer simulations.

---

The Solution: The "Constraint Effective Action"

To solve this, the authors developed a new mathematical tool called the Constraint Effective Action.

The Analogy: The "Silent Room" Experiment
Imagine you want to study how a crowd behaves, but you have a rule: The total number of happy people minus sad people must equal exactly 50.

• In a normal experiment, the crowd might naturally drift to 0, 10, or 100.
• Here, you are forcing them to stay at 50.

The authors created a mathematical "force field" (a soft constraint) that gently pushes the system to stay at that fixed number. As they crank up the force to infinity, it becomes a hard rule. This allows them to calculate the Rate Function (a fancy name for the probability curve of the system) and the Correlation Functions (how likely it is that two people far apart are feeling the same way).

Key Findings

1. Sharper Focus (The DE2 Upgrade)

The authors upgraded their tool from a "Local Potential Approximation" (LPA) to a "Second-Order Derivative Expansion" (DE2).

• LPA (The Old Lens): Like looking at a crowd from a distance and assuming everyone is a smooth, blurry blob. It missed the fine details.
• DE2 (The New Lens): Like putting on high-definition glasses. It accounts for how the "texture" of the crowd changes.
• Result: In 3D, the new lens gave a much more accurate picture, matching computer simulations (Monte Carlo) almost perfectly. In 2D, the old lens (LPA) broke down completely, but the new lens (DE2) worked, though it still had some small errors (about 10-20%).

2. The "Zero-Momentum" Quirk

One of the most interesting discoveries was about how the crowd behaves when you look at the "average" connection (zero momentum).

• The Rule: If you fix the total mood of the room, the "fluctuation" of the entire room's mood must be zero (because it's locked!).
• The Surprise: The math showed that the behavior of the crowd at this "locked" state is fundamentally different from the behavior at any other scale. It's like a drum that vibrates everywhere except for the very center point, which is glued down. The authors had to invent a new mathematical term (called $\check{\Delta}$) to describe this "glued" point, which disappears in large, infinite systems but is crucial for finite, real-world systems.

3. Checking Against Reality (Monte Carlo Simulations)

The authors didn't just do math on paper; they compared their results to massive computer simulations (Monte Carlo), which act as the "ground truth."

• In 3D: Their new method matched the computer simulations incredibly well. They could predict the shape of the probability curve and how the connections between atoms changed with distance.
• In 2D: The match was good, but not perfect. The authors noted that in 2D, the system is so sensitive that even their advanced tool struggles slightly with the extreme "tails" of the distribution (the rare, extreme moods). They also noticed some strange "wiggles" in the 2D data that they suspect are caused by "droplets" (small islands of opposite mood) forming inside the fixed magnetization.

The Conclusion

This paper is a success story for mathematical physics.

• They proved that the Functional Renormalization Group (FRG) is a robust tool, even when you add the complex constraint of a fixed magnetization.
• By upgrading the math to the second order (DE2), they fixed the failures of the old method, especially in 2D systems.
• They showed that when you lock a system's total state, the rules of how it fluctuates change in a unique way that requires special mathematical handling.

In short: They built a better telescope, pointed it at a very difficult type of star (a 2D magnet with a fixed mood), and confirmed that their new telescope sees the star much more clearly than the old one did, matching the photos taken by the best cameras (computer simulations) available.

Technical Summary

Technical Summary: Constraint Effective Action and Critical Correlation Functions at Fixed Magnetization

Problem Statement
Continuous phase transitions are characterized by universal scaling laws and non-Gaussian probability distribution functions (PDFs) of macroscopic observables, such as the order parameter. While the central limit theorem fails near criticality due to strong correlations, the PDF adopts a universal form described by a "rate function." Previous functional renormalization group (FRG) studies successfully computed these rate functions using the Local Potential Approximation (LPA). However, the LPA neglects field renormalization, imposing a vanishing anomalous dimension ($\eta = 0$), which renders it inadequate for systems with large anomalous dimensions, such as the two-dimensional Ising model ($\eta = 1/4$). Furthermore, existing FRG frameworks had not been systematically extended to compute momentum-dependent observables (correlation functions) subject to a global constraint of fixed magnetization. This paper addresses these gaps by extending the FRG framework to the second-order derivative expansion (DE2) to compute both the constraint effective action and momentum-dependent correlation functions at fixed magnetization for Ising universality classes in two and three dimensions.

Methodology
The authors employ the Functional Renormalization Group (FRG) formalism, specifically utilizing the constraint effective action ($\check{\Gamma}$). This quantity is defined as the Legendre transform of the generating function for the order parameter PDF under a fixed magnetization constraint $s$.

1. Constraint Formalism: To handle the fixed magnetization constraint $\hat{\phi}_0 = s$, the authors introduce a soft constraint via a Gaussian term with parameter $M$, taking the limit $M \to \infty$. This defines the scale-dependent constraint effective action $\Gamma_{M,k}[\phi]$.
2. Flow Equations: The exact Wetterich equation is adapted for this constrained system. The authors derive exact flow equations for the scale-dependent rate function $I_k(s)$ and the $n$-point vertex functions $\check{\Gamma}^{(n)}_k$.
3. Finite-Size Effects and Vertex Structure: A critical technical challenge identified is the behavior of vertices at zero momentum in finite-size systems. Unlike the thermodynamic limit, the constraint introduces a "finite-momentum shift" $\check{\Delta}_{n,k}(s)$ for vertices. Specifically, the two-point vertex $\check{\Gamma}^{(2)}_k(p; s)$ behaves differently at $p=0$ compared to $p \neq 0$. The authors derive that $\check{\Gamma}^{(2)}_k(p; s) = I''_k(s) + (1-\delta_{p,0})\check{\Delta}_{2,k}(s) + \check{Z}_k(s)p^2$.
4. Derivative Expansion (DE2): The hierarchy of flow equations is closed using a second-order derivative expansion. To handle the non-closure caused by the new shift functions $\check{\Delta}_{n,k}$, the authors introduce an approximation where higher-order shifts are related to derivatives of the two-point shift (e.g., $\check{\Delta}_{3,k} \approx \check{\Delta}'_{2,k}$). This allows the simultaneous numerical integration of flows for the rate function $I_k$, the field renormalization $\check{Z}_k$, and the shift $\check{\Delta}_{2,k}$.
5. Numerical Implementation: The flow equations are solved numerically for $d=2$ and $d=3$ Ising models at critical temperature $T_c$. The authors utilize an exponential regulator and apply the Principle of Minimal Sensitivity (PMS) to optimize the regulator parameter $\alpha$.

Key Contributions and Results

• Extension to DE2: The paper successfully extends the FRG calculation of critical PDFs beyond the LPA to the DE2. This is crucial for capturing the correct anomalous dimension $\eta$.
• Three-Dimensional Ising Model:
  • The DE2 calculation yields critical exponents $\eta \approx 0.0455$ and $\nu \approx 0.6275$, which are in good agreement with conformal bootstrap benchmarks ($\eta \approx 0.0363$, $\nu \approx 0.6300$).
  • The computed rate function $I(s)$ shows excellent agreement with Monte Carlo (MC) simulations up to $s \approx 1$.
  • The universal amplitude $\Delta I$ is calculated as $-0.73(6)$, consistent with MC results converging toward $-0.7878(1)$ as system size increases.
  • Momentum-dependent correlation functions (specifically the first few Fourier modes) are accurately reproduced, demonstrating the convergence of the method.
• Two-Dimensional Ising Model:
  • The LPA fails to describe the critical point in 2D (failing to find the phase transition), necessitating the DE2.
  • Using DE2, the authors obtain $\eta = 0.293$ and $\nu = 1.07$, compared to exact values of $\eta = 0.25$ and $\nu = 1.0$. The authors estimate a 10–20% error margin.
  • While the rate function minimum is somewhat overestimated compared to MC, the qualitative agreement is maintained.
  • A notable finding is the behavior of the second-momentum component of the correlation function, $A(p_2; s)$, which exhibits oscillations not seen in higher modes. The authors hypothesize these arise from droplet effects induced by the fixed magnetization, analogous to domain wall effects in 1D.
• Robustness of FRG: The study confirms that the FRG approach is robust for calculating both zero-momentum (rate function) and finite-momentum (correlation functions) observables at fixed magnetization.

Significance and Claims
The paper claims to provide a systematic, non-perturbative framework for computing universal scaling functions and correlation functions at fixed magnetization, extending previous work that was limited to the LPA. The primary significance lies in:

1. Validation of DE2: Demonstrating that the second-order derivative expansion significantly improves accuracy, particularly for 2D systems where $\eta$ is large, and provides a method to estimate error bars by comparing LPA and DE2 results.
2. Momentum Dependence: Successfully deriving and solving flow equations for momentum-dependent correlation functions under a global constraint, a task previously unaddressed in this context.
3. Finite-Size Physics: Highlighting the necessity of accounting for specific finite-size vertex structures (the $\check{\Delta}$ terms) that vanish in the thermodynamic limit but are essential for describing constrained systems.

The authors conclude that their results are in qualitative and, in many cases, quantitative agreement with Monte Carlo simulations, validating the FRG method as a precise tool for studying critical probability distributions and correlation functions in both 2D and 3D Ising universality classes. They suggest that future work could apply the Blaizot–Mendez-Galain–Wschebor (BMW) approximation to further explore the full momentum dependence of these constrained correlation functions.