Fact-checked by Grok 2 weeks ago

QCD matter

QCD matter, also known as quark matter, encompasses the various phases of strongly interacting matter where quarks and gluons serve as the fundamental degrees of freedom, governed by the theory of quantum chromodynamics (QCD).^[1] This includes the quark-gluon plasma (QGP), a deconfined state of quarks and gluons that emerges at extreme conditions of high temperature (above approximately 150–160 MeV) or high density, contrasting with the confined hadronic phase at lower temperatures and densities where quarks are bound into hadrons like protons and neutrons.^[2]^[1] In the early universe, shortly after the Big Bang, QCD matter existed as a hot QGP before cooling and expanding led to hadronization, the process by which quarks and gluons recombine into hadrons.^[2] Today, this state is recreated in laboratories through relativistic heavy-ion collisions at facilities such as the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory and the Large Hadron Collider (LHC) at CERN, where colliding heavy nuclei like gold or lead generate initial energy densities of 12–20 GeV/fm³, producing a transient QGP fireball that expands hydrodynamically before freezing out into observable hadrons.^[3]^[1] Key properties of QCD matter include deconfinement, where the strong force no longer confines quarks within hadrons, and chiral symmetry restoration, marking the transition from massive to effectively massless quarks at the critical temperature T_c.^[2] At high densities, such as those in neutron star cores, QCD matter may exhibit exotic phases like color superconductivity, where quarks pair up analogously to superconductivity in condensed matter.^[1] Theoretical studies employ lattice QCD simulations for zero baryon density, effective models for finite density, and holographic duality for strongly coupled regimes, while experimental probes include jet quenching, elliptic flow, and particle multiplicity fluctuations to map the QCD phase diagram and search for a critical endpoint.^[2]^[3]^[1]

Fundamentals of QCD

Quarks, Gluons, and Strong Interaction

Quarks are the fundamental fermions that serve as the building blocks of hadronic matter in quantum chromodynamics (QCD), carrying a non-Abelian color charge that comes in three types: red, green, or blue.^[4] These color charges are confined to the fundamental representation of the SU(3) gauge group, ensuring that physical particles, such as protons and neutrons, are color singlets formed by combinations of quarks.^[4] Gluons are the massless vector bosons responsible for mediating the strong interaction between quarks, analogous to photons in quantum electrodynamics but distinguished by carrying color charge themselves.^[4] There are eight gluons, corresponding to the adjoint representation of SU(3), which enables them to couple to both quarks and other gluons.^[4] The underlying framework of QCD is a non-Abelian gauge theory based on the SU(3)_c color symmetry group, where the subscript c denotes color.^[4] Unlike the Abelian U(1) group of electromagnetism, the non-Abelian structure of SU(3)_c introduces self-interactions among gluons through the structure constants of the group, leading to three- and four-gluon vertices that are fundamental to the dynamics of the strong force. The dynamics of quarks and gluons are described by the QCD Lagrangian density:

\mathcal{L}_{\rm QCD} = \bar{\psi} (i \gamma^\mu D_\mu - m) \psi - \frac{1}{4} G^a_{\mu\nu} G^{a \mu\nu},

where \psi represents the quark fields, m is the quark mass, D_\mu = \partial_\mu - i g_s t^a A^a_\mu is the covariant derivative incorporating the strong coupling g_s and gluon fields A^a_\mu (with t^a as the SU(3) generators), and G^a_{\mu\nu} is the gluon field strength tensor.^[4] This form captures both the kinetic and interaction terms for quarks and gluons, with the non-Abelian term -\frac{1}{4} G^a_{\mu\nu} G^{a \mu\nu} explicitly including gluon self-couplings. Quarks exist in six flavors—up (u), down (d), strange (s), charm (c), bottom (b), and top (t)—each transforming under the same SU(3)_c color group but distinguished by their masses, which span more than five orders of magnitude and influence their role in QCD matter.^[4] The lighter up, down, and strange quarks have masses on the order of a few MeV, relevant for low-energy hadronic physics, while the heavier charm, bottom, and top quarks have masses of GeV scale, with the top quark being the heaviest at approximately 173 GeV and decaying before forming hadrons.^[5]^[6]

Quark Flavor	Mass (MSbar scheme, approximate value)
Up (u)	2.2 MeV (at μ = 2 GeV)
Down (d)	4.7 MeV (at μ = 2 GeV)
Strange (s)	93 MeV (at μ = 2 GeV)
Charm (c)	1.27 GeV (at μ = m_c)
Bottom (b)	4.18 GeV (at μ = m_b)
Top (t)	172.6 GeV (pole mass)

These masses are running parameters in the modified minimal subtraction (MSbar) scheme unless noted, reflecting the scale-dependent nature of QCD.^[5]^[6]

Confinement and Asymptotic Freedom

In quantum chromodynamics (QCD), the strong interaction exhibits two fundamental and contrasting properties: asymptotic freedom at short distances and confinement at long distances. Asymptotic freedom implies that the effective strength of the interaction between quarks and gluons diminishes as the energy scale increases, allowing perturbative calculations at high energies. This behavior arises from the non-Abelian nature of the SU(3) gauge theory underlying QCD, where gluons carry color charge and self-interact. The running coupling constant \alpha_s(Q) = \frac{g_s^2}{4\pi}, where g_s is the strong coupling and Q is the momentum transfer scale, decreases logarithmically with increasing Q. This scale dependence is governed by the renormalization group beta function, whose leading-order form is \beta(g) = -\frac{g^3}{16\pi^2} \left( \frac{11}{3} N_c - \frac{2}{3} N_f \right) + \mathcal{O}(g^5), with N_c = 3 colors and N_f = 6 quark flavors. The negative sign of the first coefficient ensures that \alpha_s(Q) approaches zero in the ultraviolet limit (Q \to \infty), enabling quarks and gluons to behave as nearly free particles at sufficiently short distances, on the order of $10^{-16} m or less.^[7] The discovery of asymptotic freedom was independently reported in 1973 by David Gross and Frank Wilczek, and by Hugh David Politzer, resolving a long-standing puzzle in strong interaction physics by providing a framework for QCD as the theory of hadronic matter. Their work, which demonstrated that a non-Abelian gauge theory could be asymptotically free while remaining renormalizable, earned them the 2004 Nobel Prize in Physics. This breakthrough shifted the understanding of the strong force from a phenomenological model to a predictive quantum field theory. In stark contrast, at low energies and large distances (infrared regime), QCD displays confinement, or "infrared slavery," where the coupling \alpha_s grows, preventing the isolation of individual quarks or gluons as free particles. Quarks are eternally bound into color-neutral hadrons, such as mesons and baryons, due to the formation of a chromoelectric flux tube between color charges. This leads to a linear interquark potential V(r) \sim \sigma r, where r is the separation and \sigma \approx 1 GeV/fm is the string tension, as confirmed by lattice QCD simulations modeling the vacuum as a dual superconductor that squeezes color fields into thin tubes.^[8] The interplay between asymptotic freedom and infrared slavery dictates the dynamics of QCD matter: high-energy processes allow quark-gluon interactions to be treated perturbatively, while low-energy phenomena require non-perturbative methods, culminating in hadronization—the irreversible formation of hadrons from deconfined quarks and gluons as the system cools or expands. This dual behavior underpins the transition between free quark-gluon states and confined hadronic matter.^[9]^[10]

States of QCD Matter

Hadronic Phase

The hadronic phase of QCD matter represents the low-temperature and low-density regime where quarks are confined within color-neutral hadrons due to the strong interaction's non-perturbative effects. In this phase, observable particles are exclusively color singlets, formed by combinations of quarks and antiquarks that neutralize the SU(3) color charge. Baryons, such as the proton, consist of three quarks (qqq) in a color-antitriplet state that combines to a singlet, while mesons, like the pion, are quark-antiquark (q\bar{q}) pairs in a color octet-antioctet configuration yielding a singlet overall.^[11] At low temperatures, the hadronic phase is characterized by spontaneous chiral symmetry breaking, where the approximate SU(3)_L × SU(3)_R symmetry of massless QCD is broken to the diagonal SU(3)_V flavor symmetry. This breaking generates a nonzero quark bilinear condensate \langle \bar{q}q \rangle \approx -(250 \mathrm{MeV})^3, representing the vacuum expectation value of the scalar quark density and serving as the order parameter for chiral symmetry.^[12] The condensate arises from the dynamical generation of constituent quark masses, on the order of 300-400 MeV, far exceeding the current quark masses (a few MeV), and leads to the emergence of light pseudoscalar mesons as approximate Goldstone bosons, such as the pions with masses around 140 MeV.^[13] The QCD vacuum in the hadronic phase exhibits rich non-perturbative structure that underpins confinement, including contributions from glueballs and instantons. Glueballs are hypothetical bound states of gluons, predicted as color singlets with the lightest scalar glueball having a mass around 1.5-1.7 GeV, influencing the glue content of the vacuum and supporting the string-like flux tubes between quarks.^[14] Instantons, as topologically nontrivial gluon field configurations, contribute to the vacuum energy and induce chiral symmetry breaking by aligning quark zero modes, with the instanton liquid model describing a dilute ensemble of these objects with average size \rho \approx 1/600 MeV^{-1} and density n \approx 1 fm^{-3}. As temperature increases toward the pseudocritical scale of approximately 150-170 MeV, the hadronic phase transitions to higher-energy states through the melting of the quark condensate, where \langle \bar{q}q \rangle decreases rapidly, signaling the restoration of chiral symmetry.^[15] This melting disrupts the bound hadronic structure, allowing quarks to become less confined, though the phase remains dominated by hadronic degrees of freedom below the transition.

