Lee-Yang zeros and edge singularity in a mean-field approach

This paper investigates the analytic structure of the partition function in a finite-volume mean-field QCD model to analyze the temperature dependence of Lee-Yang zeros and edge singularities, demonstrating that while finite-size scaling methods can successfully locate the critical point, accurate determination requires a careful treatment of corrections from irrelevant operators.

The Gist

Imagine you are trying to find the exact spot on a map where a material changes from a solid to a liquid, or from a magnetized state to a non-magnetized one. In physics, this special spot is called a Critical Point (CP).

The problem is that in the real world (and in computer simulations), we can't look at an infinitely large piece of material. We are stuck looking at small, finite chunks. When you look at a small chunk, the "sharp" change at the Critical Point gets blurry and smeared out, making it very hard to pinpoint exactly where it is.

This paper is like a guidebook for finding that blurry spot using a clever mathematical trick involving "ghost numbers." Here is how the authors did it, explained simply:

1. The Problem: The "Blurry" Edge

In a perfect, infinite world, the transition at the Critical Point is sharp. But in a finite box (like a computer simulation), the transition is smooth. It's like trying to find the exact moment a sunset turns into night; on a small scale, the colors blend gradually, making it hard to say exactly when "day" ended and "night" began.

Physicists usually try to guess the location by looking at how the "sensitivity" of the material changes as they shrink or grow the box size. This is called Finite-Size Scaling.

2. The Solution: The "Ghost" Zeros

The authors used a concept called Lee-Yang Zeros. Imagine the mathematical formula that describes the material (the partition function) as a complex machine. If you plug in normal numbers, the machine works fine. But if you plug in "imaginary" or "ghost" numbers (complex numbers), the machine sometimes breaks down and outputs zero.

• The Analogy: Think of these zeros as "ghost holes" in a map. In a small box, these holes are scattered around. As you make the box bigger and bigger, these holes start to line up and form a wall.
• The Edge: The very tip of this wall of holes is called the Edge Singularity. In an infinite world, this tip touches the real map exactly at the Critical Point.

The authors' goal was to watch how these "ghost holes" move as they change the size of the box and the temperature, to see where they are heading.

3. The Method: A Better Map

The authors used a simplified model of nuclear matter (quarks and mesons) and applied a specific technique to handle the "finite size" problem.

• The Old Way: Traditional methods often assumed the material was perfectly uniform, which gave wrong answers for small boxes because it ignored tiny fluctuations.
• The New Way: The authors added a step where they "averaged out" the fluctuations of the uniform field. This kept the math simple (like a mean-field approach) but fixed the error, ensuring the math stayed smooth and accurate even for small boxes.

4. The Discovery: The "Magic" Plane

When they plotted the movement of these ghost holes, they found something interesting about the coordinate system:

• If you plot the holes on a standard map (using the chemical potential $\mu_B$), the path gets wiggly and hard to follow when the temperature gets high.
• The Trick: If you plot the holes on a map of the square of the chemical potential ($\mu_B^2$), the path becomes a straight, clean line.
• The Metaphor: It's like trying to draw a straight line on a curved piece of paper. If you flatten the paper (change the coordinate system), the line becomes perfectly straight, making it much easier to predict where it will go.

5. The Results: Finding the Spot

The team tested three different ways to find the Critical Point using these ghost holes:

1. The Ratio Method: Comparing the distance between different ghost holes.
2. The Scaled Method: Looking at the position of a single ghost hole after adjusting for size.
3. The Binder Method: A standard statistical tool used to find phase transitions.

What they found:

• All three methods worked well! They could locate the Critical Point with very high accuracy (within 1%) even when looking at relatively small boxes.
• The Catch: As they looked at larger and larger boxes, the accuracy didn't get perfectly smooth immediately. There was a tiny "bump" in the data.
• The Reason: This bump was caused by "irrelevant operators."
  • The Analogy: Imagine you are trying to hear a whisper (the main signal) in a quiet room. At first, the room is noisy (small box). As the room gets bigger, the noise fades. But then, you realize there is a very faint, high-pitched squeak (the irrelevant operator) that only becomes noticeable when the room is huge. This squeak messes up the perfect prediction if you don't account for it.

Conclusion

The paper demonstrates that by using a specific mathematical framework to track "ghost zeros" in the complex plane, physicists can accurately locate the Critical Point of nuclear matter, even when working with limited, finite-sized data. They showed that while these methods are powerful, you must be careful to account for subtle mathematical "squeaks" (corrections from irrelevant operators) to get the most precise result possible.

In short: They found a better way to draw the map of the "ghost holes" so that even with a small telescope, we can see exactly where the Critical Point is hiding.

Technical Summary

Problem Statement
Locating the second-order critical point (CP) in Quantum Chromodynamics (QCD) using thermodynamic observables in finite-volume systems is a nontrivial task. In the thermodynamic limit, the universal divergence of susceptibilities uniquely marks the critical region; however, in finite systems, these divergences are smeared, rendering all thermodynamic quantities analytic. While finite-size scaling (FSS) offers a systematic approach to locate the CP from finite-volume data, standard mean-field approaches in effective models often fail to capture the necessary analytic structure. Specifically, conventional mean-field treatments yield non-analytic phase transitions even for finite systems and cannot describe Lee-Yang zeros (LYZs), which arise from the cancellation of contributions from different field configurations. Furthermore, while functional methods can investigate finite-size rounding, they face significant difficulties in treating the behavior of LYZs in the complex parameter space at finite sizes.

