Interferometry correlations in central p+Pb collisions

We present results on interferometry correlations for pions emitted in central p+Pb collisions at sNN=5.02TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}}=5.02~\hbox {TeV}$$\end{document} in a 3+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3+1$$\end{document}-dimensional viscous hydrodynamic model with initial conditions from the Glauber Monte Carlo model. The correlation function is calculated as a function of the pion pair rapidity. The extracted interferometry radii show a weak rapidity dependence, reflecting the lack of boost invariance of the pion distribution. A cross term between the out and long directions is found to be nonzero. The results obtained in the hydrodynamic model are in fair agreement with recent data of the ATLAS Collaboration.


Introduction
In relativistic p+Pb collisions at the CERN Large Hadron Collider (LHC) a small region of large density matter is formed, which makes possible the generation of a collective flow in the expansion [1]. A number of signatures of collectivity due to final state interactions have been observed experimentally in small collision systems at relativistic energies (see e.g. [2]). We use the viscous hydrodynamic model to describe the dynamics in p+Pb interactions. We note that the initial state effect can also explain the observed two-particle correlations [3].
Quantum interferometry correlations for identical particles can serve as a measure of space-time correlations in the source [4,5]. These interferometry (also named Hanbury Brown-Twiss (HBT)) correlations have been studied in both elementary and nuclear collisions. An estimate of the size of the emission region, in the form of the HBT radii, can be extracted from a fit to the interferometry correlations. HBT radii have been measured in p+Pb collisions [6,7] and can be reproduced fairly well in hydrodynamic models posing a e-mail: piotr.bozek@fis.agh.edu.pl Glauber model or color glass condensate initial conditions [8][9][10].
In collisions of symmetric systems the correlation function in relative momentum of the pair is usually parametrized using three HBT radii [11,12]. For a source without forwardbackward symmetry, e.g. for pairs at forward/backward rapidity, an additional cross term can appear in the correlations function [13]. Such a term is predicted to be significant in the case of asymmetric collisions, d+Au or p+Pb [14].
Recently the ATLAS Collaboration has presented results on the interferometry correlations in p+Pb collisions for different rapidities of the pion pair [15]. The HBT radii show a rapidity dependence. The size of the emission region is larger on the Pb going side. The cross term is also found to be nonzero, reflecting the rapidity dependence of the charged particle density.
In this paper we present results for the interferometry radii in p+Pb collisions with 0-1% centrality calculated in a wide range of rapidities corresponding to the pioneering measurements of the ATLAS Collaboration. As in the experiment, an exponential ansatz is used that better reproduces the simulated correlation function. We show that the hydrodynamic model with Glauber model initial conditions can semiquantitatively reproduce the experimental observations. For the first time, the presented model calculations allow one to compare quantitatively the predicted magnitude of the new cross term in the correlation function with experiment. Good agreement is found. It means that the space-momentum correlations built in during the hydrodynamic expansion, including the predicted forward-backward asymmetry, are consistent with experiment.

Hydrodynamic model and HBT correlations
We describe the evolution of the matter created in p+Pb collisions using the viscous hydrodynamic model, with initial entropy density given by the nucleon Glauber Monte Carlo model. The details and the parameters of the model can be found in [8]. The model describes well the spectra and the azimuthal flow coefficients. At the end of the hydrodynamic evolution pions are emitted from the freeze-out hypersurface. In the paper we use the rapidity in the nucleon-nucleon center of mass frame, which is shifted by 0.465 rapidity units from the laboratory frame at the LHC. The Pb is going in the negative rapidity direction. The central rapidity region |y| < 0.5 in the LHC laboratory frame corresponds approximately to −1 < y < 0 in the nucleon-nucleon center of mass frame.
Pairs of same-charge pions emitted at positions x 1 and x 2 , with momenta p 1 and p 2 , are counted. The 3-dimensional correlation function is constructed by binning in relative (q = p 1 − p 2 ) and average transverse (k T = |p 1 + p 2 |/2) momentum of the pair, in the longitudinal comoving system [16]. To increase the statistics for the pion pairs, for each hydrodynamic freeze-out hypersurface 500 realistic events are generated (using a statistical emission procedure [16]) and combined together.
We use a symmetrized plane wave two-pion wave function in the definition of the correlation function. Therefore, no corrections are made for Coulomb interaction in simulated correlation function. The fitted form of the correlation function is where with the three components of the relative momentum q out , q side , and q long [11,12]. The exponential ansatz (1) is the same as the one used by the ATLAS Collaboration [15], the fitted HBT parameters R E out , R E side , and R E long and the cross term R E mix are different from the cases of the Gaussian ansatz used in Ref. [14]. To indicate that the parameters in the formula correspond to the exponential ansatz we add a superscript R E in the notation. A larger emission region gives larger HBT parameters R E in the pion correlation function. Note, however, that when using the exponential ansatz, the value of the fitted parameter cannot be directly interpreted as the size of the effective emission source, although we still use the name HBT radii for the fitted parameters.
Three cross sections of the interferometry correlation function are shown in Fig. 1, obtained in the hydrodynamic model for the pair transverse momentum 0.2 GeV < k T < 0.3 GeV and rapidity −1 < y ππ < 0. In the general case there is no reflection symmetry q out ↔ −q out or q long ↔ −q long . Both negative and positive values for the (2) one identifies terms corresponding to the three relative momenta q,