Quark-Gluon Plasma

The quark-gluon plasma (QGP) is a deconfined state of quantum chromodynamics (QCD) matter at high temperatures, where quarks and gluons propagate freely as partons over distances exceeding the typical hadron scale of approximately 1 fm. This phase emerges because asymptotic freedom renders the strong coupling weak at short distances, suppressing confinement effects and allowing the fundamental color charges to exist unbound. The QGP was first proposed theoretically by Collins and Perry in 1975, who recognized that at sufficiently high temperatures, the perturbative nature of QCD would permit a plasma of asymptotically free quarks and gluons, analogous to the ionized state of electromagnetic plasmas. Deconfinement in the QGP is quantified by the order parameter known as the Polyakov loop, a gauge-invariant Wilson line in the temporal direction whose vacuum expectation value \langle L \rangle \approx 1 signals the restoration of the Z(3) center symmetry and the absence of confinement. In this state, the active degrees of freedom consist of quarks with 2 spin states, 3 colors, and N_f flavors, alongside gluons with 2 transverse polarizations and 8 color-octet states, yielding a total degeneracy that approaches the Stefan-Boltzmann limit for thermodynamic quantities. Specifically, the energy density follows \epsilon = \frac{\pi^2}{30} g_* T^4, where the effective number of relativistic degrees of freedom g_* \approx 47.5 for N_f = 3 massless quark flavors (up, down, and strange), reflecting the bosonic contribution from gluons (16) and the fermionic contribution from quarks adjusted by the factor $7/8.^[16] The QGP manifests as a nearly ideal fluid, characterized by an extraordinarily low shear viscosity-to-entropy density ratio \eta/s \approx 0.1, approaching the universal lower bound conjectured from gauge/gravity duality and indicating minimal dissipation during collective expansion. Indirect signatures of this deconfined medium include jet quenching, the energy loss of hard partonic jets traversing the plasma via medium-induced gluon radiation and elastic scattering, which suppresses high-transverse-momentum hadron yields relative to vacuum fragmentation. Complementing this, elliptic flow quantifies the anisotropic pressure gradients in the expanding plasma, producing a second-order Fourier coefficient v_2 in particle azimuthal distributions that reflects the medium's rapid thermalization and hydrodynamic response.

Exotic Phases at High Density

At high baryon densities and low temperatures, quantum chromodynamics (QCD) predicts the emergence of exotic phases of matter distinct from the quark-gluon plasma, characterized by novel forms of quark pairing and partial deconfinement. These phases arise in the cold, dense regime of the QCD phase diagram, where the strong interaction leads to states with superconducting properties or hybrid confinement-deconfinement behaviors, potentially realized in the cores of neutron stars. Color superconductivity represents one such exotic phase, where quarks form Cooper pairs analogous to the Bardeen-Cooper-Schrieffer (BCS) mechanism in conventional superconductors, but mediated by gluon exchange. In the color-flavor-locked (CFL) phase, up, down, and strange quarks pair in a way that locks color and flavor symmetries, breaking them spontaneously to a diagonal subgroup and generating a gap \Delta \sim 10{-}100 MeV in the quark spectrum. This pairing leads to color Meissner effects, expelling magnetic fields, and renders the ground state a superfluid with both color and electromagnetic superconductivity. The CFL phase is favored at asymptotically high densities where perturbative QCD applies, but at moderate densities, mismatched Fermi surfaces due to unequal quark masses can lead to alternative patterns like the 2SC phase, where only up and down quarks pair.^[17] Quarkyonic matter proposes another paradigm for high-density QCD, featuring a deconfined quark Fermi sea at the core surrounded by a shell of confined baryonic excitations. At large baryon chemical potential \mu_B \gtrsim 1 GeV, the pressure is dominated by quark degrees of freedom near the Fermi surface, while confinement persists for excitations above it, preventing free quarks but allowing a large baryon density. This state reconciles asymptotic freedom at short distances with confinement at long distances, emerging in the large-N_c limit of QCD where N_c is the number of colors. Model calculations indicate quarkyonic matter stiffens the equation of state compared to pure hadronic matter, influencing the structure of compact stars. Recent theoretical advances, including model studies up to 2025, suggest intermediate states between the hadronic phase and full quark-gluon plasma at high density, exhibiting non-conformal behaviors such as speed-of-sound variations exceeding the conformal limit c_s^2 = 1/3. These proposals, motivated by effective models and analogies from lattice QCD in related theories, describe phases with partial chiral restoration or quarkyonic-like transitions, potentially bridging confined and deconfined regimes without sharp boundaries. Hybrid stars, incorporating quark cores amid hadronic mantles, provide an astrophysical context for these exotic phases, with observations of massive neutron stars (M \gtrsim 2 M_\odot) constraining the transition density to quark matter around 2-5 times nuclear saturation. Calculations using perturbative QCD or Nambu-Jona-Lasinio models yield stable hybrid configurations where color-superconducting quark cores occupy 10-20% of the star's radius, enhancing maximum masses and altering cooling via neutrino emission from paired quarks. These structures highlight the role of high-density QCD in explaining pulsar timing and gravitational wave signals.

Occurrence and Production

Cosmological Contexts

In the standard Big Bang model, the early universe transitioned through a quark-gluon plasma (QGP) phase shortly after the initial singularity, where quarks and gluons existed in a deconfined state due to high temperatures. This phase dominated from approximately $10^{-12} seconds, corresponding to temperatures around 100 GeV near the electroweak scale, to about $10^{-5} seconds at temperatures of roughly 150 MeV, marking the onset of the QCD phase transition.^[18] During this interval, the universe expanded rapidly while cooling, with the QGP serving as the primordial state of matter before hadronization. The shift from the electroweak era to the QCD-dominated QGP occurred as temperatures dropped from about 100 GeV to 150 MeV, a period spanning roughly $10^{-12} to $10^{-5} seconds. Lattice QCD simulations indicate that the QCD phase transition at this scale is a smooth crossover rather than a sharp first-order event, involving gradual confinement of quarks into hadrons without significant supercooling or bubble nucleation. This crossover influences baryogenesis processes, as sphaleron transitions—key to generating the observed baryon asymmetry—remain active until the electroweak scale but can be modulated by the evolving QCD dynamics, potentially affecting the efficiency of lepton-to-baryon conversion in extensions of the Standard Model. The QCD phase transition also impacts cosmological relic abundances by introducing potential density fluctuations during hadronization, which occur post-QGP at temperatures below 150 MeV. These inhomogeneities, arising from the release of latent heat and changes in the equation of state, could alter baryon distributions on small scales, influencing the subsequent big bang nucleosynthesis (BBN) around 1 MeV. For instance, even in a crossover scenario, such effects might subtly modify the predicted abundances of light elements like deuterium and helium-4 formed after full hadronization, though current observations constrain significant deviations.^[19]^[20] In inflationary cosmology, the reheating phase following the rapid expansion driven by the inflaton field can produce a primordial QGP if the reheating temperature exceeds the QCD scale of approximately 150 MeV, restoring thermal equilibrium and populating the universe with deconfined quarks and gluons. This process, occurring at temperatures potentially up to 10^{15} GeV or higher depending on the model, ensures the hot Big Bang conditions necessary for the QGP era, with non-equilibrium dynamics facilitating fast thermalization into the plasma state.^[21]

Astrophysical Environments

In the cores of neutron stars, where densities exceed several times the nuclear saturation density \rho_0 \approx 0.16 fm^{-3}, the conditions may favor the formation of quark matter, potentially transitioning from hadronic matter to deconfined quark-gluon states at densities greater than 5–10 \rho_0.^[22] This transition is theorized to stiffen the equation of state (EOS), allowing neutron stars to support masses up to approximately 2 solar masses (M_\odot) or higher, as softer purely hadronic EOS would otherwise lead to collapse. Observational evidence from pulsar timing and X-ray spectroscopy supports the existence of such stiff EOS in massive neutron stars, consistent with hybrid configurations featuring quark cores.^[23] The strange quark matter hypothesis posits that a mixture of up, down, and strange quarks could form the absolute ground state of baryonic matter, with stable strangelets—small clumps of this matter—having an energy per baryon lower than that of iron (around 930 MeV), rendering ordinary nuclei unstable against conversion. Proposed originally by Witten, this idea suggests that neutron star remnants from supernovae could convert to strange quark stars if the surface energy barrier is overcome, potentially explaining compact objects with unusual cooling or mass-radius relations.^[24] Stability analyses within chiral quark models confirm that such matter remains bound at zero pressure, unlike neutron matter.^[25] Phenomena in magnetars, highly magnetized neutron stars, such as giant flares and spin glitches, may serve as indirect probes of QCD phase transitions in their interiors. Sudden energy releases during flares could arise from magnetic field rearrangements triggering a hadron-to-quark phase transition, releasing latent heat and altering the star's rotation. Glitches, observed as abrupt spin-ups, might similarly reflect density fluctuations or superfluid responses near a first-order phase boundary, providing constraints on the EOS stiffness.^[26] Recent analyses incorporating gravitational wave data from GW170817 have tightened constraints on hybrid star models, favoring those with quark cores that match the observed tidal deformability while accommodating massive pulsars like PSR J0740+6620 (\sim 2.08 M_\odot).^[27] By 2025, multimessenger observations, including updated radius measurements from NICER, have tightened constraints on the EOS, excluding soft purely hadronic models and favoring stiffer hadronic or hybrid scenarios with quark matter components for stars above 1.4 M_\odot, though nonstrange quark cores remain viable alternatives to strange matter.^[28] Recent 2025 NICER measurements of pulsars like PSR J0437-4715 and PSR J0614-3329 have further constrained radii to ~12-13 km for 1.4 M_\odot stars, supporting stiffer EOS consistent with hybrid quark-hadron configurations.^[29] These constraints highlight the role of astrophysical environments in testing QCD at extreme densities.^[30]