Methodology
The authors investigate the analytic structure of the partition function using a minimal mean-field effective model of QCD (the quark-meson model) that incorporates finite-size effects through a specific modification of the mean-field approximation. Instead of minimizing the effective potential $U(\bar{\phi})$ to determine the order parameter, the partition function is defined by integrating over the spatially uniform mode $\bar{\phi}$:
$$Z = \int d\bar{\phi} e^{-V U(\bar{\phi})/T}$$
This approach, introduced in Ref. [59], ensures the partition function remains analytic in finite volume while recovering the conventional mean-field results in the thermodynamic limit. The study employs a Landau potential expansion near the CP to derive finite-size scaling relations. The authors analyze the distribution of LYZs and the Lee-Yang edge singularity (LYES) in the complex baryon chemical potential ($\mu_B$) plane. They compare different methods for locating the CP based on finite-size scaling, specifically:

1. Intersection analysis of Lee-Yang zero ratios (LYZR).
2. Intersection analysis of scaled LYZs ($L^{y_\eta} \text{Im} \mu_{LY}$).
3. Intersection analysis of the Binder cumulant ($B_4$).

The study utilizes two parameter sets (A and B) for the quark-meson model to examine parameter dependence and compares the results with lattice QCD findings.

Key Contributions and Results

• Analytic Structure and LYZ Trajectories: The study demonstrates that the integration over the uniform mode successfully generates analytic partition functions with discrete LYZs. The trajectories of the first LYZ, $\mu^{(1)}_{LY}(T, L)$, are shown to approach the LYES in the thermodynamic limit. The authors find that the system-size dependence of the LYZs is more naturally understood in the complex $\mu_B^2$ plane than in the complex $\mu_B$ plane. In the $\mu_B^2$ plane, the trajectories exhibit more linear behavior, and the real part of the squared LYZ, $\text{Re}(\mu_{LY}^{(1)})^2$, is insensitive to system size $L$ even at temperatures where the linear $\mu_B$ dependence shows strong $L$-dependence.
• Pseudo-critical Lines: The study defines a pseudo-critical line based on the LYZ, $\mu_{LY}^{pc}(T, L) \equiv \sqrt{\text{Re}(\mu_{LY}^{(1)}(T, L))^2}$. This definition shows good agreement with the susceptibility-based pseudo-critical line $\mu_{\chi^2}^{pc}(T)$ for system sizes $L \gtrsim 8$ fm. The authors note that the deviation of the LYES real part from the susceptibility peak at small $\mu_B$ is better described when mapping to the $\mu_B^2$ plane, consistent with charge conjugation symmetry.
• Comparison with Lattice QCD: The authors compare their results with lattice QCD simulations using Padé approximations. They observe that for fixed aspect ratios ($LT$), the first LYZ remains significantly distant from the LYES at small $LT$(e.g.,$LT=2, 4$). This suggests that identifying the CP by extrapolating the condition $\text{Im} \mu_{LY}^{(1)} = 0$ in lattice simulations may lead to an underestimation of the critical temperature $T_{CP}$ and an overestimation of the critical chemical potential $\mu_{CP}$.
• Intersection Analysis and Convergence: The intersection analysis methods (LYZR, scaled LYZ, and Binder cumulant) successfully identify the CP within 1% accuracy for moderate system sizes ($L \ge 8$ fm). However, the convergence of the intersection points toward the thermodynamic limit slows down beyond this accuracy.
• Role of Irrelevant Operators: A detailed analysis of the convergence behavior reveals that the slow convergence at large $L$ is attributed to corrections from irrelevant operators. In the specific mapping of the quark-meson model to the Landau potential, a $\bar{\phi}^5$ term emerges (due to the lack of $Z_2$ symmetry in the general system). This term scales as $L^{-3/4}$, which is a slower convergence rate than the $L^{-3/2}$ scaling expected from the linear mapping of the critical point. Consequently, the $L^{-3/4}$ term dominates the behavior of intersection points at large $L$, causing deviations from the universal values. The authors confirm this by fitting the data with a function including $L^{-3/4}$, $L^{-3/2}$, and $L^{-9/4}$ terms, which accurately describes the numerical results.

Significance
The paper establishes that the minimal extension of the mean-field approach, involving integration over the uniform mode, provides a consistent framework for studying both LYZs in finite volume and the LYES in the infinite-volume limit. The study highlights that while standard intersection methods are effective for locating the CP, a careful treatment of corrections from irrelevant operators is crucial for high-precision determination. Specifically, the presence of subleading odd-order terms (like $\bar{\phi}^5$) in systems without $Z_2$ symmetry introduces scaling violations that can significantly affect the convergence of intersection analyses. The authors argue that understanding these corrections is essential for interpreting lattice QCD results where finite-size effects are inevitable. Furthermore, the finding that the LYZ behavior is more linear in the $\mu_B^2$ plane offers a refined perspective on mapping thermodynamic variables near the CP. The work suggests that future systematic improvements should include non-homogeneous classical configurations and the $Z_3$ center symmetry to fully capture QCD-specific features like the Roberge-Weiss transition.