QCD Critical Point Research Articles

Dynamic density correlations in baryon rich fluid using Mori–Zwanzig–Nakajima projection operator method

In this work, we calculate the dynamic density correlations using Mori–Zwanzig–Nakajima projection operator method. With a judicious choice of slow variables we derive the evolution equations for these slow variables starting from generalised Langevin equation. We get the hydrodynamic form of density correlations function which consist of two acoustic peaks: also called Brillouin peaks, and one thermal peak: also called Rayleigh peak. We then estimate the dynamic density correlations near the critical point using critical exponents extracted from the statistical bootstrap model of hadronic matter. We find that the bulk viscosity contributes to the sound attenuation at leading order sim |t|^{-frac{5}{4}} while the thermal conductivity contribute at sub-leading order sim |t|^{-frac{3}{4}}. On the other hand, only the thermal conductivity contributes at leading order sim |t|^{frac{1}{4}} to the thermal diffusivity. We discuss the implications of these results in search for the QCD critical point in the heavy-ion collision experiments.

Dynamics of non-Gaussian fluctuations in model A

Motivated by the experimental search for the QCD critical point we perform simulations of a stochastic field theory with purely relaxational dynamics (model A). We verify the expected dynamic scaling of correlation functions. Using a finite size scaling analysis we obtain the dynamic critical exponent $z=2.026(56)$. We investigate time dependent correlation functions of higher moments $M^n(t)$ of the order parameter $M(t)$ for $n=1,2,3,4$. We obtain dynamic scaling with the same critical exponent $z$ for all $n$, but the relaxation constant depends on $n$. We also study the relaxation of $M^n(t)$ after a quench, where the simulation is initialized in the high temperature phase, and the dynamics is studied at the critical temperature $T_c$. We find that the evolution does not follow simple scaling with the dynamic exponent $z$, and that it involves an early time rise followed by late stage relaxation.

Uniformizing Lee-Yang singularities

Motivated by the search for the QCD critical point, we discuss how to obtain the singular behavior of a thermodynamic system near a critical point, namely the Lee-Yang singularities, from a limited amount of local data generated in a different region of the phase diagram. We show that by using a limited number of Taylor series coefficients, it is possible to reconstruct the equation of state past the radius of convergence, in particular in the critical region. Furthermore we also show that it is possible to extend this reconstruction to go from a crossover region to the first-order transition region in the phase diagram, using a uniformizing map to pass between Riemann sheets. We illustrate these ideas via the Chiral Random Matrix Model and the Ising Model.

Off-equilibrium non-Gaussian fluctuations near the QCD critical point: an effective field theory perspective

The non-Gaussian fluctuations of baryon density are sensitive to the presence of the conjectured QCD critical point. Their observational consequences are crucial for the ongoing experimental search for this critical point through the beam energy scan program at Relativistic Heavy Ion Collider (RHIC). In the expanding fireball created in a heavy-ion collision, critical fluctuations would inescapably fall out of equilibrium and require a systematic description within a dynamical framework. In this paper, we employ newly developed effective field theory (EFT) for fluctuating hydrodynamics to study the real-time critical non-Gaussian fluctuations of a conserved charge density. In particular, we derive the evolution equations for multi-point correlators of density fluctuations and obtain the closed-form solutions with arbitrary initial conditions that can readily be implemented in realistic simulations for heavy-ion collisions. We find that non-linear interactions among noise fields, which are missing in traditional stochastic hydrodynamics, could potentially contribute to the quartic (fourth-order) fluctuations in the scaling regime even at tree level in off-equilibrium situations.

QCD Critical Point and Net-Proton Number Fluctuations at RHIC-STAR

In the search of QCD phase boundary and critical point, higher-order cumulants of conserved quantities are proposed as promising observables and have been studied extensively both experimentally and theoretically. In this paper we present cumulant ratios up to 6th-order of net-proton number distributions in Au+Au collisions at √SNN = 7.7 - 200 GeV from STAR Beam Energy Scan program phase I and √S = 200 GeV p + p collisions. The results are compared with various models and Lattice QCD calculations.

Transport near the chiral critical point

The evolution of a heavy ion collision passes close to the O(4) critical point of QCD, where fluctuations of the order parameter are expected to be enhanced. Using the appropriate stochastic hydrodynamic equations in mean field near the the pseudo-critical point, we compute how these enhanced fluctuations modify the transport coefficients of QCD. Finally, we estimate the expected critical enhancement of soft pion yields, which provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations.

