Next Article in Journal
A Phenomenological Equation of State of Strongly Interacting Matter with First-Order Phase Transitions and Critical Points
Next Article in Special Issue
Geometric Aspects and Some Uses of Deformed Models of Thermostatistics
Previous Article in Journal
Vector-Interaction-Enhanced Bag Model
Previous Article in Special Issue
Exclusive Particle Production in pp and pPb Collisions at CMS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Centrality Dependent Lévy-Stable Two-Pion Bose-Einstein Correlations in \( {\sqrt{s_{NN}}} \) = 200 GeV Au+Au Collisions at the PHENIX Experiment

by
Sándor Lökös
1,2,†
1
Department of Physics, Eszterházy Károly University, Gyöngyös Mátrai út 36, H-3200 Eger, Hungary
2
Department of Atomic Physics, Eötvös University, Pázmány P. s. 1/A, H-1117 Budapest, Hungary
The author belongs to PHENIX Collaboration.
Universe 2018, 4(2), 31; https://doi.org/10.3390/universe4020031
Submission received: 4 January 2018 / Revised: 26 January 2018 / Accepted: 7 February 2018 / Published: 9 February 2018

Abstract

:
Investigation of femtoscopic correlation functions in relativistic heavy ion reactions is an important tool to access the space-time structure of particle production in the strongly interacting Quark Gluon Plasma (sQGP). The shape of the source, and thus the shape of the correlation functions, is often assumed to be Gaussian, but experimental results found evidence for heavy tails in the source distribution of pions. Recent analysis revealed that the statistically correct assumption could be the so-called Lévy distribution. The detailed investigation of correlation functions in various systems may shed light on the location of the critical endpoint on QCD (Quantum Chromodynamics) phase diagram. It could also reveal if there is partially coherent pion production or could indicate the possible in-medium mass modification of the η meson due to the (partial) restoration of the U A ( 1 ) axial symmetry. These phenomena could depend on the system size and on the collision energy. A detailed centrality-dependent analysis could explore the multiplicity dependencies of the Lévy parameters, and thus the critical and thermodynamical properties of the sQGP, and could give information about the above mentioned processes. In this paper, we present the status of the centrality dependent measurements of two-pion Lévy Bose-Einstein correlation functions s NN = 200 GeV Au+Au collisions at PHENIX.

1. 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 [1]. 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 [2] 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 power-law tails in the source distribution for pions [3]. To understand the underlying processes better one has to go beyond the Gaussian approximation to a more general distribution: the Lévy distribution [4,5,6].
Lévy-type of analyses were discussed before in Refs. [7,8]. 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 [9,10]. In this paper we present the Lévy-type of correlation functions measured by the PHENIX experiment at various centralities at s NN = 200 GeV.

2. 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’s-eye view photo can be seen in Figure 1. An exhaustive overview of the experiment can be found in Ref. [11], 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.

3. 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. [9,10]. 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 [12], 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 [13,14] 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 Equation (3) 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. [15,16,17].
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 [4,5,6]. The one-dimensional, symmetric Lévy distribution is defined by a Fourier transform:
L ( α , R , r ) = 1 ( 2 π ) 3 d 3 q e i qr e 1 2 | q R | α .
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 2. Previous measurements of two-pion Bose-Einstein correlation functions in Au+Au collisions [7] 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 [4]. 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 [4,6,18]. 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 size s 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. [16] and a Lévy-type analysis of Au+Au collisions at 200 GeV in 0–30% centrality can be found in Ref. [9].
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 [15]. Model which take into account the compositeness of the pion and its interactions exhibits similar behavior of the intercept parameter [19].
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. [20] similarly to Ref. [9] to handle this final state effect.

4. 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 .

4.1. 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 3. 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 3). We observe that all measured α values are far from the Gaussian limit of α = 2 .

4.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 4. 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 [21,22,23]) remains valid in the Lévy case too. This is visualized in Figure 5, 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 3).

4.3. 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 6. The fixed α fit gives clearer and smoother trends. We introduced the normalized λ function with a normalization factor λ max = λ ( 0.7 0.9 ) GeV / c 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.

4.4. New Scaling Parameter