Laboratory Creation

QCD matter, particularly the quark-gluon plasma (QGP), is primarily created in laboratory settings through high-energy collisions at particle accelerators, where extreme temperatures and densities briefly recreate conditions akin to the early universe. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory initiated gold-gold (Au-Au) heavy-ion collisions in 2000, producing a hot, dense medium interpreted as QGP in central collisions, with the initial overlap volume estimated at approximately 10 fm³.^[31] Similarly, the Large Hadron Collider (LHC) at CERN began lead-lead (Pb-Pb) collisions in November 2010 at a center-of-mass energy per nucleon pair of √s_{NN} = 2.76 TeV, generating even hotter and denser QGP states within comparable local volumes of about 10 fm³.^[32] These collisions involve relativistic heavy ions, where the Lorentz-contracted nuclei overlap to form a deconfined plasma that expands and cools rapidly over femtoseconds. To map the QCD phase diagram, particularly at finite baryon chemical potential (μ_B), RHIC's Beam Energy Scan (BES) program systematically varies collision energies from low values up to √s_{NN} = 200 GeV, allowing probes of increasing μ_B up to around 400 MeV in the most central Au-Au events.^[33] This multi-phase effort, launched in 2010, has collected extensive datasets across energies like 7.7, 11.5, 19.6, 27, 39, and 200 GeV, enabling studies of the transition from hadronic matter to QGP and potential critical points.^[34] In parallel, smaller systems such as proton-lead (p-Pb) and proton-proton (pp) collisions at the LHC have revealed signatures of mini-QGP droplets, where high-multiplicity events produce compact, transient deconfined regions with volumes orders of magnitude smaller than in heavy-ion collisions, yet exhibiting collective behaviors indicative of QGP formation.^[35]^[36] Future facilities will extend these investigations to higher baryon densities. The Facility for Antiproton and Ion Research (FAIR) at GSI Helmholtz Centre, with construction advancing since 2018, anticipates first heavy-ion beams for QCD experiments in 2028 using the SIS100 accelerator to probe dense matter at μ_B up to 1 GeV.^[37] Likewise, the Nuclotron-based Ion Collider fAcility (NICA) at the Joint Institute for Nuclear Research (JINR) in Dubna, under development since the early 2010s, is slated for operational heavy-ion collisions around 2025, focusing on high-density QCD probes in Au-Au interactions at √s_{NN} up to 11 GeV to explore the high-μ_B regime.^[38] These accelerators will complement RHIC and LHC by accessing longer-lived, denser QCD matter states.

Phase Diagram

Temperature and Chemical Potential Axes

The phase diagram of quantum chromodynamics (QCD) matter is conventionally mapped in the plane defined by temperature T and baryon chemical potential \mu_B, which together parameterize the thermodynamic conditions of strongly interacting systems in thermal equilibrium. The temperature T, measured in units of energy such as MeV, represents the thermal energy scale available to excitations, with relevant QCD scales spanning from near zero up to several hundred MeV, comparable to the QCD scale \Lambda_\mathrm{QCD} \approx 200 MeV. The baryon chemical potential \mu_B, also in MeV, controls the net baryon number density and extends from zero (corresponding to baryon-symmetric matter with equal numbers of baryons and antibaryons) to values around 1 GeV, the approximate scale of the nucleon mass, where dense matter relevant to astrophysical objects like neutron stars becomes accessible. For quark-level descriptions, the quark chemical potential is \mu_q = \mu_B / 3, reflecting the tripling of baryon number for three-quark constituents. Lattice QCD simulations provide the primary non-perturbative tool for exploring this diagram, particularly along the \mu_B = 0 axis, where the pseudo-critical temperature for the transition from confined to deconfined matter is determined to be T_\mathrm{pc} = 156.5 \pm 1.5 MeV for physical quark masses. At finite \mu_B, direct simulations face the severe sign problem: the fermion determinant in the path integral becomes complex for real \mu_B \neq 0, preventing efficient Monte Carlo sampling and requiring alternative approaches such as analytic continuation from imaginary chemical potentials \mu_B = i \mu_I, where the determinant remains real and positive. This workaround exploits the periodicity and analyticity of the partition function in imaginary \mu, allowing extrapolations to real values, though with increasing uncertainty at larger \mu_B / T. In the T-\mu_B plane, the low-T, low-\mu_B region describes a hadron gas of confined quarks within color-neutral mesons and baryons, while high T at moderate \mu_B yields a quark-gluon plasma of asymptotically free partons. At high \mu_B and low T, exotic phases such as color-superconducting quark matter emerge due to pairing instabilities in dense fermionic systems. The zero-\mu_B line thus serves as a reference for symmetric nuclear matter, with deviations probing net baryon densities encountered in heavy-ion collisions and compact stars.

Phase Transitions and Critical Points

In the QCD phase diagram, the transition from the hadronic phase to the quark-gluon plasma at zero baryon chemical potential \mu_B = 0 manifests as a rapid crossover rather than a sharp phase transition. Lattice QCD simulations determine the pseudocritical temperature for this chiral crossover at T_c \approx 155 MeV for physical quark masses with 2+1 flavors.^[39] This smooth behavior arises because explicit chiral symmetry breaking due to nonzero light quark masses prevents a true second-order transition, consistent with expectations from the O(4) universality class in the chiral limit where a second-order transition would occur at a lower T_c \approx 132 MeV.^[39]^[40] At higher baryon densities, the crossover is expected to evolve into a first-order phase transition line, terminating at a critical endpoint (CEP) where the transition becomes second-order. Theoretical predictions from effective models and holographic approaches place the CEP at \mu_B \sim 600 MeV and T \sim 100 MeV, marking the boundary between the crossover region and the first-order regime.^[41] This point belongs to the 3D Ising universality class, influencing observables like baryon number fluctuations near the transition.^[42] The chiral transition, associated with the restoration of approximate chiral symmetry SU(2)_L × SU(2)_R, aligns closely with the pseudocritical line at low \mu_B but weakens into a crossover due to explicit breaking by quark masses. At imaginary chemical potentials \mu = i \mu_I, QCD exhibits Roberge-Weiss periodicity with period $2\pi T / 3 in \mu_I / T, arising from the Z(3) center symmetry of pure Yang-Mills theory, leading to phase transitions between different Polyakov loop sectors that intersect with the chiral transition line.^[43]^[44] Recent lattice QCD studies in 2025, leveraging improved algorithms for continuum extrapolation and handling finite-volume effects under strangeness neutrality, provide hints on the CEP location by excluding its presence at \mu_B < 450 MeV at the 2\sigma level, suggesting it resides at higher densities consistent with model predictions.^[45] These advances refine the phase boundary mapping, enhancing constraints on the diagram's structure without direct sign problem resolution.^[45]

Thermodynamic Properties

Equation of State

The equation of state (EOS) of QCD matter describes the thermodynamic relations between pressure P, energy density \varepsilon, temperature T, and chemical potentials \mu, such as P = P(T, \mu). Lattice QCD simulations provide non-perturbative calculations of the EOS at zero or small chemical potentials, revealing a smooth crossover transition from hadronic matter to quark-gluon plasma (QGP) around the pseudocritical temperature T_c \approx 155 MeV for \mu = 0. For instance, the normalized pressure reaches P/T^4 \approx 0.2 (relative to the Stefan-Boltzmann limit) at T = 155 MeV and \mu = 0, indicating partial deconfinement with interactions suppressing the pressure below the ideal gas value.^[46] In the high-temperature limit, QCD matter approaches the conformal limit of a massless ideal gas, where \varepsilon = 3P, corresponding to the trace anomaly vanishing as \varepsilon - 3P \to 0. Deviations from this relation arise from quark masses, non-perturbative effects near T_c, and residual interactions, with lattice results showing \varepsilon/T^4 \approx 0.6 at T = 155 MeV, leading to \varepsilon \approx 2.8 P. At higher temperatures, say T > 300 MeV, the EOS closely tracks the Stefan-Boltzmann values for 2+1 flavors, with P/T^4 approaching unity in normalized units.^[46] The speed of sound squared, defined as c_s^2 = dP/d\varepsilon, provides insight into the EOS stiffness and exhibits a characteristic dip near the transition due to the softening of the medium. Lattice calculations show c_s^2 \approx 0.2 near T_c, reflecting strong interactions and the latent heat-like behavior in the crossover, before rising toward the conformal value of $1/3 in the QGP phase at T \gtrsim 200 MeV. This minimum highlights the transition's impact on hydrodynamic evolution in heavy-ion collisions.^[46] At high densities, relevant for the cores of neutron stars, the EOS is probed by effective models since direct lattice calculations are challenging due to the sign problem. These models, such as the Nambu-Jona-Lasinio (NJL) approach or perturbative QCD, predict a stiff EOS for deconfined quark matter, with c_s^2 approaching or exceeding $1/3 at baryon densities \gtrsim 5 n_0 (where n_0 is nuclear saturation density), enabling support for massive compact stars up to 2 solar masses.

