A parameter independent analysis of the QCD chiral phase transition and its universal critical behaviour

Technical Summary: Parameter-Independent Analysis of the QCD Chiral Phase Transition

Problem Statement
The nature of the chiral phase transition in Quantum Chromodynamics (QCD) at vanishing light quark masses is central to understanding the QCD phase diagram and the fate of the $U(1)_A$ axial anomaly. While lattice studies have confirmed the restoration of $SU(2)_L \times SU(2)_R$ chiral symmetry, the influence of the anomaly on the universal critical behavior remains an open question. A significant limitation in existing literature is that most numerical studies presuppose the transition belongs to the three-dimensional $O(4)$ universality class from the outset, constructing scaling ansatzes using $O(4)$ critical exponents. This approach precludes a direct, parameter-independent lattice QCD determination of the critical parameters governing the transition. The authors aim to address this by systematically determining universal critical parameters directly from lattice data without relying on a priori assumptions regarding the underlying universality class.

Methodology
The study utilizes (2+1)-flavor QCD simulations employing highly improved staggered fermions on lattices with a fixed temporal extent of $N_\tau = 8$. The analysis focuses on zero baryon density and employs the following methodological framework:

1. Improved Order Parameter: The authors utilize an improved chiral order parameter, $M$, defined as $M = M_\ell - H\chi_\ell$, where $M_\ell$ is the renormalized light quark chiral condensate, $\chi_\ell$ is the chiral susceptibility, and $H = m_\ell/m_s$ is the ratio of light-to-strange quark masses. This subtraction removes linear light quark mass dependence, ensuring the leading order dependence is cubic ($O(H^3)$) and providing a well-defined continuum limit.
2. Finite-Size Scaling (FSS): The study analyzes the finite-volume dependence of $M$ and its susceptibility across a wide range of lattice volumes ($N_\sigma$) and light quark masses ($H$). The analysis is grounded in the scaling hypothesis, where observables are expressed as sums of singular parts (described by universal scaling functions $f_G$ and $f_\chi$) and sub-leading regular terms.
3. Scaling Variables: The analysis employs scaling variables $z = t h^{-1/\Delta}$ and $z_L = l_0 h^{\nu/\Delta}/L$, where $t$ and $h$ are reduced temperature and magnetic field, and $L$ is the system size. The study investigates the behavior of the rescaled order parameter $M_{resc} = H^{-1/\delta} M$ to identify intersection points that define the critical temperature $T_c$.
4. Universality Class Testing: While the paper references $O(2)$ scaling functions (noted as a proxy for $O(4)$ in the context of the specific scaling functions used for parametrization in Ref. [12]), the primary goal is to extract critical exponents and amplitudes directly from the data to test universality rather than assuming it.

Key Results

• Intersection Point and $T_c$ Determination: Using the rescaled order parameter $M_{resc}$, the authors observe that curves for different quark mass ratios ($H$) intersect at a unique point in the chiral limit. For $N_\tau = 8$, the initial intersection analysis suggested $T_c \approx 143.7$ MeV. However, after performing a detailed finite-volume extrapolation including new data at $H = 1/240$ on $92^3 \times 8$ lattices, the estimate for the chiral phase transition temperature was refined to $T_c = 144.4(4)$ MeV.
• Finite-Volume Effects: The study quantifies finite-volume corrections to the scaling behavior. It is found that for $z_{L,b} > 1$ (where $z_{L,b} = H^{-\nu/\Delta}/L$), finite-volume effects become significant (reaching ~10% for $z_{L,b} \sim 1.5$). The approach to the infinite-volume limit is observed to be faster than $1/V$ (scaling as $z_{L,b}^3$), with leading-order dependence closer to $z_{L,b}^4$ (or $V^{-4/3}$).
• Volume Requirements: To mitigate finite-size effects and approach the chiral limit accurately, the authors determine that the lattice aspect ratio $N_\sigma/N_\tau$ must be sufficiently large. At physical quark masses, $N_\sigma/N_\tau \ge 6$ is required, increasing to approximately 11 for the smallest quark mass ratio studied ($H = 1/240$).
• Deviations from Universal Scaling: For quark mass ratios $m_\ell/m_s > 1/160$, deviations from the leading-order universal critical behavior are observed, attributed to corrections-to-scaling and regular terms. At the physical point ($H = 1/27$), these corrections remain below the 10% level.

Significance and Claims
The paper claims to provide a "parameter independent analysis" by avoiding the presumption of a specific universality class (such as $O(4)$) at the onset of the calculation. Instead, it seeks to derive critical parameters directly from the scaling behavior of lattice observables. The authors emphasize that their work is part of an ongoing research project. While they have successfully established a methodology to estimate $T_c$ and quantify finite-volume effects, they modestly state that a definitive quantitative determination of the critical exponent $\delta$—necessary to clearly distinguish between universality classes—requires further data collection, particularly at smaller quark masses ($H$) and larger volumes. The current results serve as a refined baseline ($T_c = 144.4(4)$ MeV) and a detailed characterization of finite-size scaling corrections, paving the way for a direct lattice QCD determination of the universal critical nature of the chiral phase transition.