Advances in High Energy PhysicsVolume 2013 (2013), Article ID 198928, 23 pageshttp://dx.doi.org/10.1155/2013/198928

Review Article

Femtoscopic and Nonfemtoscopic Two-Particle Correlations in and Collisions at RHIC and LHC Energies

Yu. M. Sinyukov,¹ S. V. Akkelin,¹ Iu. A. Karpenko,1,2 and V. M. Shapoval¹

¹Bogolyubov Institute for Theoretical Physics, Kiev 03680, Ukraine²FIAS, 60438 Frankfurt am Main, Germany

Received 15 April 2013; Accepted 16 June 2013

Academic Editor: Edward Sarkisyan-Grinbaum

Copyright © 2013 Yu. M. Sinyukov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The theoretical review of the last femtoscopy results for the systems created in ultrarelativistic , , and collisions is presented. The basic model, allowing to describe the interferometry data at SPS, RHIC, and LHC, is the hydrokinetic model. The model allows one to avoid the principal problem of the particlization of the medium at nonspace-like sites of transition hypersurfaces and switch to hadronic cascade at a space-like hypersurface with nonequilibrated particle input. The results for pion and kaon interferometry scales in and collisions at LHC and RHIC are presented for different centralities. The new theoretical results as for the femtoscopy of small sources with sizes of 1-2 fm or less are discussed. The uncertainty principle destroys the standard approach of completely chaotic sources: the emitters in such sources cannot radiate independently and incoherently. As a result, the observed femtoscopy scales are reduced, and the Bose-Einstein correlation function is suppressed. The results are applied for the femtoscopy analysis of collisions at  TeV LHC energy and ones at  TeV. The behavior of the corresponding interferometry volumes on multiplicity is compared with what is happening for central collisions. In addition the nonfemtoscopic two-pion correlations in proton-proton collisions at the LHC energies are considered, and a simple model that takes into account correlations induced by the conservation laws and minijets is analyzed.

1. Introduction

The two-particle correlation femtoscopy of identical particles allows one to analyze the space-time structure of a particle emission from the systems created in heavy ion, hadron and lepton collisions (for recent reviews see, e.g., [1–3]). The femtoscopy method, which is based on the Bose-Einstein or Fermi-Dirac interference of identical particles, has been proposed first in [4–7] for measurements of the geometrical sizes and shapes of the interaction region in hadronic collisions. Then it has been developed in [8–14] as a tool for a study of rapidly expanding fireballs formed in ultrarelativistic heavy ion collisions. Despite the extremely small sizes of such systems (the order of the value is around  m), they have a pronounced inhomogeneous structure. The generalized treatment of the interferometry measurements asserts that the measured scales—the interferometry radii—are associated just with the homogeneity lengths in the system [15–18]. Only in the very particular case of a finite homogeneous system such lengths correspond to the total geometrical sizes, but normally they are smaller than the latter. The interferometry scanning of femtosystems at various total momenta of pion pairs allows one to analyze the homogeneity lengths related to different space-time regions of the expanding fireball [11, 12, 17, 18]. An understanding of interferometry in terms of the homogeneity lengths, as opposed to the simple-mind geometrical picture, provides explanation to some, at first sight paradoxical, results at RHIC.

The long-term study of the peculiarities of the interferometry scales behavior in heavy ion collisions, in particular, the so-called RHIC HBT puzzle [19–21], relatively small observed radii and close to unity ratio of the two transverse femtoscales, helps much in clearing up the underlying properties of the matter created in these processes. The physical conditions explaining the RHIC HBT puzzle are [22–28] a relatively hard EoS because of a crossover transition (instead of the 1st order one) between quark-gluon and hadron phases, the presence of prethermal anisotropic transverse flow developed to thermalization time, and an “additional portion” of the transverse flow caused by the shear viscosity effect and fluctuations of the initial conditions. An account for these factors gives the possibility to describe well the hadron spectra together with the femtoscopy data within a realistic freeze-out hydrokinetic model (HKM) with a gradual decay of the fluid into observed particles [28–30].

Soon after the first Large Hadron Collider (LHC) heavy ion results were received, it became evident that the hydrodynamic picture of the collision processes, confirmed at RHIC, is clearly seen also at much higher energy. This conclusion is based on the two classes of observables. The first one is related to the azimuthal anisotropy of particle spectra expressed basically through their second harmonics or coefficients. The obtained LHC results for the transverse momentum dependence of at a given centrality bin (here and below for the quantities depending on one particle momentum means the transverse momentum of one particle, whereas for the pair quantities means the average pair transverse momentum, ) were found to be similar to the ones at RHIC, except for the higher momentum range at LHC [31]. This is the evidence of the same hydrodynamic mechanism of the anisotropy formation as at RHIC. The second type of observables deals with the direct measurements of the space-time structure of nucleus-nucleus collisions by means of the correlation femtoscopy. The hydrodynamic predictions [32] for behavior of the interferometry radii at the LHC energies were confirmed by the ALICE experiment [33]. The most impressive hydrodynamic prediction [28, 34], that the ratio of the two transverse interferometry radii, out to side, will drop in the whole -interval with increasing collision energy and will reach a value close to unity at the LHC, has been discovered experimentally [33].

However, quantitative application of the hydrodynamic approach is a nontrivial problem, because it depends on both the initial conditions for the continuous matter evolution and final state treatments for the particles production. The hydrokinetic model (HKM) [23, 29, 35] allows one to apply hydrodynamics also at the late nonequilibrated stage of gradual system decay, where it can be matched, in its hybrid version (hHKM) [30, 36], with the ultrarelativistic quantum molecular dynamics (UrQMD) cascade [37, 38]. A utilization of the UrQMD ensures an adequate description of the very rarefied stage of matter evolution and transition to particle free-streaming regime. It is especially important at the LHC energies because of relatively prolonged duration of the interacting nonequilibrated stage. It was shown that the hydrokinetic approach without such correction of later-time evolution results in overestimated effective temperature of proton spectra at RHIC energy and insufficient rise of interferometry radii and volume from top RHIC to LHC energies [36]. It is important to note that utilization of hydrokinetic model in between pure hydrodynamics and UrQMD gives the possibility to switch correctly to the UrQMD cascade at any space-like hypersurface, in particular, at isochronic one. It allows one to avoid problems that usually appear in hybrid models matching hydrodynamics with hadronic cascade at hadronisation hypersurface. The latter typically contains nonspace-like sectors that cannot be correctly accounted for in initial conditions for hadronic cascade model.

In the recent paper [39] the correlation femtoscopy analysis is going beyond the model of independent particle emitters, which is fairly good for the systems formed in heavy ion collisions but not for small systems (with sizes about 1 fm) created in , , and collisions. It is found that the uncertainty principle leads to (partial) indistinguishability of closely located emitters that fundamentally impedes their full independence and incoherence. The partial coherence of emitted particles is because of the quantum nature of particle emission and happens even if there is no specific mechanism to produce a coherent component of the source radiation. The found effect leads to reduction of the interferometry radii and suppression of the Bose-Einstein correlation functions. We review briefly the observed results and their application [40] to collisions at  TeV LHC energy.

As for the elementary particle collisions, like , there is no unambiguous interpretation of the HBT radii -dependence. It became clear [41–43] that for relatively small systems the additional two-particle correlations affect the correlation functions in the kinematic region where quantum statistical (QS) and final state interaction (FSI) correlations are usually observed. These correlations can be induced by total energy and momentum conservation laws (see, e.g., [44, 45]) and minijets [39, 46, 47]. As opposed to the QS and FSI correlations, which are familiar from the correlation femtoscopy method and so are sometimes called femtoscopy correlations, these correlations are not directly related to the spatiotemporal scales of the emitter and are therefore called nonfemtoscopic correlations. Since the latter noticeably affects correlation functions for small systems, the interferometry radii extracted from the complete correlation function in collisions depend strongly on the assumption about the so-called correlation baseline—the strength and momentum dependence of the nonfemtoscopic correlations [41–43]. It has an influence on the interpretation of the momentum dependence of the interferometry radii in collisions, where the possibility of hydrodynamic behavior of matter is questionable.

2. Escape Function Dynamics of Expanding Medium Particlization

Let us start with discussion of hydrokinetic approach to collisions and explain how it helps switch from hydrodynamical expansion to molecular dynamics of hadronic particles. It was proposed in [35] to describe the hadronic momentum spectra in collisions basing on the escape function of particles which are gradually liberated from hydrodynamically expanding systems. The escape function, first introduced in hydrodynamic framework in [48, 49] without a resort to the Boltzmann kinetics, was utilized in [35] within the Boltzmann equations in a specific approximation based on hydrodynamic approach. It was shown that such a picture corresponds to a relativistic kinetic equation with the relaxation time approximation for the collision term, where the relaxation time tends to infinity, , when , indicating a gradual transition to the free-streaming regime. It is worth noting that the hydrodynamics at fairly large times play no role in formation of locally anisotropic particle spectra.

The Boltzmann equation for the distribution function in the case of no external forces has the form The terms and are associated with number of particles which, correspondingly, came to point and left this point because of collisions. Term can easily be expressed in terms of the rate of collisions of the particle with momentum , . The term has, in general, more complicated integral structure and depends on differential cross-section.

The escape fraction, , describes the (probabilistic) distribution of the particles that reaches the hypersurface without interactions. At asymptotic hypersurface where times , the distribution function corresponds to free quanta: . In general cases, which we will need, when the space-like hypersurface is situated at the final times, we define the escape function as where is the escape probability for particle to reach freely the hypersurface at some point starting from the point . If the point belongs to this hypersurface , then , and so The escape probability is the relativistic invariant and can be expressed explicitly through the rate of collisions along the world line of free particle with momentum : where . It satisfies the differential equation

