# Baryon number diffusion with critical fluctuations

## Abstract

The description of dynamical fluctuations near the QCD critical point in heavy-ion collisions is crucial for understanding the existing and upcoming experimental data from the beam energy scan programs. In this talk we discuss the evolution of fluctuations of the net-baryon density as given by a stochastic diffusion equation. We study equilibrium as well as dynamical systems for which we can show the impact of nonequilibrium effects on the second-order moment.

###### keywords:

QCD phase diagram, QCD critical point, critical fluctuations, stochastic baryon number diffusion, real-time dynamics^{1}

00 \journalnameNuclear Physics A \runauthM. Nahrgang et al. \jidnupha \jnltitlelogoNuclear Physics A

XXVIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions

(Quark Matter 2017)

## 1 Stochastic diffusion equation

Near the conjectured critical point of QCD the net-baryon density becomes the critical mode, which exhibits diffusive dynamics (1); (2); (3). It is therefore expected that for beam energies for which the created matter reaches temperatures and baryochemical potentials in the critical region, fluctuations in the net-baryon number are enhanced (4); (5); (6); (7). Due to the fast dynamics of the expanding strongly interacting matter nonequilibrium effects become important for a quantitative description of the critical fluctuations (8); (9); (10); (11). Dynamical models of fluctuations coupled to fluid dynamical expansions have so far considered only fluctuations in the chiral order parameter, see NFD (9); (12); (13); (14); (15), which imprint traces in the cumulants of the event-by-event net-proton distributions (15); (16). In this talk we study a numerical implementation of the stochastic diffusion equation for net-baryon number near the critical point

(1) |

Here, the free energy functional

(2) |

for , where is the critical net-baryon density, contains only a Gaussian term, such that Eq. (1) is linear in . Criticality is included in the mass via the temperature dependence of the equilibrium correlation length , where as shown in Fig. 1 (left) is motivated by the 3d Ising model equation of state (17).

The stochastic current describes the internal thermal fluctuations and is given by a Gaussian white noise variable and a variance determined from the fluctuation-dissipation theorem

(3) |

The present system is a simple, linear version of fluid dynamical fluctuations capturing the essential critical phenomena up to second-order moments of the fluctuations. A full implementation of fluid dynamical fluctuations including fluctuations in the energy-momentum tensor and the coupling of Eq. (1) to energy flow is a challenge for theory and numerical implementations (18).

## 2 Static equilibrium fluctuations

We solve Eq. (1) for a one-dimensional system of length with periodic boundary conditions in an implicit scheme considering first a fixed diffusion coefficient fm expressed via and fm. Equilibrium results are independent of the transport coefficient and solely depend on temperature and the lattice spacing . The first quantity to verify our numerical algorithm against is the static structure factor . Numerical results for the discrete structure factor for two temperatures are shown in Fig. 1 (right) for different lattice spacings. Within the statistical noise at small momenta, the theoretically expected perfectly flat behavior is reproduced. The implicit scheme applied here does not introduce any additional dependences on coming from the discretization.

In Fig. 2 we show the spatial correlation function of density fluctuations for (left) and GeV (right). We compare two different system sizes with the same resolution . One observes that the fluctuations are uncorrelated over distances larger than and that the numerical results reproduce the analytical expectation . Net-baryon number conservation in the finite-size system, which is satisfied in the numerics, results in a small but finite shift from this continuum expectation, , which vanishes in the thermodynamic limit. At , we obtain the local variance of the density fluctuations, which is inversely proportional to the lattice spacing. Obviously, is enhanced at the critical temperature compared to high temperatures.

In Fig. 3 (left) this variance is shown as a function of temperature. Again, the increase around is clearly visible and analytical expectations are perfectly reproduced. We also note that higher-order moments, like the kurtosis , are vanishing as expected for purely Gaussian models.

## 3 Dynamical fluctuations

Next, we turn to dynamical fluctuations by assuming a (spatially constant but) time-dependent temperature

(4) |

