Next Article in Journal
Quantum Phenomena Inside a Black Hole: Quantization of the Scalar Field Iniside Horizon in Schwarzschild Spacetime
Next Article in Special Issue
Kaon Femtoscopy with Lévy-Stable Sources from sNN=200 GeV Au+Au Collisions at RHIC
Previous Article in Journal
Black String Solutions in Rainbow Gravity
Previous Article in Special Issue
Even and Odd Self-Similar Solutions of the Diffusion Equation for Infinite Horizon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Femtoscopic Correlation Measurement with Symmetric Lévy-Type Source at NA61/SHINE

by
Barnabás Pórfy
1,2 on behalf of the NA61/SHINE Collaboration
1
Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
2
Department of Atomic Physics, Faculty of Science, Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1111 Budapest, Hungary
Submission received: 25 May 2023 / Revised: 10 June 2023 / Accepted: 16 June 2023 / Published: 21 June 2023
(This article belongs to the Special Issue Zimányi School – Heavy Ion Physics)

Abstract

:
Measuring quantum-statistical, femtoscopic (including final state interactions) momentum correlations with final state interactions in high-energy nucleus-nucleus collisions reveal the space-time structure of the particle-emitting source created. In this paper, we report NA61/SHINE measurements of femtoscopic correlations of identified pion pairs and describe said correlations based on symmetric Lévy-type sources in Ar+Sc collisions at 150A GeV/c. We investigate the transverse mass dependence of the Lévy-type source parameters and discuss their possible interpretations.

1. Introduction

The NA61/SHINE is a fixed target experiment using a large acceptance hadron spectrometer located in the North Area H2 beam line of the CERN Super Proton Synchrotron accelerator [1]. Its main goals include the investigation and mapping of the phase diagram of strongly interacting matter, as well as measuring cross sections of processes relevant for cosmic rays and neutrino physics. In this paper, we are focusing on mapping the QCD phase diagram. In order to accomplish this, NA61/SHINE performs measurements of different collision systems at multiple energies. The experiment provides excellent tracking down to p T = 0 GeV/c. This performance is achieved by using four large Time Projection Chambers (TPC’s), which cover the full forward hemisphere. The experiment also features a modular calorimeter, called the Projectile Spectator Detector. It is located on the beam axis, after the TPC’s, and measures the forward energy which determines the collision centrality of the events. A setup of the NA61/SHINE detector system is shown in Figure 1.
The search for the critical endpoint (CEP) and investigation of the QCD phase diagram requires analysis at different temperatures and baryon-chemical potentials. To study, we need to map the phase diagram using different system sizes at various energies. NA61/SHINE investigations cover several beam momenta (13A, 20A, 30A, 40A, 75A and 150A GeV/c) and collision systems (p+p,p+Pb,Be+Be,Ar+Sc,Xe+La,Pb+Pb). In this paper, we describe the femtoscopic correlations of identical pions emitted from central Ar+Sc collisions at beam momentum of 150A GeV/c. This field is often called femtoscopy as it reveals the femtometer scale structure of particle production.

2. Femtoscopy with Lévy Shaped Sources