It follows from (1) and (5) that

The formal solution of (6) can be presented in the following form: where with corresponds to the portion of the system, which propagates without collisions until some time starting from initial hypersurface , .

The expression (7) for escape function can be explained in simple heuristic way as follows. Let us split the distribution function at each space-time point into two parts: , . The first part, , describes the fraction of the system which will continue to interact before reaching the hypersurface . The second fraction, , describes the particles that reach the hypersurface without interactions. Denote again by the space-time point where the particle at with momentum would be, if it moved freely. Consider, at each vicinity of the phase-space point , the number of particles that have escaped from the interacting system during the time interval . First, this additional portion of escaped particles can be produced only from the interacting part of the system. Second, these particles are only among particles that came to the phase-space vicinity of the point () just after the interaction during the time . Indeed, if some particles from the interacting part of the system do not interact during the given time interval , then they will interact without fail at some future time; thus, they cannot contribute to the additional portion of particles escaping during . Therefore, the additional contribution , from the time interval to the distribution function , is for and is zero for . Here is probability for any given particle at point with momentum to reach without interaction the hypersurface . The summation of such contributions is presented by (7), and the differential of this equation leads to (6).

Utilization of the escape function for the momentum spectra formation is based on (3) which can be used to describe inclusive spectra of particles, that are constructed in the standard way by means of the averages of product of creation and annihilation operators. Namely, on the hypersurface where , . Then, using the Gauss theorem, (6) and (1) analytically continued to off-mass-shell four-momenta , and taking into account that for particles on mass shell, one can get respectively, where is the probability for particle to reach the hypersurface from the space-time point , is a portion of the particles at the initial hypersurface that reach hypersurface moving freely without interactions, is the distribution function at , and is the Lorentz-invariant emission density.

The hydrokinetic approach is based on (10) where one utilizes the escaping function dynamics (6) for finding solution of the Boltzmann equation at hypersurface and calculation there the particle momentum spectra. The escaped functions and escape probabilities are calculated within the local equilibrium approximation for the term and the collision term . The latter is defined from the particle cross-sections which are similar to those in the UrQMD. The equation of state is supposed to be the same as in the Boltzmann hadron gas with changing in time chemical composition. The evolution of the parameters, such as the temperature, chemical potentials, particle concentrations, and collective velocities in the hadron-resonance gas, are defined by the equations of relativistic hydrodynamics. All the details are presented in [23, 29, 35]. The resulting distribution function at , found from the integral form (7) of the equation for escaped functions, is, of course, far from the local equilibrium one.

In the models like UrQMD, distribution function dynamics (1) is utilized for calculation of the momentum spectra in accordance with (11). Both approaches coincide in the case of free-streaming, when and therefore , then the second terms in (10), (11) are equal to zero and the first terms are coincided (note that it is not the case when . Then distribution function dynamics (1) is still as at a free-streaming, while the escaping function dynamics (6) becomes nontrivial). It is worth noting that if the initial hypersurface has nonspace-like sites, , then for some momenta we have , and so these parts give the negative contribution to the particle number density according to (10) and (11). It happens because the particles are going inside the hypersurface near the point . This negative contribution will be compensated when the particle crosses the hypersurface again in another point of the hypersurface to go out. See the possible trajectory of the particle in Figure 1. Of course, the Gauss theorem applied to the enclosed hypersurface guarantees the correct positive result for particle momentum spectra, no matter how complicated structure has and how many times the negative contributions happen. As for the external hypersurface , it can be chosen as a space-like one.

Figure 1: The cartoon of the particlization at the hypersurface (solid line) and spectra formation at (dashed line). The very initial system forms at the hypersurface (double line) and then evolves. The Gauss theorem in (10) and (11) is applied to enclosed hypersurface . In the case of large opacity for particle , also is small for particle , and so , while for the distribution function the similar value is not zero, and for it is negative if is a nonspace-like site of the hypersurface .

Typically at a modeling of collisions the initial conditions for the Boltzmann equations are selected at the hypersurface close to the hadronisation hypersurface. Such a hypersurface is not everywhere a space-like one, but at sites with particles cannot go deeply inside because of high opacity (see Figure 1). It is taken into account in the escape function formalism, where the escape probability for the particles with momentum directed into the fluid is about zero at these sites, , and so the escape function is about zero, , resulting in suppression of the negative contributions to particle momentum density in (10) from the hypersurface. Note that within the escape function formalism a transition from continuous medium to particles is described as a gradual process, because escape function and escape probability can be defined for both sides of the hypersurface (in other words, hadronisation process is treated as gradual in the escape function dynamics, and the very initial system state is defined, in fact, on a space-like hypersurface). The latter fact, accounting for (2), implies the continuity of the distribution function through in accordance with Boltzmann kinetics (1). It means that if it is the local equilibrium on the one side of the hypersurface, it will be the same on the other side, no matter whether they are space-like or nonspace-like sectors. Then in (11), dealing with the full distribution function, the negative contributions at the nonspace-like sites really exist and have to be preserved to maintain true dynamics. Summarizing, the hydrokinetic model is based on continuous behavior of the physical values. It implies a continuous behavior of the locally isotropic distribution function at any parts of the hypersurface , including “bad” nonspace-like sites. It is implemented through a continuous hydrodynamic description of the matter crossing the that is utilized for calculation of the escape probability and emission function. Therefore no discontinuity in the energy and momentum flows through arises. The hydrokinetic approach accounts also for the smooth time behavior of the escape function that forms the distribution function at some space-like hypersurface that can be used then to provide the input of the nonlocally equilibrium distribution function (the output of HKM) to the UrQMD [30].

In cascade models, which solve the Boltzmann equation in the form (11), one does not deal with the small fraction of the particles that are already free at but brings into play the distribution function itself for further numerical simulations of the collisions in the Boltzmann gas. Then one can see from (11) that near the nonspace-like site at a negative number of the particles with corresponding momenta has to be “injected” (see Figure 1). Typically, to cure this dilemma people just cut them in the full distribution function , that leads to discontinuity of the energy and momentum flows through the corresponding sites of [50]. Such a cut, introduced through the Heaviside function in the Cooper-Frye prescription [51], destroys the basic equation (11) where the negative contributions have to be preserved! Also locally, it obviously violates the continuity of the particle current . The error of such prescription is that it considers a decaying hadronic system rather as a star, practically unlimited reservoir of emitted photons/particles, while the hadronic medium formed in the heavy ion collisions is a small compact holder of emitted particles and it is rearranged when the system loses them (back reaction).

The situation can be improved if one introduces some model of the surface layer with the locally isotropic thermal distribution that includes also the particles with momenta moving towards the nonspace-like sites where . These particles have to collide with the particles that move from the inside of to outside. The attempt to substitute such real collisions inside the surface layer by the elastic reflection from the moving (with the “velocity” of the nonspace-like site) wall was performed in [52]. It results in the specific combination of the Heaviside -functions in the Cooper-Frye formula and provides the local conservation law for particles . Such a prescription is utilized in the FAST MC Freeze-out Generator [53] and is planned to be used in UrQMD [54]. However, the method [52] ensures continuity of particle and energy flows but not the momentum flow and so preserves discontinuity of momentum-energy tensor at nonspace-like segments of [52].

Because of the difficulties in building the model of the surface layer, one can use hydrokinetics, that does not deal directly with the distribution function at but only with the distributions of particles continuously escaping from the hydrodynamically expanding matter and calculate a nonequilibrated distribution function with (2) and (4) on a space-like hypersurface where one can switch to UrQMD cascade without troubles described.

Theoretical calculations, presented in the next section, are performed in hybrid hydrokinetic model (hHKM) [30]. At the initial state the Monte-Carlo Glauber model of the initial conditions is used (see [30]). After the thermalization stage of the system's evolution, the matter is supposed to be chemically and thermally equilibrated, and its expansion is described within perfect boost-invariant relativistic hydrodynamics with the lattice QCD-inspired equation of state in the quark-gluon phase [55] matched with chemically equilibrated hadron-resonance gas via crossover-type transition. The hadron-resonance gas consists of 329 well-established hadron states (according to Particle Data Group compilation [56]) made of , , and quarks, including -meson (). With such an equilibrated evolution the system reaches the chemical freeze-out isotherm with the temperature  MeV. At the second stage with , the hydrodynamically expanding hadron system gradually loses its (local) thermal and chemical equilibrium, and particles continuously escape from the system. This stage is described within the hydrokinetic approach [35] to the problem of dynamical decoupling. In hHKM model [30] the hydrokinetic stage is matching with hadron cascade UrQMD one [37, 38] at the isochronic hypersurface (with ) that guarantees the correctness of the matching (see [29, 30] for details). In the latter case of transition from the hydrokinetics to cascade, the following distribution functions are used: The different terms in (12) correspond to the following: is the collision rate of the th sort of hadrons with the rest of particles, is an income of particles into the phase-space point owing to resonance decays, and is the decay rate of a given resonance species. To calculate the collision rates, we assume meson-meson, meson-baryon, and baryon-baryon cross-sections in a way similar to the UrQMD code [37, 38].

3. Femtoscopic Correlations in Relativistic Collisions

In the processes of multiparticle production the two-particle correlation function is defined through the ratio of the one- and two-particle (semi)inclusive spectra as follows: Experimentally, the two-particle correlation function is defined as the ratio of the distribution of particle pairs from the same collision event to the distribution of pairs with particles taken from different events. In heavy ion collisions almost all the correlations between identical pions with low relative momentum are due to quantum statistics and final state interactions. As for the latter, in this review we suppose that they are already extracted from the total correlation function (the method is well known and has been proposed in [57]).