Critical fluctuations of QCD phase transitions and their related observables

We study how the dilepton production rate and the electric conductivity are affected by the soft modes inherently associated with the second-order phase transitions at the critical temperature of color superconductivity and QCD critical point. It is shown that the soft modes modify the photon self-energy significantly through the so called Aslamasov-Larkin, Maki-Thompson and density of states terms, which are known responsible for the paraconductivity in the metallic superconductivity, and they lead to an anomalous enhancement of the production rate in the low energy/momemtum region and the divergence of conductivity at the critical temperature.

Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories

The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event by event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.

Search for the critical point of strongly interacting matter (Intermittency analysis by NA61/SHINE at CERN SPS)

The existence and location of the QCD critical point is an object of both experimental and theoretical studies. The comprehensive data collected by NA61/SHINE during a two-dimensional scan in beam momentum (13A-150A GeV/c) and system size (p+p, p+Pb , Be+Be, Ar+Sc, Xe+La, Pb+Pb) allows for a systematic search for the critical point – a search for a non-monotonic dependence of various correlation and fluctuation observables on collision energy and size of colliding nuclei. In particular, fluctuations of particle number in transverse momentum space are studied. They are quantified by measuring the scaled factorial moments of multiplicity distribution.

Pileup correction on higher-order cumulants with unfolding approach

Higher-order cumulants of conserved charge distributions are sensitive observables to probe the critical fluctuations near QCD critical point in heavy-ion collisions. Due to high interaction rate, pileup event can be one of the major sources of background in the measurements of higher-order cumulants. In this paper, we studied the effects of pileup events on higher-order cumulants of proton multiplicity distributions using UrQMD model. It is found that the proposed pileup correction fails if the correction parameters are determined by the Glauber fitting of charged particle multiplicities, which is usually done in the real heavy-ion experiment. To address this, we propose a model independent unfolding approach to determine the parameters in the pileup correction. This approach can be applied in the pileup correction for the future measurement of higher-order cumulants in heavy-ion collision experiment.

Baryon transport and the QCD critical point

Fireballs created in relativistic heavy-ion collisions at different beam energies have been argued to follow different trajectories in the QCD phase diagram in which the QCD critical point serves as a landmark. Using a (1+1)-dimensional model setting with transverse homogeneity, we study the complexities introduced by the fact that the evolution history of each fireball cannot be characterized by a single trajectory but rather covers an entire swath of the phase diagram, with the finally emitted hadron spectra integrating over contributions from many different trajectories. Studying the phase diagram trajectories of fluid cells at different space-time rapidities, we explore how baryon diffusion shuffles them around, and how they are affected by critical dynamics near the QCD critical point. We find a striking insensitivity of baryon diffusion to critical effects. Its origins are analyzed and possible implications discussed.

Cumulants of the chiral order parameter in a nonequilibrium Bjorken expansion with a QCD critical point

To understand experimentally obtained net-proton number cumulants in the search for the QCD critical point, we study a dynamical model based on an effective quark-meson Lagrangian with chiral symmetry. We investigate the evolution of the expanding medium created in a heavy-ion collision using a spatially homogeneous fluid and a time-dependent order parameter, the sigma field evolved by a Langevin equation. We extract cumulants of the sigma field along a parametrized freeze-out curve and match the obtained freeze-out points to corresponding beam energies. These cumulants can be related to cumulants of the net-proton number through the sigma-proton coupling to provide a qualitative comparison to experimental data from STAR’s beam energy scan program. We demonstrate that the presence of the spinodal or mixed phase region around the first-order chiral phase transition allows for a wide interval of cumulants at the lowest beam energies.

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world.

Nucleon clustering at kinetic freezeout of heavy-ion collisions via path-integral Monte Carlo

Clustering of the four-nucleon system at kinetic freezeout conditions is studied using path-integral Monte Carlo techniques. This method seeks to improve upon previous calculations which relied on approximate semiclassical methods or few-body quantum mechanics. Estimates are given for the decay probabilities of the 4N system into various light nuclei decay channels and the strength of spatial correlations is characterized. Additionally, a simple model is presented to describe the impact of this clustering on nucleon multiplicity distributions. The effects of a possible modification of the inter-nucleon interaction due to the close critical line (and hypothetical QCD critical point) on the clustering are also studied.