The method of quantum-statistical (Bose-Einstein) correlations is based on the work of R. Hanbury Brown and R. Q. Twiss (HBT) [2], who applied it first in astrophysical intensity correlation measurements. The method was developed to determine the apparent angular diameter of stellar objects. Shortly afterwards, a similar quantum-statistical method was applied in momentum correlation measurements for proton-antiproton collisions [3,4] by Goldhaber and collaborators. Their objective was to understand pion-pion correlations and gain information on the radius, R, of the interaction volume in high-energy particle collisions. The key relationship for measuring Bose-Einstein correlations shows that the spatial momentum correlations, C ( q ) , are related to the properties of the particle emitting source, S ( x ) , that describes the probability density of particle creation for a relative coordinate x as:
C 2 ( q ) 1 + | S ˜ ( q ) | 2 ,
where S ˜ ( q ) is the Fourier transform of S ( x ) , and q is the relative momentum of the particle pair (with the dependence on the average momentum, K, of the pair suppressed and described in more detail in [5]). The usual assumption for the shape of the source based on the central limit theorem, is a Gaussian. However, such Gaussian shaped sources lead to Gaussian correlation functions. A more general assumption is the Lévy distribution [6,7]. It exhibits a power-law tail and includes a Gaussian limit, as well. Correlation functions based on this approach have been shown to describe data from different experiments, such as LEP [8], RHIC [9], and LHC [10,11] quite well. Several phenomena could explain the appearance of Lévy shaped sources. The non-Gaussianity of the source could be attributed to critical fluctuations and the emergence of spatial correlations on a large scale, which may indicate the existence of similar sources with power-law tails [12]. Further reasons include the fractal structure of QCD jets [13].
In this paper, the measured femtoscopic correlation (including final state interaction) with spherically symmetric Lévy distributions is defined as:
L ( α , R , r ) = 1 ( 2 π ) 3 d 3 ζ e i ζ r e 1 2 | ζ R | α ,
where R is the Lévy scale parameter and α defined as the Lévy stability index. In addition, ζ is the three-dimensional integration variable with dimensions of MeV/c and r is the vector of spatial coordinates. There are two special cases where the distribution can be expressed analytically. One such case is, the already mentioned, Gaussian distribution for α = 2 . Besides this, the α = 1 case leads to a Cauchy distribution. An important difference between Lévy distributions and Gaussians is the presence of a power-law tail r ( d 2 + α ) in case of α < 2 , where d represents the number of spatial dimensions. With the assumption of Lévy sources, the femtoscopic correlation functions can be expressed in the following way:
C 2 ( q ) = 1 + λ · e ( q R ) α .
C 2 ( q ) is a stretched type of exponential, where the λ intercept parameter is defined as:
C 2 ( q 0 ) = 1 + λ .
At vanishing relative momentum, the correlation function has a value of 1 + λ . This value is not accessible in the measurements and extrapolation from the region, when two tracks are experimentally resolved, is needed. However, it is commonly observed that the intercept parameter λ is less than 1. The core-halo model, explained in Refs. [14,15], can provide some insights into this parameter.
The model assumes that the source S is made up of two parts, the core and the halo ( S core and S halo ), respectively. The core contains pions created directly from hadronic freeze-out or from extremely short lived (strongly decaying) resonances. The halo consists of pions created from longer-lived resonances and the general background. It may extend to thousands of femtometers, while core part has a size of around a few femtometers. In this picture, the λ parameter turns out to be connected to the ratio of the core and the halo as:
λ = N core N core + N halo 2 .
Then, one can modify the correlation function to take the effect of the halo into account, by utilizing the Bowler-Sinyukov method [16,17] as:
C 2 ( q ) = 1 λ + λ · ( 1 + e | q R | α ) .
The halo part contributes at very small values of relative momenta, q. Therefore, it does not affect the source radii of the core part [18].
It is well known that critical points are characterized by critical exponents. One, in particular, is related to spatial correlations by the exponent denoted as η . The appearance of the parameter can be explained by the second-order phase transition at the CEP, where fluctuations will appear at all scales causing the spatial correlation function to exhibit a power-law tail with an exponent of ( d 2 + η ) , with d denoting the dimension. The Lévy exponent, α , is the exponent in the case of, previously defined, Lévy distributed sources and will also exhibit a power-law tail, with an exponent of d + 2 α [19]. Hence, α was suggested to be directly related to or being explicitly equal to, the critical exponent, η [12], in absence of other phenomena affecting the source shape. This is the basis of the idea connecting α and η . However, the Lévy-shape of the source can be attributed to several different factors besides critical phenomena, including QCD jets, anomalous diffusion, critical phenomena, and others [6,7,13,20,21]. Hence, while a non-monotonic behavior of α is expected near the critical point, a detailed understanding of the collision energy and system size dependence of α is needed to draw conclusions about the critical point.
It has been suggested that the universality class of QCD is the same as that of the 3D Ising model [22,23]. The value of η in the 3D Ising model is 0.03631 ± 0.00003 [24]. An alternative solution is to use the universality class of the 3D Ising model with a random external field, which yields an η value of 0.50 ± 0.05 [25]. The statements mentioned suggest that α would decrease to 0.50, or below, at the vicinity of the CEP. To confirm this, measurements of α are needed in different collision systems at various energies.
In this analysis, we are dealing with like-charged particles that are influenced by Coulomb repulsion. The final state Coulomb effect has been neglected in the previously defined correlation function. Thus, the correlation function in Equation (6) that lacks this effect will be denoted as C 2 0 ( q ) from now on. The correction necessitated by this effect can be done by simply taking the ratio of C 2 Coul and C 2 0 ( q ) :
K Coulomb ( q ) = C 2 Coul ( q ) C 2 0 ( q ) ,
where C 2 Coul ( q ) is the interference of solutions of the two-particle Schrödinger equation; with a Coulomb-potential [26,27]. The numerator in Equation (7) cannot be calculated analytically and requires a large numerical effort to estimate.
An approximate formula for C 2 Coul ( q ) was obtained in Ref. [10] for the case of Cauchy-shaped sources. However, a more precise treatment is required due to our assumption of Lévy-shaped sources. We are utilizing a new method in our analysis for estimating the effect of Coulomb repulsion. The treatment includes the numerical calculation presented in Refs. [27,28], the parametrization of its results, and, finally; the parametrization of the dependence of the physical parameters R, λ , and α . Thus Equation (6) is modified as:
C 2 ( q ) = N · 1 λ + λ · ( 1 + e | q R | α ) · K Coulomb ( q ) ,
where N is introduced as normalization parameter and K Coulomb ( q ) denotes the Coulomb correction.
It is important to highlight that the Coulomb correction is calculated in the pair-center-of-mass (PCMS) system, while the measurement is often done in the longitudinally co-moving system (LCMS). The assumption of Coulomb correction in the one-dimensional HBT in LCMS picture is that the shape of the source is spherical, i.e. R out = R side = R long = R R LCMS . The shape of the source, however, is spherical in the LCMS and not in the PCMS. Therefore, an approximate one-dimensional PCMS size parameter is needed. A study was done, where an average PCMS radius of
R ¯ PCMS = 1 2 3 β T 2 1 β T 2 · R
was calculated [29], with β T = K T m T .

3. Measurement Details

For the measurements this paper is based on, we analyzed Ar+Sc collisions at 150A GeV/c beam momentum in the 0–10% most central events. The available data set contains around 2.7 million events, which was reduced to around 700,000 events in the analysis. The following paragraph describes the various event, track and pair selection performed on the data. As mentioned, we have selected 0–10% of the most central events by measuring the energy contained in the projectile remnants, with the Projectile Spectator Detector (PSD). We have also selected events where no off-time beam particle was detected. These are particles which come within the drift time of the chambers which are not emptied out. Furthermore, all the events chosen have between ± 10 cm as that is the maximal distance between the main vertex z position and the center of the scandium target. Particle identification was handled by using dE/dx energy deposit in the TPC gas. The tracks were extrapolated to the interaction plane and matched with the distance against the interaction point. If the distance was ≤4 cm in the horizontal plane and ≤2 cm in the vertical plane, the track was kept. Moreover, track splitting was handled by selecting the ratio of the total number of reconstructed points on the track to the potential number of points to be between 0.5 and 1.0. Finally, to counter track merging, we have used a selection on the momentum space distance between the two tracks. It uses non-standard momentum coordinates, s x = p x / p x z , s y = p y / p x z and ρ = 1 / p x z .
We then analyzed the combination of negative pion pairs and positive pion pairs. It is important to note that this analysis was done with a one-dimensional relative momentum variable q, calculated in LCMS. These pairs were sorted into eight K T (average transverse momentum of the pair) bins in the range of 0–450 MeV/c. In each momentum bin, the relative momentum distribution of coincident pion pairs were obtained. Let us call this the actual pair momentum difference distribution, A ( q ) . A ( q ) contains quantum-statistical correlations, as well as many other residual effects related to kinematics and acceptance. The effects can be removed by constructing a combinatorial background pair distribution, denoted as B ( q ) , which is measured in the same K T or m T intervals as the A ( q ) distribution. The method we use involves randomly selecting the same number of particles as the multiplicity of the actual event. The selected particles are from other events of similar parameters and each is selected from a different event. Let us call the pair momentum difference distribution made from this method as, the background distribution, B ( q ) . This, by construction, enables us to have an uncorrelated pool of events. Then the correlation function is calculated as
C 2 ( q ) = A ( q ) B ( q ) · q 1 q 2 B ( q ) d q q 1 q 2 A ( q ) d q ,
in a [ q 1 , q 2 ] range where quantum-statistical correlations are not expected. Upon adding the contribution from the background, the previously defined Equation (8) is modified as follows:
C ( q ) = N · ( 1 + ε · q ) · 1 λ + λ · 1 + e ( q R ) α · K Coulomb ( q ) ,
with N being a normalization parameter responsible for the proper normalization of the A ( q ) / B ( q ) ratio, ε describing the linearity of background, and K Coulomb ( q ) being the Coulomb correction. We then use this formula to our measured C 2 ( q ) as shown in Figure 2. We then use this formula to our measured C 2 ( q ) as shown in Figure 2. To determine the fitting range where the effects of detector resolution do not play a significant role, we have used EPOS simulation [30] with GEANT3 for particle propagation [31]. Note that in the low-q region, the fit does not describe the data. This can be explained, according to Monte Carlo simulations of the detector response, by the limited resolution of pairs with small relative momentum.

4. Results

The three physical parameters α ,   R ,   and   λ were measured in eight bins of pair transverse momentum, K T . The three mentioned parameters were obtained through fitting the measured correlation functions with the formula shown in Equation (11). The results were investigated regarding their transverse momentum dependence. In the following, we report on the transverse mass dependence of α , λ , and R; where transverse mass is expressed as m T = m π 2 c 4 + K T 2 c 2 , with m π being the pion mass.
As explained above, the shape of the source is often assumed to be Gaussian. The Lévy stability exponent, α , can be used to extract the shape of the tail of the source. Our results, shown in Figure 3, yield values for α between 1.5 and 2.0, which imply a source closer to the Gaussian shape than the one in Be+Be collisions [5], but are still significantly lower than the α = 2 (Gaussian) case. The observed α parameter is also significantly higher than the conjectured value at the critical point ( α = 0.5 ). Altogether, these results suggest that measured correlation functions align with the assumption of a Lévy source, indicating that it is more advantageous over the Gaussian assumption. Further studies are ongoing at NA61/SHINE, using different collision energies and system sizes, in order to map the evolution of the Lévy stability index, α , as a function of collision energy and system size.
As a second parameter, let us look at the Lévy scale parameter, R, visible in Figure 4. It determines the length of homogeneity of the pion emitting source. The parameter R depends on the transverse-mass as 1 / m T . This can be derived using simple, hydrodynamical predictions for Gaussian sources [15,32]:
1 R HBT 2 = 1 R geom 2 + u T 2 · m T T 0 ,
where u T 2 is the average, transverse expansion and T 0 is the hadronisation temperature. In our case, rather surprisingly, despite the non-Gaussian nature of the source, this formula works in describing the measured femtoscopic radii. More precisely, as mentioned above, observing an R 1 / m T is particularly interesting as this type of m T dependence should rise in case of Gaussian sources ( α = 2 ) [33]. It is not entirely clear why this happens, the indicated m T dependence could form in the QGP or at a later stage. This phenomenon was also observed at RHIC [9] and in simulations at RHIC and LHC energies [20,21].
The final parameter being investigated is the intercept (also known as the correlation strength) parameter, λ , defined in Equation (5). The dependence of λ on m T is shown in Figure 5. One may observe a slight dependence on m T , but this can still can be considered constant in the investigated range. When compared to measurements from RHIC Au+Au collisions [9,34,35] and from SPS Pb+Pb interactions [36,37], an interesting phenomenon is observed. At the energies of SPS, there is no visible “hole” at lower m T values, but at RHIC energies, the “hole” appears at m T values of around a few hundred MeV. This “hole” was interpreted in Refs. [9,38] to be a sign of in-medium mass modification. The results presented in Figure 5, at the given statistical precision, do not indicate the presence of such a low- m T hole. This trend might imply that this phenomenon can be turned off at SPS. Furthermore, it can be highlighted that the λ values we obtained are significantly below unity. A possible answer can be given by the halo part of the core-halo model. It may indicate that a significant fraction of pions are the decay products of long-lived resonances.

