Lattice QCD at finite temperature and density

This paper reviews recent lattice QCD results on strongly interacting matter under extreme conditions, focusing on the finite-temperature transition at zero baryon chemical potential, the search for a critical endpoint at finite density, and thermodynamic properties under external influences such as magnetic fields, rotation, and acceleration.

The Gist

Imagine the universe as a giant, cosmic kitchen. Inside this kitchen, the fundamental ingredients of matter—quarks and gluons—are usually locked together in tight little bundles called protons and neutrons (like ingredients stuck in a jar). But under extreme conditions, like the first few microseconds after the Big Bang or inside the heart of a neutron star, these jars break open. The ingredients spill out into a hot, chaotic soup called Quark-Gluon Plasma (QGP).

This paper is a report from a team of scientists (led by Heng-Tong Ding) who are trying to cook up a perfect recipe for this soup using a supercomputer simulation called Lattice QCD. Since we can't build a real Big Bang in a lab, they use math to simulate the rules of the universe on a giant digital grid.

Here is a breakdown of their latest findings, explained with everyday analogies:

1. The Great Melting (The Phase Transition)

At normal temperatures, quarks are stuck in their jars. As you heat them up, they eventually melt into the soup.

• The Finding: The scientists confirmed that this melting isn't a sudden explosion (like ice shattering); it's a smooth transition, like butter slowly softening into a liquid.
• The Temperature: They calculated the exact "melting point" to be about 156 million degrees (156 MeV).
• The Check: They used different types of digital grids (different "recipes" for the simulation) to make sure they got the same result. It's like three different chefs using different ovens and measuring cups, all agreeing on the exact temperature water boils.

2. The Mystery of the "Ghost" Force ($U_A(1)$ Anomaly)

In the world of subatomic particles, there's a weird rule called the $U_A(1)$ anomaly. Think of it as a "ghost force" that keeps certain particles from behaving symmetrically.

• The Question: When the soup gets hot, does this ghost force disappear?
• The Finding: The scientists looked at the "spectrum" of the particles (like listening to the notes on a piano). They found that even when the soup is hot, there are still faint, ghostly notes (near-zero modes) playing. This suggests the "ghost force" doesn't vanish completely right at the melting point; it lingers a bit, which changes how the soup behaves.

3. The Treasure Hunt for the Critical Point (CEP)

Physicists believe there is a specific spot on the map of temperature and density where the soup changes from a "smooth melt" to a "sudden freeze" (a first-order phase transition). This spot is called the Critical Endpoint (CEP). Finding it is like finding the exact spot on a map where a river turns into a waterfall.

• The Problem: We can't easily simulate high-density soup because the math gets "complex" (a problem called the "sign problem"), making the computer calculations unstable.
• The Progress:
  • They have set a "ceiling" for where this point can be. It's likely not too hot (below 125 MeV).
  • They have ruled out the possibility that this point exists at low densities.
  • They are using clever tricks, like looking at the "Lee-Yang zeros" (mathematical shadows of the phase transition), to guess where the treasure is hidden. Current estimates suggest it's somewhere around 400–600 MeV in density.

4. The Magnetic Field Experiment

Imagine putting this cosmic soup inside a giant magnet.

• The Finding: Strong magnetic fields act like a giant squeeze. They change how the particles move.
• The "Magnetometer": The scientists discovered that by measuring how electric charges and baryons (protons/neutrons) wiggle together in a magnetic field, they can create a "magnetometer." This tool can tell us how strong the magnetic field is inside a heavy-ion collision (like those at the Large Hadron Collider) just by looking at the particle data.
• Charm Quarks: They also looked at "charm" quarks (heavier cousins of the usual quarks). They found that these heavy particles don't just vanish when the soup forms; they slowly turn into free-floating quarks, like heavy rocks slowly dissolving in hot tea.

5. Spinning and Accelerating the Soup

The paper also explored what happens if you spin the soup or accelerate it.

