Next Article in Journal / Special Issue
Pure-Connection Gravity and Anisotropic Singularities
Previous Article in Journal / Special Issue
Entropic Distance for Nonlinear Master Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lévy Analysis of HBT Correlation Functions in s N N = 62 GeV and 39 GeV Au + Au Collisions at PHENIX †

Faculty of Atomic Physics, Eötvös Loránd University Pázmány P. s. 1/A, H-1117 Budapest, Hungary
This paper is an extended version of our paper published in 10th Bolyai-Gauss-Lobachevsky Conference on Non-Euclidean Geometry and its Applications.
Universe 2018, 4(1), 11; https://doi.org/10.3390/universe4010011
Submission received: 17 November 2017 / Revised: 13 December 2017 / Accepted: 19 December 2017 / Published: 4 January 2018

Abstract

:
The phase diagram of strongly interacting matter can be explored by analyzing data of heavy-ion collisions at different center of mass collision energies. For investigating the space-time structure of the hadron emission source, Bose-Einstein or HBT correlation measurements are among the best tools. In this paper we present the latest results from the PHENIX experiment of the Relativistic Heavy Ion Collider (RHIC) on such measurements in s N N = 39 GeV and 62 GeV Au + Au collisions.

1. Introduction

Experimental data from heavy-ion collisions [1] are in agreement with the theoretical prediction [2] that the quark-hadron transition is cross-over at high energies (near zero baryochemical potential). At lower energies the transition is expected to be a first-order phase transition. This indicates that there could be a critical endpoint in between, where second-order phase transition takes place. Finding and characterizing this critical point is one of the most important current problems in heavy-ion physics. To study the phase diagram of strongly interacting matter, we have to perform measurements at different collision energies so that we can gain information about different regions on the phase diagram. In the past years the PHENIX experiment at the Relativistic Heavy Ion Collider of the Brookhaven National Laboratory collected a vast amount of data at various collision energies. The gold-gold beam energy scan spans from 7.7 GeV up to 200 GeV in s N N . In terms of baryochemical potential the range reaches out to approximately 400 MeV, while in freeze-out temperature this range is between 170 MeV and 140 MeV [3]. An illustration of the phase-diagram can be seen on Figure 1, while the recently extracted chemical freezeout parameters are shown on Figure 2.
There is a very practical tool to gain information about the particle-emitting source, namely the measurement of Bose-Einstein or HBT correlations (named after R. Hanbury Brown and R. Q. Twiss) of identical bosons [4,5]. In previous years, Gaussian-type sources were usually utilized to describe the correlation functions, but the latest results from the PHENIX experiment showed that we can and have to go beyond the Gaussian approximation [6,7]. In an expanding medium, on the basis of the generalized central limit theorem and the anomalous diffusion we can expect the appearance of Lévy-type sources [8,9]. The symmetric Lévy distribution is the generalization of the Gaussian distribution and it is defined by the following expression:
L ( α , R , r ) = 1 ( 2 π ) 3 d 3 q e i q r e 1 2 | q R | α .
The R parameter is called the Lévy scale parameter, and the α parameter is the index of stability or Lévy exponent. In case of α = 2 we get back the Gaussian distribution, while the α = 1 case results in the Cauchy distribution. In case of α < 2 the Lévy distribution exhibits a power-law behavior. Neglecting the final-state interactions and assuming that the source is a spherically symmetric three-dimensional Lévy distribution, furthermore using the framework of the core-halo model, the two-particle correlation function takes the following simple form (Q is a one-dimensional relative momentum variable, see details in [9] and in Section 2):
C 2 ( Q ) = 1 + λ · e ( Q R ) α .
The λ parameter in Equation (2) is called the strength of the correlation function. According to the core-halo model [10], this parameter is in connection with the core-halo ratio.
λ = ( N C / ( N C + N H ) ) 2 ,
where N C denotes the number of pions created from the freeze-out of the quark-gluon plasma (or from the decay of very short lived resonances), and N H is the number of pions created in the decays of long-lived resonances.
The Lévy-type source is also useful in the search for the critical point, as one of its parameters, the index of stability α is related to one of the critical exponents, η . In case of a second order phase transition the η exponent describes the power-law behavior of the spatial correlations at the critical point with an exponent of ( d 2 + η ) where d is the dimension [11]. In case of the three-dimensional Lévy distribution the power-law tail of the spatial correlations is described with an exponent of 1 α . From comparing these exponents we can easily see that the Lévy exponent is identical to η at the critical point [12]. The second-order QCD phase-transition is expected to be in the same universality class as the 3D Ising model [13,14], and in that case the η exponent has a definite value of 0.03631(3) at the critical point [15]. Some argue [16] that QCD more likely falls in the universality class of the random field 3D Ising model, where the value of η is 0.5 ± 0.05 [11]. If one extracts the value of the Lévy exponent α at different center-of-mass collision energies, the data may yield information on the nature of the quark-hadron phase transition, particularly it may shed light on the location of the critical endpoint (CEP) on the phase-diagram.

2. Results