The quantum-statistical enhancement of the pairs of identical pions produced with close momenta was observed first in collisions in 1959 [58]. It took more than a decade to develop the method of pion interferometry microscope based on the discovered phenomenon. This was done at the beginning of 70 s by Kopylov and Podgoretsky [4–6]. Their theoretical analysis assumed the radiating source to consist of independent incoherent emitters. In fact, such a representation is used for a long time for the analysis of the space-time structure of particle sources. For such chaotic sources the four-point average in (8) can be decomposed in the following sum of pair products: Then using (8) and (10) one can express the correlation function (13) for chaotic sources where , , is the emission function; if , then as it follows from (10) (in general case it includes in addition the contribution from ).

When calculating these femtoscopic correlation functions in a quasiclassical particle production model (or event generator) that produces particles without any quantum correlations, the output of event generator is the list of particle positions (at the point of their last interaction) and their momenta, and is equal to unity there. So one has to construct in addition a numerical procedure to calculate quantum-statistical correlations in accordance with (15). It is done usually in a way that is similar to the final state interaction method: one takes outcome of a given classical event generator and then constructs a numerical approximation of (15) based on smoothness conditions [59]. This can be done using the binning technique, also used in several event generators [53, 60] including hHKM. One takes outcome of a given distribution of the particle pairs in the bins according to their relative momentum and the average momentum of the particle pair . If one calculates the correlations arising only due to the Bose-Einstein enhancement, for example,for pion pairs, the corresponding numerical equivalent of (15) at each transverse bin looks like where if and 0 otherwise, with being the bin size in histograms. We decompose the relative momentum into projections and perform analysis in the longitudinal center of mass system (LCMS), where the mean longitudinal momentum of the pair vanishes. Evaluation according to (16) can be done with the help of 3D histograms, implemented in ROOT library classes [61], and in hHKM two separate histograms are used to calculate the numerator and the denominator of (16), which are divided one by another to get the correlation function.

Some remarks to such a receipt are in order here. First, note that this procedure does not change single-particle momentum spectra, while one can expect that they will be changed if the quantum statistics were taken into account in the event generators. Even if hydrodynamic evolution accounts for the quantum statistics through corresponding EoS and gives Bose-Einstein and Fermi-Dirac heat distributions as the input for event generator like UrQMD, the subsequent quasiclassical UrQMD hadronic cascade destroys the true quantum-statistical picture. Unfortunately, an account for the quantum statistics is still not realized for realistic event generators because, in particular, any direct account of quantum statistics is very time consuming (for current developments and recent attempts to overcome this problem see [62, 63] and references therein). Therefore, the theoretical analysis of the quantum effects at multiparticle production, that goes beyond the simple prescription (16), is still very important. Some significant results have been already obtained. One of them concerns quantum corrections to spectra and correlations in the case when homogeneity length in the system is less than the particle wavelength, [15, 16, 64, 65]—this drastically changes the form of the spectra and BE correlations. The other is the multibosonic effects when particle number is close to the Bose-Einstein condensation [65–67]—this leads to decrease of the interferometry radii. Also, the effects of coherence for femtoscopic correlations of charged particles were considered in detail in [68] based on the formalism of generalized coherent states—it gives the tool for the correlation search of a coherence component by measuring the correlations between like and oppositely charged pions. One more quantum effect, which also cannot be taken into account by prescription (16), is connected with the uncertainty principle and is analyzed very recently in [39, 40]. We will discuss this effect in Sections 4 and 5.

The resulting correlation function obtained with (16) is fitted with the Bertsch-Pratt parameterization

Next, we show some results for the correlation radii from hHKM model. Following the experimental cuts (which are somewhat different for STAR, PHENIX, and ALICE Collaborations), we consider pions in central pseudorapidity region . Owing to longitudinal boost invariance and approximate azimuthal symmetry for the most central collisions which we consider for the present HBT studies, the cross-terms , , and are neglected.

The comparison of interferometry radii, calculated in hHKM with the experimental data from  GeV collisions at RHIC, is shown in Figure 2. The parameters of the model are chosen to reproduce the basic set of observables: charged hadron density at midrapidity as a function of collision centrality [30]. Note that PHENIX presented its results for the 0–30% centrality bin, which corresponds to a smaller average multiplicity than the 0–5% STAR bin; thus PHENIX radii lie slightly below the ones calculated by STAR; in our model we observe the same tendency with the initial conditions set for 0–30% centrality.

(a)

(b)

(c)

(d)

Figure 2: HBT radii of pairs at 200 A GeV RHIC energy, calculated in different models and compared to experimental data from the STAR [69] and PHENIX [70] collaborations. Dashed lines correspond to the hydrokinetic procedure of matching (hHKM), while solid lines stand for the “hybrid” model case. The results for the “hybrid-isochronic” model are presented for central 0–5% events with a gray solid line. The ICs used for hydrodynamic evolution:  fm/c with zero initial transverse flow and the MC-Glauber profiles for initial entropy density.

From Figure 2 one can conclude that both the hHKM and the “hybrid” cases describe the data quite satisfactory, except for HBT radii for 30–50% centrality, which are clearly underestimated (and here we note that the hHKM results for 40–50% centrality, which describe well the hadron spectra and flow [30], seem to correspond to smaller effective system size). The shear viscosity should also reduce the ratio because it enhances the transverse flow, and the corresponding influence on HBT radii is shown in [71] (however the model presented in [71] does not include cascade stage; thus only resonance decays are considered after the hydrodynamic freezeout at isothermal hypersurface). As for the “hybrid-isochrone” model, it fails to describe the shape of pion, kaon, and proton transverse spectra, , and long-, side- and out-interferometry radii. The main difference between the first two (hHKM, hybrid) and “hybrid-isochrone” scenarios is that the first two matching procedures do not use the local equilibrium particle distribution functions as input for UrQMD cascade at the space-time regions where the system should be far from equilibrium. The peripheral regions at isochronic hypersurface are spatially and temporally distant from the freeze-out isotherm and have rather small temperatures, and in this transition area the finite and rapidly expanding system cannot be described hydrodynamically: the free-streaming regime of particle propagation already starts there.

Next, recent results from ALICE Collaboration show considerable rise of both and (and the corresponding rise of interferometric volume) with increase of collision energy from 200 A GeV RHIC to 2.76 TeV LHC. As one can see from Figure 3, this behavior is well reproduced in hHKM (see also [36]), and it is found to be caused by the protracted cascade stage at LHC energy. We keep unchanged the model parameters when passing from RHIC to LHC energies, except for a general normalization of initial entropy (or energy) density for increased , contribution from binary collisions, and the baryonic chemical potentials at freeze-out; we also find a decent reproduction of basic set of observables with hHKM at LHC. One can conclude that this supports the same physical picture of bulk matter evolution at both top RHIC and LHC energies.

(a)

(b)

(c)

(d)

Figure 3: HBT radii of pairs for the most central events, calculated in the hHKM model and compared to experimental data from ALICE [33]. Dashed lines correspond to hydrokinetic procedure of matching, while solid lines stand for the “hybrid” model case. Corresponding HBT radii for top RHIC energy are shown for comparison purposes.

In Figure 4 (taken from recent STAR publication [72]) we show HBT radii of kaon pairs measured recently by STAR collaboration, compared to the older results by PHENIX and to the ones calculated in hHKM. The figure demonstrates a good reproduction of kaon femtoscopy in hHKM and Buda-Lund models. This is connected with the fact that in hHKM we see no exact scaling between kaon and pion radii. Generally, due to different cross-sections with the hadron mixture, pions and kaons decouple (i.e., suffer their last interactions) from the system at different times and have different contributions from strong resonance decays. It was demonstrated in the pure HKM [29] that, despite the smaller cross-section, the emission duration of kaons with the same as for pions is slightly larger than for pions, at least for intermediate , because the same corresponds to smaller for kaons (and the duration of emission process is highly dependent), which results in somewhat larger values of HBT radii for kaons. Generally, HBT measurements for different particle species seem to be a valuable tool to study the details of the particle liberation process and discriminate between models of particle production.

(a)

(b)

(c)

(d)

Figure 4: Transverse mass dependence of Gaussian radii (a) , (b) , and (c) for midrapidity kaon pairs from the 30% most central collisions at  GeV. The results from STAR [72] and PHENIX [73] are compared to the calculations in hHKM model (dotted lines).

4. Correlation Femtoscopy for Small Systems

