Kinetic freeze out from an anisotropic fluid in high-energy heavy-ion collisions: particle spectra, Hanbury Brown-Twiss radii, and anisotropic flow

Dissipative relativistic fluid-dynamical descriptions of the extended fireball formed in high-energy heavy-ion collisions are quite successful, yet require a prescription for converting the fluid into particles. We present arguments in favour of using a locally anisotropic momentum distribution for the particles emitted from the fluid, so as to smooth out discontinuities introduced by the usual conversion prescriptions. Building on this ansatz, we investigate the effect of the asymmetry on several observables of heavy ion physics.


I. INTRODUCTION
A large amount of the dynamical properties of the fireball created in high energy collisions of heavy nucleibe it at the Brookhaven Relativistic Heavy Ion Collider (RHIC) or at the CERN Large Hadron Collider (LHC)can be described to a good approximation within the framework of relativistic fluid dynamics (see Ref. [1] for a critical review). The relevant equations of motion have to be supplemented with appropriate initial conditions for the continuous medium and with a recipe for the end of its evolution, namely the conversion of the fluid into particles [2]. Our focus in this work will be on the latter point and on how adopting a new ansatz for the transition can help mitigate a few issues in the usual approaches.
After their emission from the fluid-in which the mean free path is assumed to be very small-, the particles can be left to propagate freely, i.e. they at once acquire an infinitely large mean free path: the conversion step is the so-called (kinetic) freeze out, after which the particle momenta no longer evolve. Alternatively, the emitted particles can be fed into an "afterburner" that implements a set of transport equations for the various particle species and thereby ensures a more gradual change of the mean free path (see Ref. [3] for a recent overview). The fluid-particle transition is then a switch between model descriptions, rather than a physical phenomenon.
Irrespective of the subsequent fate of the particles, their emission from the fluid is often modelled in a similar way. For that reason, we shall generally for the sake of simplicity refer to the conversion process as "freeze out", although the actual decoupling occurs later when the particles are evolved with a transport code.
Most existing studies follow some version of the Cooper-Frye prescription [4]: the fluid breaks up, more or less suddenly, when reaching a three-dimensional hypersurface Σ-sometimes replaced by a thin fourdimensional shell, to mitigate the inherent abruptness of * borghini@physik.uni-bielefeld.de † s.feld@physik.uni-bielefeld.de ‡ chlang@physik.uni-bielefeld.de the description-defined by some a priori criterion, like a constant temperature or energy density. At each point on the freeze-out hypersurface, particles are emitted with a given phase space distribution f (x, p). 1 Integrating over the whole hypersurface, the resulting invariant spectrum of the emitted particles of type i reads In this expression, we have taken into account the fact that the phase space distribution depends on the particle species, especially its bosonic or fermionic nature, and we explicitly factorized out the particle degeneracy factor g i . The phase space occupation factor in the Cooper-Frye formula is chosen so as to ensure the conservation of energy, momentum, and charges-if any-across the freezeout hypersurface. Accordingly, f i (x, p) is usually taken to be the equilibrium thermal distribution-which is appropriate for a perfect fluid-or a near-equilibrium distribution including "correction terms" that match the stress energy tensor of a dissipative fluid. Determining these corrections either from pure theory [5][6][7][8][9][10][11][12][13] or within more phenomenological data-driven approaches [14,15] is an ongoing effort. In any case, it is always implicitly assumed that the decoupling medium is not far from local equilibrium, so that dissipative effects remain small.
Thus, the occupation factors f i at freeze out considered in the literature are, up to small corrections, isotropic in the fluid local rest frame, reflecting the assumption of (near) local thermal equilibrium. Accordingly, f i depends on position only through the corresponding dependence of thermodynamic variables, namely the flow velocity u(x) and its gradients, the freeze-out temperature T f.o. (x), and possibly the chemical potential µ i (x).
In this work, we shall depart from this local isotropy of f i and assume instead a locally asymmetric momentum distribution at decoupling.
Before proceeding any further, let us mention that the existence of some local momentum anisotropy at freeze out was already considered in Ref. [16]. As will become clear in the following section, the anisotropy we are interested in is of a different kind, reflecting the dissimilar underlying motivation. Nevertheless, some of the findings of Ref. [16] naturally translate into similar results in our case.

