Dynamics of causal hydrodynamic fluctuations in an expanding system

We develop a framework of causal hydrodynamic fluctuations in one-dimensional expanding system performing linearisation of the hydrodynamic equations around the boost invariant solution. Through the description of space-time evolution of thermodynamic variables and flow velocity, we find a novel phenomenon that the structure of thermodynamic variables is almost frozen. We also show that two-particle correlation functions of final hadrons after freezeout are closely related with the mass of hadrons and properties of the medium such as viscosity, relaxation time and equation of state.


Introduction
Nowadays, relativistic dissipative hydrodynamics including shear and/or bulk viscosities has been used to extract transport properties of the QGP from experimental data.According to the fluctuation-dissipation relations (FDR), hydrodynamic fluctuations and dissipations are, however, always accompanied with each other.Since phenomena induced by the hydrodynamic fluctuations include the information of transport coefficients through FDRs, these enable us to analyse transport properties of the QGP from a viewpoint of hydrodynamic fluctuations.Therefore, FDRs require hydrodynamic fluctuations in the dynamical framework of relativistic heavy ion collisions.
The first application of hydrodynamic fluctuations to the phenomenology of relativistic heavy ion collisions was made in Ref. [1] and later extended to causal framework in Ref. [2].They linearised the hydrodynamic equations around the boost invariant solution in one-dimensional expanding system [3] and regarded the linearised equations as equations of motion (EoM) for fluctuations.Solving these equations, they obtained correlation of pion yield fluctuations as a function of rapidity gap.Although the basic idea of our work is almost the same as Ref. [2], we focus more on analysis of properties of the QGP, identified hadron spectra and event-by-event phenomena induced by hydrodynamic fluctuations.
Although our formalism is almost the same as the one obtained in Ref. [2], we regard the linearised equations as stochastic differential equations and solve them numerically on an event-by-event basis.We demonstrate the event-by-event space-time evolution of fluctuations of thermodynamic variables and analyse the two-point correlation functions of them.To see the effects of hydrodynamic fluctuations on experimental observables, we also analyse twoparticle correlation functions of final hadrons.
To perform linearisation, we first assume small deviations of four flow velocity, u µ , from the boost invariant solution [3], u µ = (cosh (η s + δy(τ, η s )) , 0, 0, sinh (η s + δy(τ, η s ))) , where η s ≡ tanh −1 (z/t) and δy are space-time rapidity and deviation of flow rapidity from boost invariant solution, respectively.Correspondingly, all thermodynamic variables in the energy-momentum tensor can be expanded as, e.g., e ≈ e 0 (τ) + δe(τ, η s ).All variables with subscripts 0 denote variables of background and these with δ denote fluctuations.Under these assumptions, we obtain EoMs of energy density fluctuations and flow rapidity fluctuations: So far, we have not assumed any specific forms of constitutive equations for shear pressure π ≡ π 00 − π 33 and bulk pressure Π in the course of derivation.Therefore, we employ the simplest causal constitutive equations [4,5] including noise terms ξ π and ξ Π and linearise these equations following the same prescription as explained above: where transport coefficients η and τ π are shear viscosity and relaxation time for shear pressure, respectively.Note that the specific form of the equation of bulk pressure Π is quite similar to that of shear pressure π.The fluctuations of hydrostatic pressure δp and transport coefficients δη and δτ π can be converted from energy density fluctuations once we assume models of EoS and transport coefficients.We next set the power of noises and its probability distribution.The noise term ξ π obeys the following FDR [6] in the Milne coordinate neglecting possible corrections due to the dynamical evolution of backgrounds [7]: where the delta functions in transverse plane and in η s direction are replaced with 1/∆x∆y and Gaussian function, G(η s − η ′ s ), with standard deviation σ η s , respectively.We employ two models of the EoS: the conformal EoS, p = 1 3 e, with the degrees of freedom, d = 47.5, for massless N f = 3 QCD and a parametrization of lattice EoS results [8].In this study we neglect the bulk pressure Π for simplicity, which is justified when we use the conformal EoS.As for the transport coefficients, we employ the specific shear viscosity η/s = 1/4π [9] and relaxation time for shear pressure τ π = (2 − ln2)/2πT [10].
To see how correlations of fluctuations of thermodynamic variables and flow rapidity are inherited by two-particle correlations, we calculate momentum distribution of hadrons via the Cooper-Frye formula [11] assuming the Boltzmann distribution with viscous correction [12,13].Figure 2 shows normalised two-particle correlations dN dY 1 dN dY 2 / dN dY 1 dN dY 2 as functions of rapidity gap ∆Y ≡ |Y 1 − Y 2 | with various settings, here subscripts 1 and 2 are labels of particle 1 2, respectively.As seen in Fig. 2 (a), the pattern of the correlations is seen more clearly for heavier hadrons.It indicates that the heavier hadrons are better probes of two-particle correlations.In Fig. 2 (b), one sees viscosity enhances the correlations.Note that relaxation time turns out to suppress the correlations (not shown).The magnitude of two-particle correlations is highly sensitive to the transport coefficients.Nevertheless the positions of dips and bumps do not depend on them.We also find that the shape of two-particle correlations is different from each other according to the choice of the EoS.It indicates that two-particle correlations include information of not only transport coefficients but also the EoS.

