Searching for the QCD critical point through constant entropy density contours

Basic Information

• Paper ID: 2510.13001
• Title: Searching for the QCD critical point through constant entropy density contours
• Authors: Hitansh Shah, Mauricio Hippert, Jorge Noronha, Claudia Ratti, Volodymyr Vovchenko
• Classification: hep-ph (High Energy Physics - Phenomenology), nucl-th (Nuclear Theory)
• Publication Date: October 14, 2024
• Paper Link: https://arxiv.org/abs/2510.13001

Abstract

This paper proposes a novel method for locating the QCD critical point by constructing constant entropy density contours. Through two independent analyses of lattice QCD data at zero baryon chemical potential, the critical point is found to be located at $T_c = 114 \pm 7$ MeV and $\mu_{Bc} = 602 \pm 62$ MeV in an expansion truncated at $\mu_B^2$ order. The method remains consistent with the latest lattice QCD results within the range $\mu_B/T \leq 3$. More precise determination of the required expansion coefficients from lattice simulations is crucial for reliably establishing the location of the QCD critical endpoint.

Research Background and Motivation

Core Problem

This research addresses the fundamental problem of locating the critical point in the QCD phase diagram of strongly interacting matter. The QCD critical point (CP) is the endpoint of the first-order phase transition line connecting the quark-gluon plasma phase and the hadronic gas phase. Its precise location is essential for understanding the phase structure of strongly interacting matter.

Problem Significance

1. Experimental Relevance: One of the primary objectives of collider experiments such as RHIC and FAIR is to search for the QCD critical point, which is crucial for understanding the properties of baryon-rich matter
2. Theoretical Challenge: Direct lattice QCD simulations at non-zero net baryon density are limited by the fermion sign problem
3. Phase Diagram Completeness: Determining the critical point location is a key step in providing a complete description of the QCD phase diagram

Limitations of Existing Methods

• Lattice QCD simulations at zero baryon chemical potential show an analytic crossover transition from hadronic gas to quark-gluon plasma
• Direct lattice calculations at finite real net baryon density face technical difficulties
• Existing theoretical predictions and model estimates exhibit substantial uncertainties

Research Motivation

This paper proposes a novel method utilizing constant entropy density contours to extrapolate from $\mu_B = 0$ to finite $\mu_B$, based on the physical property that entropy becomes multivalued in the first-order phase transition region to locate the critical point.

Core Contributions

1. Proposed a novel method based on constant entropy density contours for locating the QCD critical point
2. Established a systematic mathematical framework for extrapolation from zero chemical potential to finite chemical potential through $\mu_B^2$ order expansion
3. Obtained specific coordinates of the critical point: $(T_c, \mu_{Bc}) = (114 \pm 7, 602 \pm 62)$ MeV
4. Verified method consistency with lattice QCD results showing excellent agreement within the range $\mu_B/T \leq 3$
5. Provided systematic error analysis by evaluating result reliability through two independent data processing methods

Methodology Details

Theoretical Foundation

The method is based on the physical property that entropy density becomes multivalued during a first-order phase transition. Near the critical point, the crossing of entropy density contours can be used to determine the first-order phase transition line and the critical endpoint.

Mathematical Framework

Constant Entropy Density Contour Expansion

Define temperature $T_s$ such that the constant entropy density condition is satisfied: $s[T_s(\mu_B; T_0), \mu_B] = s(T_0, \mu_B = 0)$

Expand $T_s$ in powers of $\mu_B$: $T_s(\mu_B; T_0) \approx T_0 + \alpha_2(T_0)\frac{\mu_B^2}{2} + O(\mu_B^4)$

Expansion Coefficient Calculation