(d)
The formula used by the ATLAS Collaboration deviates from the standard Gaussian ansatz in two respects. First, the exponential form is used instead of a Gaussian. It is justified by the fact that in the range of the fit the exponential ansatz works better (compare Fig. 1 from this work and Fig. 1 from [8]). Second, the terms in (3)  gives the scale of the term q out q long we use the symbol R E mix to denote the cross-term parameter in the metric (2).

Transverse momentum dependent interferometry correlations
In Fig. 2 are shown the three HBT radii and the cross term R 2 mix as functions of the average transverse momentum k T , for −1 < y ππ < 0. The values of the fitted parameters depend slightly on the fit range for q. We have varied the range of the relative momentum in the fit q low < |q| < 0.3 GeV, q low ∈ [0.02, 0.03] GeV. The vertical size of the boxes in Figs. 2 and 3 represents an estimate of the uncertainty in the fitted parameters obtained by varying the fit range, for the experimental data the vertical size of the boxes represents the experimental error.
The size of the effective emission region decreases with k T . This indicates the existence of strong correlations between the momentum and the emission point for the emitted particles. Due to the collective flow pion pairs of high transverse momentum are effectively emitted from a smaller region of the source [17]. The model describes well the experimental data for R E out , while R E side is slightly underestimated. The pion distribution in rapidity is not flat, which makes the correlations asymmetric in the forward-backward direction. We find a nonzero cross term R E mix , the magnitude of this term decrease with k T in a similar way as R E out and R E long . The asymmetry of the interferometry correlation function can be related to the forward-backward asymmetry of the average emission times of pions from the source (Fig. 4 in [14]). Pions are emitted earlier on the proton going side and the correlation coefficient between the longitudinal position and time becomes nonzero.

Rapidity dependent interferometry correlations
In Fig. 3 presents the rapidity dependence of the HBT parameters and the R E mix cross term for pion pairs with 0.2 GeV < Fig. 3 Interferometry radii in central p+Pb collisions as a function of average pair rapidity y ππ . Results of 3-dimensional hydrodynamic calculations are compared to data from the ATLAS Collaboration [15]. R E out , R E side , R E long , and R E mix are shown in (a) through (d) [fm] Both the experimental data and the simulation results show some dependence on the pion pair rapidity y ππ . However, the experimentally observed rapidity dependence is slightly stronger. We notice that the model correctly describes the experimental data for the cross term R E mix (panel (d)) in Fig. 3. Both the magnitude and the sign of the cross term are fairly well reproduced.
The rapidity distribution of pions emitted in central p+Pb collisions is far from boost invariant (Fig. 4). The asymmetry results from the large difference in the number of participants between the two projectiles, one for the proton and 19 on average for the Pb nucleus. The change in the number of emitted pions with rapidity determines the change in the effective size of the emission region. The sign of the cross term R E mix follows approximately the reverse sign of the slope of the rapidity distribution of pions.

Summary
We analyze the interferometry correlations for pions emitted in central p+Pb collisions at the LHC. The correlations functions are calculated in the viscous hydrodynamic model with Glauber Monte Carlo initial conditions. The correlation function is constructed from same-charged pion pairs. Three interferometry radii are extracted from an exponential fit to the correlation function. The lack of forwardbackward symmetry in the emission region leads to an additional cross-term coupling the q out and q long directions. The corresponding cross-term parameter R E mix is found to be nonzero.
The interferometry parameters R E out , R E side , R E long , and the cross term R mix are calculated as functions of the average transverse momentum k T and rapidity y ππ of the pion pair. The model reproduces the recent ATLAS Collaboration data to 10% accuracy for the k T dependence of the HBT parameters. We find that the HBT radii as functions of the pair rapidity y ππ are smaller on proton going side. The rapidity dependence of the cross term follows the slope of the rapidity distribution of emitted pions and is in good agreement with experimental data. The calculation confirms that the hydrodynamic model can reproduce qualitatively the space-time features of the emission source produced in high energy p+Pb collisions.