The empirically found new scaling parameter [9,10] 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 7.
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.

5. Conclusions

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.

Acknowledgments

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

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Brown, R.H.; Twiss, R.Q. Correlation between Photons in two Coherent Beams of Light. Nature 1956, 177, 27–29. [Google Scholar] [CrossRef]
  2. 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. [Google Scholar] [CrossRef]
  3. Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Amirikas, R.; Aphecetche, L.; Aronson, S.H.; Averbeck, R.; et al. Evidence for a long-range component in the pion emission source in Au + Au collisions at s N N = 200 GeV. Phys. Rev. Lett. 2007, 98, 132301. [Google Scholar] [CrossRef] [PubMed]
  4. Csörgő, T.; Hegyi, S.; Zajc, W.A. Bose-Einstein correlations for Levy stable source distributions. Eur. Phys. J. C 2004, 36, 67–78. [Google Scholar] [CrossRef]
  5. 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. [Google Scholar] [CrossRef]
  6. Csanád, M.; Csörgő, T.; Nagy, M. Anomalous diffusion of pions at RHIC. Braz. J. Phys. 2007, 37, 1002–1013. [Google Scholar] [CrossRef]
  7. Csanád, M. Measurement and analysis of two- and three-particle correlations. Nucl. Phys. A 2006, 774, 611–614. [Google Scholar] [CrossRef]
  8. Kincses, D. PHENIX results on Lévy analysis of Bose-Einstein correlation functions. In Proceedings of the 10th International Workshop on Critical Point and Onset of Deconfinement (CPOD 2016), Wrocław, Poland, 30 May–4 June 2016. [Google Scholar]
  9. Adare, A.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Alexander, J.; Apadula, N.; Alfred, M.; Al-Ta’ani, H.; Angerami, A.; et al. Lévy-stable two-pion Bose-Einstein correlations in s N N = 200 GeV Au+Au collisions. arXiv, 2017; arXiv:1709.05649. [Google Scholar]
  10. Csanád, M. Lévy femtoscopy with PHENIX at RHIC. Universe 2017, 3, 85. [Google Scholar] [CrossRef]
  11. Adcox, K.; Bennett, M.J.; Bobrek, M.; Boissevain, J.B.; Boose, S.; Bosze, E.; Conway, R.; Aphecetche, L.; Aronson, S.H.; Zhou, S.; et al. PHENIX detector overview. Nucl. Instrum. Meth. A 2003, 499, 469–479. [Google Scholar] [CrossRef]
  12. Yano, F.B.; Koonin, S.E. Determining Pion Source Parameters in Relativistic Heavy Ion Collisions. Phys. Lett. B 1978, 78, 556–559. [Google Scholar] [CrossRef]
  13. Bolz, J.; Ornik, U.; Plumer, M.; Schlei, B.R.; Weiner, R.M. Resonance decays and partial coherence in Bose-Einstein correlations. Phys. Rev. D 1993, 47, 3860–3870. [Google Scholar] [CrossRef]
  14. Csörgő, T.; Lorstad, B.; Zimanyi, J. Bose-Einstein correlations for systems with large halo. Z. Phys. C Part. Fields 1996, 71, 491–497. [Google Scholar] [CrossRef] [Green Version]
  15. Csörgő, T. Particle interferometry from 40-MeV to 40-TeV. Heavy Ion Phys. 2002, 15, 1–80. [Google Scholar] [CrossRef]
  16. Vértesi, R.; Csörgő, T.; Sziklai, J. Significant in-medium η’ mass reduction in s N N = 200 GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider. Phys. Rev. C 2011, 83, 054903. [Google Scholar] [CrossRef]
  17. Csörgő, T.; Vértesi, R.; Sziklai, J. UA(1) Symmetry Restoration from an In-Medium η ’ Mass Reduction in s N N = 200 GeV Au + Au Collisions, Quantum chromodynamics and beyond: Gribov-80 memorial volume. In Proceedings of the Memorial Workshop Devoted to the 80th Birthday of V N Gribov, Trieste, Italy, 26–28 May 2010; pp. 307–318. [Google Scholar]
  18. 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. [Google Scholar]
  19. 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. B 2015, 891, 466–481. [Google Scholar] [CrossRef]
  20. Sinyukov, Yu.; Lednicky, R.; Akkelin, S.V.; Pluta, J.; Erazmus, B. Coulomb corrections for interferometry analysis of expanding hadron systems. Phys. Lett. B 1998, 432, 248–257. [Google Scholar] [CrossRef]
  21. Makhlin, A.N.; Sinyukov, Y.M. The hydrodynamics of hadron matter under a pion interferometric microscope. Z. Phys. C Part. Fields 1988, 39, 69–73. [Google Scholar] [CrossRef]
  22. Csörgő, T.; Lörstad, B. Bose-Einstein correlations for three-dimensionally expanding, cylindrically symmetric, finite systems. Phys. Rev. C 1996, 54, 1390–1403. [Google Scholar] [CrossRef]
  23. 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. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A schematic beam view drawing of the detector system of the accelerator complex at the PHENIX experiment and the RHIC facilities from a bird’s-eye view.