Transport Coefficients

Transport coefficients describe the dissipative response of QCD matter, such as the quark-gluon plasma (QGP), to external gradients, playing a key role in its hydrodynamic evolution. In the QGP phase, these include shear viscosity η, which governs momentum diffusion; bulk viscosity ζ, related to volume changes; electrical conductivity σ, characterizing charge transport; and baryon diffusion coefficient D_B, describing net baryon number propagation. These properties arise from interactions among quarks and gluons, with values determined theoretically via lattice QCD, effective models, and dualities like AdS/CFT, and constrained by collective flow in heavy-ion collisions. Shear viscosity η quantifies the fluid's resistance to shear deformations, often normalized as the ratio η/s to the entropy density s for scale invariance. The AdS/CFT correspondence predicts a universal lower bound η/s ≥ 1/(4π) ≈ 0.08 for strongly coupled relativistic fluids, derived from gravitational perturbations in anti-de Sitter spacetime dual to conformal field theories. In the QGP, this bound is nearly saturated, with hydrodynamic analyses of heavy-ion collision data yielding η/s ≈ 0.1–0.5, indicating a low-viscosity, nearly perfect fluid behavior distinct from the ideal limit where η = 0 as per the equation of state.^[47] Bulk viscosity ζ measures resistance to uniform compression or expansion and vanishes in conformally invariant theories but is nonzero in QCD due to scale symmetry breaking from the running coupling and quark masses. Near the pseudocritical temperature T_c ≈ 155 MeV, where conformal invariance is most violated during the hadron-QGP transition, ζ exhibits a pronounced peak, with lattice QCD computations showing ζ/s reaching ~1, far exceeding shear viscosity in this regime.^[48] This enhancement reflects rapid changes in the trace anomaly and equation of state across the crossover. Electrical conductivity σ and baryon diffusion D_B are derived from linear response theory via Kubo formulas, relating them to low-frequency limits of retarded correlation functions of conserved currents. For σ, the Kubo relation σ = (1/3) lim_{ω→0} (1/ω) Im G^R_{ii}(ω,0), where G^R is the electromagnetic current correlator, yields lattice QCD values of σ/T ≈ 0.1–0.4 in the QGP at temperatures above T_c, scaling with the strong coupling. Similarly, D_B emerges from the baryon current correlator, with recent lattice results indicating D_B T ≈ 1–2 at high temperatures, decreasing near finite baryon density due to enhanced scattering. Recent 2025 hydrodynamic studies of heavy-ion collisions provide stringent lower bounds on these coefficients, confirming η/s ≳ 0.08 from flow suppression patterns and aligning with AdS/CFT predictions, while constraining ζ/s < 0.5 away from T_c to match expansion dynamics.

Theoretical Approaches

Lattice QCD Simulations

Lattice QCD offers a rigorous, non-perturbative approach to investigating the properties of QCD matter by discretizing the theory on a hypercubic lattice in four-dimensional Euclidean spacetime, allowing numerical simulations via Monte Carlo methods. The continuum QCD action S = \int d^4x \, \mathcal{L}, where \mathcal{L} includes gluon and quark kinetic terms, is replaced by a discrete sum over lattice sites separated by spacing a, with the continuum limit recovered as a \to 0. Gluons are represented by SU(3) link variables U_\mu(x) on the links, while quarks are described using fermion discretizations such as Wilson fermions or staggered fermions. Wilson fermions incorporate a non-derivative Wilson term to eliminate the 15 unphysical doubler modes that arise in naive discretizations, ensuring a single continuum fermion species per flavor, though at the cost of explicit chiral symmetry breaking that requires careful tuning. Staggered fermions, by contrast, preserve a remnant chiral symmetry and reduce doublers to four "tastes" per flavor, which must be taken to the fourth root in the action to match the continuum with N_f degenerate flavors, enabling efficient simulations of light quarks. These formulations allow computation of thermodynamic observables, such as the pressure and energy density, from the partition function Z = \int \mathcal{D}U \, \det M \, e^{-S_g}, where S_g is the pure gauge action (e.g., the plaquette or improved actions) and M the fermion matrix. Phase transitions in QCD matter are identified through lattice observables sensitive to symmetry changes. Deconfinement, associated with the breaking of Z(3) center symmetry in the pure gauge limit, is probed via the renormalized Polyakov loop \langle L \rangle, whose non-zero value above the transition signals quark liberation; the corresponding susceptibility \chi_L = \partial^2 \ln Z / \partial \beta^2 (with \beta = 6/g^2) exhibits a peak marking the pseudocritical temperature.^[49] Chiral symmetry restoration is monitored by the quark chiral condensate \langle \bar{\psi} \psi \rangle, which acts as an order parameter in the chiral limit and decreases rapidly near the transition, with its susceptibility providing another indicator of the crossover. A major challenge in lattice simulations of QCD matter arises at finite baryon chemical potential \mu_B \neq 0, where the fermion determinant \det M(\mu) becomes complex due to the imaginary part from the chemical potential term, leading to the infamous sign problem that renders standard Monte Carlo importance sampling inefficient as the phase factor oscillates wildly.^[50] This issue is addressed through indirect methods, such as Taylor expansion of observables (e.g., pressure) in powers of \mu_B / T around \mu_B = 0, where coefficients are computed directly, or reweighting techniques that use configurations at \mu_B = 0 to estimate ratios at small finite \mu_B, though both approaches have limitations in convergence radius and computational cost.^[51] Recent advances in lattice QCD, including highly improved staggered quark actions, multi-level integration algorithms, and exascale computing resources, have facilitated precise continuum extrapolations on finer lattices with physical quark masses. These efforts have refined the pseudocritical temperature for the chiral-deconfinement crossover at zero chemical potential to T_c = 156.5 \pm 1.5 MeV, consistent with features of the QCD phase diagram such as the absence of a critical point at low \mu_B.

Perturbative QCD and Weak Coupling

At temperatures T \gg \Lambda_\mathrm{QCD} \approx 200 MeV, where \Lambda_\mathrm{QCD} is the intrinsic QCD scale setting the onset of non-perturbative effects, the running strong coupling g(T) becomes weak due to asymptotic freedom, enabling a perturbative expansion of QCD matter properties using Feynman diagrams and weak-coupling techniques. This regime corresponds to the quark-gluon plasma phase at asymptotically high temperatures, where the deconfined quarks and gluons behave as a weakly interacting gas, but collective plasma effects necessitate resummations to handle infrared sensitivities arising from long-range interactions.^[52] Naive perturbative calculations encounter infrared divergences from soft modes with momenta of order gT, which are resolved through the hard thermal loop (HTL) resummation scheme. HTL perturbation theory reorganizes the expansion by resummed self-energies into effective propagators and vertices that capture the leading plasma physics, such as Debye screening of static electric fields. The resulting Debye screening mass is m_D \sim gT, specifically m_D = gT \sqrt{N_c/3} for pure gluodynamics with N_c = 3, which exponentially suppresses the Yukawa potential between color charges at distances larger than $1/(gT).90508-B) This resummation systematically improves convergence, allowing computations of thermodynamic and transport properties beyond leading order. A key transport observable in this weak-coupling framework is the jet quenching parameter \hat{q}, which measures the average transverse momentum squared acquired per unit length by a high-energy parton propagating through the plasma due to elastic and inelastic scatterings. Perturbative evaluations yield \hat{q} \sim g^4 T^3 \ln(1/g), with the leading logarithmic term originating from multiple soft gluon exchanges screened at the Debye scale, while the overall scaling reflects the density of scattering centers proportional to T^3. This parameter encodes medium-induced energy loss and has been computed to next-to-leading order, highlighting the importance of HTL resummation for the soft sector. To address non-perturbative magnetic screening at even longer distances \sim 1/(g^2 T), dimensional reduction provides a systematic effective field theory approach. The hard Matsubara modes with frequencies \sim T are first integrated out, yielding electrostatic QCD (EQCD), a dimensionally reduced 3D SU(N_c) gauge theory coupled to an adjoint Higgs field representing the temporal gluon component A_0, with 3D coupling g_3^2 \sim g^2 T and Higgs mass m_3 \sim gT matched perturbatively from 4D QCD.00864-M) Further integrating out the electrostatic modes produces magnetostatic QCD (MQCD), a pure 3D Yang-Mills theory with coupling g_M^2 \sim g^4 T, capturing the non-perturbative magnetic sector where the spatial string tension sets the confinement scale.00263-7) These perturbative and effective theory methods enable precise calculations of the equation of state (EOS), quantifying the pressure P, energy density \epsilon, and related thermodynamic quantities as expansions in \alpha_s = g^2/(4\pi). The pressure for a gluon plasma is P = \frac{\pi^2}{45} (N_c^2 - 1) T^4 \left[ 1 - \frac{15}{4} \left( \frac{\alpha_s}{\pi} \right) + \cdots \right], extended up to next-to-next-to-leading order (NNLO), or O(\alpha_s^2) and including O(g^6 \ln(1/g)) terms via EQCD matching and 3D lattice inputs for non-perturbative contributions.129) Such NNLO results, incorporating quark flavors, agree with lattice QCD simulations for T \gtrsim 4 T_c (where T_c is the transition temperature), establishing the scale where weak-coupling approximations become reliable.

