Next Article in Journal / Special Issue
Scaling Properties of Spectra in New Exact Solutions of Rotating, Multi-Component Fireball Hydrodynamics
Previous Article in Journal / Special Issue
Investigating the Poor Match among Different Precessing Gravitational Waveforms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

PHENIX Results of Three-Particle Bose-Einstein Correlations in \({\sqrt{s_{NN}}}\) = 200 GeV Au+Au Collisions

Eszterháy Károly University Károly Róbert Campus, H-3200 Gyöngyös, Mátrai út 36, Hungary
Universe 2018, 4(3), 57; https://doi.org/10.3390/universe4030057
Submission received: 17 January 2018 / Revised: 13 February 2018 / Accepted: 19 February 2018 / Published: 19 March 2018

Abstract

:
Bose-Einstein correlations (BECs) of identical hadrons reveal information about hadron creation from the strongly interacting matter formed in ultrarelativistic heavy-ion collisions. The measurement of three-particle correlations may in particular shed light on hadron creation mechanisms beyond thermal/chaotic emission. In this paper, we show the status of PHENIX measurements of three-pion correlations as a function of momentum differences within the triplets. We analyze the shape of the correlation functions through the assumption of Lévy sources and a proper treatment of the Coulomb interaction within the triplets. We measure the three-particle correlation strength ( λ 3 ), which, together with the two-particle correlation strength λ 2 , encodes information about hadron creation mechanisms. From a consistent analysis of two- and three-particle correlation strengths, we establish a new experimental measure of thermalization and coherence in the source.

1. Introduction