The second-order expansion coefficient $\alpha_2$ is calculated at $\mu_B = 0$: $\alpha_2(T_0)\Big|_{\mu_B=0} = -\frac{2T_0\chi_2^B(T_0) + T_0^2\chi_2^{B'}(T_0)}{s'(T_0)}$

where:

• $\chi_2^B = \left[\frac{\partial^2(p/T^4)}{\partial(\mu_B/T)^2}\right]_T$ is the second-order baryon number susceptibility
• $s$ is the entropy density
• Prime denotes temperature derivative

Critical Point Condition

The critical point is an inflection point satisfying $\partial^2s/\partial^2T$ divergence. Combined with the constant entropy density condition, the determining equation for the critical point location is obtained:

$\alpha_2''(T_{0,c}) = 0$$\mu_{B,c} = \sqrt{-\frac{2}{\alpha_2'(T_{0,c})}}$$T_c = T_{0,c} + \alpha_2(T_{0,c})\frac{\mu_{B,c}^2}{2}$

Technical Innovations

1. Clear Physical Motivation: Utilizes the multivalued property of entropy density in phase transitions with a clear physical picture
2. Systematic Extrapolation: Establishes a systematic extrapolation framework from zero to finite chemical potential
3. Error Propagation: Implements lattice QCD uncertainty propagation through Monte Carlo sampling
4. Method Verification: Validates method effectiveness within known lattice QCD result ranges

Experimental Setup

Data Sources

Uses continuum extrapolated lattice QCD results from the Wuppertal Budapest collaboration:

• Baryon number susceptibility $\chi_2^B$ data from reference 8
• Entropy density $s$ data from reference 9
• Data covers temperature ranges at zero baryon chemical potential

Data Processing Methods

Employs two independent data analysis approaches:

1. Parametrization Method: Parametric fitting of lattice QCD results
2. Spline Interpolation Method: Smooth spline interpolation of raw data

Error Analysis

• Propagates Gaussian errors from lattice QCD inputs through Monte Carlo sampling
• Calculates 68% confidence intervals
• Does not include systematic assessment of truncation errors

Verification Strategy

Compares and verifies against existing lattice QCD calculations within the range $\mu_B/T \leq 3.5$

Experimental Results

Main Results

Critical Point Location

Using the parametrization method:

• $T_{0,c} = 140.9 \pm 2.0$ MeV
• Critical point: $(T_c, \mu_{Bc}) = (114.3 \pm 6.9, 602.1 \pm 62.1)$ MeV

Using the spline interpolation method:

• Critical point: $(T_c, \mu_{Bc}) = (119.5 \pm 5.4, 556.5 \pm 49.8)$ MeV

Method Consistency Verification

Within the range $\mu_B/T \leq 3.5$, the extrapolation results demonstrate excellent consistency with lattice QCD data, validating the method's reliability in known regions.

Phase Transition Characteristics

At $\mu_B \gtrsim 602$ MeV, entropy density exhibits pronounced S-shaped multivalued behavior, a characteristic signature of first-order phase transitions.

Systematic Error Assessment

Results from the two independent analysis methods are consistent within one standard deviation, indicating method robustness. The confidence intervals from parametrization and spline interpolation methods partially overlap, enhancing result credibility.

Physical Picture

• At $\mu_B = 450$ MeV: crossover transition
• At $\mu_B = 602$ MeV: critical behavior
• At $\mu_B = 750$ MeV: mixed phase

Related Work

Theoretical Methods

Existing methods for locating the QCD critical point include:

• Model-based estimates
• Direct extrapolation from zero and imaginary chemical potential lattice data
• Functional renormalization group methods
• Dyson-Schwinger equation approaches

Advantages of This Work

1. Directly Based on Lattice QCD Data: Avoids model dependence
2. Clear Physical Motivation: Utilizes entropy multivalued property with explicit physical foundation
3. Systematic Error Control: Validates result robustness through multiple methods
4. Experimental Relevance: Results correlate with energy ranges of heavy-ion collision experiments

Conclusions and Discussion

Main Conclusions

1. The proposed constant entropy density contour method successfully predicts the QCD critical point location
2. The critical point is located at $(T_c, \mu_{Bc}) = (114 \pm 7, 602 \pm 62)$ MeV
3. Results indicate that heavy-ion collisions at $\sqrt{s_{NN}} \approx 4-6$ GeV are closest to the critical point region

Limitations

1. Truncation Error: Current analysis is limited to $\mu_B^2$ order expansion; contributions from higher-order terms are unknown
2. Extrapolation Range: Extrapolations beyond lattice QCD verification ranges contain uncertainties
3. Systematic Error: Truncation error impact on results is not fully quantified

Future Directions

1. Higher-Order Corrections: Include $\mu_B^4$ and higher-order terms
2. Imaginary Chemical Potential Analytic Continuation: Improve extrapolation by combining imaginary chemical potential lattice data
3. Direct Simulation: Develop direct lattice QCD methods at non-zero $\mu_B$
4. Experimental Verification: Search for critical point signals within the predicted energy range

In-Depth Evaluation

Strengths

1. Strong Method Innovation: The entropy density contour-based method possesses clear physical motivation and mathematical foundation
2. High Result Reliability: Consistency between two independent analysis methods enhances result credibility
3. Experimental Relevance: Predicted critical point location correlates with current heavy-ion experiment energy ranges
4. Strong Systematicity: Establishes a complete framework from theory to numerical implementation

Weaknesses

1. Theoretical Limitations: $\mu_B^2$ order truncation may be insufficient for accurately describing physics in large chemical potential regions
2. Incomplete Error Estimation: Lacks systematic assessment of truncation errors
3. Extrapolation Assumptions: The validity of constant entropy density contour extrapolation in large $\mu_B$ regions lacks sufficient verification
4. Lattice Data Dependence: Result quality entirely depends on precision of input lattice QCD data

Impact

1. Theoretical Contribution: Provides new theoretical tools for QCD critical point searches
2. Experimental Guidance: Offers specific energy targets for heavy-ion collision experiments
3. Method Generalization: The method can be extended to phase transition studies in other systems
4. Field Advancement: Promotes integration of lattice QCD with heavy-ion physics

Application Scenarios

1. Heavy-Ion Collision Experiments: Guides optimization of operating parameters for RHIC, FAIR and other experimental facilities
2. Theoretical Calculations: Provides benchmark comparisons for other QCD critical point search methods
3. Astrophysics: May provide insights for phase transition studies in neutron star matter
4. Condensed Matter Physics: Method can be generalized to critical phenomena research in other strongly correlated systems

References

This paper cites important literature in the QCD phase diagram research field, including:

• Benchmark lattice QCD data from the Wuppertal Budapest collaboration
• Major recent works on QCD critical point theoretical predictions
• Freeze-out parameter measurements from heavy-ion collision experiments
• Related theoretical methods and numerical technique developments

Overall Assessment: This is a high-quality paper with significant contributions to QCD phase diagram research. The method is novel, results are practically meaningful, and it provides valuable theoretical guidance for experimental searches of the QCD critical point. Despite some theoretical limitations, its innovation and practical utility make it an important advance in the field.