Effective Models and Theories

Effective models in quantum chromodynamics (QCD) provide phenomenological approximations to describe the complex dynamics of quark-gluon matter at intermediate energy scales, where full QCD calculations are computationally intensive. These models simplify the strong interaction by incorporating key symmetries like chiral symmetry and confinement through effective Lagrangians, enabling studies of phase transitions and thermodynamic properties without relying on perturbative expansions. They are particularly useful for modeling the chiral restoration and deconfinement transitions in hot and dense QCD matter, often tuned to reproduce lattice QCD results at zero baryon density. The Nambu-Jona-Lasinio (NJL) model captures chiral dynamics through a four-fermion interaction, originally proposed to describe spontaneous chiral symmetry breaking in quark matter. The model's Lagrangian includes a term \mathcal{L}_{\text{int}} = G [(\bar{q} q)^2 + (\bar{q} i \gamma_5 \vec{\tau} q)^2], where q represents quark fields, G is the coupling constant, and \vec{\tau} are Pauli matrices for flavor degrees of freedom; this interaction leads to dynamical mass generation for quarks via a gap equation, mimicking the constituent quark mass in the vacuum. In the context of QCD matter, the NJL model is extended to finite temperature and chemical potential using the Matsubara formalism, allowing computation of the chiral condensate and susceptibility, which signal the chiral phase transition around 150-200 MeV in temperature. The model has been widely applied to predict the equation of state and phase diagram, though it neglects explicit confinement effects. To incorporate quark confinement, the Polyakov-Nambu-Jona-Lasinio (PNJL) model augments the NJL framework with a background field for the Polyakov loop, a gauge-invariant order parameter for the deconfinement transition. The Polyakov loop \Phi = \frac{1}{3} \Tr_c \exp(i \beta A_4), where A_4 is the temporal gluon field and \beta = 1/T, modifies the quark distribution functions in the thermodynamic potential, effectively suppressing quark contributions below the critical temperature. This extension yields a more realistic phase structure with separate chiral and deconfinement transitions, often coinciding near T_c \approx 170 MeV, and has been fitted to lattice data for the pressure and energy density. The PNJL model also predicts a critical endpoint in the QCD phase diagram at finite baryon chemical potential, around \mu_B / T_c \approx 2-3. The quark-meson (QM) model and the linear sigma model (LSM) offer alternative effective descriptions focused on meson-mediated interactions and phase transitions in QCD matter. In the QM model, quarks couple linearly to scalar and pseudoscalar meson fields (sigma and pion), with the meson potential derived from chiral symmetry considerations, leading to a mean-field treatment where the sigma field generates the quark mass. This setup reproduces the chiral phase transition as a crossover at zero chemical potential, with critical temperature T_c \approx 150 MeV, and extends to finite density for exploring the first-order transition line. The LSM, primarily for the light quark sector, uses a chiral-invariant potential V(\sigma, \vec{\pi}) = \lambda (\sigma^2 + \vec{\pi}^2 - v^2)^2 - H \sigma, where explicit breaking is included via the term H \sigma; it effectively models the restoration of chiral symmetry at high temperature through the vanishing of the sigma vev. Both models are instrumental in studying meson properties in medium and have been used to compute susceptibilities matching experimental dilepton spectra. Hydrodynamic models, particularly viscous relativistic hydrodynamics, describe the collective evolution of QCD matter after an initial thermalization stage, treating it as a near-perfect fluid. Based on the conservation of the energy-momentum tensor \partial_\mu T^{\mu\nu} = 0, with viscous corrections from Israel-Stewart theory, these models incorporate shear and bulk viscosities to simulate expansion and flow patterns in heavy-ion collisions. Seminal applications, such as boost-invariant Bjorken flow, predict elliptic flow coefficients v_2 \approx 0.1-0.2 at RHIC energies, capturing the transition from ideal to viscous hydrodynamics as entropy density decreases. While effective for spatiotemporal evolution, these models rely on inputs like the equation of state from microscopic theories and briefly connect to transport coefficients, such as shear viscosity over entropy density \eta/s \approx 0.1-0.2 near T_c.

Other Methods

In addition to lattice simulations, perturbative methods, and effective field theories, several alternative theoretical frameworks provide valuable insights into the non-perturbative dynamics of QCD matter through formal expansions and dualities. These approaches leverage symmetries, limits, or mappings to gravity to probe confinement, chiral symmetry breaking, and transport properties in regimes inaccessible to direct computation. The 1/N_c expansion, introduced by 't Hooft, considers the limit of large number of colors N_c while keeping the 't Hooft coupling g^2 N_c fixed, transforming QCD into a solvable theory dominated by planar diagrams.90474-1) In this limit, quark loops are suppressed, gluons drive the dynamics leading to confinement, and the spectrum consists of an infinite tower of narrow mesons and glueballs with interactions scaling as 1/N_c, making meson scattering perturbative.90216-4) Chiral symmetry breaking persists, with the light pseudoscalar mesons emerging as nearly stable Goldstone bosons, while heavier mesons acquire widths of order 1/N_c.90262-1) This framework has been applied to understand the hadron spectrum and thermodynamics, revealing that the deconfinement transition sharpens as N_c increases.^[53] Holographic duality, inspired by the AdS/CFT correspondence, maps strongly coupled QCD-like theories to weakly coupled gravity in anti-de Sitter (AdS) space, offering a tool to study confinement and real-time dynamics. In AdS/QCD models, such as the hard-wall construction, an infrared cutoff in the extra dimension mimics confinement, reproducing linear Regge trajectories for mesons and the quark-gluon plasma's equation of state. These duals predict universal transport coefficients; notably, the shear viscosity to entropy density ratio η/s equals 1/(4π) in the infinite coupling limit, saturating the conjectured bound and aligning with heavy-ion collision data for near-ideal hydrodynamics. Soft-wall variants further incorporate dilaton profiles to model chiral symmetry breaking and light hadron masses without fine-tuning. Supersymmetric QCD (SQCD) extends QCD by adding supersymmetric partners, enabling exact non-perturbative solutions via dualities that illuminate the infrared behavior of the non-supersymmetric theory.00048-6) In SQCD with N_f flavors, Seiberg duality relates the electric theory (SU(N_c) gauge group) to a magnetic dual (SU(N_f - N_c)) for N_f > N_c + 1, where confinement and chiral symmetry breaking manifest as a smooth IR fixed point or Affleck-Dine-Seiberg superpotential.00484-I) These exact results, derived from holomorphy and anomalies, provide benchmarks for confinement mechanisms and the quark condensate in QCD, particularly in the conformal window where N_f ≈ (9/2) N_c. The functional renormalization group (FRG) method integrates out quantum fluctuations via Wetterich's flow equation for the effective average action Γ_k, yielding non-perturbative evolution of effective potentials from ultraviolet to infrared scales.90762-M) In QCD applications, FRG solves coupled flow equations for quark, gluon, and meson propagators in truncated schemes like the quark-meson model, capturing the chiral phase transition and critical endpoints with fluctuations beyond mean field. This approach reveals tricritical scaling near μ_B ≈ 0 and first-order transitions at high density, consistent with lattice benchmarks, while avoiding sign problems in the complex plane.

Experimental Investigations

Facilities and Collisions

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) in the United States is a dedicated facility for studying QCD matter through symmetric collisions of heavy ions, such as gold nuclei (Au-Au), as well as polarized proton-proton (p-p) collisions.^[54] Operating since 2000 and entering its final run in 2025, RHIC accelerates ions to energies up to \sqrt{s_{NN}} = 200 GeV per nucleon pair for heavy-ion runs and protons to 255 GeV, enabling the creation of hot, dense matter conditions.^[55] A key feature is the Beam Energy Scan (BES) program, which systematically varies collision energies from 3 GeV to 200 GeV to map the QCD phase diagram, including searches for a critical point in the transition to quark-gluon plasma (QGP).^[56] The Large Hadron Collider (LHC) at CERN in Switzerland provides higher-energy collisions for QCD matter investigations, including symmetric lead-lead (Pb-Pb) runs at \sqrt{s_{NN}} = 5.36 TeV and asymmetric proton-lead (p-Pb) collisions at similar energies per nucleon pair, alongside proton-proton (p-p) reference data up to 5.02 TeV in the heavy-ion program context, as of the 2024-2025 runs.^[57] These collisions, conducted since 2010, produce larger and hotter QGP volumes than at RHIC due to the increased center-of-mass energy.^[57] The primary detectors for heavy-ion physics are ALICE, optimized for tracking and particle identification in the forward region; ATLAS, focusing on high-momentum particles; and CMS, emphasizing calorimetry and muon detection, all of which record Pb-Pb and p-Pb events. In heavy-ion collisions at these facilities, the initial geometric overlap of the colliding nuclei forms an almond-shaped region in the transverse plane for non-central impact parameters, due to the offset between the ion centers.^[58] This asymmetry drives the subsequent hydrodynamic expansion of the produced QGP fireball, which evolves from an initial size on the order of the nuclear radius (~7 fm for lead) to a roughly spherical volume with a radius of approximately 10 fm at thermal freeze-out.^[31] The fireball lifetime, from formation to hadronization, is typically around 10 fm/c, allowing collective expansion before particles decouple.^[59] Looking ahead, the Electron-Ion Collider (EIC), under construction at BNL with installation beginning after RHIC operations conclude in 2025, will enable electron-proton and electron-nucleus collisions using polarized beams to probe the three-dimensional spin and momentum structure of quarks and gluons within QCD matter.^[60] Expected to start operations in the early 2030s, the EIC will complement hadron collider studies by providing high-resolution imaging of nucleon substructure without the complexities of strong initial-state interactions.^[61]

Observables and Signatures

Observables in heavy-ion collisions provide signatures of quark-gluon plasma (QGP) formation by probing its collective behavior, energy loss mechanisms, and thermodynamic fluctuations. These measurements are extracted from particle yields, angular distributions, and event-by-event variations in the final-state hadrons, photons, and leptons produced in collisions. Key observables include anisotropic flow patterns, suppression of high-momentum particles, multiplicity fluctuations, and electromagnetic radiation, each sensitive to distinct aspects of the QGP's properties such as its viscosity and initial conditions.^[62] Anisotropic flow harmonics v_n quantify the azimuthal anisotropy in the transverse momentum distribution of produced particles relative to the reaction plane, with the second-order harmonic v_2, known as elliptic flow, being particularly prominent. Elliptic flow arises from the initial spatial eccentricity of the collision overlap region, which drives pressure gradients that convert geometric asymmetry into momentum space anisotropy through hydrodynamic expansion of the QGP. Higher-order harmonics v_3 and beyond reflect more complex initial-state fluctuations and non-linear responses in the medium. Measurements of v_n are influenced by transport coefficients like shear viscosity, which dampen flow at higher orders.^[62]^[63] Jet quenching manifests as the suppression of high-transverse momentum (p_T) hadrons and jets traversing the QGP, quantified by the nuclear modification factor R_{AA}, defined as the ratio of yield in heavy-ion collisions to that in proton-proton collisions scaled by the number of binary nucleon-nucleon interactions. In central collisions, R_{AA} < 1 indicates significant energy loss of partons due to interactions with the dense medium, where gluons and quarks radiate or collide inelastically, losing energy proportional to the medium's density and path length. This suppression is stronger for light hadrons than for heavy quarks, providing a tomographic probe of the QGP's opacity.^[63] Event-by-event fluctuations in conserved quantities, such as net baryon number, serve as probes for the QCD critical endpoint (CEP) by enhancing higher-order cumulants near a phase transition. The kurtosis of net-proton distributions, a proxy for net-baryon number, measures the non-Gaussianity of multiplicity fluctuations, with deviations from Poisson statistics signaling critical correlations that diverge at the CEP. Higher moments \chi_n of the baryon number susceptibility, related to cumulants via \kappa_4 \sigma^2 = \chi_4 / \chi_2, are expected to show sign changes or peaks in the vicinity of the critical region, allowing searches across beam energies. Dilepton and real photon yields offer penetrating probes of the QGP's early thermal evolution, as these electromagnetic particles interact weakly and escape without rescattering. Low-mass dileptons (virtual photons) and direct photons are dominantly produced from quark-antiquark annihilation and bremsstrahlung in the hot medium, with their invariant mass or transverse momentum spectra encoding the emission temperature. The yield's exponential dependence on temperature allows extraction of an effective early QGP temperature, typically around 200-300 MeV, distinguishing thermal radiation from hadronic backgrounds.^[64]

Evidence and Recent Discoveries

Heavy-Ion Collision Results

Heavy-ion collision experiments at facilities like the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) have provided key experimental evidence for the formation and properties of quark-gluon plasma (QGP). In 2005, the STAR and PHENIX collaborations at RHIC reported strong elliptic flow, quantified by the second-order flow harmonic v_2, in Au+Au collisions at \sqrt{s_{NN}} = 200 GeV. These measurements, scaling with the number of constituent quarks, aligned closely with predictions from relativistic ideal hydrodynamics, confirming the nearly frictionless, collective expansion of a hot, dense QGP medium created in the collisions.^[65] Subsequent data from the LHC's ALICE, ATLAS, and CMS experiments in 2010 further characterized the QGP's fluidity. Viscous hydrodynamic modeling of Pb+Pb collisions at \sqrt{s_{NN}} = 2.76 TeV reproduced observed flow anisotropies with a shear viscosity to entropy density ratio \eta/s \approx 1/(4\pi), the minimal value allowed by quantum field theory bounds, indicating the QGP behaves as the most perfect fluid observed in nature.^[66] This low viscosity enabled efficient conversion of initial spatial asymmetries into final-state momentum anisotropies, supporting the interpretation of a strongly coupled QGP phase. Recent 2025 analyses from ALICE have illuminated the pre-equilibrium dynamics preceding full QGP thermalization. These studies reveal a highly non-equilibrium initial state, characterized by glasma fields and classical Yang-Mills evolution, before transitioning to the equilibrated QGP. The temperature evolves from an initial value of approximately 500 MeV, reflecting extreme early densities, to a kinetic freeze-out temperature around 250 MeV, where particle interactions cease and the system hadronizes.^[67]^[68] The RHIC Beam Energy Scan (BES) program, spanning energies from \sqrt{s_{NN}} = 200 GeV down to 3 GeV, aims to map the QCD phase diagram and locate the critical endpoint (CEP) separating first-order and crossover transitions. As of 2025, no definitive CEP signature has been observed in the scanned region, though enhanced net-proton number fluctuations and higher-order cumulants at lower energies (\sqrt{s_{NN}} \lesssim 20 GeV) suggest proximity to a critical region with growing correlations. September 2025 results from BES-II by the STAR collaboration show signs of phase-change turbulence, further constraining the CEP to higher baryon densities beyond current reach.^[69] These results underscore the need for future BES-III data.

Neutron Star Observations

Neutron star observations offer a unique window into the properties of QCD matter under extreme densities and low temperatures, where the equation of state (EOS) governs the structure of these compact objects and potentially reveals phase transitions to deconfined quark phases. Unlike high-temperature probes from heavy-ion collisions, neutron stars probe cold, ultra-dense matter, with core densities exceeding several times nuclear saturation density, where perturbative QCD suggests the possible onset of quark deconfinement or exotic hadronic phases. In 2019, NASA's Neutron Star Interior Composition Explorer (NICER) provided the first precise mass-radius measurement for the millisecond pulsar PSR J0030+0451 through Bayesian modeling of its X-ray pulse profiles. The analysis inferred an equatorial radius of $12.71^{+1.14}_{-1.19} km for a mass of approximately 1.34 M_\odot, with credible intervals at 68% confidence, implying a stiff EOS at densities around 2-3 times nuclear saturation to support the observed compactness. This constraint rules out overly soft EOS models that would predict smaller radii and favors those compatible with a gradual transition to quark matter, though it does not require deconfined phases below about 1.4 M_\odot.^[70] A 2020 analysis integrated the tidal deformability from the LIGO/Virgo gravitational wave event GW170817 with mass measurements of heavy pulsars exceeding 2 M_\odot, such as PSR J0348+0432, to probe the high-density EOS. Employing a Bayesian framework with a parametrized EOS that allows for phase transitions, the study concluded that massive neutron stars (M \gtrsim 1.4 M_\odot) likely harbor sizable quark-matter cores, as purely hadronic EOS fail to simultaneously match the low tidal deformability (indicating softness at merger densities) and the high maximum masses. The posterior distributions indicate a greater than 80% probability for a first-order transition to two-flavor quark matter in the cores of the most massive observed neutron stars.^[71] Recent multimessenger observations from gravitational wave events, including updated analyses of GW170817 and subsequent detections up to 2025, have further refined constraints on hybrid star configurations, where a quark core is surrounded by hadronic matter. Incorporating pulsar timing, X-ray radius measurements, and gravitational wave signals, these studies indicate ongoing constraints on the prevalence of hybrid configurations.^[72] The hypothetical stability of strange quark matter, a candidate for the ground state of QCD matter, is tested through searches for strangelets—small nuggets of strange matter—in cosmic rays, with no detections reported to date. The PAMELA experiment set stringent upper limits on the strangelet flux at < 10^{-7} particles per square meter per second per steradian above 10^{10} eV.^[73]

Cosmological Implications