Figure 1. A schematic beam view drawing of the detector system of the accelerator complex at the PHENIX experiment and the RHIC facilities from a bird’s-eye view.
Universe 04 00031 g001
Figure 2. Comparison of the Gaussian and the Lévy distribution with a given α , R and λ . It can be seen that the Lévy distribution has a power-law-like tail.
Figure 2. Comparison of the Gaussian and the Lévy distribution with a given α , R and λ . It can be seen that the Lévy distribution has a power-law-like tail.
Universe 04 00031 g002
Figure 3. The α Lévy stability parameter as a function of m T (left); The α as a function of N part (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Figure 3. The α Lévy stability parameter as a function of m T (left); The α as a function of N part (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Universe 04 00031 g003
Figure 4. The Lévy scale parameter as a function of m T with free α (left) and in the case when α = α is fixed during the fit (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Figure 4. The Lévy scale parameter as a function of m T with free α (left) and in the case when α = α is fixed during the fit (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Universe 04 00031 g004
Figure 5. The hydro-type behavior of the Lévy scale parameter as a function of m T with free α (left) and in the case when α = α is fixed during the fit (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Figure 5. The hydro-type behavior of the Lévy scale parameter as a function of m T with free α (left) and in the case when α = α is fixed during the fit (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Universe 04 00031 g005
Figure 6. The λ ( m T ) function in the different centrality ranges (left); The normalized λ ( m T ) function with fixed α = α value (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Figure 6. The λ ( m T ) function in the different centrality ranges (left); The normalized λ ( m T ) function with fixed α = α value (right). The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Universe 04 00031 g006
Figure 7. The new, empirically found scaling variable does not depend on the α fixation. The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Figure 7. The new, empirically found scaling variable does not depend on the α fixation. The statistical errors are shown with errorbars and the systematical uncertainties boxes.
Universe 04 00031 g007

Share and Cite

MDPI and ACS Style

Lökös, S. Centrality Dependent Lévy-Stable Two-Pion Bose-Einstein Correlations in \( {\sqrt{s_{NN}}} \) = 200 GeV Au+Au Collisions at the PHENIX Experiment. Universe 2018, 4, 31. https://doi.org/10.3390/universe4020031

AMA Style

Lökös S. Centrality Dependent Lévy-Stable Two-Pion Bose-Einstein Correlations in \( {\sqrt{s_{NN}}} \) = 200 GeV Au+Au Collisions at the PHENIX Experiment. Universe. 2018; 4(2):31. https://doi.org/10.3390/universe4020031

Chicago/Turabian Style

Lökös, Sándor. 2018. "Centrality Dependent Lévy-Stable Two-Pion Bose-Einstein Correlations in \( {\sqrt{s_{NN}}} \) = 200 GeV Au+Au Collisions at the PHENIX Experiment" Universe 4, no. 2: 31. https://doi.org/10.3390/universe4020031

APA Style

Lökös, S. (2018). Centrality Dependent Lévy-Stable Two-Pion Bose-Einstein Correlations in \( {\sqrt{s_{NN}}} \) = 200 GeV Au+Au Collisions at the PHENIX Experiment. Universe, 4(2), 31. https://doi.org/10.3390/universe4020031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop