Two scales in Bose-Einstein correlations

We argue that the secondaries produced in high energy hadron collisions are emitted by small size sources distributed over a much larger area in impact parameter space occupied by the interaction amplitude. That is, Bose-Einstein correlation of two emitted identical particles should be described by a `two-radii' parametrization ansatz. We discuss the expected energy, charged multiplicity and transverse momentum of the pair (that is, $\sqrt{s},~N_{\rm ch}, k_t$) behaviour of both the small and large size components.


Introduction
An effective tool to study the space-time structure of the production amplitude is to measure the Bose-Einstein correlations (BEC) between two identical particles produced in the inclusive hadron interaction, see, for example, [1]− [4]. Consider the situation where we have one pion with momentum p 1 emitted at point r 1 and another identical pion with p 2 and r 2 . The inclusive cross section for the two identical particles takes the form where M is the production amplitude, and where 4-vectors Q = p 2 −p 1 and r = r 1 −r 2 . The ... denote the averaging over r 1 and r 2 , The e irQ term is due to the permutation of the identical pions; that is, it allows for the pion with p 2 to be emitted from the point r 1 and simultaneously for p 1 from r 2 . As a rule the Q dependence of the amplitude M is relatively flat in comparison with the Q dependence of e irQ . Thus we are able to evaluate the size of the pion production domain by studying the Q dependence of the whole cross section d 2 σ.
To extract the effect we compare the measured Q spectrum with a similar one but without BEC. To be precise we form the ratio where dN/dQ is the two pion distribution integrated over all the variables except Q, and dN ref /dQ is the distribution expected in a world without BEC. There are different ways to choose dN ref /dQ. We may measure the π + π − distribution for non-identical pions; or we may change the sign of the three momentum of the second pion p 2 → − p 2 ; and so on. None of these approaches compensates for the Q dependence of M completely; but for the conventional 'one-radius' fit the different values of the mean radius 1 ,r, extracted from the data are close to each other. Such analyses of high energy proton-proton interactions at the LHC have been performed by ATLAS [5], CMS [6] and ALICE [7].
The problem is that the value of λ turns out to be less than 1. In particular, CMS claim the λ = 0.62 ± 0.01 [6]. On the other hand from (2) we expect R(Q = 0) = 1. Moreover, it is clear from Fig. 1 of [6], and the analogous plots of other groups, that the fit does not describe the very low Q data points. This indicates that there should be another component of R(Q) with a larger radius populating the region of small Q.
In the present paper we argue that the expected structure of the pion emission domain is highly inhomogeneous. We should consider small size pion sources distributed in the much larger area of the proton-proton interaction. That is, we are led to parametrize R(Q) by two different mean radii. We explain the physical origin of this situation below.

Mechanisms for multiparticle production
It was shown long ago that to describe a high energy (say, proton-proton) interaction it is convenient to first select the subset of diagrams which provides the interaction across a large rapidity gap, and whose contribution does not decrease when the rapidity separation increases [8]. The resulting ladder-like set of diagrams forms Reggeon exchange. In terms of the hadronic degrees of freedom the corresponding subset of diagrams -called multiperipheral ladder diagrams ( Fig. 1(a)) -was studied first by Amati el al. [9]. In terms of QCD 2 they form the BFKL pomeron [12]. In a general purpose Monte Carlo these are the diagrams for the DGLAP evolution amplitude ( Fig. 1(b)). Recall that it is not sufficient to consider only one ladder. To describe multiparticle production we have to consider the possibility of the exchange of a few ladders (Reggeons), see Fig. 1

(c). In a Monte Carlo this is called the Multiple Interaction (MI) option.
Already at this stage we observe two quite different scales. The slope, B el , of elastic protonproton scattering is usually parametrized in the form where the constant B 0 is driven by the size of an incoming proton. On the other hand, the value of α P reflects the internal transverse size of the ladder. Based on the pre-LHC data, typical numbers are B 0 ∼ 10 GeV −2 and α P ∼ 0.25 [13].
Strictly speaking, (4) corresponds to one-ladder (Reggeon) exchange. At very high energies the speed of the shrinkage of the diffraction cone increases due to the stronger screening of the amplitude in the centre of the disc, at small impact parameters. Asymptotically, in the black disc limit, where the total cross section grows as ln 2 s, the slope B el (s) also increases as ln 2 s. It was shown [14] that indeed the high energy LHC data indicate the presence of a ln 2 s component in the elastic slope which is consistent with the growth of the total cross section. However if we consider an individual ladder then the effective value of α P in the ladder even decreases with energy due to the larger available k t space for the intermediate partons. It is known that α P → 0 in the BFKL case. Another example is the Monte Carlo description of multiparticle production. In order to tune the generator to the high energy data one has to introduce an infrared cutoff, k min t whose value grows as s 0.12 [15] -that is, the transverse size of the ladder decreases.
Thus the interaction of two high energy protons should be described by a diagram like Fig.1(c) in which the size of each individual ladder is rather small (as seen from the small value of α P ). Yet the separation between the ladders is of the order of the radius of the interaction amplitude which should be correlated with the total value of B el (s). With increasing energy we expect the transverse size of an individual ladder (measured in the central rapidity interval) will decrease (as indicated by the behaviour of k min t (s)). On the other hand the separation between the ladders is expected to increase, as indicated by the behaviour of B el (s).

Two components in Bose-Einstein correlations
Having the above picture in mind, we expect in BEC to observe a new object -a small-size pion source. In other words, BEC should be described by two different radii. One radius corresponds to the case when both pions are emitted from the same ladder -this will measure the size of an individual 'pomeron'. Since the pion is not a point-like object the radius will be smeared out by fluctuations in the process of the formation of the pions. The second radius will correspond to the pions being emitted from two different ladders -it is a measure of the separation between the ladders. Therefore we propose to fit the observed correlation R(Q) by a formula with two different mean radii 3 which better reflects the complicated structure of the pion emission domain. In the ideal case we expect the low multiplicity events to be produced via a diagram with only one 'cut' pomeron exchange ( Fig. 1(b)). In general there may be more pomerons in the whole amplitude, but only one ladder radiates the secondaries. As the multiplicity becomes larger the secondaries are mainly emitted from a few different ladders, and the probability to find the two identical pions originating from the same pomeron decreases. That is, we expect the relative contribution λ of the large component (described by, say,r 1 ) to increase with N ch , while on the other hand, the strength of the small size component (1 − λ) decreases. Unfortunately we cannot predict that λ → 0 as N ch → 0, and that λ → 1 for very large N ch . The situation is complicated by the strong fluctuations of the number of secondaries in each ladder (or pomeron). Recall that actually we do not measure the total charged multiplicity of an event, but rather the number of secondaries in a limited central rapidity interval (like |η| < 2.5 in the case of ATLAS and CMS). Then the multiplicity corresponding to one pomeron is relatively small (N ch ∼ 4) and the fluctuations strongly wash out the relation between the measured values and the number of cut pomerons. Moreover for very low N ch we may sample contributions from diffractive dissociation which have a qualitatively different structure. Nevertheless at large N ch the multipomeron contribution dominates; that is BEC are driven mainly by the component with the largest radius,r 1 . Indeed it is seen in Fig. 3(b) of [5] that, in the one radius fit, (3), the radius increases with multiplicity reaching saturation ofr 2fm for N ch > ∼ 50.
Such saturation was predicted in [16], where it was explained that the radius measured in a one-radius fit, is driven, not by the initial energy, but mainly by the number of cut pomerons, n P ∝ N ch . Indeed, the radius of the individual pomeron depends weakly on energy (α P const.), while the number of cut pomerons observed in the event (that is the separated pion sources) is proportional to N ch . So the probability to have two identical pions from two different sources increases with N ch . When n p = 1 (at low N ch ) we observe one pomeron and measure its radius. On the other hand, for large N ch we study the separation between the pomerons and the value ofr is saturated at the radius of the interaction amplituder ∝ B el (s). It was shown in [16] that B el ∼ 20 GeV −2 [17] corresponds in a 'one-radius' fit tor = 2.2fm, which is in good agreement with Fig. 3(b) of the ATLAS fit [5].
We emphasize that the separation between the pomerons (the pion sources) is not equal to the incoming hadron (proton) radius. First note, the radius corresponding to the interaction amplitude is larger. It increases with energy. Recall that in each successive step of the ladder in Fig. 1 the impact parameter changes by ∆b t . This leads to a diffusion in b t , which results in the second (or α p ) term in the equation for the elastic slope, (4). Next, the pion is not a point-like particle and its formation also occupies some volume. Finally, there may be an interaction between the final state secondaries.
Strictly speaking the picture that we describe above corresponds to the initial stage of the interaction, and does not allow for possible final state rescattering. If there are final state interactions (either in terms of a hadron gas or a quark-gluon plasma) then BEC will measure the radius given by the point of the last interaction -that is, the domain occupied by secondaries is extended up to the stage where the particle density becomes so low that further interaction is very unlikely, and the size of the domain at which the last interaction occurs will be seen in BEC. The probability of final state rescattering increases with particle density N ch . If this contribution were to dominate, then it is natural to expectr ∝ (N ch ) 1/3 . However, the data of Fig. 3(b) of [5] do not show any evidence of such a behaviour for large N ch .
So what are the expectations of the 'two-radii' fit? The contribution of small raduis dominates at low N ch and decreases with increasing N ch [16]. The value of the small radius,r 2 , is almost independent of energy -there is a small tendency to decrease due to a larger k t in the ladder (see, for example [18]). On the other hand, the strength of ther 1 component increases with N ch -the value of the radius,r 1 , correlated with B el (that is, the radius of interaction of the incoming protons) slowly increases with energy. At asymptotic energies we expect B el ∝ ln 2 s. A tendency already seen in LHC data [14]. Therefore for s → ∞ we expect r 1 ∝ c lns. However, the coefficient c is numerically quite small.
For very large N ch the value orr 1 may additionally increase due to the final state rescattering, as was discussed above. However it is not seen in the 7 TeV pp collision data.

The k t dependence of BEC
Here we discuss the dependence of the radii,r 1 ,r 2 , on the transverse momentum of the identical pair, k t = (p 1 + p 2 ) t /2, First, we make a trivial remark -for larger k t we have better space resolution in implementing the BEC method. At low k t it may be hard to distinguish between r 1 andr 2 components, since the radiir 1 andr 2 will be larger than that measured at large k t due to uncertainty principle smearing, and thus closer to each other.
Recall that the probability to produce two large k t pions from differents ladders (pomerons) decreases with increasing k t as the single particle inclusive cross section decreases steeply with k t . That is, two identical pions, each with large k t , should be produced from the same large k t jet. So with increasing k t of the pair we expect a large contribution from ther 2 component with the value ofr 2 decreasing, reaching saturation corresponding to the jet size. To be more precise, we mean the size of pion formation due to the hadronization of large k t jets. This tendency ofr to decrease with increasing k t was indeed observed in a 'one-radii' fit of LHC data, see Figs. 5,6 of [5]. However, the value of k t was not sufficiently large to see saturation in these plots.
Another effect at very large k t , which may give an important contribution to large N ch events, concerns the multiplicity of jets which increase as exp(−c √ lnE T ), see, for example, [19]. However the high E T jet cross section is too small to identify this effect in the present data.

Conclusions
We emphasize that the dynamics of high energy hadron interactions is based on subamplitudes of small transverse size which are distributed over the whole domain occupied by the full interaction amplitude. Thus, in BEC we have to observe a small size object corresponding to the emission of both pions (or both kaons etc.) from a single subamplitude and a larger radius caused by events where the pions are produced from different subamplitudes. At large N ch the relative contribution of the large radius component increase. Enlarging the k t of the identical pion pair we improve the space resolution of the BEC analyzer. This allows a better separation of the contributions from the small and large radii components. When k t becomes too large (say, k t > ∼ 1 GeV) the probability to produce such large k t pions from two different subamplitudes becomes small. In this case BEC measure the radius of the 'jet' which emits this high k t pair of identical pions. That is, we expect the radius to decrease with increasing k t , reaching saturation for k t > ∼ 1 GeV.