with , . The system is first equilibrated at high temperature GeV, before the cooling begins at fm/c. The critical temperature is reached at fm/c. In Fig. 3 (right) we compare the results for three different initial diffusion coefficients fm. For the fastest diffusion process the dynamical values of the variance come close to the equilibrium values around . With decreasing diffusion strength nonequilibrium effects become stronger and lead to a decrease of the maximally achieved variance. A slight shift of the maximum to later times is also observed. This retardation effect is, however, rather weak in the present purely Gaussian model.

## 4 Discussion

We have benchmarked a numerical realization of the Gaussian stochastic diffusion equation versus analytical results for static and equilibrated systems. Because of the uncorrelated nature of the fluctuations in the net-baryon density, we find that the local variance depends on the lattice spacing. This unphysical situation will be cured in future work by the introduction of a surface tension term. We have observed an increase of the local variance near the critical temperature. We have additionally verified that higher-order moments, like the skewness and kurtosis, are vanishing, hence that propagated net-baryon number fluctuations are again purely Gaussian. Finally, we investigated systems in which the (spatially constant) temperature changes with time. In this case the dynamics is subject to nonequilibrium effects and becomes relaxational. We observe a reduction of the maximal variance and a shift of this maximum to later times corresponding to temperatures below . Depending on the transport coefficient these nonequilibrium and retardation effects are more or less pronounced.

In ongoing work we include non-Gaussian couplings in the free energy functional motivated by the d Ising universality class and investigate in particular the real-time dynamics of non-Gaussian cumulants (19).

## Acknowledgments

M.N. acknowledges support from the TOGETHER Project Région Pays de la Loire (France). The work of M.B. is funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska Curie grant agreement No 665778 via the National Science Center, Poland, under grant Polonez UMO-2016/21/P/ST2/04035. This work was supported in parts by the U.S. Department of Energy under grants DE-FG02-03ER41260 and DE-FG02-05ER41367. The authors acknowledge fruitful discussions within the Beam Energy Scan Theory (BEST) Topical Collaboration.

### Footnotes

- volume:

### References

- P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49 (1977) 435.
- D. T. Son and M. A. Stephanov, Phys. Rev. D 70 (2004) 056001.
- H. Fujii and M. Ohtani, Phys. Rev. D 70 (2004) 014016.
- M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81 (1998) 4816.
- M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60 (1999) 114028.
- M. A. Stephanov, Phys. Rev. Lett. 102 (2009) 032301.
- M. Asakawa, S. Ejiri and M. Kitazawa, Phys. Rev. Lett. 103 (2009) 262301.
- B. Berdnikov and K. Rajagopal, Phys. Rev. D 61 (2000) 105017.
- M. Nahrgang, S. Leupold, C. Herold and M. Bleicher, Phys. Rev. C 84 (2011) 024912.
- M. Kitazawa, M. Asakawa and H. Ono, Phys. Lett. B 728 (2014) 386.
- S. Mukherjee, R. Venugopalan and Y. Yin, Phys. Rev. C 92 (2015) no.3, 034912.
- M. Nahrgang, S. Leupold and M. Bleicher, Phys. Lett. B 711 (2012) 109.
- M. Nahrgang, C. Herold, S. Leupold, I. Mishustin and M. Bleicher, J. Phys. G 40 (2013) 055108.
- C. Herold, M. Nahrgang, I. Mishustin and M. Bleicher, Phys. Rev. C 87 (2013) no.1, 014907.
- C. Herold, M. Nahrgang, Y. Yan and C. Kobdaj, Phys. Rev. C 93 (2016) no.2, 021902.
- L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 112 (2014) 3, 032302.
- M. Bluhm, M. Nahrgang, S. A. Bass and T. Schaefer, Eur. Phys. J. C 77 (2017) no. 4, 210.
- M. Nahrgang, M. Bluhm, T. Schäfer and S. Bass, arXiv:1704.03553 [nucl-th].
- M. Nahrgang, M. Bluhm, T. Schäfer, S. Bass to be published soon.