II. MOTIVATION
The sudden-decoupling scenario embodied in the Cooper-Frye formula (1) aims at gluing together two rather different descriptions. The mismatch of the models is obvious if the fluid freezes out into free-streaming particles, as exemplified by the jump of the Knudsen number from very small to very large values. Even when the Cooper-Frye prescription is used to switch from a dissipative fluid to a collection of interacting hadrons, there remain issues [2,3]. An often mentioned problem is the existence of sectors of the conversion hypersurface Σ where dσ µ (x) dσ µ (x) < 0, which can locally lead to negative contributions to the Cooper-Frye integral. Cures to this issue have been proposed (see e.g. Ref. [17] and references therein), which themselves remain incomplete since they introduce discontinuities across Σ either in the stress energy tensor or in the velocity. These shocks are however artefacts of the modelling, not physical ones.
Another issue of the usual sudden freeze out recipe is the sensitivity of the observables computed with the emitted particles, in particular their spectra, to the parameters in the Cooper-Frye formula. This is in our eyes quite a crucial point: it means that the matching between a "microscopic" approach and a long-wavelength effective theory thereof, namely the kinetic modelling in terms of particles and the fluid-dynamical description, depends significantly on the parameter that separates them, which makes the whole procedure questionable.
A possible way out of the problem is to drop the assumption of a sudden freeze out in favour of a continuous one [18,19]. However, in the current implementations of this approach, the particles decoupling from the fluid do not reinteract with each other afterwards. This again implies for each particle a sudden transition from a vanishingly small to an infinitely large mean free path-where the latter is viewed somewhat abusively as the average length that a given particle is likely to travel in its next step-, which is again unsatisfactory, even though this does not happen at once for the whole fluid.
Despite its deficiencies, the "naive" Cooper-Frye formula remains attractive because of its simplicity, which makes it easier to test novel ideas. In order to ensure a better transition between the fluid and particle description, it seems desirable to "twist" one of the models or both, so as to bring them closer to each other. In this spirit, we suggest that anisotropic hydrodynamics [20][21][22] can improve the smoothness of the transition between the continuous and particle frameworks. As we shall demonstrate in next section, this ansatz helps alleviating the sensitivity to the freeze-out temperature T f.o. : Introducing new control parameters, namely those governing the anisotropy of the phase-space distribution at decoupling, widens the possible range of values for T f.o. . In the remainder of this section, we wish to list a few arguments in favour of distorting the particle distribution at freeze out.
First, in the context of heavy-ion physics there is an obvious analogy with the advocated use of anisotropic hydrodynamics at early stages of the medium evolution, to ease the transition from the locally asymmetric energymomentum tensor of the fields left by the colliding nuclei to the almost isotropic tensor needed to apply usual hydrodynamics consistently. In the early evolution stage, the phase-space distribution is deformed along the axis of the nucleus-nucleus collision (z-axis), while in the case we are interested here we do not expect such a global direction for the anisotropy.
As a matter of fact, our second incentive to resort to a possibly strongly anisotropic freeze-out distribution is the observation of a similar asymmetry, parametrized as two different translation temperatures along the streamlines and perpendicular to them, in hypersonic nonrelativistic flows [23]. These findings help us specify the kind of anisotropy we want to consider hereafter. Let us for simplicity focus on particle emission around midrapidity, so as to discard any anisotropy along the z-direction on symmetry grounds. Far from the fluid, each particle will tend to fly away radially, as implied by the simultaneous conservation of angular momentum and (kinetic) energy. The dispersion of the momentum components transverse to the radial direction will thus be much smaller than that of the radial component.
Eventually, a third argument for assuming a deformed particle distribution is that such an anisotropy was actually found for post-freeze-out distributions arising from the decoupling through time-like portions of freeze-out hypersurfaces [24,25].
Accordingly, we argue that it would be helpful to assume in the Cooper-Frye picture a freeze-out distribution which is already deformed, with a larger mean squared momentum along the radial direction. That is, adopting the Cartesian (out, side, long) system of femtoscopic studies, we assume a larger pressure along the local "out"-direction than in the sidewards and longitudinal directions.
As such, the anisotropy we consider hereafter differs from that considered in Ref. [16], in which the distribution is assumed to be distorted along the z-axis, as a remnant of the distortion along that direction in the initial state of the nucleus-nucleus collision. Both deformations can naturally be present at once, yet our purpose here is to examine the influence of a larger radial-momentum dispersion, so that we keep the pressures in the side-and long-directions equal.
The actual functional form of the phase-space distri-bution at freeze out, as well as the size of the parameters measuring the anisotropy, should emerge from a detailed kinetic description of the decoupling process, which we postpone to a future study [26]. In next section, we shall postulate such a form and examine the change induced by the momentum-space asymmetry on various observables of heavy-ion collisions.