In our analysis we used s N N = 62 GeV and s N N = 39 GeV Au + Au data recorded by PHENIX during the 2010 running period. We measured two-particle HBT correlation functions for π + π + and π π pairs. For the experimental construction of the correlation functions we used the usual event mixing method as described in Section III.B. in Reference [7]. It is well known that the in the longitudinal co-moving system (LCMS) the Bertsch-Pratt radii ( R o u t , R s i d e , R l o n g ) are of similar magnitude in s N N = 200 GeV Au + Au collisions at RHIC [17]. This means that the Bose-Einstein correlation functions are approximately spherically symmetric in the LCMS frame. This was one of the reasons why we used the absolute value of the three-momentum difference in the LCMS frame for our one-dimensional relative momentum variable. In case of the 62 GeV dataset we used 8 average transverse mass ( m T ) bins and 4 centrality ranges, while in case of the 39 GeV dataset we used 6 m T bins and 2 centrality ranges. We checked that the correlation functions of π + π + and π π pairs are consistent, so we combined the data. The aforementioned formula for the correlation function (Equation (2)) does not take into account the final-state Coulomb interaction. In our fits we incorporated that effect and also introduced a linear background. In all m T bins the fits described the measured correlation functions in a statistically acceptable manner.
In Figure 3 we present the resulting fit parameters (R, λ , α ) versus pair m T , as well as an empirically found scaling parameter R ^ = R / ( λ · ( 1 + α ) ) . We observe that most features of the 200 GeV results [18] are resembled by these low energy data. Just like in the 200 GeV case, we can observe a decrease in λ values at small m T at both 62 and 39 GeV which may point to the increase of the halo fraction in that region. In case of a partial chiral U A ( 1 ) restoration, the in-medium mass of the η meson could be modified [19] which then could lead to an increase in its yield. This effect would result in an increased number of halo pions in the small- m T region, and hence it could provide one explanation for our observation [20]. Of course there could be other effects [4] that would result in similar behavior of λ so a precise measurement is important.
Moving on to the Lévy scale parameter we can see a decreasing behavior, and we can clearly observe a geometrical centrality dependence. In more central collisions, the R parameter (which is related to the physical size of the system) takes larger values. In case of the Lévy exponent α we can see that for all centrality and m T ranges the values are between 1 and 2 for both 39 and 62 GeV. The measured values are far from the Gaussian α = 2 case and also far from the conjectured α 0.5 at the critical point. Although we have not seen an indication for critical behavior in neither the 62 GeV nor the 39 GeV measurements, we can not exclude these ranges from the critical point search because finite-size and finite-time effects can modify the picture immensely as discussed in Reference [21].
We also show an empirically found scaling parameter, R ^ . We observed that the Levy parameters ( α , R, λ ) are strongly correlated and we tried to search for less correlated combinations of these parameters. We indeed found such a parameter, R ^ . If this parameter is used as a fit parameter instead of R, the correlation coefficients for ( R ^ , λ ) and ( R ^ , α ) are reduced substantially. Due to the reduction of the correlation the uncertainty of R ^ is also significantly reduced compared to that of R. What is even more interesting is the linear scaling of 1 / R ^ with m T . The possible physical interpretation of this remarkable affine linear dependence of 1 / R ^ on m T is still unknown at the time of writing this paper.
In the near future we would like to explore further the energy and centrality dependence in s N N = 27 GeV, 19 GeV and 15 GeV Au + Au collisions. We will also conduct a three-dimensional analysis to investigate the shape of the particle-emitting source in more detail.

Acknowledgments

This research was supported by the funding agencies listed in Reference [7]. In addition, D.K. was supported by the ÚNKP-17-3 New National Excellence program of the Hungarian Ministry of Human Capacities as well as the NKFIH grant FK-123842.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Csanád, M.; Csörgő, T.; Lörstad, B.; Ster, A. Indication of quark deconfinement and evidence for a Hubble flow in 130-GeV and 200-GeV Au + Au collisions. J. Phys. G 2004, 30, S1079–S1082. [Google Scholar] [CrossRef]
  2. Aoki, Y.; Endrődi, G.; Fodor, Z.; Katz, S.D.; Szabó, K.K. The order of the quantum chromodynamics transition predicted by the standard model of particle physics. Nature 2006, 443, 675–678. [Google Scholar] [CrossRef] [PubMed]
  3. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Ajitanand, N.; Alekseev, I.; Anderson, D.M.; Aoyama, R.; Aparin, A.; et al. Bulk Properties of the Medium Produced in Relativistic Heavy-Ion Collisions from the Beam Energy Scan Program. Phys. Rev. C 2017, 96, 044904. [Google Scholar] [CrossRef]
  4. Csörgő, T. Particle interferometry from 40-MeV to 40-TeV. Heavy Ion Phys. 2002, 15, 1–80. [Google Scholar] [CrossRef]
  5. Lisa, M.A.; Pratt, S.; Soltz, R.; Wiedemann, U. Femtoscopy in relativistic heavy ion collisions. Ann. Rev. Nucl. Part. Sci. 2005, 55, 357–402. [Google Scholar] [CrossRef]
  6. Csanad, M. Measurement and analysis of two- and three-particle correlations. Nucl. Phys. A 2006, 774, 611–614. [Google Scholar] [CrossRef]
  7. Adare, A.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Alexander, J.; Alfred, M.; Al-Ta’ani, H.; Angerami, A.; Aoki, K.; et al. Lévy-stable two-pion Bose-Einstein correlations in s N N = 200 GeV Au + Au collisions. arXiv, 2017; arXiv:nucl-ex/1709.05649. [Google Scholar]
  8. Csanád, M.; Csörgő, T.; Nagy, M. Anomalous diffusion of pions at RHIC. Braz. J. Phys. 2007, 37, 1002–1013. [Google Scholar] [CrossRef]
  9. 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]
  10. 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]
  11. Rieger, H. Critical behavior of the three-dimensional random-field Ising model: Two-exponent scaling and discontinuous transition. Phys. Rev. B 1995, 52, 6659–6667. [Google Scholar] [CrossRef]
  12. Csörgő, T. Correlation Probes of a QCD Critical Point. In Proceedings of the PoS (HIGH-PTLHC 08), Tokaj, Hungary, 16–19 March 2008. [Google Scholar]
  13. Halasz, A.M.; Jackson, A.D.; Shrock, R.E.; Stephanov, M.A.; Verbaarschot, J.J.M. On the phase diagram of QCD. Phys. Rev. D 1998, 58, 096007. [Google Scholar] [CrossRef]
  14. Stephanov, M.A.; Rajagopal, K.; Shuryak, E.V. Signatures of the tricritical point in QCD. Phys. Rev. Lett. 1998, 81, 4816–4819. [Google Scholar] [CrossRef]
  15. El-Showk, S.; Paulos, M.F.; Poland, D.; Rychkov, S.; Simmons-Duffin, D.; Vichi, A. Solving the 3d Ising Model with the Conformal Bootstrap II. c-Minimization and Precise Critical Exponents. J. Stat. Phys. 2014, 157. [Google Scholar] [CrossRef]
  16. 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]
  17. Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N. Bose-Einstein correlations of charged pion pairs in Au + Au collisions at s(NN)**(1/2) = 200-GeV. Phys. Rev. Lett. 2004, 93, 152302. [Google Scholar] [CrossRef] [PubMed]
  18. Csanad, M. Lévy femtoscopy with PHENIX at RHIC. arXiv, 2017; arXiv:nucl-ex/1711.05575. [Google Scholar] [CrossRef]
  19. Kapusta, J.I.; Kharzeev, D.; McLerran, L.D. The Return of the prodigal Goldstone boson. Phys. Rev. D 1996, 53, 5028–5033. [Google Scholar] [CrossRef]
  20. Vance, S.E.; Csörgő, T.; Kharzeev, D. Partial U(A)(1) restoration from Bose-Einstein correlations. Phys. Rev. Lett. 1998, 81, 2205–2208. [Google Scholar] [CrossRef]
  21. Lacey, R.A. Indications for a critical point in the phase diagram for hot and dense nuclear matter. Nucl. Phys. A 2016, 956, 348–351. [Google Scholar] [CrossRef]
Figure 1. The phase diagram of strongly interacting matter.
Figure 1. The phase diagram of strongly interacting matter.
Universe 04 00011 g001
Figure 2. Extracted chemical freezeout parameters from Reference [3].
Figure 2. Extracted chemical freezeout parameters from Reference [3].
Universe 04 00011 g002
Figure 3. Centrality and m T dependence of Lévy source parameters ( R , λ , α , R ^ ) in s N N = 62 GeV (top) and s N N = 39 GeV (bottom) Au + Au collisions. The different colors and marker styles are denoting the different centrality ranges. The auxiliary figures at the bottom show the relative systematic uncertainties.
Figure 3. Centrality and m T dependence of Lévy source parameters ( R , λ , α , R ^ ) in s N N = 62 GeV (top) and s N N = 39 GeV (bottom) Au + Au collisions. The different colors and marker styles are denoting the different centrality ranges. The auxiliary figures at the bottom show the relative systematic uncertainties.
Universe 04 00011 g003

Share and Cite

MDPI and ACS Style

Kincses, D. Lévy Analysis of HBT Correlation Functions in s N N = 62 GeV and 39 GeV Au + Au Collisions at PHENIX. Universe 2018, 4, 11. https://doi.org/10.3390/universe4010011

AMA Style

Kincses D. Lévy Analysis of HBT Correlation Functions in s N N = 62 GeV and 39 GeV Au + Au Collisions at PHENIX. Universe. 2018; 4(1):11. https://doi.org/10.3390/universe4010011

Chicago/Turabian Style

Kincses, Dániel. 2018. "Lévy Analysis of HBT Correlation Functions in s N N = 62 GeV and 39 GeV Au + Au Collisions at PHENIX" Universe 4, no. 1: 11. https://doi.org/10.3390/universe4010011

APA Style

Kincses, D. (2018). Lévy Analysis of HBT Correlation Functions in s N N = 62 GeV and 39 GeV Au + Au Collisions at PHENIX. Universe, 4(1), 11. https://doi.org/10.3390/universe4010011

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