Critical Probability Distributions of the order parameter at two loops II: $O(n)$ universality class

This paper presents a two-loop perturbative method to compute the family of probability distributions for the $O(n)$ model's order parameter, indexed by the ratio of system size to correlation length, and validates these results against Monte Carlo simulations and Functional Renormalization Group data.

The Gist

Imagine you are trying to predict how a massive crowd of people will behave at a music festival. Will they all stay spread out, or will they suddenly clump together in one spot? In physics, we call this "clumping" an order parameter, and the way that clumping happens is described by a probability distribution.

This paper is a highly technical "instruction manual" for calculating exactly how that clumping happens in a specific family of physical systems called the O(n) models.

Here is the breakdown of the paper using everyday analogies.

---

1. The Core Concept: The "Clumping" Rule

In physics, many things—from magnets to fluids—undergo a "phase transition" (like water turning to ice). At the exact moment this happens (the critical point), the system is in a state of chaos and order at the same time.

The researchers want to know the PDF (Probability Distribution Function).

• The Analogy: Think of the PDF as a "Weather Forecast for Chaos." Instead of saying "it will rain," the PDF says, "There is a 10% chance the crowd will be perfectly even, a 70% chance they will form small groups, and a 20% chance they will form one giant mosh pit."

2. The "O(n)" Family: Different Flavors of Chaos

The paper focuses on the O(n) universality class. In physics, "universality" means that very different things (like a liquid and a magnet) actually follow the same rules. The "$n$" is just a setting that changes the "flavor" of the system:

• $n=1$ (Ising Model): Like a crowd where people can only face North or South.
• $n=2$ (XY Model): Like a crowd where people can spin around in circles.
• $n=3$ (Heisenberg Model): Like a crowd where people can move in any direction in 3D space.

The author is taking a formula that worked for the simple "North/South" crowd and upgrading it to work for the "Spinning" and "3D" crowds.

3. The "Two-Loop" Upgrade: Sharpening the Lens

The paper mentions "two-loop order of perturbation theory."

• The Analogy: Imagine you are looking at a distant mountain through a pair of binoculars.
  • Zero-loop is like looking with your naked eye; you see the mountain, but it's a blurry blob.
  • One-loop is like using basic binoculars; you see the shape and the trees.
  • Two-loop is like using a high-powered telescope; you can suddenly see the individual rocks and cracks in the cliffside.

By going to "two loops," the scientist is making the mathematical prediction much more precise and detailed.

4. The "$\zeta$" Factor: The Size of the Room