Dynamical spectral structure of density fluctuation near the QCD critical point

The expression for the dynamical spectral structure of the density fluctuation near the QCD critical point has been derived using linear response theory within the purview of Israel–Stewart relativistic viscous hydrodynamics. The change in the spectral structure of the system as it approaches the critical point has been studied. The effects of the critical point have been introduced in the system through a realistic equation of state and the scaling behaviour of various transport coefficients and thermodynamic response functions. We have found that the Rayleigh and Brillouin peaks are distinctly visible when the system is away from the critical point but the peaks tend to merge near the critical point. The sensitivity of the structure of the spectral function on wave vector (k) of the sound wave has been demonstrated. It has been shown that the Brillouin peaks get merged with the Rayleigh peak because of the absorption of sound waves in the vicinity of the critical point.

Intermittency analysis of proton numbers in heavy-ion collisions at energies available at the BNL Relativistic Heavy Ion Collider

Local density fluctuations near the QCD critical point has been suggested to exhibit a power-law behavior which can be probed by an intermittency analysis on scaled factorial moment (SFM) in relativistic heavy-ion collisions. The collision energy and centrality dependence of the second-order SFMs are systematically investigated in Au $+$ Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27, 39, 62.4, and 200 GeV within the UrQMD model. We estimate the noncritical background in the measurement of intermittency and propose a cumulative variable method to effectively remove the background contributions. We further study the effect of particle detection efficiency by implementing the RHIC (STAR) experimental tracking efficiencies in the UrQMD events. A cell-by-cell method is proposed for experimental application of efficiency corrections on SFM. This work can provide a guidance of background subtraction and efficiency correction for the experimental measurement of intermittency in the search of the QCD critical point in heavy-ion collisions.

Critical exponents and transport properties near the QCD critical endpoint from the statistical bootstrap model

We present an estimate of the behavior of the shear and bulk viscosity coefficients when the QCD critical point is approached from the hadronic side, describing hadronic matter within the statistical bootstrap model of strong interactions. The bootstrap model shows critical behavior near the quark-hadron transition temperature if the parameter characterizing the degeneracy of Hagedorn states is properly chosen. We calculate the critical exponents and amplitudes of relevant thermodynamic quantities near the QCD critical point and combine them with an Ansatz for the shear and bulk viscosity coefficients to derive the behavior of these coefficients near the critical point. The shear viscosity to entropy density ratio is found to decrease when the temperature is increased, and to approach the Kovtun–Son–Starinets bound 1/(4pi ) faster near the critical point, while the bulk viscosity coefficient is found to rise very rapidly.

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $O(4)$ critical point of QCD, and thus the fluctuations of the order parameter $(\sigma, \vec{\pi})$ are expected to be enhanced. Our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudo-critical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $O(4)$ critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudo-critical point, and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real time properties of thermal QCD near the pseudo-critical point.

Many-body forces and nucleon clustering near the QCD critical point

It has been proposed that one can look for the QCD critical point (CP) by the Beam Energy Scan (BES) accurately monitoring event-by-event fluctuations. This experimental program is under way at the BNL RHIC collider. Separately, it has been studied how clustering of nucleons at freezeout affects proton multiplicity distribution and light nuclei production. It was found that even a minor increase of the range of nuclear forces dramatically increases clustering, while large correlation length $\xi$ near CP makes attraction due to binary forces unrealistically large. In this paper we show that repulsive many-body forces near CP should overcome the binary ones and effectively suppress clustering. We also discuss current experimental data and point out locations at which a certain drop in clustering may already be observed.

The fate of nonlinear perturbations near the QCD critical point

The impact of the QCD critical point on the propagation of nonlinear waves has been studied. The effects have been investigated within the scope of second-order causal dissipative hydrodynamics by incorporating the critical point into the equation of state, and the scaling behavior of transport coefficients and of thermodynamic response functions. Near the critical point, the nonlinear waves are found to be significantly damped which may result in the disappearance of the Mach cone effects of the away side jet. Such damping may lead to enhancement in the fluctuations of elliptic and higher flow coefficients. Therefore, the disappearance of Mach cone effects and the enhancement of fluctuations in flow harmonics in the event-by-event analysis may be considered as signals of the critical endpoint.