5. Conclusions

In the report above, we discussed the NA61/SHINE measurement of one-dimensional, identified, two-pion, femtoscopic correlation functions; in the 0–10% most central Ar+Sc collisions at 150A GeV/c. We discussed the transverse mass dependencies of the Lévy source parameters. Results on the Lévy scale parameter, α , showed a significant deviation from Gaussian sources and are not in the vicinity of the conjectured value at the critical point. The Lévy scale parameter, R, shows a visible decrease with m T . The correlation strength parameter, λ , does not show any significant m T dependence, but, maps different patterns at RHIC and similar trends at SPS energies. With these results at hand, we plan to measure Bose-Einstein correlations in larger systems, as well as at smaller energies, to continue mapping the phase diagram of the strongly interacting matter.

Funding

This research was supported by the ÚNKP-22-3 New National Excellence Program of the Ministry for Culture and Innovation from the source of the National Research, Development and Innovation Fund, and the NKFIH OTKA K-138136 grant.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available.

Acknowledgments

The author would like to thank the NA61/SHINE collaboration.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
QCDquantum chromodynamics
CERNConseil européen pour la recherche nucléaire
SPSSuper Proton Synchrotron
HBTHanbury Brown and Twiss
BEBose-Einstein
CEPcritical endpoint
NA61/SHINENorth Area 61 / SPS Heavy Ion and Neutrino Experiment
LCMSLongitudinally Co-Moving System