Quark-gluon plasma (QGP) formed a significant portion of the early universe's matter content during the Big Bang, persisting from approximately 10^{-6} seconds to 10^{-5} seconds after the initial singularity, when temperatures ranged from about 150 MeV to 2 GeV. This phase of deconfined quarks and gluons contributed to the universe's entropy production, which influences the power spectrum of cosmic microwave background (CMB) anisotropies by altering the initial conditions for structure formation. Specifically, the entropy density released during the QCD phase transition helps set the baryon-to-photon ratio, impacting the damping of acoustic oscillations observed in CMB temperature fluctuations. Big Bang nucleosynthesis (BBN) provides stringent constraints on the QCD critical temperature T_c, requiring it to be below approximately 200 MeV to avoid deviations from the standard model's predicted light element abundances, such as helium-4 and deuterium. Observations of primordial abundances from quasar absorption spectra confirm that the QCD transition did not introduce significant non-equilibrium effects or extra entropy that would disrupt BBN yields, aligning with lattice QCD estimates of T_c \approx 155 MeV. These constraints underscore the smooth integration of QCD matter dynamics into the electroweak epoch without altering the expansion rate beyond observational tolerances. Analyses of COBE/FIRAS data have constrained spectral distortions in the CMB arising from energy injections during the QCD transition, particularly those modulated by the plasma's shear viscosity. These distortions, manifesting as deviations from a perfect blackbody spectrum at the \mu- and y-type levels, offer indirect evidence for the viscous dissipation in the expanding QGP, with upper limits on the distortion parameters (|\mu| < 9 \times 10^{-5}, |y| < 1.5 \times 10^{-5}) constraining the transport coefficients derived from holographic models of QCD matter.^[74] Such insights refine our understanding of how the QCD epoch's hydrodynamics influenced the thermal history prior to recombination. QCD matter in the early plasma also serves as a portal for dark matter interactions, where light mediators coupled to strong sector fields facilitate freeze-out or freeze-in mechanisms for dark matter production. These QCD axion-like particles or hidden sector models predict feeble couplings that preserve the standard cosmological timeline while contributing to the relic density, consistent with CMB and large-scale structure data. Observations from future CMB experiments, such as the Simons Observatory, are expected to further test these interactions by searching for spectral imprints in the early universe's thermal bath.

References

[1]
https://arxiv.org/pdf/1108.2939.pdf
[2]
Hot and Dense QCD: The Physics of the Quark-Gluon-Plasma
The basic constituents of QCD are quarks which interact through the exchange of gluons. It is believed that shortly after the creation of the universe in the ...
[3]
How Do Quark-Gluon-Plasma Fireballs Explode into Hadrons?
Dec 20, 2023 · Quark gluon plasma (QGP) is an exciting state of matter that scientists create in a laboratory by colliding two heavy nuclei.
[4]
[PDF] 9. Quantum Chromodynamics - Particle Data Group
May 31, 2024 · Hadrons are color-singlet (i.e. color- neutral) combinations of quarks, anti-quarks, and gluons. Ab-initio predictive methods for QCD ...
[5]
[PDF] 60. Quark Masses - Particle Data Group
May 31, 2024 · The physical val- ues of the quark masses are determined by computing as many physical quantities as there are. 31st May, 2024. Page 3. 3. 60.
[6]
[PDF] 61. Top Quark - Particle Data Group
May 31, 2024 · The top quark is a left-handed quark with a large mass, unique semi-weak decay, and a short lifetime, making it behave as a free quark.
[7]
The discovery of asymptotic freedom and the emergence of QCD
In this Nobel lecture, I shall describe the turn of events that led to the discovery of asymptotic freedom, which in turn led to the formulation of QCD.
[8]
Effective confining potentials for QCD | Phys. Rev. D
Oct 15, 2014 · In this paper we will compare the physical descriptions, their effective potentials, and the mass scales controlling quark confinement obtained from lattice, ...
[9]
Infrared slavery and quark confinement - ScienceDirect.com
We pose the question of whether the so-called infrared slavery mechanism as, e.g., being manifest in non-Abelian gauge theories, necessarily confines quarks ...Missing: hadronization | Show results with:hadronization
[10]
[PDF] QCD Phenomenology at High Energy
Beta Function. ○ Running of of the QCD coupling αS is determined by the β function, which has the expansion β(αS) = −bα. 2. S. (1 + b. ′. αS) + O(α. 4. S. ) b ...
[11]
[PDF] 9. Quantum Chromodynamics - Particle Data Group
Jun 1, 2020 · Hadrons are color-singlet (i.e. color- neutral) combinations of quarks, anti-quarks, and gluons. Ab-initio predictive methods for QCD ...
[12]
[PDF] arXiv:1401.5355v2 [hep-ph] 6 Mar 2014
Mar 6, 2014 · We show that for the lower branch of the quark condensate and values higher than approximately. −(250 MeV)3 the chiral critical end-point in ...
[13]
[PDF] qq condensate for light quarks beyond the chiral limit - arXiv
Dec 6, 2006 · The value of the qq condensate for the u and d quarks now well determined to be −(235 ± 15 MeV)3 by experiment [7] is the result in the chiral ...
[14]
The Glueball; The Fundamental Particle of Non-Perturbative QCD
Aug 7, 1996 · Theoretical ideas related to the existence of glueballs in QCD are reviewed. These include non-perturbative phenomena such as confinement, ...
[15]
Colloquium: Phase diagram of strongly interacting matter
Jul 17, 2009 · Historically this was one of the great triumphs of relativistic field theories. In QCD, quarks and gluons are the elementary degrees of freedom.<|separator|>
[16]
[PDF] 3. Thermodynamics of the QGP
quark-gluon plasma: gluons as massless bosons with degeneracy g g. = 2(spin) x 8(color) = 16 quarks massless fermions with degeneracy g q. = N f x 2(spin) x 3 ...
[17]
[0812.2831] Color superconductivity and dense quark matter - arXiv
Dec 15, 2008 · We review some of the basic concepts of color superconductivity and some of the aspects of this phase of matter which are relevant for ...
[18]
Quark-Gluon Plasma: from accelerator experiments to early Universe
Oct 14, 2015 · In the Big Bang scenario, the early Universe is characterized by the {\it particle era}, i.e. a Universe made of particles. This period connects ...Missing: timeline | Show results with:timeline
[19]
Cosmic separation of phases, density inhomogeneities and ...
We investigate the role of the cosmic QCD phase transition on the relative abundances of light elements in primordial nucleo-synthesis, with particular ...
[20]
[hep-ph/9901308] Relics of the Cosmological QCD Phase Transition
Jan 14, 1999 · Abstract: The abundance and size distribution of quark nuggets (QN), formed a few microseconds after the big bang due to first order QCD ...Missing: light elements
[21]
The cosmological heavy ion collider: Fast thermalization after ...
Oct 10, 2020 · We apply this to the reheating epoch of inflationary cosmology, and the formation of a cosmological quark gluon plasma (QGP). We compute the ...
[22]
Constraining the equation of state in neutron-star cores via the long ...
Feb 3, 2025 · Evidence for quark-matter cores in massive neutron stars. ... How perturbative QCD constrains the equation of state at neutron-star densities.Missing: rho0 stiffening
[23]
THE DENSE MATTER EQUATION OF STATE FROM NEUTRON ...
We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence.
[24]
[PDF] Strange Quark Matter and Compact Stars.
Theoretically the energy per baryon of strange quark matter may be below 930 MeV, which would render such matter more stable than nuclear matter. (surface ...
[25]
[PDF] Conditions for the existence of stable strange quark matter
Abstract. We discuss the possible existence of absolutely stable strange quark matter within three different types of chiral models.
[26]
Could spinning soup start starquakes and explain a magnetar ...
May 8, 2023 · This rapid shift in radius would cause a starquake called a glitch, where the magnetar suddenly increases its spin thanks to conservation of angular momentum.
[27]
Comprehensive analysis of constructing hybrid stars with a ...
This paper investigates hybrid stars using a hadronic model and a Nambu-Jona-Lasino model for quark matter, including color superconductivity, to study phase ...
[28]
Study of Hybrid Stars with Nonstrange Quark Matter Cores
Feb 19, 2025 · It is found that pure nonstrange quark matter cores can exist in hybrid stars, possessing 0.026–0.04 M ⊙.
[29]
The mixed phase quark core in massive hybrid stars
Mar 29, 2025 · We investigate the properties of hybrid star and the mixed phase core to explore the radius ratio of the mixed phase in hybrid star.
[30]
Strongly coupled quark-gluon plasma in heavy ion collisions
Jul 19, 2017 · Important QGP signals suggested in the 1980s were the jet quenching ( Bjorken, 1982 ) and subsequent melting of charmonium and bottonium ...
[31]
[PDF] Advances in Quark Gluon Plasma - arXiv
In the third part, the heavy-ion collisions at ultra-relativistic energies will be proposed as a unique experimental method to study QGP in the laboratory, as.
[32]
The RHIC Beam Energy Scan Phase II: Physics and Upgrades - arXiv
Oct 10, 2018 · Phase 2 of the BES at RHIC is scheduled to start in 2019 and will explore with precision measurements the intermediate-to-high \mu_B region of the QCD phase ...
[33]
RHIC Beam Energy Scan Program: Phase I and II - Inspire HEP
In 2010, the Relativistic Heavy Ion Collider (RHIC) launched a multi-step experimental program to investigate the phase diagram of strongly interacting ...Missing: μ_B √ s_NN
[34]
Fresh, Direct Evidence for Tiny Drops of Quark-Gluon Plasma
Jan 15, 2025 · A new analysis of PHENIX data provides fresh evidence that collisions of small particles with large nuclei create small specks of quark-gluon plasma.
[35]
Possible Formation of QGP-droplets in Proton-Proton Collisions at ...
Aug 28, 2019 · Abstract page for arXiv paper 1908.10566: Possible Formation of QGP-droplets in Proton-Proton Collisions at the CERN Large Hadron Collider.Missing: mini- sources
[36]
[PDF] Status of the CBM experiment at FAIR - EPJ Web of Conferences
The CBM collaboration prepares for first beams with the CBM detector at the FAIR SIS 100 accelerator in 2028, allowing for a detailed investigation of dense ...<|separator|>
[37]
Status of the NICA project at JINR - ResearchGate
fabrication and tests of the magnets for the NICA facility were started in 2016 [7]. of 290 tons. Work is in a good progress.
[38]
[PDF] T µB ? - arXiv
temperature Tpc ≃ 155 MeV at zero baryon density [3]. Our simulation results for Tpc,lat are ambiguous, as expected for a crossover, see below. We ...Missing: T_c | Show results with:T_c
[39]
[PDF] QCD phase diagram in a magnetized medium from the chiral ... - arXiv
Jun 21, 2021 · The result is that a crossover occurs at a pseudocritical temper- ature Tc(µB = 0) ≃ 155 MeV [25]. Unfortunately, LQCD calculations cannot ...
[40]
[PDF] QCD critical point: Recent developments - EPJ Web of Conferences
It is, therefore, notable, that they all point to the ex- istence of the critical point in the same region of the QCD phase diagram: Tc ∼ 100 − 110. MeV and µc ...
[41]
Variations of the crossover and first-order phase transition curve in ...
Aug 4, 2025 · It is generally expected that the critical point of QCD belongs to the 3D Ising universality class [9] .
[42]
Toward the chiral phase transition in the Roberge-Weiss plane
We discuss the interplay between chiral and center sector phase transitions that occur in QCD with an imaginary quark chemical potential.
[43]
Topological feature and phase structure of QCD at complex ...
Nov 12, 2015 · 2. QCD with imaginary chemical potential · μ I , QCD has a special periodicity so-called the RW periodicity [4]. The RW periodicity can effect ...
[44]
[PDF] arXiv:2502.10267v1 [hep-lat] 14 Feb 2025
Feb 14, 2025 · This paper uses lattice simulations to search for the QCD critical point, excluding it at µB < 450 MeV at the 2σ level.Missing: T_c | Show results with:T_c
[45]
[1407.6387] The equation of state in (2+1)-flavor QCD - arXiv
Jul 23, 2014 · Abstract:We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with ...
[46]
“QGP Signatures” revisited | The European Physical Journal C
Mar 8, 2024 · The data are in good agreement with hybrid model calculations that use values \eta /s \sim 0.1-0.2 in the QGP phase. ... The specific shear ...
[47]
[0705.4280] Bulk viscosity of QCD matter near the critical temperature
May 29, 2007 · We find that close to the deconfinement temperature bulk viscosity becomes large, with viscosity-to-entropy ratio zeta/s about 1.Missing: Tc conformal breaking
[48]
Probing deconfinement with Polyakov loop susceptibilities
Jul 9, 2013 · We show that the ratios of these susceptibilities are excellent probes of the deconfinement transition, independent of the renormalization of the Polyakov loop.
[49]
Finite-density lattice QCD and sign problem - ScienceDirect.com
In this article, we summarize the past development and current status of the field of finite-density lattice QCD.Missing: mu | Show results with:mu
[50]
New Way to Determine Lattice QCD Equation of State
Sep 15, 2022 · In this work, we present a cumulant expansion procedure that allows to investigate and regulate these biased estimates at different orders in \mu.
[51]
[2002.10188] Perturbative Thermal QCD: Formalism and Applications
Feb 24, 2020 · This paper discusses perturbation theory for studying QCD matter, covering thermal field theory, electromagnetic rates, transport coefficients, ...Missing: NNLO alpha_s Laine
[52]
[0905.1061] Large-N_c QCD - arXiv
May 7, 2009 · Abstract: A brief review of large-N_c QCD and the 1/N_c expansion is given. Important results for large-N_c mesons and baryons are highlighted.Missing: Nc limit gluons Goldstone
[53]
BNL | Relativistic Heavy Ion Collider (RHIC)
RHIC is the world's only machine capable of colliding beams of polarized protons to investigate the 'missing' spin of the proton.RHIC Physics · Electron-Ion Collider · Accelerator Complex · sPHENIX Detector
[54]
Relativistic Heavy Ion Collider (RHIC) - DOE Office of Science
RHIC collides all ion beam species from protons to Uranium at energies up to 100 GeV (Giga Electron Volts) per nucleon (250 GeV for protons). Two concentric ...
[55]
Beam Energy Scan Theory (BEST) Collaboration Project | Physics
The BEST Collaboration, funded by the US Department of Energy, provides a theoretical framework for interpreting the Beam Energy Scan program results.
[56]
First lead-ion collisions in the LHC at record energy - CERN
Nov 23, 2022 · On Friday, 18 November, a test using collisions of lead ions was carried out in the LHC and provided an opportunity for the experiments to ...
[57]
Heavy-ion run at the LHC begins - CERN
Nov 6, 2024 · A new heavy-ion run began at the LHC, smashing together lead ions, containing 208 nucleons (82 protons and 126 neutrons), at an energy of 5.36 TeV per nucleon ...Missing: sqrt( s_NN) PbPb
[58]
Mapping the hydrodynamic response to the initial geometry in heavy ...
Feb 24, 2012 · For midcentral collisions, elliptic flow is driven by the almond-shaped overlap area; therefore, fluctuations are smaller and the estimate is ...
[59]
Jet measurements in heavy ion physics | Rev. Mod. Phys.
Jun 12, 2018 · This process, called jet quenching, is studied with experimental measurements of high momentum hadrons, two particle correlations, and jet ...
[60]
Electron-Ion Collider (EIC)
- **Electron-Ion Collider (EIC):**
[61]
Electron-Ion Collider Plans Advance — with Global Support
Jul 28, 2025 · ... EIC becomes operational in the mid-2030s. “Participation in the EIC is viewed as a major scientific and strategic opportunity,” said Jana ...
[62]
[1102.3010] Elliptic Flow: A Brief Review - arXiv
Feb 15, 2011 · In this review I will give a brief introduction to the physics of ultra-relativistic heavy-ion collisions and discuss the current status of elliptic flow ...
[63]
[0908.2306] Jet Quenching in Heavy Ion Collisions - arXiv
Aug 17, 2009 · This review article was prepared for the Landolt-Boernstein volume on Relativisitc Heavy Ion Physics.
[64]
Phenomenology of photon and dilepton production in relativistic ...
We discuss the latest theoretical results on direct photon and dilepton production from relativistic heavy-ion collisions. While the dilepton spectra at low ...
[65]
[nucl-ex/0510029] The role of elliptic flow correlations in the ... - arXiv
Oct 10, 2005 · Abstract: Flow measurements are reviewed with particular emphasis on the hydrodynamic character of elliptic flow at RHIC.
[66]
[1108.5323] The viscosity of quark-gluon plasma at RHIC and the LHC
Aug 26, 2011 · The paper extracts the shear viscosity of quark-gluon plasma (QGP) using elliptic flow data, finding the same value at RHIC and LHC energies.
[67]
[2510.04661] Overview of the latest developments in understanding ...
Oct 6, 2025 · We will focus on the latest developments regarding the features and the numerical description of the non-equilibrium pre-QGP matter, as well as ...Missing: ALICE LHC temperature evolution
[68]
[PDF] HASCO Summer School 2024 Klaus Reygers, Heidelberg University
Jul 30, 2024 · Expected initial temperatures of 500 MeV or higher. Expansion and cooling: T ~ τ–1/3. Transition QGP to hadrons at Tpc = 150–160 MeV. Chemical ...
[69]
More Signs of Phase-change 'Turbulence' in Nuclear Matter
Sep 29, 2025 · STAR Collaboration unveils most precise data yet in search for a 'critical point' in nuclear phase transitions.Missing: end | Show results with:end
[70]
A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter ...
Dec 12, 2019 · We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451, conditional on pulse-profile modeling.
[71]
Evidence for quark-matter cores in massive neutron stars - Nature
Jun 1, 2020 · (1). where μCET is the baryon chemical potential corresponding to the density nCET, that is, nCET ≡ n(μCET). This result is then further ...
[72]
New Upper Limit on Strange Quark Matter Abundance in Cosmic ...
Sep 8, 2015 · In this work we present results of a direct search for strange quark matter (SQM) in cosmic rays with the PAMELA space spectrometer.Missing: bounds none