In the work [39] it is shown that for small systems formed in particle collisions (e.g., , ), where the observed interferometry radii are about 1-2 fm or smaller, the uncertainty principle does not allow one to distinguish completely between individual emission points. Also the phases of closely emitted wave packets are mutually coherent. All that is taken into account in the formalism of partially coherent phases in the amplitudes of closely spaced individual emitters. The measure of distinguishability and partial coherence is then the overlap integral of the two emitted wave packets. In thermal systems the role of corresponding coherent length is played by the thermal de Broglie wavelength that defines also the size of a single emitter. The Monte-Carlo method (16) cannot account for such effects since it deals with classical particles and point-like emitters (points of the particle's last collision). The classical probabilities are summarized according to the event generator method (16), while in the quantum approach a superposition of partially coherent amplitudes, associated with different possible emission points, serves as the input for further calculations [39]. Such approach leads to the reduction of the interferometry radii as compared to (16). In addition, the ascription of the factor to the weight of the pion pair in (16) is not correct for very closely located points and , because there is no Bose-Einstein enhancement if the two identical bosons are emitted from the same point [39, 74]. The effect is small for large systems with large number of independent emitters, but for small systems it can be significant, and one has to exclude the excessive contributions (“double counting” [39]) in the two-particle emission amplitude. Such corrections lead to a suppression of the Bose-Einstein correlations that is manifested as a reduction of the observed correlation function intercept as compared to one in the standard method (16).

The results of [39] are presented in the nonrelativistic approximation related to the rest frame of the source moving with four-velocity . In hydrodynamic/hydrokinetic approach the role of such a source at given pair's half-momentum bin near some value is played by the fluid element or piece of the matter with the sizes equal to the homogeneity lengths [15, 16]. These lengths are extracted from the HKM simulations, namely, from the interferometry radii defined by the Gaussian fits to the correlation functions obtained in HKM. All the pairs in procedure (1) are considered in the longitudinally comoving system (LCMS) in the boost-invariant approximation automatically selects the longitudinal rest frame of the source and longitudinal homogeneity length in this frame (it is Lorentz dilated as compared to one in the global system [13]). The femtoscopy analysis is typically related to some bin, and so one needs also to determine the transverse source size in the transverse rest frame, marked by the asterisk. The corresponding Lorentz transformations do not change the side-homogeneity length; as for the out-direction, we proceed in the way proposed first in [13]. The corresponding transformations then are [13, 40] Here , is the emission duration, is a rapidity of the source in transverse direction, and is half sum of transverse rapidities of the particles forming the pair. Note that , and if rapidity of the pair is equal to the rapidity of the source, , then in this particular case the radius in the rest frame is Lorentz-dilated by the factor . Generally, the reference system where the pair's momentum is zero does not coincide with the rest frame of the source that emits the pair. Therefore the direct application of these formulas is not an easy task for rather complicated emission structure in hypothetical hydrodynamic/hydrokinetic model of collisions. Of course, the details of the transformation will be different for the string event generator; therefore one can estimate the theoretical uncertainties using the radii transformation just in the two limiting cases and .

In [40] the quantum corrections are calculated at each -bin in the rest frame of the corresponding source using (18), and then a transition is made again to the global reference system. In what follows the asterisk mark is omitted and all values are related to the source rest frame. To account that due to the uncertainty principle the emitters (strictly speaking emitted wave packets) have the finite sizes ( is the momentum variance of the particle radiation) when defining the lengths of coherence, one should at first consider the amplitude of the radiation processes and only then make statistical averaging over phases of the wave packets using the overlap integral as the coherence measure [39].

Following to [39] we present the quantum state corresponding to a boson with mass emitted at the time from the point as a wave packet with momentum variance which then propagates freely: where is some phase and defines the primary momentum spectrum that we take in the Gaussian form, with the variance . The effective temperature of particle emission in the local rest frames in HKM, , is close to the chemical freeze-out temperature .

The amplitude of the single-particle radiation from some 4-volume is supposed to be a superposition of the wave functions with the Gaussian coefficients , where is the probability distribution in the case of random phases: Then in the rest frame of the source where is the normalization constant.

The single- and two-particle spectra, averaged over the ensemble of emission events with partially correlated phases are

The phase averages are associated with corresponding overlap integrals [39] where are the wave functions of single bosonic states in coordinate representation.

Then the correlation function can be expressed through the homogeneity lengths in the local rest frame , , and which in its turn are expressed through the HBT radii obtained from the Gaussian fit (17) of the HKM correlation functions and transformation law (18) as described above: where , parameter is defined from the model numerically (it is the order of unity for  fm and tends to zero for the large sources—see [39] for details), and the subtracted term responds for the elimination of the double counting.

Now we can see that the apparent interferometry radii extracted from the Gaussian fits to the correlation function (26) are reduced as compared to those obtained in the standard approach.

Particularly, if we neglect the double counting effects, truncate the subtracted term in (26), and fit the correlation function with the Gaussian form (17), we obtain the femtoscopic radii , , and related to the standard ones , , and as follows: where according to the nonrelativistic approximation. At large source sizes, for example, when the homogeneity lengths correspond to -collisions, ,  , all these ratios tend to unity.

The mean emission duration is supposed to be proportional to the average system's size, , that leads to the quadratic equation expressing (and ) through . The latter are connected with ones taken in the global reference system according to transformation laws (18). The value is a free model parameter. Then we put these extracted values into the expression (26) for the correlation function and perform its fitting with the Gaussian (17). This gives us finally the interferometry radii , , and in view of the uncertainty principle. The radii are presented then in the global system using the transformations inverse to (18).

The correlation function is the ratio of the two- and one-particle spectra. It is found [39] that quantum corrections to this ratio are not so sensitive to different forms of the wave packets as the spectra itself. In particular, the effective temperature of the corrected transverse spectra depends on whether the parameter of mean particle momentum is included or not into the wave packet formalism. If yes, the corrected effective temperature for small sources  fm is equal or even higher than that the individual emitters have, , while for the wave packets in the form (19) it is lower [39]. Besides this, in the nonrelativistic approximation one can describe only very soft part of the spectra. That is why we focus in the paper on the corrections to the Bose-Einstein correlation functions where in the rest frame of the source the total and relative momenta of the boson pairs are fairly small.

5. Femtoscopy of Collisions in View of the Uncertainty Principle

This section of the review is based on the results of [40].

The attempt of the systematic theoretical analysis of the pion femtoscopy of collisions at the top RHIC and  TeV LHC energies was made in [75] within Quark-Gluon String Model (QGSM). It was found that for a satisfactory description of the interferometry radii one needs to reduce significantly the formation time by increasing the string tension value relative to the one fixed by the QGSM description of the spectra and multiplicity. Otherwise, the radii obtained within QGSM are too large as compared to the measured femtoscopy scales. Hypothetically one can hope to reduce the predicted radii suggesting the other approach—the hydrodynamic mechanism of the bulk matter production in collisions, at least, for high multiplicity events. Although the applicability of hydrodynamics to collisions is still questionable, there are certain arguments in its favor. For instance, in [76, 77] multiparticle production in nuclear collisions is related to that in hadronic ones within the model based on dissipating energy of participants and their types, which includes Landau relativistic hydrodynamics and constituent quark picture. Then to reproduce high multiplicity, the initially very small system has to be initially superdense. This leads to very large collective velocity gradients, and so the homogeneity lengths should be fairly small. However, as we will demonstrate, even at the maximally possible velocity gradients at the given multiplicity, one gets again an overestimate of the interferometry radii in collisions. Therefore, one can conclude that the problem of theoretical description of the interferometry radii in collisions may be a general one for different types of event generators associated with various particle production mechanisms. Here we try to correct the femtoscopic results of event generators using for this aim the quantum effects accounting for partial indistinguishability and mutual coherence described in the previous section.

We will employ the hydrokinetic model (HKM) [23, 29, 35] in its hybrid form [30] where the UrQMD hadronic cascade is considered as the semiclassical event generator at postfreeze-out (“afterburn”) stage of the hydrodynamic/hydrokinetic evolution. The analysis provided in [30] shows a fairly small difference of the one- and two-particle spectra obtained in hHKM and in the case of the direct matching of hydrodynamics and UrQMD cascade at the chemical freeze-out hypersurface. Thus, in this paper we utilize just the latter simplified “hybrid” variant for the afterburn stage.

Let us try to apply the above hydrokinetic picture to the LHC collisions at  TeV aiming to get the minimal interferometry radii/volume at the given multiplicity bin. As it is known [23] the maximal average velocity gradient and so the minimal homogeneity lengths can be reached at the Gaussian-like initial energy density profile. For the same aim we use the minimal transverse scale in ultrahigh energy collision, close to the size of gluon spots [78, 79] in proton moving with a speed . In detail, the initial boost-invariant tube for collisions has the energy density Gaussian distribution in transverse plane with the width (rms)  fm [78, 79], and, following [30], we attribute it to the initial proper time  fm/c. At this time there is no initial transverse collective flow. The maximal initial energy density is defined by all charged particle multiplicity bin.

In Figures 5 and 6, taken from [40], for the two multiplicity classes and we present the three curves for interferometry radii as a function of : the experimental one, the one taken just from the HKM simulations, and the other one obtained after application of the quantum corrections. The basic parameters used correspond to the limiting case, when the homogeneity lengths in the source rest frame coincide with the ones taken from the model, . For this case ,  GeV/c, and  GeV/c. The parameter values linearly increase with from 1.15 to 1.35 for the case and from 1.02 to 1.1 for . As one can see, the quantum corrected -dependency of the radii gets closer to the experimental values, but for large the corrections are insufficient to fully describe the observable femtoscopy scales behavior. This fact may indicate that sources of particles with large cannot be described in hydrodynamical approximation. Note that just for such large the nontrivial baseline corrections, already provided in presented experimental data, are very essential.

(a)

(b)

(c)

(d)

Figure 5: Interferometry radii dependency on , . The data are taken from [42].

(a)

(b)

(c)

(d)

Figure 6: Interferometry radii dependency on , . The data are taken from [42].

Considering the multiplicity dependence of femtoscopy scales in and collisions we cannot bypass the scaling hypothesis issue [44, 45], that suggests a universal linear dependence of the HBT volume on the particle multiplicity. It means that the observed interferometry volume depends roughly only on the multiplicity of particles produced in collision but not on the geometrical characteristics of the collision process. At the same time, as it was found in the theoretical analysis in [88, 89], the interferometry volume should depend not only on the multiplicity but also on the initial size of colliding systems. In more detail, the intensity of the transverse flow depends on the initial geometrical size of the system: roughly, if the pressure is , then the transverse acceleration . The interferometry radii , that are associated with the homogeneity lengths, depend on the velocity gradient and geometrical size and for nonrelativistic transverse expansion can be approximately expressed through , the averaged transverse velocity , and inverse of the temperature at some final moment [17, 18, 90, 91]:

The result (29) for the HBT radii depends obviously on and, despite its roughness, demonstrates the possible mechanism of compensation of the growing geometrical radii of an expanding fireball in the femtoscopy measurements. For some dynamical models of expanding fireballs the interferometry radii, measured at the final time of system's decoupling, are coincided with the initial geometrical ones, no matter how large the multiplicity is [35, 92]. The reason for such a behavior is clear from (11): if there is no dissipation in the expanding system, , then the spectra and correlation functions are coincided with the initial ones. The detail study of hydrodynamically expanding systems is provided in [88, 89]. It is found that at the boost-invariant isentropic and chemically frozen evolution the interferometry volume, if it were possible to measure at some evolution time , is approximately constant: where is the averaged phase-space density [93, 94] which is found to be approximately conserved during the hydrodynamic evolution under the above conditions as well as [88, 89]. As for the effective temperature of the hadron spectra, , one can see that, when the system's temperature drops, the mean increases; therefore does not change much during the evolution (it slightly decreases with time for pions and increases for protons). Hence , if it has been measured at some evolution time , will also approximately conserve. Of course, the real evolution is neither isentropic nor chemically frozen and includes also QGP stage, but significant dependence of the femtoscopy scales on the initial system size preserves anyway.

Figure 7 shows the dependency for considered case of collisions at the LHC,  TeV, and for the most central (only!) collisions of nuclei having similar sizes, and , at the SPS, RHIC, and LHC. We have also put on the plot our prediction for the interferometry volume of -system, that has initially larger size than -system. As one can see, the different groups of points corresponding to , , and events cannot be fitted by the same straight line. This apparently confirms the result obtained in [88, 89] that the interferometry volume is a function of both variables: the multiplicity and the initial size of colliding system. The latter depends on the atomic number of colliding objects and the collision centrality .

Figure 7: The interferometry volume dependency on charged particles multiplicity. The curve fragments in the middle correspond to the HKM prediction for collision at the LHC energy  GeV. The upper one is related to the initial transverse system size  fm and for the two lower ones  fm. The curves at the left and at the right represent the HKM results for and central collisions respectively, compared to the experimental data at AGS, SPS, RHIC, and LHC, taken from papers [42, 80–87]. The volumes are calculated as a product of respective experimental radii. The blue lines correspond to pure HKM results, whereas the quantum corrections to them are presented by the red lines. Figure is taken from [40].

One can conclude that quantum corrections to the pion interferometry radii in collisions at the LHC can significantly improve the (semiclassical) event generator results that typically give an overestimate of the experimental interferometry radii and volumes. The corrections account for the basic (partial) indistinguishability and mutual coherence of the closely located emitters because of the uncertainty principle [39]. The additional suppression of the Bose-Einstein correlation function also appears. The effects become important for small sources, 1-2 fm or smaller. Such systems cannot be completely random and so require the modification of the standard theoretical approach for the correlation femtoscopy. The predicted femtoscopic scales for collisions need some small corrections only for its minimal values corresponding to the initial transverse size of system 0.9 fm.

6. Nonfemtoscopic Two-Pion Correlations in Small Systems

The interest to the nonfemtoscopic correlations is motivated in particular by the fact that for relatively small systems they appreciably affect the complete two-particle correlation function, forming the so-called correlation baseline. It has an influence on the interpretation of the interferometry radii momentum dependence in collisions with high multiplicities, where the different mechanisms of spectra formation are under discussion. Therefore, for successful and unambiguous application of the correlation femtoscopy method to the case of elementary particle collisions, one needs to know the mechanisms of nonfemtoscopic correlations to separate the femtoscopic and nonfemtoscopic ones.

The two-particle correlation function is defined as where is the probability of observing two particles with momenta and , while and designate the single-particle probabilities. In heavy ion collisions where , , and denotes the femtoscopic correlation function that takes into account quantum statistics and final state interactions. In elementary particle collisions additional (nonfemtoscopic) correlations, like those arising from jet/string fragmentation and from energy and momentum conservation (see, e.g., [41–43]), can also give a significant contribution. The important problem in this regard is whether we can factorize out the part corresponding to the nonfemtoscopic correlations, , from the total correlation function:

It [39], such a factorization, was demonstrated numerically within the simple models of three- and two-particle emission accounting for minijets fragmentation and event-by-event initial state fluctuation induced nonfemtoscopic correlations, and below we will assume this factorization property.

Recently, the ALICE Collaboration utilized some event generators, which do not include effects of quantum statistics, for an estimate of the correlation baseline (i.e., nonfemtoscopic correlation function of identical pions) under the Bose-Einstein peak at LHC energies [41, 42]. It was motivated by a reasonable agreement of the corresponding event generator simulations with the experimental data for correlation functions of oppositely charged pions in collisions at the same energy [41, 42]. The calculated correlation baseline has been utilized by the ALICE Collaboration to extract femtoscopic correlations from measured identical pion two-particle correlation functions [41]. Because the utilized event generators account for energy-momentum conservation and emission of minijets, it was conjectured in [41, 42] that some specific peculiarities of the unlike-sign pion correlations as well as like-sign nonfemtoscopic pion correlations can be caused by the jet-like and energy-momentum conservation induced correlations. In [46], we support this conjecture. There, to describe the nonfemtoscopic pion correlations in a simple analytical model, we assume that pions are produced with momenta in multiparticle production events and consider pions as distinguishable, yet equivalent particles with symmetrical probability density functions. A distinguishability of equivalent particles means that there is no quantum interference between possibilities that correspond to all permutations of the particle momenta . Then single-particle probability, , and the two-particle probability, , can be written as where . The nonsymmetrized -pion probability density in such events reads as where is nonsymmetrized -particle production amplitude, and are 4-momenta of colliding particles (protons), and is the normalization factor, Expression (36) for is rather complicated because, in particular, it depends on particles that are produced in addition to pions. The latter means also that one can hardly expect that total energy or momentum of the pion subsystems are constants in the system’s center of mass; instead, one can expect that they fluctuate in event-by-event basis. Here we assume that the total transverse momentum of pions is equal to zero in the system's center of mass (keeping, however, in mind that this constraint is, in fact, too strong and can be weakened if necessary) and neglect the constrains conditioned by the conservation of energy and longitudinal momentum supposing that the system under consideration is barely -pion subsystem in a small midrapidity region of the total system. Then, motivated by (36), we assume that a nonsymmetrized -pion probability density can be written as where is a nonsymmetrized function of pionic momenta, denotes average constraints on the -pion states that appear due to energy and momentum conservations in multiparticle production events, and we assume that where are transverse components of the momenta of the particles. Then the normalization factor is

Now, to describe the nonfemtoscopic pion correlations in a simple analytical model, we assume that there are no other correlations in the production of -pion states except for the correlations induced by transverse momentum conservation and cluster (minijet) structures in momentum space. For the sake of simplicity we assume here that the only two-particle clusters appear. Then one can write for fairly large where denotes the jet-like correlations between momenta and ; existence of such correlations means that cannot be expressed as a product of one-particle distributions. Then, utilizing the integral representation of the -function by means of the Fourier transformation, , and accounting for (34), (38), (39), and (41), the single-particle probability reads as where A possibility of different cluster configurations of particles means, in particular, that registered particles with momenta and can belong either to different minijets or to the same minijet. Then, taking into account (35), (38), (39), and (41), we get where Here . The first term in the right hand side of (44) is associated with events where the two registered particles belong to the same minijet, and the second term corresponds to events where the particles are from different minijets. Evidently, the former is relatively rare; however, notice that the first term can be significant for small systems with not very large .

It was demonstrated in [46] that this model can reproduce, with reasonable parameters, the correlation functions of unlike-sign pions measured by the ALICE Collaboration [41] and nonfemtoscopic correlations of like-sign pions that are generated in the PHOJET simulations and utilized as the correlation baseline by the ALICE Collaboration [41]. To calculate nonfemtoscopic correlations, analytical parameterizations of the functions in interest were used, namely, where , , and are some parameters, and we assume that . In accordance with ALICE baseline obtained from the PHOJET event generator simulations, we assume that only is measured for each bin. Assuming that longitudinal components of the registered particles are equal to zero, , we approximate as where denotes unregistered angle between and , . Then and, taking into account (44), we get where

It is well known (see, e.g., [44, 45]) that the influence of exact conservation laws on single-particle and two-particle momentum probability densities at the -particle production process depends on the value of and disappears at . Since one considers a subsystem of pions but not the total system, to weaken the influence of the total transverse momentum conservation on pions we will consider and with instead of and in (53). This is the simplest way to account for a weakened conservation law in our model. At the same time, the factor in (53) remains the same since it is associated with the combinatorics of the distribution of particles between clusters in momentum space (“minijets”), which happens whether or not one weakens the total momentum conservation law. Also, for more exact fitting of the data points in each average transverse momentum bin, we utilize the auxiliary factors ; when we compared results of our calculations with ALICE two-pion correlation and simulation data, these proportionality factors differ slightly from unit in our calculations (nearly ). Then (53) gets the form

The results of our calculations of the nonfemtoscopic correlation functions are shown in Figures 8 and 9 in comparison with correlation functions reported by the ALICE Collaboration [41] for different transverse momenta of pion pairs (actually, we performed calculations for the mean value in each bin). The data for unlike-sign pion correlations measured by the ALICE Collaboration as well as for the PHOJET simulations of like-sign two-pion nonfemtoscopic correlation functions at midrapidity for the total charged multiplicity bin in collisions at  GeV are taken from [41, 95]. Note that correlations of nonidentical pions measured by the ALICE Collaboration, as well as the PHOJET simulations of identical two-pion correlation functions, demonstrate Coulomb FSI correlations at the lowest bin and peaks coming from resonance decays. These Coulomb FSI and contributions from resonance production are not taken into account and so are not reproduced in our model. The presented results are obtained for ,  GeV (to minimize the number of fit parameters, we fixed for all calculations. Note that with these parameter values the mean transverse momentum is about  GeV), and the fitted values of are different for like-sign and unlike-sign pion pairs, namely, for the former and for the latter. The relatively high value of can be interpreted as a residual effect on the pion subsystem of total energy-momentum conservation in a multiparticle production process. The relation between fitted values means that the magnitude of the correlations induced by a minijet for unlike-sign pion pairs is higher than that for like-sign ones. This happens because in the former there is no local charge conservation constraint for the production of oppositely charged pion pairs, and, therefore, one can expect less identically charged pion pairs from the fragmenting minijets than oppositely charged ones. One can see from the figures that the behavior of the nonfemtoscopic correlation functions of pions, , is reproduced well despite the simplicity of our model. This is a result of the competition of the two trends: an increase of the correlation function with because of momentum conservation and its decrease due to fragmentation of one minijet into the registered pion pair. Figure 9 also demonstrates the relative contribution of the first and second terms in (56) to the nonfemtoscopic correlation functions.

(a)

(b)

Figure 8: The nonfemtoscopic correlation functions of like-sign (a) and unlike-sign (b) pions in the  GeV bin from a simulation using PHOJET [41, 95] (solid dots) and those calculated from the analytical model: minijets + momentum conservation (solid line). See the text for details.

(a)

(b)

Figure 9: The nonfemtoscopic correlation functions of like-sign pions (a) from a simulation using PHOJET and unlike-sign pions (b) measured by the ALICE Collaboration in the  GeV bin from [41, 95] (solid dots) and those calculated from the analytical model (solid line). The contributions to the nonfemtoscopic correlation function from the first term of (56) (dotted line) and from the second one (dashed line) are also presented.

Another type of correlations not induced by the QS effects is the correlations connected with existence of subensembles of events with different emission functions that leads to the corresponding fluctuations in single-particle and two-particle momentum spectra. In hydrodynamical models of nucleus-nucleus and proton-proton collisions these fluctuations can be caused by asymmetrically fluctuating initial densities used for the hydrostage of the model. Let us consider the effect of such correlations on the resulting correlation function on the example of a simple analytical model [46]. Disregarding the QS correlations (for model of two-particle emission that accounts for the QS see [39]) we suppose that the -particle probability density is defined as where is the -particle probability density for some type of the initial conditions, and denotes the distribution over initial conditions, . To analyze the possible effect of fluctuating initial conditions, here we neglect conservation law constraints and the production of minijets. Because we assume uncorrelated particle emissions for each specific initial condition, one can write where we normalize as , and then ; see (38) and (40). The two-particle nonfemtoscopic correlation function then reads as

Evidently, the different type of fluctuation, that is, the form of the distribution , leads to a different behavior of the nonfemtoscopic correlations. To illustrate that fluctuations can lead to the nonfemtoscopic correlation functions that are similar to the ones induced by minijets, let us consider the toy model where and normalization is chosen in such a way that and . The main feature of such a model is that event-by-event single-particle transverse momentum spectra have a maximum for event-by-event fluctuating values. Such momentum spectrum fluctuations could take place, for example, in hydrodynamics with a highly inhomogeneous initial energy density profile without cylindrical or elliptic symmetry. One can easily see that in this case decreases with , and this means (after taking into account (51) and (52)) that decreases with too, which is similar to the behavior of if the nonfemtoscopic correlations are induced by minijets. At the same time, unlike the latter, the hydrodynamical fluctuations lead to similar correlations for like-sign and unlike-sign pion pairs. Then, our analysis suggests that, up to different resonance yields, the value of the slope of the correlation baseline at relatively low can be somewhere between pure hydrodynamic (i.e., the same as for nonidentical pion pairs) and pure minijet (i.e., lower than for nonidentical pion pairs) scenarios.

7. Summary

It is often said that the femtoscopic measurements allow one to restore the space-time structure of the particle/nucleus collision processes. It is really possible but only through the theoretical constructions allowing to interpret these processes and describe the wide range of the corresponding experimental data. If such an advanced model is based on some kind of the space-time picture, then the crucial point for the model is a recreation of the femtoscopic data. If it is successful, then, first, this model claims to be true, and, second, the same is related to recover by the model the space-time structure of a collision process.

In this review we analyzed hydrokinetic model (HKM) for collisions as the reliable theoretical construction for description of the momentum spectra of pions, kaons, protons/antiprotons, all charged particles, and elliptic flows, in wide range of different centralities [29, 30] at RHIC and LHC. We demonstrate here that it describes well the available femtoscopic data for and collisions.

A particular interest is the theoretical basis of HKM which is grounded on the escape functions (not on the distribution ones) when the particle liberation process from the expanding medium is described within the Boltzmann equation. This allows one to avoid the principal problem of the particlization of a medium at nonspace-like sites of transition hypersurfaces, where medium evolution converts into hadronic cascade, like the UrQMD one. So, using the hydrokinetics, as the “buffer” zone, one can switch to hadronic cascade at some space-like hypersurface with a nonequilibrated particle input from HKM. In review we describe just this hybrid variant of the HKM + UrQMD (hHKM).

The HKM, UrQMD, and other existing models/event generators are quasiclassical constructions, and quantum elements, like Bose-Einstein correlations, are inputted typically at the very late stage as the external weights of the boson pair states. In its turn such a procedure is based on the model of independent chaotic sources [4–6]. It is found recently [39] that such a model is inadequate for fairly small sources where the uncertainty principle leads to (partial) indistinguishability of closely located emitters that fundamentally impedes their full independence and incoherence. We review the recent results as for the correlation femtoscopy method that is going beyond the standard approach of independent/random particle emission. The partial coherence of emitted particles is because of the quantum nature of particle emission and happens even if there is no specific mechanism to produce a coherent component of the source radiation. The measure of distinguishability/indistinguishability and mutual coherence of the two emitted wave packets is associated with their overlap integral. In thermal systems the role of corresponding coherent length is played by the thermal de Broglie wavelength.

The application of the new femtoscopy method is demonstrated using the results of [40], where it was found that quantum corrections to the pion interferometry radii in collisions at the LHC can significantly improve the (semiclassical) event generator results that typically give an overestimate of the experimental interferometry radii and volumes. The effects become important for small sources, 1-2 fm or smaller. Such systems cannot be completely random and so require the modification of the standard theoretical approach for the correlation femtoscopy.

More sophisticated result of this study is a good applicability of the hydrodynamics/hydrokinetics with the quantum corrections for description of the femtoscopy scales not only in collisions but also, at least for large multiplicities, in events. These scales are well reproduced for not large . Whether it means the validity of the hydrodynamic mechanism for the bulk matter production in the LHC collisions is still an open question. It is also related to the problem of early thermalization in the processes of heavy ion collisions; the nature of such phenomenon is still a fundamental theoretical issue.

For small systems the specific nonfemtoscopic correlations can appear in the same kinematic region as the femtoscopic ones. They could be a result of the cluster (minijet) structures in final momentum space of produced particles and the global energy-momentum conservation constraints. The latter typically result in an increase with for fairly high of the nonfemtoscopic two-pion correlation functions of small systems, whereas the former mostly determines a decrease of the ones at relatively low . We discuss the simple analytical model [46] that takes into account correlations induced by the total transverse momentum conservation as well as minijets and show that the model gives reasonable description of the two-pion nonfemtoscopic correlations of identical and nonidentical pions in proton-proton collision events at  GeV reported by the ALICE Collaboration [41]. The important issue is that femtoscopic and nonfemtoscopic correlations are factorized [39].

There can be different types of multiparticle production mechanisms in collisions, and some of them could result in qualitatively similar nonfemtoscopic correlation functions. We discuss heuristic arguments [46] that the two-pion nonfemtoscopic correlation functions calculated in hydrodynamics with event-by-event fluctuating initial conditions can be qualitatively similar at relatively low to the ones calculated in the PHOJET-like generators, where the nonfemtoscopic correlations for low are mainly caused by minijets. It is worth noting an important difference between the nonfemtoscopic correlations induced by minijets and hydrodynamical fluctuations: while the former leads to a higher magnitude of the nonfemtoscopic correlations for unlike-sign pion pairs as compared to like-sign pions, the latter results in a similar (up to the resonance contributions) strength of the nonfemtoscopic correlations for identical and nonidentical pions. Then, if the applicability of hydrodynamics to collisions will be justified, such an analysis allows one to estimate the correlation baseline and, so, to extract the femtoscopic scales in these collisions by means of tuning the hydrokinetic model to reproduce the experimental unlike-sign pion correlations.

Acknowledgments

The research was carried out within the scope of the EUREA: European Ultra Relativistic Energies Agreement (European Research Group: “Heavy ions at ultrarelativistic energies”) and is supported by the National Academy of Sciences of Ukraine (Agreement F4-2013) and by the State fund for Fundamental Researches of Ukraine (Agreement F33/24-2013). The work was supported by the Program of Fundamental Research of the Department of Physics and Astronomy of NAS Ukraine.

References

M. A. Lisa, S. Pratt, R. Soltz, and U. Wiedemann, “Femtoscopy in relativistic heavy ion collisions: two decades of progress,” Annual Review of Nuclear and Particle Science, vol. 55, pp. 357–402, 2005.

M. Lisa and S. Pratt, “Femtoscopically probing the freeze-out configuration in heavy ion collisions,” in The Final Expansion State of Relativistic Nuclear Collisions, R. Stock, Ed., vol. 23 of Relativistic Heavy Ion Physics, pp. 653–685, 2010.

Z. Chajȩcki, “Femtoscopy in hadron and lepton collisions: world systematics,” Acta Physica Polonica B, vol. 40, no. 4, pp. 1119–1136, 2009.

G. I. Kopylov and M. I. Podgoretsky, “Correlations of identical particles emitted by highly excited nuclei,” Soviet Journal of Nuclear Physics, vol. 15, p. 219, 1972.

G. I. Kopylov and M. I. Podgoretsky, “Multiple production and interference of particles emitted by moving sources,” Soviet Journal of Nuclear Physics, vol. 18, p. 336, 1974.

G. I. Kopylov and M. I. Podgoretsky, “Mutual coherence functions of elementary particles and multiple production,” Soviet Journal of Nuclear Physics, vol. 19, p. 215, 1974.

G. Cocconi, “Second-order interference as a tool for the determination of hadron fireball dimensions,” Physics Letters B, vol. 49, no. 5, pp. 459–461, 1974.

S. Pratt, “Pion interferometry for exploding sources,” Physical Review Letters, vol. 53, no. 13, pp. 1219–1221, 1984.

S. Pratt, “Pion interferometry of quark-gluon plasma,” Physical Review D, vol. 33, no. 5, pp. 1314–1327, 1986.

A. N. Makhlin and Yu. M. Sinyukov, “The pion interferometry theory for the hydrodynamic stage of multiple processes,” Soviet Journal of Nuclear Physics, vol. 46, p. 345, 1987.

V. A. Averchenkov, A. N. Makhlin, and Yu. M. Sinyukov, “Study of collective motion in hadronic matter by a pion interferometry method,” Soviet Journal of Nuclear Physics, vol. 46, p. 905, 1987.

A. N. Makhlin and Yu. M. Sinyukov, “The hydrodynamics of hadron matter under a pion interferometric microscope,” Zeitschrift für Physik C, vol. 39, no. 1, pp. 69–73, 1988.

Yu. M. Sinyukov, “Hanbury Brown/Twiss correlations for expanding hadron and quark-gluon matter,” Nuclear Physics A, vol. 498, no. C, pp. 151–159, 1989.

Y. Hama and S. S. Padula, “Bose-Einstein correlation of particles produced by expanding sources,” Physical Review D, vol. 37, no. 11, pp. 3237–3245, 1988.

Yu. M. Sinyukov, “Spectra and correlations in locally equilibrium hadron and quark-gluon systems,” Nuclear Physics A, vol. 566, no. C, pp. 589–592, 1994.

Yu. M. Sinyukov, “Boson spectra and correlations in small thermalized systems,” in Hot Hadronic Matter: Theory and Experiment, J. Letessier, H. H. Gutbrod, and J. Rafelski, Eds., pp. 309–322, Plenum, New York, NY, USA, 1995.

S. V. Akkelin and Yu. M. Sinyukov, “The HBT-interferometry of expanding sources,” Physics Letters B, vol. 356, no. 4, pp. 525–530, 1995.

S. V. Akkelin and Yu. M. Sinyukov, “The HBT-interferometry of expanding inhomogeneous sources,” Zeitschrift fur Physik C, vol. 72, no. 3, pp. 501–507, 1996.

U. Heinz, “The quark-gluon plasma at RHIC,” Nuclear Physics A, vol. 721, pp. C30–C39, 2003.

S. Pratt, “Correlations and fluctuations, a summary of Quark Matter 2002,” Nuclear Physics A, vol. 715, pp. 389c–398c, 2003.

S. Soff, S. A. Bass, D. H. Hardtke, and S. Y. Panitkin, “Particle correlations at RHIC—scrutiny of a puzzle,” Nuclear Physics A, vol. 715, pp. 801c–804c, 2003.

Yu. M. Sinyukov, “Matter evolution and soft physics in collisions,” Acta Physica Polonica B, vol. 37, no. 12, pp. 3343–3369, 2006.

S. V. Akkelin, Y. Hama, Iu. A. Karpenko, and Yu. M. Sinyukov, “Hydro-kinetic approach to relativistic heavy ion collisions,” Physical Review C, vol. 78, no. 3, Article ID 034906, 2008.

Yu. M. Sinyukov, A. V. Nazarenko, and Iu. A. Karpenko, “Is early thermalization really needed in collisions?” Acta Physica Polonica B, vol. 40, no. 4, pp. 1109–1117, 2009.

M. Gyulassy, Iu. Karpenko, A. V. Nazarenko, and Yu. M. Sinyukov, “HBT and initial conditions for hydrodynamic expansion in collisions,” Brazilian Journal of Physics, vol. 37, no. 3A, pp. 1031–1038, 2007.

W. Broniowski, W. Florkowski, M. Chojnacki, and A. Kisiel, “Free-streaming approximation in early dynamics of relativistic heavy-ion collisions,” Physical Review C, vol. 80, no. 3, Article ID 034902, 2009.

S. Pratt, “The long slow death of the HBT puzzle,” Nuclear Physics A, vol. 830, no. 1–4, pp. 51c–57c, 2009.

Iu. A. Karpenko and Yu. M. Sinyukov, “Energy dependence of pion interferometry scales in ultra-relativistic heavy ion collisions,” Physics Letters B, vol. 688, no. 1, pp. 50–54, 2010.

Iu. A. Karpenko and Yu. M. Sinyukov, “Kaon and pion femtoscopy at the highest energies available at the BNL Relativistic Heavy Ion Collider (RHIC) in a hydrokinetic model,” Physical Review C, vol. 81, no. 5, Article ID 054903, 2010.

Iu. A. Karpenko, Yu. M. Sinyukov, and K. Werner, “description of bulk observables in the hydrokinetic model of collisions at the BNL Relativistic Heavy Ion Collider and the CERN Large Hadron Collider,” Physical Review C, vol. 87, Article ID 024914, 2013.

K. Aamodt, B. Abelev, A. A. Quintana, et al., “Elliptic flow of charged particles in Pb-Pb collisions at TeV,” Physical Review Letters, vol. 105, Article ID 252302, 11 pages, 2010.

N. Armesto, N. Borghini, S. Jeon, et al., Eds., “Heavy-ion collisions at the LHC—last call for predictions,” Journal of Physics G, vol. 35, no. 5, Article ID 054001, 2008.

K. Aamodt, A. A. Quintana, D. Adamová, et al., “Two-pion Bose-Einstein correlations in central Pb-Pb collisions at TeV,” Physics Letters B, vol. 696, no. 4, pp. 328–337, 2011.

M. Chojnacki, W. Florkowski, W. Broniowski, and A. Kisiel, “Soft heavy-ion physics from hydrodynamics with statistical hadronization: Predictions for collisions at TeV,” Physical Review C, vol. 78, no. 1, Article ID 014905, 13 pages, 2008.

Yu. M. Sinyukov, S. V. Akkelin, and Y. Hama, “Freeze-out problem in hydrokinetic approach to collisions,” Physical Review Letters, vol. 89, no. 5, pp. 052301/1–052301/4, 2002.

Iu. A. Karpenko and Yu. M. Sinyukov, “Femtoscopic scales in central collisions at RHIC and LHC energies in a hydrokinetic model,” Journal of Physics G, vol. 38, no. 12, Article ID 124059, 2011.

S. A. Bass, M. Belkacem, M. Bleicher, et al., “Microscopic models for ultrarelativistic heavy ion collisions,” Progress in Particle and Nuclear Physics, vol. 41, pp. 255–369, 1998.

M. Bleicher, E. Zabrodin, C. Spieles, et al., “Relativistic hadron-hadron collisions in the ultra-relativistic quantum molecular dynamics model,” Journal of Physics G, vol. 25, no. 9, Article ID 1859, 1999.

Yu. M. Sinyukov and V. M. Shapoval, “Correlation femtoscopy of small systems,” Physical Review D, vol. 87, no. 9, Article ID 094024, 15 pages, 2013.

V. M. Shapoval, P. Braun-Munzinger, Iu. A. Karpenko, and Yu. M. Sinyukov, “Femtoscopic scales in and collisions in view of the uncertainty principle,” Physics Letters B, vol. 725, no. 1–3, pp. 139–147, 2013.

K. Aamodt, N. Abel, U. Abeysekara, et al., “Two-pion Bose-Einstein correlations in pp collisions at GeV,” Physical Review D, vol. 82, no. 5, Article ID 052001, 14 pages, 2010.

K. Aamodt, A. A. Quintana, D. Adamová, et al., “Femtoscopy of pp collisions at and 7 TeV at the LHC with two-pion Bose-Einstein correlations,” Physical Review D, vol. 84, no. 11, Article ID 112004, 22 pages, 2011.

M. M. Aggarwal, Z. Ahammed, A. V. Alakhverdyants, et al., “Pion femtoscopy in collisions at GeV,” Physical Review C, vol. 83, Article ID 064905, 17 pages, 2011.

Z. Chajȩcki and M. Lisa, “Global conservation laws and femtoscopy of small systems,” Physical Review C, vol. 78, no. 6, Article ID 064903, 2008.

Z. Chajecki and M. Lisa, “Conservation laws and multiplicity evolution of spectra at energies available at the BNL Relativistic Heavy Ion Collider,” Physical Review C, vol. 79, no. 3, p. 034908, 2009.

S. V. Akkelin and Yu. M. Sinyukov, “Deciphering nonfemtoscopic two-pion correlations in collisions with simple analytical models,” Physical Review D, vol. 85, no. 7, Article ID 074023, 2012.

M. Janik, “ angular correlations in pp collisions at the LHC registered by the ALICE experiment,” http://arxiv.org/abs/1203.2844

F. Grassi, Y. Hama, and T. Kodama, “Continuous particle emission: a probe of thermalized matter evolution?” Physics Letters B, vol. 355, no. 1-2, pp. 9–14, 1995.

F. Grassi, Y. Hama, and T. Kodama, “Particle emission in the hydrodynamical description of relativistic nuclear collisions,” Zeitschrift für Physik C, vol. 73, pp. 153–160, 1996.

K. A. Bugaev, “Shock-like freeze-out in relativistic hydrodynamics,” Nuclear Physics A, vol. 606, no. 3-4, pp. 559–567, 1996.

F. Cooper and G. Frye, “Single-particle distribution in the hydrodynamic and statistical thermodynamic models of multiparticle production,” Physical Review D, vol. 10, no. 1, pp. 186–189, 1974.

Yu. M. Sinyukov, “Direct conversion of mixed phase into free hadronic gas and transverse pion momenta in ultra-relativistic nuclear collisions,” Zeitschrift für Physik C, vol. 43, pp. 401–409, 1989.

N. S. Amelin, R. Lednicky, T. A. Pocheptsov et al., “Fast hadron freeze-out generator,” Physical Review C, vol. 74, no. 6, Article ID 064901, 2006.

P. Huovinen and H. Petersen, “Particlization in hybrid models,” The European Physical Journal A, vol. 48, article 171, 2012.

M. Laine and Y. Schroder, “Quark mass thresholds in QCD thermodynamics,” Physical Review D, vol. 73, Article ID 085009, 13 pages, 2006.

K. Nakamura, J. Beringer, J.-F. Arguin, et al., “Review of particle physics,” Journal of Physics G, vol. 37, Article ID 075021, 2010.

Yu. M. Sinyukov, R. Lednicky, S. V. Akkelin, J. Pluta, and B. Erazmus, “Coulomb corrections for interferometry analysis of expanding hadron systems,” Physics Letters B, vol. 432, no. 3-4, pp. 248–257, 1998.

G. Goldhaber, S. Goldhaber, W. Lee, and A. Pais, “Influence of Bose-Einstein statistics on the antiproton-proton annihilation process,” Physical Review, vol. 120, no. 1, pp. 300–312, 1960.

S. Pratt, “Validity of the smoothness assumption for calculating two-boson correlations in high-energy collisions,” Physical Review C, vol. 56, no. 2, pp. 1095–1098, 1997.

A. Kisiel, W. Florkowski, W. Broniowski, and J. Pluta, “Femtoscopy in hydrodynamics-inspired models with resonances,” Physical Review C, vol. 73, no. 6, Article ID 064902, 2006.

“ROOT Data Analysis Framework,” http://root.cern.ch/drupal/

O. Utyuzh, G. Wilk, and Z. Włodarczyk, “Modeling Bose-Einstein correlations via elementary emitting cells,” Physical Review D, vol. 75, no. 7, Article ID 074030, 2007.

G. A. Kozlov, O. V. Utyuzh, G. Wilk, and Z. Włodarczyk, “Some forgotten features of the Bose-Einstein correlations,” Physics of Atomic Nuclei, vol. 71, no. 9, pp. 1502–1504, 2008.

Yu. M. Sinyukov, “Boson spectra and correlations for thermal locally equilibrium systems,” Acta Physica Hungarica Series A, vol. 10, pp. 113–136, 1999.

Yu. M. Sinyukov and B. Lorstad, “On intensity interferometry at high multiplicities,” Zeitschrift für Physik C, vol. 61, pp. 587–592, 1994.

Yu. M. Sinyukov, S. V. Akkelin, and R. Lednicky, “Thermodynamics of multi-boson phenomena,” in Proceedings of the 8th International Workshop on Multiparticle Production in Matrahaza, T. Csorgo, S. Hegyi, G. Hwa, et al., Eds., p. 66, World Scientific.

R. Lednicky, V. Lyuboshitz, K. Mikhailov, Yu. M. Sinyukov, A. Stavinsky, and B. Erazmus, “Multiboson effects in multiparticle production,” Physical Review C, vol. 61, no. 3, pp. 349011–3490119, 2000.

S. V. Akkelin, R. Lednicky, and Yu. M. Sinyukov, “Correlation search for coherent pion emission in heavy ion collisions,” Physical Review C, vol. 65, no. 6, pp. 649041–6490416, 2002.

STAR Collaboration, “Pion interferometry in collisions at GeV,” Physical Review C, vol. 71, no. 4, Article ID 044906, 21 pages, 2005.

S. S. Adler, S. Afanasiev, C. Aidala, et al., “Bose-Einstein correlations of charged pion pairs in collisions at GeV,” Physical Review Letters, vol. 93, no. 15, Article ID 152302, 6 pages, 2004.

P. Bozek, “Flow and interferometry in -dimensional viscous hydrodynamics,” Physical Review C, vol. 85, no. 3, Article ID 034901, 9 pages, 2012.

L. Adamczyk, J. K. Adkins, G. Agakishiev, et al., “Freeze-out dynamics via charged kaon femtoscopy in GeV central collisions,” http://arxiv.org/abs/1302.3168

S. Afanasiev, C. Aidala, N. N. Ajitanand, et al., “Charged kaon interferometric probes of space-time evolution in collisions at   GeV,” Physical Review Letters, vol. 103, no. 14, Article ID 142301, 6 pages, 2009.

R. Lednicky, V. L. Lyuboshitz, and M. I. Podgoretsky, “Interference correlations of identical particles in models with closely spaced sources,” Yad. Fiz., vol. 38, p. 251, 1983, Soviet Journal of Nuclear Physics, vol. 38, pp. 147, 1983

M. S. Nilsson, L. V. Bravina, E. E. Zabrodin, L. V. Malinina, and J. Bleibel, “Study of π π correlations at LHC and RHIC energies in pp collisions within the quark-gluon string model,” Physical Review D, vol. 84, no. 5, Article ID 054006, 2011.

E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” The European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010.

E. K. G. Sarkisyan and A. S. Sakharov, “Multihadron production features in different reactions,” AIP Conference Proceedings, vol. 828, pp. 35–41, 2006.

B. Kopeliovich, A. Schäfer, and A. Tarasov, “Nonperturbative effects in gluon radiation and photoproduction of quark pairs,” Physical Review D, vol. 62, no. 5, Article ID 054022, pp. 1–31, 2000.

E. V. Shuryak and I. Zahed, “Understanding nonperturbative deep-inelastic scattering: instanton-induced inelastic dipole-dipole cross section,” Physical Review D, vol. 69, no. 1, Article ID 014011, 2004.

K. Aamodt, A. A. Quintanab, D. Adamová, et al., “Two-pion Bose–Einstein correlations in central Pb-Pb collisions at TeV,” Physics Letters B, vol. 696, no. 4, pp. 328–337, 2011.

D. Antończyk, “Two pion correlations at SPS energies,” Acta Physica Polonica B, vol. 40, no. 4, pp. 1137–1144, 2009.

S. V. Afanasiev, T. Anticic, D. Barna, et al., “Energy dependence of pion and kaon production in central collisions,” Physical Review C, vol. 66, no. 5, Article ID 054902, 9 pages, 2002.

C. Alt, T. Anticic, B. Baatar, et al., “Bose-Einstein correlations of π-π pairs in central collisions at 20A,30A,40A,80A, and 158A GeV,” Physical Review C, vol. 77, no. 6, Article ID 064908, 20 pages, 2008.

J. Adams, C. Adler, M. M. Aggarwal, et al., “Identified particle distributions in pp and collisions at   GeV,” Physical Review Letters, vol. 92, no. 11, Article ID 112301, 6 pages, 2004.

J. Adams, M. M. Aggarwal, Z. Ahammed, et al., “Pion interferometry in collisions at GeV,” Physical Review C, vol. 71, no. 4, Article ID 044906, 21 pages, 2004.

S. S. Adler, S. Afanasiev, C. Aidala, et al., “Identified charged particle spectra and yields in collisions at GeV,” Physical Review C, vol. 69, no. 3, Article ID 034909, 32 pages, 2004.

S. S. Adler, S. Afanasiev, C. Aidala, et al., “Bose-Einstein correlations of charged pion pairs in collisions at   GeV,” Physical Review Letters, vol. 93, no. 15, Article ID 152302, 6 pages, 2004.

S. V. Akkelin and Yu. M. Sinyukov, “Phase-space densities and effects of resonance decays in a hydrodynamic approach to heavy ion collisions,” Physical Review C, vol. 70, no. 6, Article ID 064901, 17 pages, 2004.

S. V. Akkelin and Yu. M. Sinyukov, “Hanbury-Brown-Twiss search for new states of matter in collisions,” Physical Review C, vol. 73, no. 3, Article ID 034908, 2006.

T. Csorgo and B. Lorstad, “Bose-Einstein correlations for three-dimensionally expanding, cylindrically symmetric, finite systems,” Physical Review C, vol. 54, no. 3, pp. 1390–1403, 1996.

S. V. Akkelin, P. Braun-Munzinger, and Yu. M. Sinyukov, “Reconstruction of hadronization stage in collisions at 158 A GeV/c,” Nuclear Physics A, vol. 710, no. 3-4, pp. 439–465, 2002.

P. Csizmadia, T. Csorgo, and B. Lukacs, “New analytic solutions of the non-relativistic hydrodynamical equations,” Physics Letters B, vol. 443, no. 1–4, pp. 21–25, 1998.

G. F. Bertsch, “Meson phase space density in heavy ion collisions from interferometry,” Physical Review Letters, vol. 72, no. 15, pp. 2349–2350, 1994.

G. F. Bertsch, “Meson phase space density in heavy ion collisions from interferometry,” Physical Review Letters, vol. 77, pp. 789–789, 1996.

“The Durham HepData Project,” http://hepdata.cedar.ac.uk/