References

  1. Abgrall, N. et al. [NA61/SHINE Collaboration] NA61/SHINE facility at the CERN SPS: Beams and detector system. JINST 2014, 9, P06005. [Google Scholar] [CrossRef]
  2. Hanbury Brown, R.; Twiss, R.Q. A Test of a new type of stellar interferometer on Sirius. Nature 1956, 178, 1046. [Google Scholar] [CrossRef]
  3. Goldhaber, G.; Fowler, W.B.; Goldhaber, S.; Hoang, T.F. Pion-pion correlations in antiproton annihilation events. Phys. Rev. Lett. 1959, 3, 181. [Google Scholar] [CrossRef]
  4. Goldhaber, G.; Goldhaber, S.; Lee, W.Y.; Pais, A. Influence of Bose-Einstein statistics on the anti-proton proton annihilation process. Phys. Rev. 1960, 120, 300. [Google Scholar] [CrossRef]
  5. Adhikary, H. et al. [NA61/SHINE Collaboration] Measurements of two-pion HBT correlations in Be + Be collisions at 150A GeV/c beam momentum, at the NA61/SHINE experiment at CERN. arXiv 2023, arXiv:2302.04593. [Google Scholar]
  6. 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]
  7. 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] [Green Version]
  8. Achard, P. et al. [L3 Collaboration] Test of the τ-Model of Bose-Einstein Correlations and Reconstruction of the Source Function in Hadronic Z-boson Decay at LEP. Eur. Phys. J. C 2011, 71, 1648. [Google Scholar] [CrossRef] [Green Version]
  9. Adare, A. et al. [PHENIX Collaboration] Lévy-stable two-pion Bose-Einstein correlations in s N N = 200 GeV Au + Au collisions. Phys. Rev. C 2018, 97, 064911. [Google Scholar] [CrossRef] [Green Version]
  10. Sirunyan, A.M. et al. [CMS Collaboration] Bose-Einstein correlations in pp,pPb, and PbPb collisions at s N N = 0.9 − 7 TeV. Phys. Rev. C 2018, 97, 064912. [Google Scholar] [CrossRef] [Green Version]
  11. Measurement of Two-Particle Bose-Einstein Momentum Correlations and Their Levy Parameters at s NN = 5.02 TeV PbPb Collisions; Technical Report CMS-PAS-HIN-21-011; CERN: Geneva, Switzerland, 2022.
  12. 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] [CrossRef]
  13. Csorgo, T.; Hegyi, S.; Novak, T.; Zajc, W.A. Bose-Einstein or HBT correlations and the anomalous dimension of QCD. Acta Phys. Polon. B 2005, 36, 329–337. [Google Scholar]
  14. Csörgő, T. Particle interferometry from 40-MeV to 40-TeV. Acta Phys. Hung. A 2002, 15, 1–80. [Google Scholar] [CrossRef]
  15. 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] [PubMed]
  16. Sinyukov, Y.; 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] [Green Version]
  17. Bowler, M.G. Coulomb corrections to Bose-Einstein correlations have been greatly exaggerated. Phys. Lett. B 1991, 270, 69–74. [Google Scholar] [CrossRef]
  18. Maj, R.; Mrowczynski, S. Coulomb Effects in Femtoscopy. Phys. Rev. C 2009, 80, 034907. [Google Scholar] [CrossRef]
  19. Csörgő, T. Correlation Probes of a QCD Critical Point. arXiv 2009, arXiv:0903.0669. [Google Scholar] [CrossRef] [Green Version]
  20. Kincses, D.; Stefaniak, M.; Csanád, M. Event-by-Event Investigation of the Two-Particle Source Function in Heavy-Ion Collisions with EPOS. Entropy 2022, 24, 308. [Google Scholar] [CrossRef]
  21. Kórodi, B.; Kincses, D.; Csanád, M. Event-by-event investigation of the two-particle source function in s N N = 2.76 TeV PbPb collisions with EPOS. arXiv 2022, arXiv:2212.02980. [Google Scholar]
  22. 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]
  23. 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] [Green Version]
  24. 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, 869. [Google Scholar] [CrossRef] [Green Version]
  25. Rieger, H. Critical behavior of the three-dimensional random-field Ising model: Two-exponent scaling and discontinuous transition. Phys. Rev. B 1995, 52, 6659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kincses, D.; Nagy, M.I.; Csanád, M. Coulomb and strong interactions in the final state of Hanbury-Brown–Twiss correlations for Lévy-type source functions. Phys. Rev. C 2020, 102, 064912. [Google Scholar] [CrossRef]
  27. Csanád, M.; Lökös, S.; Nagy, M. Expanded empirical formula for Coulomb final state interaction in the presence of Lévy sources. Phys. Part. Nucl. 2020, 51, 238–242. [Google Scholar] [CrossRef]
  28. Csanád, M.; Lökös, S.; Nagy, M. Coulomb final state interaction in heavy ion collisions for Lévy sources. Universe 2019, 5, 133. [Google Scholar] [CrossRef] [Green Version]
  29. Kurgyis, B.; Kincses, D.; Nagy, M.; Csanád, M. Coulomb interaction for Lévy sources. arXiv 2020, arXiv:2007.10173. [Google Scholar]
  30. Pierog, T.; Werner, K. EPOS Model and Ultra High Energy Cosmic Rays. Nucl. Phys. Proc. Suppl. 2009, 196, 102–105. [Google Scholar] [CrossRef] [Green Version]
  31. Brun, R.C.F. GEANT Detector Description and Simulation Tool, CERN Program Library Long Writeup W5013. 1993. Available online: http://wwwasdoc.web.cern.ch/wwwasdoc/geant/geantall.html (accessed on 1 May 2023).
  32. Csanád, M.; Vargyas, M. Observables from a solution of 1+3 dimensional relativistic hydrodynamics. Eur. Phys. J. A 2010, 44, 473–478. [Google Scholar] [CrossRef] [Green Version]
  33. Sinyukov, Y.M. Spectra and correlations in locally equilibrium hadron and quark-gluon systems. Nucl. Phys. A 1994, 566, 589C–592C. [Google Scholar] [CrossRef]
  34. 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] [Green Version]
  35. Abelev, B.I. et al. [STAR Collaboration] Pion Interferometry in Au+Au and Cu+Cu Collisions at RHIC. Phys. Rev. C 2009, 80, 024905. [Google Scholar] [CrossRef]
  36. Beker, H. et al. [NA44 Collaboration] m(T) dependence of boson interferometry in heavy ion collisions at the CERN SPS. Phys. Rev. Lett. 1995, 74, 3340–3343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Alt, C. et al. [NA49 Collaboration] Bose-Einstein correlations of pi-pi-pairs in central Pb+Pb collisions at A-20, A-30, A-40, A-80, and A-158 GeV. Phys. Rev. C 2008, 77, 064908. [Google Scholar] [CrossRef] [Green Version]
  38. Vance, S.E.; Csörgő, T.; Kharzeev, D. Partial UA1 restoration from Bose-Einstein correlations. Phys. Rev. Lett. 1998, 81, 2205–2208. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The setup of the NA61/SHINE detector system during the run of Ar+Sc.