• Rotation: If you spin the soup, the center becomes a liquid (deconfined) while the edges stay solid (confined). It's like a spinning pizza dough where the center stretches out but the crust stays tight.
• Acceleration: If you push the soup hard, it creates a similar split between hot and cold regions, governed by a law of physics called the Tolman-Ehrenfest law (which relates heat to gravity/acceleration).

The Big Picture

This paper is essentially a quality control report for our understanding of the early universe.

1. We know the melting point very precisely now.
2. We know the "ghost force" is still active near the melting point.
3. We are narrowing down the search for the mysterious Critical Endpoint.
4. We are learning how to read the "signatures" of magnetic fields and rotation in the data coming from real-world experiments.

By refining these digital recipes, scientists are getting closer to understanding exactly how the universe evolved from a hot, chaotic soup into the structured matter (stars, planets, and us) we see today.

Technical Summary

1. Problem Statement

Understanding the phase structure of Quantum Chromodynamics (QCD) matter under extreme conditions is a central goal of nuclear and particle physics. Key challenges include:

• The QCD Phase Diagram: Determining the nature of the transition from hadronic matter to the Quark-Gluon Plasma (QGP) at zero baryon chemical potential ($\mu_B = 0$) and mapping the phase boundary at $\mu_B > 0$.
• The Chiral Limit and $U_A(1)$ Anomaly: Clarifying the fate of the axial $U_A(1)$ symmetry at finite temperature, which dictates the universality class of the chiral phase transition and the existence of a critical endpoint (CEP).
• The Sign Problem: The inability to perform standard Monte Carlo simulations at finite baryon density due to the complex fermion determinant.
• External Conditions: Understanding how QCD thermodynamics responds to extreme external fields (magnetic, rotational, acceleration) relevant to heavy-ion collisions and compact stars.

2. Methodology

The paper reviews results primarily derived from Lattice QCD, a non-perturbative approach using discretized spacetime. Key methodological aspects include:

• Discretization Schemes: Comparing results from Staggered fermions (HISQ, 4HEX) and Chirally symmetric fermions (Möbius Domain-Wall Fermions - MDWF, Overlap fermions) to control lattice artifacts and chiral symmetry restoration.
• Continuum Extrapolation: Using multiple lattice spacings ($N_\tau$) to extrapolate results to the continuum limit ($a \to 0$).
• Finite Density Strategies:
  • Taylor Expansion: Expanding thermodynamic quantities around $\mu_B = 0$.
  • Imaginary Chemical Potential: Simulating at imaginary $\mu_B$ and using analytic continuation (e.g., Lee-Yang edge singularities) to infer real $\mu_B$ behavior.
  • Complex Langevin (CL): Stochastic simulations in complexified configuration space, utilizing drift-distribution diagnostics to check for correctness.
  • Diffusion Models & Worldvolume HMC: Emerging machine-learning and geometric approaches to sample complex measures.
• Spectral Analysis: Utilizing the Dirac eigenvalue spectrum ($\rho(\lambda)$) and mass-derivative methods to probe microscopic mechanisms of symmetry breaking.

3. Key Contributions and Results

A. QCD Thermodynamics at $\mu_B = 0$

• Crossover Confirmation: Confirmed that the transition at physical quark masses is a smooth crossover with a pseudocritical temperature $T_{pc} \approx 156\text{--}159$ MeV. Results from MDWF and Staggered fermions are consistent within uncertainties.
• $U_A(1)$ Anomaly Fate:
  • Macroscopic Probes: Investigated via correlator degeneracies ($\pi$ vs. $\sigma$, vector vs. axial-vector) and topological susceptibility ($\chi_t$). Results suggest $U_A(1)$ remains effectively broken at $T_{pc}$ in the chiral limit, though the anomaly weakens at higher temperatures.
  • Microscopic Mechanism: Analysis of the Dirac spectrum reveals a near-zero-mode contribution scaling as $\rho(\lambda \to 0) \propto m_\ell^2$. This supports a scenario where $U_A(1)$ breaking persists at $T_{pc}$, implying the chiral transition for $N_f=2+1$ QCD belongs to the O(4) universality class (second order) in the chiral limit.