III. EFFECT OF THE LOCAL ANISOTROPY ON OBSERVABLES
Let us assume for the phase-space distribution at decoupling of a particle species with mass m a Romatschke-Strickland-like profile [27], namely f an.
where p ′ out denotes the out component of the particle momentum p ′ with respect to the local rest frame (lrf) of the fluid at position x. Λ, which generalizes the temperature, characterizes the scale over which the particle momentum takes significant values. As hinted at by the notations, both Λ and the anisotropy parameter ξ depend a priori on position and the particle type. Hereafter they will be treated as parameters, and for simplicity taken as constant over the freeze-out hypersurface Σ.
The anisotropy parameter ξ must be larger than −1, to ensure the positivity of the expression under the square root. In order to obtain a larger pressure along the radial direction than perpendicular to it, ξ should be negative.
To test the influence of the momentum anisotropy in Eq. (2), we assume some specific freeze-out flow profile and hypersurface Σ. We thus let the fluid decouple at a constant proper time τ f.o. on an longitudinally infinite, azimuthally symmetric cylinder of radius R. Taking as coordinates in the laboratory frame the proper time τ , space-time rapidity ς, and cylindrical coordinates r, φ, we assume for the fluid velocity on Σ a generalised blast wave-like profile, namely [28,29] for the radial coordinate, u φ = u ς = 0 in the azimuthal and ς directions, and eventually u τ = 1 + (u r ) 2 . With this choice, the phase-space occupation factor (2) reads, when expressed in the laboratory frame Under these assumptions, we can numerically compute the Cooper-Frye integral, from which we can obtain the transverse momentum spectrum, Hanbury-Brown-Twiss (HBT) radii (R out , R side , R long ) [30,31], and the anisotropic flow coefficients v n . We shall focus on pions (m = 140 MeV) produced at midrapidity.
We first present results obtained with fixed values of the "effective temperature" Λ = 150 MeV and of the parameters of the blast wave profile: τ f.o. = 7.5 fm/c, R = 10 fm,ū max = 1, V 2 = V 3 = 0.05, except for HBT radii for which all V n vanish. In contrast, we let the anisotropy parameter ξ vary, giving it values from −0.5 to 0 in steps of 0.1, together with 0.15 and 0.3. According to our argumentation in the previous section, these positive values, which lead to smaller pressure in the radial direction as perpendicular to it, should not be relevant for freeze out; yet we included them for reference sake. Figure 1 shows the resulting transverse momentum distributions. As is to be expected, non-zero values of ξ lead to deviations from the almost exponential shape valid in the isotropic case. More precisely, the spectrum becomes harder when ξ goes to increasingly negative values. This clearly reflects the growing radial pressure-or equivalently effective radial temperature Λ/ √ 1 + ξ-obtained by assuming ξ < 0.
In figure 2, we display the various HBT radii, together with the ratio R out /R side , as functions of the pair transverse momentum K T . To be more precise, the radii R 2 side and R 2 long are the f an. -weighted averages over the freezeout hypersurface of y 2 = r 2 sin 2 φ and z 2 = τ 2 sinh 2 ς, re-  spectively, while R 2 out is the average of x 2 − (K T t/E K ) 2 , where x = r cos φ and t = τ cosh ς.
As was just mentioned, negative values of ξ amount to a larger "radial temperature", and thus to higher thermal velocities in the outwards direction. Since at the same time the emission duration barely changes, this naturally leads to a larger R out , as observed in the upper left panel, as well as to a larger ratio R out /R side (lower right panel) In turn, the longitudinal radius R long shown in the lower left panel is to a large extent unaffected by ξ; this could be anticipated since the longitudinal part of the occupation factor remains unchanged. On the other hand, the behaviour of the sidewards radius R side with varying ξ seen in the upper right panel of figure 2 is more involved, and we did not find a satisfactory explanation describing all its details.
The transverse-momentum dependence of elliptic flow v 2 for various ξ values is shown in figure 3; triangular flow v 3 follows exactly the same trend, so that we do not show it. Thus, anisotropic flow decreases when ξ becomes more negative, that is, as the radial temperature grows. This behaviour reflects the fact that an increase in random thermal motion tends to dilute the effect of directed collective behaviour encoded in the flow velocity and its anisotropies, i.e., it diminishes the v n values, as seen here. Before going any further, let us note that in a more complete approach, the local anisotropy parametrized in this work by ξ should not be uniform, but rather positiondependent. In particular, ξ (or similar parameters) would normally be function of the azimuthal angle φ, paralleling the corresponding dependence of the velocity profile, as we now argue. 2 The fluid-particle conversion, whose modelling ξ is supposed to facilitate, roughly happens when the fluid expansion rate ∇ µ u µ (x) becomes comparable to that of elastic scatterings. Since the flow velocity varies with φ, so does the expansion rate, which motivates an azimuthal dependence of ξ. On the other hand, the scattering rate depends on particle density, obtained by integrating the occupancy factor over momentum, and on the relative velocity of particles. As follows from a straightforward change of integration variable [16], the density is inversely proportional to 1 + ξ(x), thus a priori φ-dependent. In turn, the typical relative velocity is controlled by the (effective) temperature(s) of the decoupling medium, thus function of φ as well. . . All in all, every relevant physical quantity depends on azimuth, so it is non-trivial-and within the scope of this paper rather academic-to determine the actual dependence of ξ. In any case, there will be such a dependence, which will affect the anisotropic flow coefficients v n . The results shown in figures 3 and 6 are thus to be taken with a grain of salt, since they neglect this ingredient.
After having investigated the influence of ξ when all other parameters are fixed, we now want to illustrate the degeneracy introduced by this new parameter, showing that very similar values of the observables can be obtained with different pairs (Λ, ξ). Note that we did not attempt to optimize the results we now report by fine tuning the parameters, as will be made apparent by the values of the latter.
In figure 4, we display the transverse momentum spectra for four sets of values of (Λ, ξ), with Λ varying between 130 and 160 MeV and ξ ranging from −0.5 to 0.3. In all four cases, the values of all other parameters are the same as above, in particularū max = 1. All four curves are barely distinguishable below p T = 1.5 GeV, above which that with (Λ = 130 MeV, ξ = −0.5) starts curving up. The spectrum for (Λ = 140 MeV, ξ = −0.25) only starts differing from those with larger Λ from about 2 GeV onwards, while the remaining two stay very close up to at least 3 GeV. In addition, we show in the same figure the spectrum for (Λ = 130 MeV, ξ = −0.5) and a different flow velocity, namely withū max = 0.8. The change inū max makes the spectrum almost collapse on that for (Λ = 150 MeV, ξ = 0), with at most a 15% relative difference over the whole momentum range.
The HBT radii R out and R side and the elliptic flow v 2 for the same sets of parameters as in figure 4 are respectively shown in figures 5 and 6. As in the case of the transverse spectra, the values of R out or v 2 for all four pairs (Λ, ξ) in the caseū max = 1 are very close to each other, with (Λ = 130 MeV, ξ = −0.5) being most apart  from the other three. We also include the result of the computation withū max = 0.8 which gives a good approximation to the p T -distribution: for v 2 , it basically makes no difference with respect to the caseū max = 1, whereas the departure is more marked for R out . All in all, the results for transverse-momentum distributions, R out , and v 2 support our claim that introducing an extra parameter opens a much wider range for the "freeze-out temperature", here Λ, without affecting drastically the values of the observables.
In contrast, the sidewards HBT radius R side displayed in the bottom panel of figure 5 is much more sensitive to the choice of decoupling parameters (Λ, ξ). This is actually somewhat reassuring, since femtoscopic measurements are precisely designed to probe the space-time configuration at decoupling [32].

