Centrality dependent Lévy-stable two-pion Bose-Einstein correlations at the PHENIX experiment
Sándor Lökös
nucl-ex hep-ex




§ INTRODUCTION


Intensity correlation measurements represent a widely used technique in high energy physics. The technique was discovered in astrophysics by R. Hanbury Brown and Q. R. Twiss in correlation measurements, that were performed in radio and optical astronomy to measure the angular diameters of stars <cit.>. Hanbury Brown and Twiss are considered the experimental founders of the HBT effect. Independently, the intensity correlations of identical pions were observed in proton-antiproton annihilation. These correlations were explained by G. Goldhaber, S. Goldhaber, W-Y. Lee and A. Pais <cit.> on the basis of the Bose-Einstein symmetrization of the wave-function of identical pion pairs. 

Because the Bose-Einstein correlations are related to the Fourier transform of the source, correlation measurements can provide an insight into the shape and the dynamics of the source. The shape of the source is usually assumed to be Gaussian universally across collision energy, particle type and centrality. PHENIX found evidence for and power-law tails in the source distribution for pions <cit.>. To understand the underlying processes better one has to go beyond the Gaussian approximation to a more general distribution: the Lévy distribution <cit.>.

Lévy-type of analyses were discussed before in Refs. <cit.>. Furthermore, a detailed analysis 0-30% Au+Au data at √(s_ NN)=200 GeV has been performed, where further details of Lévy femtoscopy are also given <cit.>. In this paper we present the Lévy-type of correlation functions measured by the PHENIX experiment at various centralities at 200 GeV.




§ PHENIX EXPERIMENT


The PHENIX experiment at RHIC was designed to study various type of particles produced in heavy ion collisions. A schematic beam view drawing of the detector system and a bird view photo can be seen in Figure <ref>. An exhaustive overview of the experiment can be found in Ref. <cit.>, here we restrict ourselves to mention the detectors that played role in the present analysis.



This analysis uses the beam-beam counters (BBC) for event characterization. Its two arms cover 2π in azimuth angle. The charge sum in both BBC arms is used to determine the centrality of the event. The tracking and the determination of the longitudinal position of the vertex were done by Drift Chamber (DC) and Pad Chambers (PC1). The charged pion identification was done by lead scintillator Electromagnetic Calorimeter (PbSc) as well as the high resolution time-of-flight detectors (TOF East and TOF West). We identified charged pions in 0.2 GeV/c ≤ p_ T≤ 0.85 GeV/c transverse momentum range.



§ TWO-PARTICLE CORRELATION FUNCTIONS AND THE LÉVY DISTRIBUTION



The detailed theoretical description of the basis of the Lévy-type of analysis can be found in Ref. <cit.>. Here we just give a short introduction to the topic. In general, the two-particle correlation function is defined as


    C_2(p_1,p_2)=N_2(p_1,p_2)/N_1(p_1)N_1(p_2),

where N_k is the k-particle invariant momentum distribution function. With the S(x,p) source function <cit.>, the following formula can be obtained for the two-particle correlation function


    C_2^(0)(q,K)=1+|S̃(q,K)|^2/|S̃(0,K)|^2,

where the four-momentum difference of the pair q=p_1-p_2 and the average momentum of the pair K=(p_1+p_2)/2 were introduced as new variables. The ^(0) index indicates that none of the final state effects were taken into account. Without Coulomb interaction in the final state the correlation function takes the value of 2 at zero relative momentum. However, the experimental limits do not allow us to resolve the relative momentum below 4-8 MeV/c. Hence it is usually observed that the experimentally measured correlation functions, when extrapolated to zero relative momentum, yield a value different from 2. This can be quantified in the following form:


    λ(K) = lim_q→ 0 C_2(q,K)-1 ≠ 1.

This λ(K) is called the intercept parameter or the strength of the correlation function. This observation can be interpreted easily with the core-halo model. In this model of the Bose-Einstein correlations <cit.> the source is treated as a composite one with two parts: a core part which is a hydrodynamically behaving, fireball-type source and a halo part, which contains long-lived resonances.


Several long-lived resonances with a decay width smaller than the resolvable momentum difference (e.g. η, η', K_S^0 or in an appropriate experiment ω) decay to pions, which contribute to the halo region. The pions coming from the core do not correlate with the pions coming from the decays of the resonances of the halo part, but both type of pions will be detected. Therefore the long-lived resonances could reduce the correlation strength in a momentum range given by the decay kinematics. So the correlation function of Eq. (<ref>) is modified into the


    C_2(q,K)=1 + (N_c(K)/N_c(K)+N_h(K))^2 |S̃(q,K)|^2/|S̃(0,K)|^2 = 1 + λ|S̃(q,K)|^2/|S̃(0,K)|^2

form. If there is no halo part and the source is chaotic and fully thermal then λ=1. If λ < 1 it is the sign of the presence of the decay mesons. For more detailed calculations, discussions and conclusions see Ref. <cit.>.


Now let us discuss the shape of the correlation functions. In an expandig system, the mean free path is increasing as the system expands. In such case a generalized form of the central limit theorem leads to the appearance of the Lévy distribution as source function <cit.>. The one-dimensional, symmetric Lévy distribution is defined by a Fourier transform:


    ℒ(α,R,r)=1/(2π)^3∫ d^3 q e^iqr e^-1/2 |qR|^α.

Here α is called the Lévy-index of stability or the shape parameter, while R is the Lévy scale parameter. In the case of α=2 the Gaussian case is restored, the case of α=1 corresponds to the Cauchy distribution. From the Lévy-type of source the following correlation function can be obtained:


    C_2^(0)(q,K)=1+λ(K) e^-(R(K)q)^α(K).

The effect of the parameters can be seen on Figure <ref>. Previous measurements of two-pion Bose-Einstein correlation functions in Au+Au collisions <cit.> indicate the presence of a long-range power-law-like component which could be understood in terms of Lévy distributions.

The α shape parameter or Lévy-index indicates the deviation from the Gaussian case. Moreover, it could be associated to one of the critical exponents namely to the critical exponent of the spatial correlation in 3 dimension <cit.>. Because of this possible relation between the two exponents, the measurements of α could help to locate the critical point on the phase diagram of the QCD if there is any <cit.>. The critical exponent of the spatial correlation could vary near to a critical point. If α could be associated to this critical exponent, it should show changes near to the critical point. This association of the two exponent gives us the motivation to preform precise HBT measurements in various systems with different sizes or created with different energies.

In a hot, dense matter like the sQGP the U_ A(1) symmetry might be restored. If it is, then the η' meson should have smaller mass than with broken U_ A(1) which means enhanced production of the η' meson. This particle could decay into five pions and because of the dynamics the pions will have p_ T≈ 150-200 MeV/c transverse momentum. These pions contribute to the halo part of the source which means that the correlation strength will be supressed by the effect around these p_ T range. So if we measured that the λ has a decreasing trend at lower p_ T range then it might be caused by the (partially) restored U_ A(1) axial symmetry. A detailed explanation can be found in Ref. <cit.> and a Lévy-type analysis of Au+Au collisions at 200 GeV in 0-30% centrality can be found in Ref. <cit.>.


There are other possible explanation for the non-constant behavior of the λ(m_ T). Partial coherence could effect the λ(m_ T) but the investigation of this phenomena requires 3-particle correlation measurements <cit.>. Model which take into account the compositeness of the pion and its interactions exhibits similar behavior of the intercept parameter <cit.>.


Because charged pions are measured, the Coulomb interaction in the final state has an effect that should be taken into account. We utilized a generalization of the Sinyukov method detailed in Ref. <cit.> similarly to Ref. <cit.> to handle this final state effect.



§ RESULTS OF THE CENTRALITY DEPENDENT ANALYSIS


We studied the Lévy-parameters in 10% wide centrality bins in 0-60% range and in 18 transverse mass (m_ T) bins at √(s_ NN) = 200 GeV collision energy in Au+Au system. We determine the m_ T and N_ part dependencies of the Lévy parameters. The transverse mass is defined with m_ T = √(m^2 + p_ T^2).



 §.§ The Lévy shape parameter α


First, we measured the α(m_ T) dependencies in the above mentioned centrality ranges. The results can be seen in Figure <ref>. On the left, it can be observed that the value of α does not depend strongly on m_ T, compared to the systematic uncertainties. Hence the m_ T average ⟨α⟩ is meaningful and describes well the exponent for each centrality bin, we could perform fix α = ⟨α⟩ fits.


If we study the centrality dependence of α(m_ T) points or the averages it turns out that there is a clear non-monotonic behavior as a function of N_ part (as shown in the right panel of Figure <ref>.). We observe that all measured α values are far from the Gaussian limit of α=2.





 §.§ The Lévy scale parameter R


In a Lévy case, the source has no second moment, hence the Lévy scale is not equivalent to the Gaussian radius or RMS value. Nevertheless, the Lévy scale exhibits similar features, for example a decreasing trend as a function of m_ T as it can be seen in Figure <ref>. in case of free and fix α. It turns out that the 1/R^2 ∝ a + b · m_ T kind of proportionality (deduced from hydrodynamical calculations <cit.>) remains valid in the Lévy case too. This is visualized in Figure <ref>., where we show 1/R^2 vs m_ T, also in case of free and fix α. The trends are smoother and clearer if the α parameter are fixed to the averaged value corresponding to the given centrality bin (for the averaged value, see Figure <ref>.). 







 §.§ The Lévy strength λ


As it was mentioned above the λ parameter could be related to the core-halo ratio. λ(m_ T) is shown in Figure <ref>. The fixed α fit gives clearer and smoother trends. We introduced the normalized λ function with a normalization factor λ_ max= ⟨λ⟩_(0.7-0.9)/^2. If we plot the λ/λ_ max as a function of m_ T it shows that the decreasing trends, or the characteristic of the “hole” does not depend strongly on the centrality.





 §.§ New scaling parameter


The empirically found new scaling parameter <cit.> was observed in the centrality dependent case as well. If one introduces the combination of the three Lévy parameters:


    R̂ = R/λ(1+α)

a surprising scaling property of this variable rises up in all centrality ranges and in the whole investigated m_ T range as it can be seen in Figure <ref>.

A very remarkable property of this scaling is that it is not sensitive to the α fixation, even though all the parameters above get clearer and smoother trends with fixed α=⟨α⟩ fits. In the case of this variable the α could be free or fixed, the trends are the same within statistical errors. The precise meaning of this variable is not clear yet and as we know there are not any theoretical prediction or explanation of it.



§ CONCLUSION


We measured correlation functions at 200 GeV at PHENIX experiment in 18 m_ T and 6 centrality range. We fit the measured data with Coulomb corrected  Lévy distribution, yield the Lévy parameters and determine their m_ T and centrality dependencies. We can find that the Lévy-stability index α has non-monotonic behavior as a function of centrality. We also observed that the α is far from the Gaussian case corresponds to α=2. The Lévy scale parameter R has hydrodynamical trends in all observed centrality bin as a function of m_ T. The strength of the correlation function λ(m_ T) decreases at lower m_ T. We also observed that the normalized λ / λ_ max function has no strong centrality dependence. Further detailed and dedicated analysis is needed to clarify these observed phenomena.

The author is thankful for the support of EFOP-3.6.1-16-2016-00001 and NKFIH grant FK-123842

——-


[Brown and Twiss(1956)]HBT_orig
Brown, R.H.; Twiss, R.Q.
Correlation between Photons in two Coherent Beams of Light.
Nature 1956, 177, 27–29.

[Goldhaber et al.(1960)Goldhaber, Goldhaber, Lee, and
  Pais]PhysRev.120.300
Goldhaber, G.; Goldhaber, S.; Lee, W.; Pais, A.
Influence of Bose-Einstein Statistics on the Antiproton-Proton
  Annihilation Process.
Phys. Rev. 1960, 120, 300–312.

[Adler et al.(2007)Adler et al.]Adler:2006as
Adler, S.S.; others.
Evidence for a long-range component in the pion emission source in
  Au + Au collisions at s(NN)**(1/2) = 200-GeV.
Phys. Rev. Lett. 2007, 98, 132301,
  http://xxx.lanl.gov/abs/nucl-ex/0605032[arXiv:nucl-ex/nucl-ex/0605032].

[Csörgő et al.(2004)Csörgő, Hegyi, and
  Zajc]Csorgo:2003uv
Csörgő, T.; Hegyi, S.; Zajc, W.A.
Bose-Einstein correlations for Levy stable source distributions.
Eur. Phys. J. 2004, C36, 67–78,
  http://xxx.lanl.gov/abs/nucl-th/0310042[arXiv:nucl-th/nucl-th/0310042].

[Metzler et al.(1999)Metzler, Barkai, and Klafter]Metzler:1999zz
Metzler, R.; Barkai, E.; Klafter, J.
Anomalous Diffusion and Relaxation Close to Thermal Equilibrium: A
  Fractional Fokker-Planck Equation Approach.
Phys. Rev. Lett. 1999, 82, 3563–3567.

[Csanád et al.(2007)Csanád, Csörgő, and
  Nagy]Csanad:2007fr
Csanád, M.; Csörgő, T.; Nagy, M.
Anomalous diffusion of pions at RHIC.
Braz. J. Phys. 2007, 37, 1002–1013,
  http://xxx.lanl.gov/abs/hep-ph/0702032[arXiv:hep-ph/hep-ph/0702032].

[Csanád(2006)]Csanad:2005nr
Csanád, M.
Measurement and analysis of two- and three-particle correlations.
Nucl. Phys. 2006, A774, 611–614,
  http://xxx.lanl.gov/abs/nucl-ex/0509042[arXiv:nucl-ex/nucl-ex/0509042].

[Kincses(2016)]Kincses:2016jsr
Kincses, D.
PHENIX results on Lévy analysis of Bose-Einstein correlation
  functions.
10th International Workshop on Critical Point and Onset of
  Deconfinement (CPOD 2016) Wrocław, Poland, May 30-June 4, 2016,  2016,
  http://xxx.lanl.gov/abs/1610.05025[arXiv:nucl-ex/1610.05025].

[Adare et al.(2017)Adare et al.]Adare:2017vig
Adare, A.; others.
Lévy-stable two-pion Bose-Einstein correlations in
  √(s__NN)=200 GeV Au+Au collisions 2017.
http://xxx.lanl.gov/abs/1709.05649[arXiv:nucl-ex/1709.05649].

[Csanád()]MateBGL
Csanád, M.
Lévy femtoscopy with PHENIX at RHIC.
http://xxx.lanl.gov/abs/1711.05575[1711.05575].

[Adcox et al.(2003)Adcox et al.]Adcox:2003zm
Adcox, K.; others.
PHENIX detector overview.
Nucl. Instrum. Meth. 2003, A499, 469–479.

[Yano and Koonin(1978)]Yano:1978gk
Yano, F.B.; Koonin, S.E.
Determining Pion Source Parameters in Relativistic Heavy Ion
  Collisions.
Phys. Lett. 1978, 78B, 556–559.

[Bolz et al.(1993)Bolz, Ornik, Plumer, Schlei, and
  Weiner]Bolz:1992hc
Bolz, J.; Ornik, U.; Plumer, M.; Schlei, B.R.; Weiner, R.M.
Resonance decays and partial coherence in Bose-Einstein
  correlations.
Phys. Rev. 1993, D47, 3860–3870.

[Csörgő et al.(1996)Csörgő, Lörstad, and
  Zimányi]Csorgo:1994in
Csörgő, T.; Lörstad, B.; Zimányi, J.
Bose-Einstein correlations for systems with large halo.
Z. Phys. 1996, C71, 491–497,
  http://xxx.lanl.gov/abs/hep-ph/9411307[arXiv:hep-ph/hep-ph/9411307].

[Csörgő(2002)]Csorgo:1999sj
Csörgő, T.
Particle interferometry from 40-MeV to 40-TeV.
Heavy Ion Phys. 2002, 15, 1–80,
  http://xxx.lanl.gov/abs/hep-ph/0001233[arXiv:hep-ph/hep-ph/0001233].

[Vértesi et al.(2011)Vértesi, Csörgő, and
  Sziklai]Vertesi:2009wf
Vértesi, R.; Csörgő, T.; Sziklai, J.
Significant in-medium η' mass reduction in √(s_NN)=200
  GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider.
Phys. Rev. 2011, C83, 054903,
  http://xxx.lanl.gov/abs/0912.0258[arXiv:nucl-ex/0912.0258].

[Csörgő et al.(2011)Csörgő, Vértesi, and
  Sziklai]Csorgo:2010hj
Csörgő, T.; Vértesi, R.; Sziklai, J.
U_A(1) Symmetry Restoration from an In-Medium η ' Mass
  Reduction in √(s_NN)=200 GeV Au+Au Collisions.
Quantum chromodynamics and beyond: Gribov-80 memorial volume.
  Proceedings, Memorial Workshop devoted to the 80th birthday of V.N. Gribov,
  Trieste, Italy, May 26-28, 2010,  2011, pp. 307–318,
  http://xxx.lanl.gov/abs/1012.5058[arXiv:nucl-ex/1012.5058].

[Csörgő et al.(2006)Csörgő, Hegyi, Novák,
  and Zajc]Csorgo:2005it
Csörgő, T.; Hegyi, S.; Novák, T.; Zajc, W.A.
Bose-Einstein or HBT correlation signature of a second order QCD
  phase transition.
AIP Conf. Proc. 2006, 828, 525–532,
  http://xxx.lanl.gov/abs/nucl-th/0512060[arXiv:nucl-th/nucl-th/0512060].

[Gavrilik and Mishchenko(2015)]Gavrilik:2014pxa
Gavrilik, A.M.; Mishchenko, Yu.A.
Correlation function intercepts for μ̃,q-deformed Bose
  gas model implying effective accounting for interaction and compositeness of
  particles.
Nucl. Phys. 2015, B891, 466–481,
  http://xxx.lanl.gov/abs/1411.5955[arXiv:hep-ph/1411.5955].

[Sinyukov et al.(1998)Sinyukov, Lednicky, Akkelin, Pluta, and
  Erazmus]Sinyukov:1998fc
Sinyukov, Yu.; Lednicky, R.; Akkelin, S.V.; Pluta, J.; Erazmus, B.
Coulomb corrections for interferometry analysis of expanding hadron
  systems.
Phys. Lett. 1998, B432, 248–257.

[Makhlin and Sinyukov(1988)]Makhlin1988
Makhlin, A.N.; Sinyukov, Y.M.
The hydrodynamics of hadron matter under a pion interferometric
  microscope.
Zeitschrift für Physik C Particles and Fields 1988,
  39, 69–73.

[Csörgő and Lörstad(1996)]Csorgo:1995bi
Csörgő, T.; Lörstad, B.
Bose-Einstein correlations for three-dimensionally expanding,
  cylindrically symmetric, finite systems.
Phys. Rev. 1996, C54, 1390–1403,
  http://xxx.lanl.gov/abs/hep-ph/9509213[arXiv:hep-ph/hep-ph/9509213].

[Chapman et al.(1995)Chapman, Scotto, and Heinz]Chapman:1994yv
Chapman, S.; Scotto, P.; Heinz, U.W.
A New cross term in the two particle HBT correlation function.
Phys. Rev. Lett. 1995, 74, 4400–4403,
  http://xxx.lanl.gov/abs/hep-ph/9408207[arXiv:hep-ph/hep-ph/9408207].