• Universality and Scaling: Developed improved order parameters and generalized Banks-Casher relations for higher-order cumulants to better isolate singular scaling behavior near the critical point.

B. Finite Density: CEP and Phase Boundary

• Critical Endpoint (CEP) Constraints:
  • Established an upper bound on the CEP temperature: $T_{CEP} \lesssim 125$ MeV based on the curvature of the chiral critical line.
  • Wuppertal-Budapest Results: Using constant-entropy contours and Lee-Yang zero analysis, they exclude a CEP for $\mu_B < 450$ MeV at the $2\sigma$ level.
  • Estimates: Current lattice-informed estimates place the CEP in the range $T_{CEP} \sim 80\text{--}120$ MeV and $\mu_{B, CEP} \sim 400\text{--}650$ MeV.
• Phase Boundary:
  • Staggered Approach: Used strangeness-neutral proxies (constant $\chi_S^2$ or $\mu_S/\mu_B$) to map the crossover line up to $\mu_B \approx 500$ MeV, finding agreement with Hadron Resonance Gas (HRG) models.
  • Wilson Approach: Used vector-axial-vector degeneracy to determine curvature, yielding results consistent with staggered formulations.

C. Methodological Advances for Complex Measures

• Validated Complex Langevin simulations using drift-distribution fall-off criteria, identifying validity windows and breakdown points near transitions.
• Highlighted emerging techniques combining diffusion models (generative AI) and worldvolume HMC (Lefschetz thimbles) to tackle the sign problem.

D. External Conditions and Thermodynamics

• Strong Magnetic Fields ($eB$):
  • Pions: Neutral pion mass decreases monotonically with $eB$; charged pion mass shows non-monotonic behavior (rise then saturation/turnover) due to internal structure effects.
  • Magnetometer: Proposed $\chi_{BQ}^{11}/\chi_{QS}^{11}$ as a sensitive probe for magnetic fields, showing significant enhancement near the crossover.
  • Topology: Strong fields enhance topological susceptibility at low $T$ but suppress it near/above the crossover (Inverse Magnetic Catalysis).
  • CME: Equilibrium global CME conductivity vanishes; however, non-uniform fields induce localized currents that integrate to zero globally.
• Isospin Chemical Potential ($\mu_I$): Confirmed the transition to a pion-condensed BEC phase at $\mu_I \gtrsim m_\pi/2$ and mapped the Equation of State (EoS) stiffening.
• Rotation and Acceleration: Simulations suggest spatially inhomogeneous confinement-deconfinement structures (e.g., deconfined cores in rotating systems, boundaries following the Tolman-Ehrenfest law in accelerated frames).

4. Significance

• Precision Benchmarking: The work solidifies the crossover picture at $\mu_B=0$ and provides high-precision continuum-extrapolated data for the Equation of State, crucial for interpreting heavy-ion collision data (e.g., STAR, ALICE).
• Theoretical Clarity: The Dirac spectrum analysis provides a microscopic justification for the persistence of $U_A(1)$ breaking, refining our understanding of the chiral phase transition's universality class.
• Phenomenological Bridge: By calculating conserved charge correlations (cumulants) and chemical potentials ($\mu_S, \mu_Q$) under realistic constraints (strangeness neutrality), the paper creates a direct quantitative link between lattice predictions and experimental observables in relativistic heavy-ion collisions.
• New Physics Frontiers: The exploration of QCD under rotation, acceleration, and strong magnetic fields opens new avenues for understanding matter in neutron stars and the early universe, while the development of methods for complex measures (CL, AI) offers a path toward solving the sign problem at high density.

In summary, this review highlights that Lattice QCD has evolved from a qualitative tool to a precision instrument capable of mapping the QCD phase diagram, resolving symmetry restoration mechanisms, and predicting the behavior of matter under extreme external deformations.