IV. DISCUSSION
We have argued that resorting to anisotropic hydrodynamics might smooth out the transition from usual dissipative fluid dynamics to a particle description at the end of the evolution of the fireball created in ultrarelativistic heavy-ion collisions. In particular, this could help diminish the sensitivity of computed observables on the parameters introduced by the decoupling prescription.
As a matter of fact, our findings for transverse spectra, R out , and v 2 (figures 4-6) support the idea that introducing an extra parameter, which governs the local momentum anisotropy at decoupling, opens a much wider range for the "freeze-out temperature", here Λ, without changing significantly the values of the observables. It thus seems possible to find a whole region of parameters to which the "early time" signals like anisotropic flowwhich carry information on the properties of the fireball along its whole evolution [33], rather than on decoupling itself-are to a large extent insensitive. On the other hand, some sensitivity remains for the observables which are governed by the freeze-out process.
In the present exploratory study, we postulated the asymmetric form of the occupation factor at decoupling f an. , and investigated some of the consequences within a toy model. The actual form of f an. , together with that of the associated hydrodynamical quantities, still has to be calculated in a more microscopic approach [26]. This involves at the same time a discussion of the freeze-out hypersurface Σ, whose position in space-time obviously depends on the amount of momentum anisotropy in the phase-space distribution.
Once this is done, it will be necessary to study how the "improved" prescription can be implemented in practice, i.e. how numerical simulations of dissipative fluid dynamics, anisotropic hydrodynamics, and particle transport can be glued together in a satisfactory manner. An important point will be to check what the shortcomings of the sudden freeze-out scenario, in particular the backflow of particles through Σ [34], become in the new approach: if there is more freedom in choosing the decoupling hypersurface, some choices may be more convenient than others. Eventually, it will be interesting to investigate the possible relation of the new prescription, which in essence still assumes a sudden fluid-particle conversion, with continuous emission [18,19]. For instance, one may wonder if it is possible to mimic the latter within the former, or whether one has to formulate a continuous version of the "anisotropic decoupling" scenario.