Figure 1. The setup of the NA61/SHINE detector system during the run of Ar+Sc.
Universe 09 00298 g001
Figure 2. Example fit with Bose–Einstein correlation function at K T = 0.25 0.30 GeV/c for the sum π + + π + + π + π . Blue points with error bars represent the data, the green dash-dotted line shows the fitted function with Coulomb correction given by Equation (11) within the range of 0.0525 GeV/c to 0.2 GeV/c and red dashed line represent normalization to background 0.2 GeV/c to 0.4 GeV/c. In the low-q region, the black dotted line indicates the extrapolated function outside of the fit range due to prominent detector resolution effects, mentioned in the text.
Figure 2. Example fit with Bose–Einstein correlation function at K T = 0.25 0.30 GeV/c for the sum π + + π + + π + π . Blue points with error bars represent the data, the green dash-dotted line shows the fitted function with Coulomb correction given by Equation (11) within the range of 0.0525 GeV/c to 0.2 GeV/c and red dashed line represent normalization to background 0.2 GeV/c to 0.4 GeV/c. In the low-q region, the black dotted line indicates the extrapolated function outside of the fit range due to prominent detector resolution effects, mentioned in the text.
Universe 09 00298 g002
Figure 3. The Lévy stability index, α , for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Special cases corresponding to a Gaussian ( α = 2 ) or a Cauchy ( α = 1 ) source are shown, as well as α = 0.5 , the conjectured value corresponding to the critical endpoint. Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Figure 3. The Lévy stability index, α , for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Special cases corresponding to a Gaussian ( α = 2 ) or a Cauchy ( α = 1 ) source are shown, as well as α = 0.5 , the conjectured value corresponding to the critical endpoint. Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Universe 09 00298 g003
Figure 4. The radial scale parameter, R, for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Figure 4. The radial scale parameter, R, for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Universe 09 00298 g004
Figure 5. The correlation strength parameter, λ , for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Figure 5. The correlation strength parameter, λ , for 0–10% central Ar+Sc at 150A GeV/c, as a function of m T . Boxes denote systematic uncertainties, bars represent statistical uncertainties.
Universe 09 00298 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pórfy, B., on behalf of the NA61/SHINE Collaboration. Femtoscopic Correlation Measurement with Symmetric Lévy-Type Source at NA61/SHINE. Universe 2023, 9, 298. https://0-doi-org.brum.beds.ac.uk/10.3390/universe9070298

AMA Style

Pórfy B on behalf of the NA61/SHINE Collaboration. Femtoscopic Correlation Measurement with Symmetric Lévy-Type Source at NA61/SHINE. Universe. 2023; 9(7):298. https://0-doi-org.brum.beds.ac.uk/10.3390/universe9070298

Chicago/Turabian Style

Pórfy, Barnabás on behalf of the NA61/SHINE Collaboration. 2023. "Femtoscopic Correlation Measurement with Symmetric Lévy-Type Source at NA61/SHINE" Universe 9, no. 7: 298. https://0-doi-org.brum.beds.ac.uk/10.3390/universe9070298

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