In particle and nuclear physics, intensity interferometry provides a direct experimental method for the determination of sizes, shapes and lifetimes of particle-emitting sources (for reviews, see [1,2,3,4,5]). In particular, boson interferometry provides a powerful tool for the investigation of the spacetime structure of particle production processes, as Bose-Einstein correlations (BECs) of two or three identical bosons reflect both geometrical and dynamical properties of the particle-radiating source.
The size (radius) of the source in heavy-ion collisions has been found to decrease with increasing transverse momentum, p t , or transverse mass, m t = m 2 + p t 2 , of the bosons. This effect can also be explained by hydrodynamical models [6,7]. The main purpose of the present paper is to determine the three-particle correlation strength ( λ 3 ) as a function of transverse momentum. Recently [8], the two-particle correlation strength ( λ 2 ) was determined. These measurements of λ 2 , when combined with this analysis, may test the limits of the core-halo model [9,10,11] with a thermalized core. For this purpose, we introduce a new parameter, κ 3 , which is a function of λ 3 and λ 2 . This new parameter is not equal to 1 when there are extra effects in the core, for example, a core not fully thermalized, or when there is partial coherence in the core. The main purpose of the present work is to investigate whether κ 3 indicates extra effects or not.
We start with the definition of the three-particle correlation function:
C 3 ( k 1 , k 2 , k 3 ) = N 3 ( k 1 , k 2 , k 3 ) N 1 ( k 1 ) N 1 ( k 2 ) N 1 ( k 3 )
where N 3 is the three-particle invariant momentum distribution, defined by
N 3 ( k 1 , k 2 , k 3 ) = S ( r 1 , k 1 ) S ( r 2 , k 2 ) S ( r 3 , k 3 ) | Ψ k 1 , k 2 , k 3 ( r 1 , r 2 , r 3 ) | 2 Π i = 0 3 d 4 r i
and N 1 is the single-particle invariant momentum distribution, defined as
N 1 ( k ) = S ( r 1 , k 1 ) | Ψ k ( r 1 ) | 2 d 4 r
In the above equations, Ψ k 1 , k 2 , k 3 is the three-particle wave function, and Ψ k is the single-particle wave function. Furthermore, S ( r , k ) is the source distribution, which describes the probability density of particle creation at the spacetime point r with momentum k .
Our assumption for the source function is the symmetric Lévy distribution [12,13,14,15,16], which is defined by
L ( α , R , r ) = 1 ( 2 π ) 3 d 3 q e i q r e 1 2 | q R | α
where α is the Lévy index, and R is the Lévy scale. Then α = 2 gives back the Gaussian case, and α = 1 yields a Cauchy distribution. For the Lévy stable source distribution, 0 < α 2 in general.
Assuming properly symmetrized plane-waves for the wave functions and the Lévy distribution as the source function, the three-particle correlation function is given as
C 3 ( 0 ) ( k 12 , k 13 , k 23 ) = 1 + 3 e 0 . 5 ( | 2 k 12 R | α + | 2 k 13 R | α + | 2 k 23 R | α ) + 2 e | 2 k 12 R | α + e | 2 k 13 R | α + e | 2 k 23 R | α
where k i j = | k i k j | / 2 , that is, the half-momentum difference of the ith and jth particles. Furthermore, the 2 parameter is the two-particle correlation strength parameter in the three-particle correlation function, and 3 is the three-particle strength parameter.
To be able to fit the model to the data, the correlation function has to be completed with a background that parameterizes possible long-range correlations. In this analysis, we use a simple linear background in every direction, with the same slope. Thus the correlation function is the following:
C 3 , fit ( 0 ) ( k 12 , k 13 , k 23 ) = N ( 1 + ϵ k 12 ) ( 1 + ϵ k 13 ) ( 1 + ϵ k 23 ) C 3 ( 0 ) ( k 12 , k 13 , k 23 )
We also have to take into account the Coulomb interaction, as this decreases the number of particle pairs at low-momentum differences. To add Coulomb correction to the model, we define a Coulomb correction factor:
K 3 ( k 12 , k 23 , k 31 ) = d 3 r 1 d 3 r 2 d 3 r 3 S ( r 1 ) S ( r 2 ) S ( r 3 ) | Ψ k 12 , k 23 , k 31 0 ( r 1 , r 2 , r 3 ) | 2 d 3 r 1 d 3 r 2 d 3 r 3 S ( r 1 ) S ( r 2 ) S ( r 3 ) | Ψ k 12 , k 23 , k 31 C ( r 1 , r 2 , r 3 ) | 2
where Ψ 0 is the free-particle wave function (symmetrized plane wave), and Ψ C is the (symmetrized) solution of the full three-particle Coulomb problem. Instead of solving the three-particle Coulomb problem, we use a general approach to handle the three-particle Coulomb interaction, which is called the “Generalized Riverside” method [17]:
K 3 ( k 12 , k 13 , k 23 ) K 2 ( k 12 ) K 2 ( k 13 ) K 2 ( k 23 ) , with K 2 ( k i j ) = d 3 r 1 d 3 r 2 S ( r 1 ) S ( r 2 ) | Ψ k 12 0 ( r 1 , r 2 ) | 2 d 3 r 1 d 3 r 2 S ( r 1 ) S ( r 2 ) | Ψ k 12 C ( r 1 , r 2 ) | 2
This approximates the three-particle Coulomb correction well if certain conditions [17] are satisfied, and it does not depend on the directions of the momentum differences, only their magnitude. With this Coulomb correction factor, the full BEC function is
C 3 ( k 12 , k 13 , k 23 ) = C 3 , fit ( 0 ) ( k 12 , k 13 , k 23 ) K C o u l o m b ( k 12 , k 13 , k 23 )
By using this approximation, we are able to reuse the Coulomb correction calculation method developed in [18]. After some consistency checks, we conclude that using this approximation can be considered, as the relative systematic error due to this approximation is less than 2%.
The two- and three-particle correlation strengths are defined as the extrapolated intercept of the given correlation function:
λ 2 C 2 ( k 12 0 ) 1
λ 3 C 3 ( k 12 = k 13 = k 23 0 ) 1
With respect to these, it is important to note that not all pions are created directly from the strongly interacting matter. A significant fraction of pions are secondary, coming from decays. Hence the source will have two components: a core of primordial pions, and a halo consisting of the decay products of long-lived resonances (such as η , η , K s 0 , and ω ):
S = S c o r e + S h a l o
In the core-halo model, the fraction of the core is defined by
f c = N core N core + N halo
where N core is the number of pions in the core, and N halo is the number of pions in the halo. If there is partial coherence in the core, a coherence parameter can be introduced [19,20], examining the fraction of coherently produced pions in the core:
p c = N coherent N coherent + N incoherent
where N coherent + N incoherent = N core . With these, λ 2 and λ 3 can be expressed as functions of f c and p c :
λ 2 = f C 2 ( 1 p c ) 2 + 2 p c ( 1 p c )
λ 3 = 2 f C 3 ( 1 p c ) 3 + 3 p c ( 1 p c ) 2 + 3 f C 2 ( 1 p c ) 2 + 2 p c ( 1 p c )
To be able to investigate the partial coherence within the limits of this model, we have to take Equations (15)–(16) and use the measured λ 2 and λ 3 parameters as input in these equations; we also have to look for f c and p c values that solve both equations. To determine whether a solution exists, we can plot f c as a function of p c with given λ 2 and λ 3 . In this way, we obtain two curves (one from the equation with λ 2 , and the other from the equation with λ 3 ), and if these curves overlap, then the second- and third-order BEC functions are consistent with each other, and the degree of partial coherence can be determined from the region of the overlap.

2. Results

The data used in the analyses were collected by the PHENIX detector in 2010; 7.3 billion minimum bias Au+Au collisions at s N N = 200 GeV were recorded. We measured three-particle correlation functions of π π π and π + π + π + triplets in 31 m T bins ranging from 228 to 871 MeV, with event selection, single track cuts and pair cuts identical to the analysis of [8]. The visualization of three-dimensional fits is difficult; thus instead of showing two-dimensional projections of these three-dimensional fits, we have decided to show, as an illustration, various one-plus-one dimensional slices of them. These slices were obtained by requiring that the two-particle relative momenta k i , j between the pairs ( i , j ) = ( 1 , 2 ) , ( 2 , 3 ) and ( 3 , 1 ) were equal. The χ 2 / N D F values and the confidence levels were statistically acceptable. Figure 1 indicates the PHENIX preliminary three-particle correlation data and the fitted curve as a univariate slice of a multi-variate function. We note that no significant long-range correlation was observed: ϵ was close to zero. In these fits, the values of the Lévy scale R and the Lévy index of stability α were taken from PHENIX preliminary two-particle correlation data in the same transverse-momentum bins and in the same selection of minimum-bias s N N = 200 GeV Au+Au collisions. (See these preliminary R and α data in [21].) Within one standard deviation, the same parameter values were obtained if we released the values of R and α in these fits, indicating an internal consistency between the two- and three-particle correlation measurements. Given this internal consistency and our focus on the comparison of the strength or intercept parameters λ 2 and λ 3 of the two- and three-particle BEC functions, here we report these values using the R and α values as determined from the two-particle BEC function measurements in minimum-bias s N N = 200 GeV Au+Au collisions [21].

2.1. Three-Particle Correlation Strength

The three-particle correlation strength, derived from 2 and 3 fit parameters as λ 3 = 2 + 3 , can be seen in Figure 2. According to the core-halo model, λ 3 had to be in the range [ 0 , 5 ] , as indicated by a green box in Figure 2. One can see that within errors, this parameter fell in the allowed range. The detailed description of the systematic uncertainty sources is given in [8].

2.2. core-halo Independent Parameter

In this subsection, the possibility of partial coherence is investigated. We introduce a new parameter, a combination of λ 3 and λ 2 , as follows:
κ 3 = λ 3 3 λ 2 2 λ 2 3
This parameter does not depend on the value f c , and its value is unity independently of the transverse mass, if the particle emission has no coherent component or if p c = 0 in the core-halo picture of n-particle BECs.
The measured κ 3 parameter is shown in Figure 3 and was calculated by using λ 3 from the fits described in this paper, and λ 2 from [21]. Our errors were still preliminary, and a full systematic analysis is in progress. Figure 3 shows the κ 3 parameter for the average m T of the triplets.
We observe that κ 3 was different from unity at m T 450 MeV. The deviation of κ 3 from 1 may indicate some kind of partial coherence; however, the preliminary status of these data points and their error bars prevent us from a more definitive conclusion. Further investigations are ongoing and will be reported elsewhere. Figure 4 shows f c as a function of p c for different m T regions. Here we took Equations (15) and (16) and the measured λ 2 and λ 3 parameters and determined the allowed ( f c , p c ) region, as shown in Figure 4.

3. Conclusions

In this paper, we simultaneously investigated two- and three-particle femtoscopic correlation functions in PHENIX Au+Au collisions at 200 GeV/nucleon collision energy. Lévy distributions yielded correlation functions that provide a statistically acceptable description of the measurements. In this paper, we used the source parameters R and α from two-particle fits and obtained the three-particle correlation strength from the analysis of three-particle correlations. We have investigated the results in view of the core-halo model and derived a core-halo independent parameter κ 3 . These preliminary results have been compared to predictions based on chaotic emission, and very small deviations were observed. This, however, needs a detailed final analysis, which we plan to carry out in the near future.

Acknowledgments

T.N. was supported by the EFOP-3.6.1-16-2016-00001 grant. This work was also supported by the NKFIH Grants FK-123942 and FK-123959.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gyulassy, M.; Kauffmann, S.K.; Wilson, L.W. Pion interferometry of nuclear collisions. Phys. Rev. C 1979, 20, 2267–2292. [Google Scholar] [CrossRef]
  2. Boal, D.H.; Gelbke, C.K.; Jennings, B.K. Intensity interferometry in subatomic physics. Rev. Mod. Phys. 1990, 62, 553–602. [Google Scholar] [CrossRef]
  3. Baym, G. The physics of Hanbury Brown-Twiss intensity interferometry: From stars to nuclear collisions. Acta Phys. Pol. B 1998, 29, 1839–1884. [Google Scholar]
  4. Kittel, W. Bose-Einstein Correlations in Z Fragmentation and Other Reactions. Acta Phys. Pol. B 2001, 32, 3927–3972. [Google Scholar]
  5. Csörgő, T. Particle Interferometry from 40 MeV to 40 TeV. Heavy Ion Phys. 2002, 15, 1–80. [Google Scholar] [CrossRef]
  6. Wiedemann, U.; Scotto, P.; Heinz, U. Transverse momentum dependence of Hanbury-Brown-Twiss correlation radii. Phys. Rev. C 1996, 53, 918–931. [Google Scholar] [CrossRef]
  7. 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]
  8. 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]
  9. 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]
  10. Csörgő, T.; Lörstad, B.; Zimányi, J. Bose-Einstein correlations for systems with large halo. Z. Phys. C 1996, 71, 491–497. [Google Scholar] [Green Version]
  11. Csanád, M. Measurement and analysis of two- and three-particle correlations. Nucl. Phys. A 2006, 774, 611–614. [Google Scholar] [CrossRef]
  12. Nolan, J.P. Univariate Stable Distributions. In Stable Distributions—Models for Heavy Tailed Data; Publishing House: Boston, MA, USA, 2018. [Google Scholar]
  13. 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]
  14. 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]
  15. Csanád, M.; Csörgő, T.; Nagy, M. Anomalous diffusion of pions at RHIC. Braz. J. Phys. 2007, 37, 1002–1013. [Google Scholar] [CrossRef]
  16. Csanád, M. Lévy femtoscopy with PHENIX at RHIC. Universe 2017, 3, 85. [Google Scholar] [CrossRef]
  17. Gangadharan, D. Techniques for multiboson interferometry. Phys. Rev. C 2015, 92, 014902. [Google Scholar] [CrossRef]
  18. Alt, E.O.; Csörgo, T.; Lörstad, B.; Schmidt-Sorensen, J. Coulomb corrections to the three-body correlation function in high-energy heavy ion reactions. Phys. Lett. B 1999, 458, 407–414. [Google Scholar] [CrossRef]
  19. Boggild, H.; Boissevainb, J.; Doddc, J.; Esumid, S.; Fabjane, C.W.; Ferencef, D.; Franzg, A.; Hardtkeh, D.; van Heckeb, H.; Humanich, T.J.; et al. Three pion correlations in sulphur lead collisions at the CERN SPS. Phys. Lett. B 1999, 455, 77–83. [Google Scholar] [CrossRef]
  20. Csörgo, T.; Lörstad, B.; Schmid-Sørensen, J.; Ster, A. Partial coherence in the core-halo picture of Bose– Einstein n-particle correlations. Eur. Phys. J. C 1999, 9, 275–281. [Google Scholar] [CrossRef]
  21. Kincses, D.; PHENIX Collaboration. PHENIX results on Lévy analysis of Bose-Einstein correlation functions. Acta Phys. Polon. Supp. 2017, 10, 627. [Google Scholar] [CrossRef]
Figure 1. Diagonal visualization of the three-particle Bose-Einstein correlation functions C 3 with the results of fits.
Figure 1. Diagonal visualization of the three-particle Bose-Einstein correlation functions C 3 with the results of fits.
Universe 04 00057 g001
Figure 2. The three-particle correlation strength, λ 3 , as a function of m T . Statistical and systematic uncertainties are shown as bars and boxes, respectively.
Figure 2. The three-particle correlation strength, λ 3 , as a function of m T . Statistical and systematic uncertainties are shown as bars and boxes, respectively.
Universe 04 00057 g002
Figure 3. Core-halo independent parameter as a function of m T . Statistical and preliminary systematic uncertainties are shown as bars and boxes, respectively.
Figure 3. Core-halo independent parameter as a function of m T . Statistical and preliminary systematic uncertainties are shown as bars and boxes, respectively.
Universe 04 00057 g003
Figure 4. The fraction of the core f c parameter as a function of coherence parameter p c for two different m T values.
Figure 4. The fraction of the core f c parameter as a function of coherence parameter p c for two different m T values.
Universe 04 00057 g004

Share and Cite

MDPI and ACS Style

Novák, T. PHENIX Results of Three-Particle Bose-Einstein Correlations in \({\sqrt{s_{NN}}}\) = 200 GeV Au+Au Collisions. Universe 2018, 4, 57. https://doi.org/10.3390/universe4030057

AMA Style

Novák T. PHENIX Results of Three-Particle Bose-Einstein Correlations in \({\sqrt{s_{NN}}}\) = 200 GeV Au+Au Collisions. Universe. 2018; 4(3):57. https://doi.org/10.3390/universe4030057

Chicago/Turabian Style

Novák, Tamás. 2018. "PHENIX Results of Three-Particle Bose-Einstein Correlations in \({\sqrt{s_{NN}}}\) = 200 GeV Au+Au Collisions" Universe 4, no. 3: 57. https://doi.org/10.3390/universe4030057

APA Style

Novák, T. (2018). PHENIX Results of Three-Particle Bose-Einstein Correlations in \({\sqrt{s_{NN}}}\) = 200 GeV Au+Au Collisions. Universe, 4(3), 57. https://doi.org/10.3390/universe4030057

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