Summary
We developed a framework of causal hydrodynamic fluctuations in one-dimensional expanding system by means of linearisation.Through the description of space-time evolution of thermodynamic variables, we found a frozen structure of thermodynamic variables.From the analysis of the structure, it could be possible to extract the information of the early stage of hydrodynamic evolution.We found that two-particle correlations as functions of rapidity gap are more enhanced for heavier hadrons and that the magnitude of two-particle correlations is sensitive to the properties of the medium such as viscosity, relaxation time and EoS.These results provide us with an opportunity for multidimensional analysis of the QGP properties created in the relativistic heavy ion collisions.

Figure 1 .
Figure 1.(Left) Space-time evolution of total energy density distribution from one sampled event.(Right) Time evolution of correlation of energy density fluctuations as a function of space-time rapidity gap ∆η s ≡ |η s − η ′ s | at τ = 2, 3 and 4 fm.The total 10,000 events are averaged to obtain the correlation.The conformal EoS is assumed for both results.

Figure 1 (
Figure 1 (left) shows space-time evolution of total energy density distribution from one sampled event.The background energy density decreases rapidly from the initial value due to the rapid expansion of the system.Remarkably, streak-like structure appears through the time evolution and is kept until the final time τ = 11 fm.It means that the pattern of the energy density distribution is almost frozen and could carry the information of the early stage.One possible reason of such a structure formation is that the interplay between the diffusion of fluctuations and the effect of stretching fluctuations due to the rapid expansion of the system.Figure 1 (right) shows time evolution of two-point correlation of normalised energy density fluctuations, δẽ ≡ δe/e 0 .Correlation of normalised energy density fluctuations rapidly grows up around the origin ∆η s ∼ 0 and a dip appears at ∆η s ∼ 0.5.This behavior is plausible from the viewpoint of the conservation law.It is clearly seen that the information propagates in ∆η s direction as time passes.

Figure 1 (
Figure 1 (left) shows space-time evolution of total energy density distribution from one sampled event.The background energy density decreases rapidly from the initial value due to the rapid expansion of the system.Remarkably, streak-like structure appears through the time evolution and is kept until the final time τ = 11 fm.It means that the pattern of the energy density distribution is almost frozen and could carry the information of the early stage.One possible reason of such a structure formation is that the interplay between the diffusion of fluctuations and the effect of stretching fluctuations due to the rapid expansion of the system.Figure 1 (right) shows time evolution of two-point correlation of normalised energy density fluctuations, δẽ ≡ δe/e 0 .Correlation of normalised energy density fluctuations rapidly grows up around the origin ∆η s ∼ 0 and a dip appears at ∆η s ∼ 0.5.This behavior is plausible from the viewpoint of the conservation law.It is clearly seen that the information propagates in ∆η s direction as time passes.

Figure 2 .
Figure 2. Normalised two-particle correlation functions as functions of rapidity gap ∆Y ≡ |Y 1 − Y 2 |.(a) Hadron species dependence (massless particles, charged pions, charged kaons and protons).(b) Shear viscosity dependence of proton two-particle correlation functions (default, 3 times larger η and 5 times larger η).Isochronous freezeout at τ = 10.0 fm with lattice EoS[8] is performed in both results.The total 10,000 events are averaged to obtain the correlation functions.