The paper introduces a variable called $\zeta$ (zeta). This represents the ratio between the size of the system (the room) and the "correlation length" (how much one person's movement affects their neighbor).

• The Analogy: If you are in a tiny elevator, everyone bumps into each other immediately. If you are in a massive stadium, you can move without affecting someone on the other side. $\zeta$ tells the math whether we are talking about the "Elevator Scenario" or the "Stadium Scenario."

5. The Results: Checking the Math against Reality

After doing pages of incredibly dense math (the "Appendix" sections), the author compares their "theoretical forecast" against Monte-Carlo simulations (which are essentially massive, high-speed computer experiments).

The Verdict: The "Two-Loop Telescope" was much more accurate than the "One-Loop Binoculars." The math predicted the computer simulations very closely, especially for small movements. However, for "extreme" movements (the "large field behavior"), the math still struggles a bit—it's like the telescope is great for the mountain, but starts to blur when you try to look at a tiny bird flying far away.

Summary for a Non-Scientist

"I have created a much more precise mathematical formula to predict how particles or atoms will cluster together during a major physical change. I tested my formula against computer simulations, and it proved to be significantly more accurate than previous versions, especially when dealing with complex, multi-directional systems."

Technical Summary

Technical Summary: Critical Probability Distributions of the Order Parameter at Two Loops II: $O(n)$ Universality Class

1. Problem Statement

The paper addresses the computation of the Probability Distribution Functions (PDFs) of the order parameter (specifically the normalized total spin) for systems in the $O(n)$ universality class at criticality. While previous research has successfully computed these PDFs at the one-loop level for both Ising and $O(n)$ models, a systematic two-loop perturbative approach for the $O(n)$ model was lacking.

A central complexity in these computations is that the PDF is not a single function but a family of functions indexed by a parameter $\zeta = L/\xi_\infty$, where $L$ is the system size and $\xi_\infty$ is the bulk correlation length. This parameter arises because the double limit of infinite volume ($L \to \infty$) and criticality ($t \to 0$) is non-unique and depends on how the ratio of these two scales is approached.

2. Methodology

The author employs a field-theoretic perturbative approach using the $\epsilon = 4-d$ expansion. The key methodological steps include:

• Modified Gibbs Free Energy: The PDF is interpreted via a modified partition function $Z_{M,s}[h]$ using a sharply peaked Gaussian delta function (with mass $M \to \infty$). The problem is then mapped to calculating the modified Gibbs free energy $\Gamma_M[\phi]$ and the associated rate function $I(s)$.
• Two-Loop Expansion: The author generalizes the method used for the Ising model to the $O(n)$ model. This involves expanding the logarithm of the partition function around the mean-field configuration and computing 1-Particle Irreducible (1PI) diagrams up to two-loop order.
• Renormalization: The theory is parameterized using renormalized temperature $t$, coupling constant $g$, and fields $\phi_i$ defined at an infrared scale $\mu$. The author utilizes the Minimal Subtraction (MS) scheme and dimensional regularization to handle divergences.
• Scale Selection: To capture all relevant fluctuations, the renormalization scale is chosen as $\mu = L^{-1}$. This allows the rate function to be expressed in terms of dimensionless variables: $\tilde{s}$ (scaled spin), $\tilde{g}$ (scaled coupling), and $\tilde{t}$ (scaled temperature).
• Decomposition: The rate function is decomposed into an infinite-volume fixed-point potential $I_{\zeta,n,\text{inf}}(\tilde{s})$ and finite-size corrections $I_{\zeta,n,\text{fin}}(\tilde{s})$.

3. Key Contributions

• Generalization to $O(n)$: The primary contribution is the extension of the two-loop perturbative framework from the $n=1$ (Ising) case to the general $O(n)$ model.
• Analytical Expressions: The paper provides explicit, highly complex analytical expressions for the rate function up to second order in $\epsilon$ (Eq. 27 and 28). These include terms involving Jacobi-theta functions and various integral functions ($\Delta, \theta, I, K$) defined in the appendices.
• Large-Field Behavior Prediction: The author analytically derives the sub-leading universal large-field behavior of the rate function, showing it follows the form $I_{\zeta,n}(x) \sim x^{\delta+1} - n(\delta-1)^2 \log x$.
• Numerical Framework: The work provides the necessary mathematical machinery (discrete sum computations via Poisson summation) to evaluate the non-trivial integrals required for these distributions.

4. Results

• Improved Accuracy: Comparison with Monte-Carlo (MC) simulations for the $O(2)$ (XY model) and $O(3)$ (Heisenberg model) shows that the two-loop results provide a significantly better fit to the data than the one-loop results, particularly in the small-field region.
• Minimum of the Rate Function: The paper provides improved determinations for the minimum of the rate function ($I_{\min}$) for $O(2)$ and $O(3)$ models, showing closer agreement with MC data than previous perturbative orders.
• Comparison with FRG: The results are compared against Functional Renormalization Group (FRG) calculations, showing qualitative and quantitative consistency.
• Large-Field Discrepancy: The author notes that while the small-field behavior is well-reproduced, the large-field behavior deviates from MC data. This is attributed to the fact that the two-loop expansion provides an $\epsilon$-expansion of the exponent $\delta$ rather than the true non-perturbative exponent.

5. Significance

This work advances the understanding of finite-size scaling and critical phenomena in $O(n)$ symmetric systems. By providing a systematic way to compute the entire family of PDFs indexed by $\zeta$, it offers a powerful tool for studying how boundary conditions and system geometry influence critical fluctuations. Furthermore, the paper identifies the need for RG improvement to correct large-field behavior and sets the stage for future investigations into the "triviality" of field theories at $d=4$ and the application of these methods to out-of-equilibrium systems.