# Femtoscopic and Nonfemtoscopic Two-Particle Correlations inA+A

Authors: Yu. M. Sinyukov; S. V. Akkelin; Iu. A. Karpenko; V. M. Shapoval
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

The theoretical review of the last femtoscopy results for the systems created in ultrarelativisticA+A, p+p, and p+Pb 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 Pb+Pb and Au+Au 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 p+p collisions at s=7 TeV LHC energy and p+Pb  ones at s=5.02 TeV. The behavior of the corresponding interferometry volumes on multiplicity is compared with what is happening for central A+A 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.


---

## Body

## 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 10-14 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 orv2 coefficients. The obtained LHC results for the transverse momentum dependence of v2(pT) at a given centrality bin (here and below for the quantities depending on one particle momentum pT means the transverse momentum of one particle, whereas for the pair quantities pT means the average pair transverse momentum, pT=(p1T+p2T)/2) 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 pT 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 pT-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 p+p, p+Pb, and e++e- 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 p+p collisions at s=7 TeV LHC energy.As for the elementary particle collisions, likep+p, there is no unambiguous interpretation of the HBT radii pT-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 p+p 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 p+p 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 toA+A 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 A+A 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, τrel→∞, when t→∞, 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 functionf(x,p) in the case of no external forces has the form
(1)pip0∂f(x,p)∂xi=Fgain(x,p)-Floss(x,p).
The terms Fgain and Floss are associated with number of particles which, correspondingly, came to point (x,p) and left this point because of collisions. Term Floss(x,p)=R(x,p)f(x,p) can easily be expressed in terms of the rate of collisions of the particle with momentum p, R(x,p)=〈σvrel〉n(x). The term Fgain has, in general, more complicated integral structure and depends on differential cross-section.The escape fraction,fescσ(x,p), describes the (probabilistic) distribution of the particles that reaches the hypersurface σ without interactions. At asymptotic hypersurface σ where times t→∞, the distribution function corresponds to free quanta: fescσ(x,p)→f(x,p). In general cases, which we will need, when the space-like hypersurface σ is situated at the final times, we define the escape function as
(2)fescσ(x,p)=𝒫σ(x,p)f(x,p),
where 𝒫σ(x,p) is the escape probability for particle to reach freely the hypersurface σ at some point (tσ,xσ) starting from the point (t,x). If the point x belongs to this hypersurface σ, then 𝒫σ(x,p)=1, and so
(3)fescσ(x,p)|σ=f(x,p)|σ.
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 p:
(4)𝒫σ(x,p)=exp(-∫ttσdt¯R(x¯,p)),
where x¯≡(t¯,xσ+(p/p0)(t¯-tσ)). It satisfies the differential equation
(5)1𝒫σ(x,p)pμp0∂∂xμ𝒫σ(x,p)=R(x,p)=Floss(x,p)f(x,p).It follows from (1) and (5) that
(6)pμ𝒫0∂∂xμfescσ(x,p)=𝒫σ(x,p)Fgain(x,p).The formal solution of (6) can be presented in the following form:
(7)fescσ(x,p)=fescσ(x0,p)+∫tσ0(x0)tdt′𝒫σ(x′,p)Fgain(x′,p),
where x=(t,x),σ:t=const,fescσ(x0,p) with x0≡(tσ0(x0),x0=x+(p/p0)(tσ0(x0)-t)) corresponds to the portion of the system, which propagates without collisions until some time t starting from initial hypersurface σ0, x′≡(t′,x+(p/p0)(t′-t)).The expression (7) for escape function fescσ(x,p) can be explained in simple heuristic way as follows. Let us split the distribution function at each space-time point into two parts: f(x,p)=fintσ(x,p)+fescσ(x,p), x=(t,x). The first part, fintσ(x,p), describes the fraction of the system which will continue to interact before reaching the hypersurface σ. The second fraction, fescσ(x,p), describes the particles that reach the hypersurface σ without interactions. Denote again by x′≡(t′,x+(p/p0)(t′-t)) the space-time point where the particle at x with momentum p would be, if it moved freely. Consider, at each vicinity of the phase-space point (x,p), the number of particles that have escaped from the interacting system during the time interval (t′,t′+dt′). 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 (x′,p) just after the interaction during the time dt′. Indeed, if some particles from the interacting part of the system fint(x′,p) do not interact during the given time interval (t′,t′+dt′), then they will interact without fail at some future time; thus, they cannot contribute to the additional portion of particles escaping during dt′. Therefore, the additional contribution Δfescσ(x,p;t′), from the time interval (t′,t′+dt′) to the distribution function fescσ(x,p), is Δfescσ(x,p;t′)=𝒫σ(x′,p)Fgain(x′,p)dt′ for t′<t and is zero for t′>t. Here 𝒫σ(x′,p) is probability for any given particle at point x′ with momentum p 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,
(8)p0dNdp=〈ap+ap〉,p10p20dNdp1dp2=〈ap1+ap2+ap1ap2〉,…,
that are constructed in the standard way by means of the averages of product of creation and annihilation operators. Namely, on the hypersurface σout(9)〈ap1+ap2〉=∫σout‍dσμ(x)pμexp(iqx)fesc(x,p),
where p=(p1+p2)/2, q=p1-p2. Then, using the Gauss theorem, (6) and (1) analytically continued to off-mass-shell four-momenta p, and taking into account that ∂μ(pμexp(iqx))=0 for particles on mass shell, one can get
(10)〈ap1+ap2〉=p0∫σ0‍d3σμ(x0)pμfescσ(x0,p)e-iqx0+p0∫σ0σ‍d4x𝒫σ(x,p)Fgain(x,p)e-iqx,(11)〈ap1+ap2〉=p0∫σ0‍d3σμ(x0)pμf(x0,p)e-iqx0+p0∫σ0σ‍d4x×[Fgain(x,p)-Floss(x,p)]e-iqx,
respectively, where 𝒫σout(x,p) is the probability for particle to reach the hypersurface σ from the space-time point x, fescσ(x0,p)=f(x0,p)𝒫σ(x0,p) is a portion of the particles at the initial hypersurface σ0 that reach hypersurface σ moving freely without interactions, f(x0,p) is the distribution function at σ0, and p0𝒫Fgain 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 Fgain and the collision term R. 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 Fgain(x,p)=Floss(x,p)≡0 and therefore 𝒫σ(x0,p)≡1, 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 Fgain(x,p)=Floss(x,p)≠0. 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 σ0 has nonspace-like sites, Δσμ(x1), then for some momenta p1 we have p1μΔσμ(x1)<0, 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 σ0 near the point x1. This negative contribution will be compensated when the particle crosses the hypersurface again in another point x2 of the hypersurface to go out. See the possible trajectory of the particle 1 in Figure 1. Of course, the Gauss theorem applied to the enclosed hypersurface σ0⊕σ guarantees the correct positive result for particle momentum spectra, no matter how complicated structure σ0 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σ0 (solid line) and spectra formation at σ (dashed line). The very initial system forms at the hypersurface σin (double line) and then evolves. The Gauss theorem in (10) and (11) is applied to enclosed hypersurface σ0⊕σ. In the case of large opacity fescσ(x1,p1)≈0 for particle 1, also is small for particle 2, and so pμΔσ0μ(x1)fescσ(x1,p)≈0, while for the distribution function f(x1,p) the similar value pμΔσ0μ(x1)f(x1,p) is not zero, and for p=p1 it is negative if Δσ0 is a nonspace-like site of the hypersurface σ0.Typically at a modeling ofA+A collisions the initial conditions for the Boltzmann equations are selected at the hypersurface σ0 close to the hadronisation hypersurface. Such a hypersurface is not everywhere a space-like one, but at sites with pμΔσμ<0 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, 𝒫σ(x0,p)≈0, and so the escape function is about zero, fescσ(x0,p)≈0, resulting in suppression of the negative contributions to particle momentum density in (10) from the σ0 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 σ0 (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 σ0 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 σ0, including “bad” nonspace-like sites. It is implemented through a continuous hydrodynamic description of the matter crossing the σ0 that is utilized for calculation of the escape probability and emission function. Therefore no discontinuity in the energy and momentum flows through σ0 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 σ0 but brings into play the distribution function f(x,p) 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 x1 a negative number of the particles with corresponding momenta p1has to be “injected” (see Figure 1). Typically, to cure this dilemma people just cut them in the full distribution function f(x1,p), that leads to discontinuity of the energy and momentum flows through the corresponding sites of σ0 [50]. Such a cut, introduced through the Heaviside function Θ(dσμpμ) 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 nμ=nuμ:∂μnμ(x)|Δσ0(x1)≠0. 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 momentap1 moving towards the nonspace-like sites where p1μΔσ0(x1)<0. These particles have to collide with the particles that move from the inside of Δσ0(x1) 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 ∂μnμ(x)|Δσ0(x1)=0. 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 σ0 [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σ0 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 (2+1)D 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 u, d, and s quarks, including σ-meson (f0(600)). With such an equilibrated evolution the system reaches the chemical freeze-out isotherm with the temperature Tch=165 MeV. At the second stage with T<Tch, 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 σ:t=const (with Tσ(r=0)=Tch) 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:
(12)fi(τ,x→,p→)=fil.eq.(τ0,x→(τ0)(τ),p→)×exp(-∫τ0τR~i(s,x→(s)(τ),p→)ds)+∫τ0τdλ[fil.eq.(λ,x→(λ)(τ),p→)×R~i(λ,x→(λ)(τ),p→)+G~idecay(λ,x→(λ)(τ),p→)]×exp(-∫λτ[R~i(s,x→(s)(τ),p→)+D~i(s,x→(s)(τ),p→)]ds∫τ0τdλ).
The different terms in (12) correspond to the following:R~i(λ,x→,p→) is the collision rate of the ith sort of hadrons with the rest of particles, G~idecay(λ,x→,p→) is an income of particles into the phase-space point owing to resonance decays, and D~i(λ,x→,p→) 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 RelativisticA+A 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:(13)C(p1,p2)=p10p20(dN/dp1dp2)p10(dN/dp1)p20(dN/dp2).
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 inp-+p 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:
(14)〈ap1+ap2+ap1ap2〉=〈ap1+ap1〉〈ap2+ap2〉+〈ap1+ap2〉〈ap2+ap1〉.
Then using (8) and (10) one can express the correlation function (13) for chaotic sources
(15)C(p,q)=1+∫d4x1d4x2S(x1,p)S(x2,p)eiq(x1-x2)(∫‍d4x1S(x1,p1))(∫‍d4x2S(x2,p2)),
where p=(p1+p2)/2, q=p1-p2, S(x,p) is the emission function; if fesc|σ0σ≈0, then S(x,p)=p0Fgain(x,p)𝒫σ(x,p) as it follows from (10) (in general case it includes in addition the contribution from σ0).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, andC(p,q) 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 q→ and the average momentum of the particle pair p→. 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 pT bin looks like
(16)C(q→)=∑i∑j≠iδΔ(q→-pi→+pj→)(1+cos((pj-pi)(xj-xi)))∑i∑jδΔ(q→-pi→+pj→),
where δΔ(x)=1 if |x|<Δp/2 and 0 otherwise, with Δp being the bin size in histograms. We decompose the relative momentum q→ into (qout,qside,qlong) 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, λ<1/p0 [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
(17)C(q)=1+λ·exp(-Rout2qout2-Rside2qside2-Rlong2qlong2-Ros2qoutqside-Rsl2qsideqlong-Rol2qoutqlong).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|η|<0.5. 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 Ros2, Rsl2, and Rol2 are neglected.The comparison of interferometry radii, calculated in hHKM with the experimental data froms=200 GeV Au+Au 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.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 c= 0–5% events with a gray solid line. The ICs used for hydrodynamic evolution: τ0=0.1 fm/c with zero initial transverse flow and the MC-Glauber profiles for initial entropy density.
(a)
(b)
(c)
(d)From Figure2 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 Rout/Rside 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, v2, 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 bothRside and Rlong (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 dN/dy, 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.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.
(a)
(b)
(c)
(d)In Figure4 (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 mT as for pions is slightly larger than for pions, at least for intermediate mT, because the same mT corresponds to smaller pT for kaons (and the duration of emission process is highly pT 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.Transverse mass dependence of Gaussian radii (a)Rout, (b) Rside, and (c) Rlong for midrapidity kaon pairs from the 30% most central Au+Au collisions at sNN=200 GeV. The results from STAR [72] and PHENIX [73] are compared to the calculations in hHKM model (dotted lines).
(a)
(b)
(c)
(d)


## 4. Correlation Femtoscopy for Small Systems

In the work [39] it is shown that for small systems formed in particle collisions (e.g., pp, e+e-), 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 1+cos((x1-x2)(p1-p2)) to the weight of the pion pair in (16) is not correct for very closely located points x1 and x2, 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 uμ. In hydrodynamic/hydrokinetic approach the role of such a source at given pair's half-momentum bin near some value p is played by the fluid element or piece of the matter with the sizes equal to the homogeneity lengths λ(p) [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 pT 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]
(18)Rout*(p)=Rout(p)coshyTcosh(yT-ηT),Rside*=Rside=λside,Rlong*=Rlong=λlong,λout*=λoutcoshyTcosh(yT-ηT),pout*p0*T*=TsinhyTcosh(yT-ηT).
Here Rout*2(p)=λout*2+(pout*2/p0*2)T*2, T is the emission duration, ηT is a rapidity of the source in transverse direction, and yT=(y1T+y2T)/2 is half sum of transverse rapidities of the particles forming the pair. Note that yT*=yT-ηT, and if rapidity of the pair is equal to the rapidity of the source, yT*=0, 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 p+p 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 Rout*=Rout and Rout*=γRout(γ=coshyT).In [40] the quantum corrections are calculated at each pT-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 〈(Δx)2〉~k2 (k 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 ψxi(p,t) corresponding to a boson with mass m emitted at the time ti from the point xi as a wave packet with momentum variance k which then propagates freely:
(19)ψxi(p,t)=eipxi-iEteiφ(xi)f~(p),
where φ(xi) is some phase and f~ defines the primary momentum spectrum f(p) that we take in the Gaussian form,
(20)f(p)=f~2(p)=1(2πk2)3/2e-(p2/2k2),
with the variance k2=mT. The effective temperature of particle emission in the local rest frames in HKM, T, is close to the chemical freeze-out temperature Tch.The amplitude of the single-particle radiation from some 4-volume is supposed to be a superposition of the wave functionsψx(p) with the Gaussian coefficients ρ^(x)=ρ(x), where ρ(x) is the probability distribution in the case of random phases:
(21)ρ(x,t)∝e-∑ixi*2/2λi*2(p)-t*2/2T*2(p).
Then in the rest frame of the source
(22)A(p,t)=c∫‍d4xψx(p,t)ρ^(x),
where c is the normalization constant.The single- and two-particle spectra, averaged over the ensemble of emission events with partially correlated phasesφ(x) are
(23)W(p)¯=c2∫‍d4xd4x′eip(x-x′)ρ^(x)ρ^(x′)mmmmi×f(p)〈ei(φ(x)-φ(x′))〉,W(p1,p2)¯=c4∫‍d4x1d4x2d4x1′d4x2′ei(p1x1+p2x2-p1x1′-p2x2′)·f(p1)f(p2)ρ^(x1)ρ^(x2)ρ^(x1′)ρ^(x2′)×〈ei(φ(x1)+φ(x2)-φ(x1′)-φ(x2′))〉.The phase averages are associated with corresponding overlap integrals [39]
(24)〈ei(φ(x)-φ(x′))〉=Gxx′=Ixx′=|∫‍d3rψx(t,r)ψx′*(t,r)|,(25)〈ei(φ(x1)+φ(x2)-φ(x1′)-φ(x2′))〉=Gx1x1′Gx2x2′+Gx1x2′Gx2x1′-Gx1x2′Gx2x1′Gx1x2,
where ψxi(t,r)=(1/(2π)3/2)∫f(p)e-ip(r-xi)e-i(p2/2m)(ti-t)d3p are the wave functions of single bosonic states in coordinate representation.Then the correlation functionC(p,q) can be expressed through the homogeneity lengths in the local rest frame RL≡λlong*(p), RS≡λside*(p), and RO≡λout*(p) 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:
(26)C(p,q)=W(p1,p2)¯W(p1)¯W(p2)¯=1+e-qO2RO2(4k02RO2/(1+4k02RO2))-qS2RS2(4k02RS2/(1+4k02RS2))×e-qL2RL2(4k02RL2/(1+4k02RL2))×e-((q·p)2T2/m2)(4k2T2/(1+4k2T2))-Cd(p,q),
where k02=k2/(1+αk4T2/m2), parameter α(k2R2) is defined from the model numerically (it is the order of unity for R~1 fm and tends to zero for the large sources—see [39] for details), and the subtracted term
(27)Cd(p,q)=e-2qO2k02RO4(1+8k02RO2)/(1+4k02RO2)(1+8k02RO2+8k04RO4)×e-2qS2k02RS4(1+8k02RS2)/(1+4k02RS2)(1+8k02RS2+8k04RS4)×e-2qL2k02RL4(1+8k02RL2)/(1+4k02RL2)(1+8k02RL2+8k04RL4)×e-2k2T4(p·q)2(1+8k2T2)/m2(1+4p2T2)(1+8k2T2+8k4T4)×F(k02Ri2,k2T2),F(k02Ri2,k2T2)=(k0k)3/2×((1+4k2T2)(1+8k2T2+8k4T4)(1+4k02RO2)(1+8k02RO2+8k04RO4)·(1+4k02RS2)(1+8k02RS2+8k04RS4)(1+4k02RL2)(1+8k02RL2+8k04RL4))1/2
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 termCd(p,q) in (26), and fit the correlation function with the Gaussian form (17), we obtain the femtoscopic radii Rout, Rside, and Rlong related to the standard ones Rout,st, Rside,st, and Rlong,st as follows:
(28)Rout2Rout,st2=(RO24k02RO21+4k02RO2+T2vout24k2T21+4k2T2)×(RO2+T2vout2)-1,Rside2Rside,st2=4k02RS21+4k02RS2,Rlong2Rlong,st2=4k02RL21+4k02RL2,
where vout=pout*/p0*≪1 according to the nonrelativistic approximation. At large source sizes, for example, when the homogeneity lengths correspond to A+A-collisions, k02R2≫1,  k2T2≫1, all these ratios tend to unity.The mean emission duration is supposed to be proportional to the average system's size,T=a(RO+RS+RL)/3, that leads to the quadratic equation expressing RO (and T) through Ri,st. The latter are connected with ones taken in the global reference system according to transformation laws (18). The value a 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 Rout, Rside, and Rlong 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 R~1 fm is equal or even higher than that the individual emitters have, T=k2/m, 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 ofp+p 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 ofp+p collisions at the top RHIC and s=0.9 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 p+p collisions, at least, for high multiplicity events. Although the applicability of   hydrodynamics to p+p 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 p+p 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 p+p collisions. Therefore, one can conclude that the problem of theoretical description of the interferometry radii in p+p 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 LHCp+p collisions at s=7 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 p+p collision, close to the size of gluon spots [78, 79] in proton moving with a speed v≈c. In detail, the initial boost-invariant tube for p+p collisions has the energy density Gaussian distribution in transverse plane with the width (rms) R=0.3 fm [78, 79], and, following [30], we attribute it to the initial proper time τ0=0.1 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 Figures5 and 6, taken from [40], for the two multiplicity classes 〈dNch/dη〉=9.2 and 〈dNch/dη〉=17.9 we present the three curves for interferometry radii as a function of pT: 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, γ=1. For this case a≈1, k=0.18 GeV/c, and pT*=0.13 GeV/c. The α parameter values linearly increase with pT from 1.15 to 1.35 for the 〈dNch/dη〉=9.2 case and from 1.02 to 1.1 for 〈dNch/dη〉=17.9. As one can see, the quantum corrected pT-dependency of the radii gets closer to the experimental values, but for large pT the corrections are insufficient to fully describe the observable femtoscopy scales behavior. This fact may indicate that sources of particles with large pT cannot be described in hydrodynamical approximation. Note that just for such large pT the nontrivial baseline corrections, already provided in presented experimental data, are very essential.Interferometry radii dependency onpT, 〈dNch/dη〉=9.2. The data are taken from [42].
(a)
(b)
(c)
(d)Interferometry radii dependency onpT, 〈dNch/dη〉=17.9. The data are taken from [42].
(a)
(b)
(c)
(d)Considering the multiplicity dependence of femtoscopy scales inp+p and p+Pb 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 R0g of the system: roughly, if the pressure is p=c02ϵ, then the transverse acceleration a=∇xTp/ϵ∝p(xT=0)/(R0gϵ)=c02/R0g. The interferometry radii RT, 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 R0g, the averaged transverse velocity 〈|vT|〉, and inverse of the temperature β at some final moment τ [17, 18, 90, 91]:
(29)RT=Rg(τ)1+(2/π)〈|vT|〉2βmT≈R0g(1+τ2c022(R0g)2-βmTτ2c02π2(R0g)2).The result (29) for the HBT radii depends obviously on R0g 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, Fgain(t,x)=Floss(t,x), 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:
(30)Vint(τ)≃C(s)dN/dy(τ)〈f〉τTeff3(τ),
where 〈f〉 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 dN/dy [88, 89]. As for the effective temperature of the hadron spectra, Teff(τ)=T(τ)+m(〈vT(τ)〉2/2), one can see that, when the system's temperature T drops, the mean vT2 increases; therefore Teff does not change much during the evolution (it slightly decreases with time for pions and increases for protons). Hence Vint, 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.Figure7 shows the dependency Vint(〈dNch/dη〉) for considered case of p+p collisions at the LHC, s=7 TeV, and for the most central (only!) collisions of nuclei having similar sizes, Pb+Pb and Au+Au, at the SPS, RHIC, and LHC. We have also put on the plot our prediction for the interferometry volume of pPb-system, that has initially larger size than pp-system. As one can see, the different groups of points corresponding to p+p, p+Pb, and A+A 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 A of colliding objects and the collision centrality c.Figure 7
The interferometry volume dependency on charged particles multiplicity. The curve fragments in the middle correspond to the HKM prediction forp+Pb collision at the LHC energy s=5.02 GeV. The upper one is related to the initial transverse system size R=1.5 fm and for the two lower ones R=0.9 fm. The curves at the left and at the right represent the HKM results for p+p and A+A central collisions respectively, compared to the experimental data at AGS, SPS, RHIC, and LHC, taken from papers [42, 80–87]. The pp volumes are calculated as a product RoutRsideRlong 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 inp+p 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 p+Pb collisions need some small corrections only for its minimal values corresponding to the initial transverse size of pPb 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 inp+p 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(31)C(p1,p2)=P(p1,p2)P(p1)P(p2),
where P(p1,p2) is the probability of observing two particles with momenta p1 and p2, while P(p1) and P(p2) designate the single-particle probabilities. In heavy ion collisions
(32)C(p1,p2)≈CF(p,q),
where p=(p1+p2)/2, q=p2-p1, and CF 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, CNF, from the total correlation function:
(33)C(p1,p2)=CF(p,q)CNF(p,q).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 p+p 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 N pions are produced with momenta p1,…,pN in (N+X) 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 N! permutations of the particle momenta pi. Then single-particle probability, PN(p1), and the two-particle probability, PN(p1,p2), can be written as
(34)PN(p1)=1N∑i=1N‍∫‍dΩp*Ei*δ(3)(p1-pi*)P^N(p1*,…,pN*),(35)PN(p1,p2)=1N(N-1)∑i≠j=1N‍∫‍dΩp*Ei*δ(3)(p1-pi*)Ej*×δ(3)(p2-pj*)P^N(p1*,…,pN*),
where dΩp=(d3p1)/E1⋯(d3pN)/EN. The nonsymmetrized N-pion probability density in such events reads as
(36)P^N(p1,…,pN)=1K∑X‍∫‍dΩkδ(4)(pa+pb-∑i=1N‍pi-∑j=1X‍kj)×|MN+X(p1,…,kX)|2,
where MN+X(p1,…,kX) is nonsymmetrized (N+X)-particle production amplitude, pa and pb are 4-momenta of colliding particles (protons), and K is the normalization factor,
(37)K=∑X‍∫‍dΩkdΩpδ(4)(pa+pb-∑i=1N‍pi-∑j=1X‍kj)×|MN+X(p1,…,kX)|2.
Expression (36) for P^N(p1,…,pN) is rather complicated because, in particular, it depends on X particles that are produced in addition to N 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 N 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 N-pion subsystem in a small midrapidity region of the total system. Then, motivated by (36), we assume that a nonsymmetrized N-pion probability density can be written as
(38)P^N(p1,…,pN)=1Kδ(p1,…,pN)FN(p1,…,pN),
where FN(p1,…,pN) is a nonsymmetrized function of pionic momenta, δ(p1,…,pN) denotes average constraints on the N-pion states that appear due to energy and momentum conservations in multiparticle production events, and we assume that
(39)δ(p1,…,pN)=δ(2)(pT1+pT2+⋯+pTN),
where pT1,pT2,…,pTN are transverse components of the momenta of the N particles. Then the normalization factor is
(40)K=∫‍dΩpδ(p1,…,pN)FN(p1,…,pN).Now, to describe the nonfemtoscopic pion correlations in a simple analytical model, we assume that there are no other correlations in the production ofN-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 N≫1(41)FN(p1,…,pN)=f(p1)⋯f(pN)Q(p1,p2)⋯Q(pN-1,pN),
where Q(pi,pj) denotes the jet-like correlations between momenta pi and pj; existence of such correlations means that FN cannot be expressed as a product of one-particle distributions. Then, utilizing the integral representation of the δ-function by means of the Fourier transformation, δ(2)(pT)=(2π)-2∫‍d2rTexp(irTpT), and accounting for (34), (38), (39), and (41), the single-particle probability reads as
(42)PN(p1)=1(2π)2K∫‍d2rTGN(p1,rT),
where
(43)GN(p1,rT)=∫dΩp*E1*δ(3)(p1-p1*)×eirT(pT1*+⋯+pTN*)×FN(p1*,…,pN*).
A possibility of different cluster configurations of particles means, in particular, that registered particles with momenta p1 and p2 can belong either to different minijets or to the same minijet. Then, taking into account (35), (38), (39), and (41), we get
(44)PN(p1,p2)=NN(N-1)PN1jet(p1,p2)+N(N-1)-NN(N-1)PN2jet(p1,p2),
where
(45)PN1jet(p1,p2)=1(2π)2K∫‍d2rTGN1jet(p1,p2,rT),(46)PN2jet(p1,p2)=1(2π)2K∫‍d2rTGN2jet(p1,p2,rT),(47)GN1jet(p1,p2,rT)=∫‍dΩp*Ei*δ(3)(p1-p1*)Ej*×δ(3)(p2-p2*)eirT(pT1*+⋯+pTN*)FN,(48)GN2jet(p1,p2,rT)=∫‍dΩp*Ei*δ(3)(p1-p1*)Ej*×δ(3)(p2-p3*)eirT(pT1*+⋯+pTN*)FN.
Here FN≡FN(p1*,…,pN*). 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 N.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,
(49)f(pi)=Eiexp(-pi,T2TT2)exp(-pi,L2TL2),(50)Q(pi,pj)=exp(-(pi-pj)2α2),
where TT, TL, and α are some parameters, and we assume that TL≫TT. In accordance with ALICE baseline obtained from the PHOJET event generator simulations, we assume that only qinv=(p2-p1)2-(E2-E1)2 is measured for each pT bin. Assuming that longitudinal components of the registered particles are equal to zero, p1L=p2L=0, we approximate qinv2 as
(51)qinv2≈qT2(m2+pT2sin2ϕm2+pT2),
where ϕ denotes unregistered angle between pT and qT, pTqT=|pT||qT|cosϕ. Then
(52)CNF(|pT|,qinv)=∫02π‍dϕPN(p1,p2)∫02π‍dϕPN(p1)PN(p2),
and, taking into account (44), we get
(53)CNF(|pT|,qinv)=N-2N-1×(1N-2CN2jet(|pT|,qinv)+1N-2CN1jet(|pT|,qinv)),
where
(54)CN2jet(|pT|,qinv)=∫02π‍dϕPN2jet(p1,p2)∫02π‍dϕPN(p1)PN(p2),(55)CN1jet(|pT|,qinv)=∫02π‍dϕPN1jet(p1,p2)∫02π‍dϕPN(p1)PN(p2).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 N-particle production process depends on the value of N and disappears at N→∞. Since one considers a subsystem of N pions but not the total system, to weaken the influence of the total transverse momentum conservation on pions we will consider CM1jet and CM2jet with M>N instead of CN1jet and CN2jet in (53). This is the simplest way to account for a weakened conservation law in our model. At the same time, the factor 1/(N-2) 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 0.9). Then (53) gets the form
(56)CNF(|pT|,qinv)=Λ(pT)×(1N-2CM2jet(|pT|,qinv)+1N-2CM1jet(|pT|,qinv)).The results of our calculations of the nonfemtoscopic correlation functionsCNF 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 Nch≥12 bin in p+p collisions at s=900 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 qinv 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 M=50, TT=α=0.65 GeV (to minimize the number of fit parameters, we fixed TT=α for all calculations. Note that with these parameter values the mean transverse momentum 〈pT〉 is about 0.58 GeV), and the fitted values of N are different for like-sign and unlike-sign pion pairs, namely, N±±=20 for the former and N+-=11 for the latter. The relatively high value of M can be interpreted as a residual effect on the pion subsystem of total energy-momentum conservation in a multiparticle production process. The relation N+-<N±± between fitted N 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, CNF, 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 qinv 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.The nonfemtoscopic correlation functions of like-sign (a) and unlike-sign (b) pions in the0.1<pT<0.25 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)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 the0.7<pT<1.0 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.
(a)
(b)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 N-particle probability density is defined as
(57)PN(p1,p2,…,pN)=∑i‍w(ui)PN(p1,p2,…,pN;ui),
where PN(p1,p2,…,pN;ui) is the N-particle probability density for some ui type of the initial conditions, and w(ui) denotes the distribution over initial conditions, ∑i‍w(ui)=1. 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
(58)PN(p1,p2,…,pN;ui)=f(p1;ui)f(p2;ui)⋯f(pN-1;ui)f(pN;ui),
where we normalize f(p;ui) as ∫‍(d3p/E)f(p;ui)=1, and then K=1; see (38) and (40). The two-particle nonfemtoscopic correlation function CNF then reads as
(59)CNF(p1,p2)=∑iw(ui)f(p1;ui)f(p2;ui)∑iw(ui)f(p1;ui)∑jw(uj)f(p2;uj).Evidently, the different type of fluctuation, that is, the form of the distributionw(ui), 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
(60)w(uT)=α2πexp(-uT2α2),(61)f(p;uT)=β2γπ3/2Eexp(-(pT-uT)2β2)exp(-pL2γ2),
and normalization is chosen in such a way that ∫‍d2uTw(uT)=1 and ∫‍(d3p/E)f(p;uT)=1. 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 pT 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 CNF decreases with qT2,
(62)CNF(p,q)~exp(-β42(α2+β2)qT2),
and this means (after taking into account (51) and (52)) that CNF decreases with qinv2 too, which is similar to the behavior of CNF 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 qinv 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) forA+A 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 Au+Au and Pb+Pb 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 p+p 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 inA+A collisions but also, at least for large multiplicities, in p+p events. These scales are well reproduced for not large pT. Whether it means the validity of the hydrodynamic mechanism for the bulk matter production in the LHC p+p 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 withqinv for fairly high qinv of the nonfemtoscopic two-pion correlation functions of small systems, whereas the former mostly determines a decrease of the ones at relatively low qinv. 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 s=900 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 inp+p 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 qinv to the ones calculated in the PHOJET-like generators, where the nonfemtoscopic correlations for low qinv 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 p+p 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.



---

*Source: 198928-2013-09-01.xml*