Next Article in Journal
Cosmic Microwave Background from Effective Field Theory
Next Article in Special Issue
Higgs and BSM Studies at the LHC
Previous Article in Journal / Special Issue
Estimating the Variation of Neutron Star Observables by Dense Symmetric Nuclear Matter Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lévy HBT Results at NA61/SHINE

by
Barnabás Pórfy
1,2
1
Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
2
Department of Atomic Physics, Institute of Physics, Faculty of Science, Eötvös Loránd University, H-1111 Budapest, Hungary
Universe 2019, 5(6), 154; https://doi.org/10.3390/universe5060154
Submission received: 14 May 2019 / Revised: 4 June 2019 / Accepted: 14 June 2019 / Published: 16 June 2019
(This article belongs to the Special Issue The Zimányi School and Analytic Hydrodynamics in High Energy Physics)

Abstract

:
Bose–Einstein (or Hanbury–Brown and Twiss (HBT)) momentum correlations reveal the space–time structure of the particle emitting source created in high energy nucleus–nucleus collisions. In this paper we present the latest NA61/SHINE measurements of Bose–Einstein correlations of identified pion pairs and their description based on Lévy distributed sources in Be + Be collisions at 150A GeV/c. We investigate the transverse mass dependence of the Lévy source parameters and discuss their possible interpretations.

1. Introduction

NA61/SHINE is a fixed target experiment at the CERN Super Proton Synchrotron (SPS). One of its main aims is to study the phase diagram of QCD. In order to accomplish that, different collision systems at multiple energies are investigated. The NA61/SHINE detector is equipped with four large time projection chambers (TPC) [1], these are covering the full forward hemisphere providing excellent tracking down to a transverse momentum of 0 GeV/c. The experiment also features a modular calorimeter, located on the beam axis after the TPCs. This detector is called the projectile spectator detector, and it 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.
In order to study the QCD phase diagram and search for the critical end point (CEP), vastly different collision systems (p + p, p + Pb, Be + Be, Ar + Sc, Xe + La, Pb + Pb) are investigated at various beam momenta (13A, 20A, 30A, 40A, 75A and 150A GeV/c). There are many observables to accomplish this goal. In the analysis described in this paper we measure Bose–Einstein (or Hanbury–Brown and Twiss (HBT)) correlations of identical pions in Be + Be collisions at 150A GeV/c. These, based on the principles of quantum-statistical correlations, reveal the femtometer scale structure of pion production, hence this field is often called femtoscopy.

2. Femtoscopy, Lévy Sources and the Critical End Point

The method of femtoscopy is based on the work of R. Hanbury Brown and R. Q. Twiss [2] as well as Goldhaber and collaborators [3]. The key relationship of this method shows that the spatial momentum correlations ( C ( q ) ) are related to the properties of the particle emitting source ( S ( x ) , describing the probability density of particle creation) in the following way:
C ( q ) 1 + | S ˜ ( q ) | 2 ,
where S ˜ ( q ) is the Fourier transform of S ( x ) , and q is the momentum difference of the pair (dependence on the average momentum K is suppressed here). See more details e.g., in ref. [4]. The usual assumption for the shape of the source is, based on the central limit theorem, Gaussian. A generalization of this assumption is to assume Lévy distributed sources. A possible reason is that due to the expanding medium, the mean free path may increase and thus anomalous diffusion and Lévy distributed sources may appear [5,6]. Alternatively, due to critical fluctuations and the appearance of large scale spatial correlations, similar power-law tailed sources may be present [7]. Another reason for Lévy distributed sources may be the fractal structure of QCD jets, as discussed in ref. [8]. Here we restrict our investigation to symmetric Lévy distributions, as they have proven to provide a suitable description of Bose–Einstein correlations in nucleus-nucleus collisions [4]. Furthermore, we restrict ourselves to describe the spatial part of the source, and the time dependence is absorbed through the connection of momentum difference q and average momentum K in case of identical particles:
q K = q 0 K 0 .
Then the symmetric Lévy distribution is characterized by two parameters: Lévy scale parameter R and the Lévy exponent α . The distribution is defined as follows:
L ( α , R , r ) = 1 ( 2 π ) 3 d 3 q e i q r e 1 2 | q R | α .
This distribution can be expressed analytically in two special cases. One is the already mentioned Gaussian distribution for α = 2 ; furthermore, α = 1 leads to a Cauchy distribution. An important difference between Lévy distributions and Gaussians is the presence of a power-law tail in case of α < 2 , i.e., for large distances (r), the following holds:
L ( α , R , r ) r ( d 2 + α ) ,
where d represents the number of spatial dimensions. With Lévy sources, the Bose–Einstein or HBT correlation functions can be expressed in the following way:
C ( q ) = 1 + λ · e ( q R ) α ,
where the λ intercept parameter was introduced, which is defined as
λ = lim q 0 C ( q ) ,
where the q 0 extrapolation is done in experimentally available q regions (limited from below by the two-track resolution in the given measurement). The Core-Halo model [9,10] may be utilized to understand this λ intercept parameter. The core–halo model splits the source into two pieces. The core part contains the primordially created pions (directly from hadronic freeze-out) or from short lived (strongly decaying) resonances. The halo consists of pions created from longer lived (compared to the usual source size of a few femtometers) resonances and general background. In this picture, the λ parameter turns out to be connected to the ratio of the Core and the Halo as follows:
λ = N core N core + N halo 2 .
Finally, let us come back to the above mentioned point of connecting Lévy sources to the search for the Critical End Point of QCD. Critical points are characterized by critical exponents, one of which is the exponent of spatial correlations. This appears because due to the second order phase transition at the CEP, the spatial correlation functions becomes a power-law with an exponent of ( d 2 + η ) (where d is the dimension and η is the critical exponent of spatial correlations). We can see that the Lévy exponent α , given in Equation (4), defines a similar power-law, and hence α may be regarded as identical to the critical exponent η . See further discussions in Ref. [7]. Critical exponents are universal in the sense that they take the same values in case of physical systems belonging to the same universality class. It has been shown [11] that the universality class of QCD is that of the 3D Ising model. The value of the critical exponent η has been calculated to be 0.03631 ( 3 ) [12]. Alternatively, one may rely on the universality class of the 3D Ising model with a random external field, in which case an η value of 0.5 ± 0.05 was calculated [13]. Considering the previous statements, if we “scan” the phase diagram with different energies and systems and measure the values of the α exponent, we might be able to gain more information on the location and characteristics of the CEP.

3. Measurement Details

In this measurement we analyzed the 0–20% most central Be + Be collisions at 150A GeV/c. This dataset consists of about three million events, which after various event and track quality selections was reduced to around 300,000 events. The track acceptance in this analysis was as follows. The rapidity region of analyzed particles is 0.85 < η < 4.85 (corresponding to | η | < 2 in the center-of-mass frame), the azimuthal coverage is 2 π ; n this track sample, we identified pions based on their deposited energy d E / d x in the TPC gas and charge obtained from the curvature of their trajectories in the magnetic field. We then analyzed negative pion pairs and positive pion pairs, as well as the combination of these two (i.e., created a dataset of identically charged pion pairs). These pairs were sorted into four K T (average pair transverse momentum) bins in the range of 0–600 MeV/c. In each momentum bin, we measured the pair distribution of pairs from the same event, let us call this the A ( q ) actual pair distribution. This contains quantumstatistical correlations, as well as many other residual effects related to kinematics and acceptance. To remove this undesirable effects, we created a mixed event for each actual event, by randomly selecting particles from other events of similar parameters, and making sure particles are each selected from a different event. Let us call the pair distribution from this sample B ( q ) , the background distribution. Then the correlation function is calculated as C ( q ) = A ( q ) / B ( q ) , provided a proper normalization is done in a q range where quantumstatistical correlations are not expected. Let us mention here, that our analysis was done with a one-dimensional momentum difference variable q, calculated in the longitudinally co-moving system (LCMS), as in this frame, an approximately spherically symmetric source can be expected, furthermore, the extrapolation of q 0 is equivalent with the three dimensional case.
Final state effects are still present in the C ( q ) correlation function. Among these, for pion pairs, the most important is the Coulomb effect, responsible for the repulsion of same charged pairs. It is usually handled by a so-called Coulomb correction as follows. The Coulomb correction for Lévy type sources is a complicated numerical integral to calculate and fit, as discussed in Ref. [4]. However it can be observed that the Coulomb correction does not strongly depend on the Lévy exponent α . Hence we can use then the approximate formula published by the CMS Collaboration in ref. [14], valid for α = 1 :
K Coulomb ( q , R ) = Gamow ( q ) · 1 + π η ( q ) q R c 1.26 + q R c , where
Gamow ( q ) = 2 π η ( q ) e 2 π η ( q ) 1   and   η ( q ) = α QED · π q ,
where η ( q ) is the so-called Sommerfeld parameter, and α QED is the fine-structure constant (neither should be confused with the above discussed exponents η and α ). Utilizing the usual Bowler–Sinyukov method (i.e., Coulomb correcting only for the Core part of the source) as indicated in Refs. [15,16], one obtains:
C ( q ) = N · 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. This is the final fit function we are using to describe our data.
As also found in ref. [4], the Lévy parameters ( α , R , λ ) are highly correlated, and especially in a low statistics dataset, it is hard to determine them precisely. We might be able to reduce this correlation and the statistical uncertainty of the parameters, if we fix one of the three parameters to a well motivated value. The resulting statistical uncertainties and free parameter values are modified due to the additional physical assumptions used to fix one of the parameters. From a statistical point of view, a bootstrap type of method may also be used. However, our main aim with this is to see a more clear trend of the m T dependence of the parameters, with additional physical assumptions. One assumption is that α (i.e., the shape of the pion emitting source) is independent of m T , with that we may fix α to a weighted average of the four α values obtained in free parameter fits performed in each K T bin. The other option is fixing R with the following equation motivated by hydrodynamical predictions of the particle emission homogeneity length (essentially the HBT radii) in case of expanding fireballs [10]. In this case, we fit the following equation to the m T dependence (where m T = m 2 + K T 2 , i.e., the average transverse mass of the pair) of the R Lévy scale:
R ( m T ) = A 1 + m T / B .
Previous results with free parameter fits were shown in Ref. [17], hence here we concentrate on the results of the above mentioned fixed parameter fits. We again note that our aim with fixing one of the parameters to a physically motivated value is to show the trend of the m T dependence of the parameters.

4. Results

First, the measured correlation functions were fitted with the above mentioned (Equation (10)) function with three free parameters ( α ,   λ   and   R ), as shown in Ref. [17]. Using the results from the free parameter fit, we fitted a constant function to the α values for all m T bins, as well as the formula of Equation (11) to the R values in each bin. Then we analyzed first with one parameter fixed. All results are shown in Figure 2, Figure 3 and Figure 4. Let us note here, that all the measurement settings (event selection, track selection, pair cuts, fitting interval) were varied systematically to obtain an estimate of systematic uncertainties. These, along with the statistical uncertainty of the parameters (obtained by the Minos algorithm) are shown in Figure 2, Figure 3 and Figure 4.
The Lévy stability exponent α determines the source shape, and a value of 2 corresponds to a Gaussian source, a value of 1 to a Cauchy source, and 0.5 is the conjectured value at the critical point. Our results, along with these special cases (dotted yellow lines), are shown in Figure 2. It is clear from the figure that the statistical uncertainties of α from the fixed R fit are reduced by a factor 4–5; however, the values of α are similar in both cases. These values are far from the Gaussian case as well as the conjectured CEP value, motivating us to perform this measurement in different systems and at different energies as well.
The Lévy scale R determines the correlation length of pion pairs from the given system. From a simple hydrodynamical picture one obtains a R 1 / m T type of affine linear dependence, as already mentioned above. This, i.e., Equation (11), describes the free parameter fit data points of Figure 3 well. The fixed α fits are within uncertainties compatible with the free parameter results, however, they suggest a more or less constant trend of R versus m T in the higher m T bins. This motivated us to perform the same measurement in collision systems with larger multiplicities (where statistical uncertainties are expected to be reduced proportionally to the square of the mean multiplicity).
The last parameter to study was the correlation strength parameter λ , given in Equation (7). The transverse mass dependence of λ is shown in Figure 4. Comparing the three different fits (free parameter fit, fixed α fit, fixed R fit), it was visible that λ in free parameter fit was compatible (within statistical uncertainties) with both other cases. All fits showed a roughly constant λ ( m T ) trend. This was in contrast to the findings at RHIC, see e.g., the compilation in ref. [18]. This finding was, however, compatible with previous SPS measurements (in different systems), see e.g., ref. [19].

5. Conclusions

We reported above on the NA61/SHINE measurement of one-dimensional identified two-pion Bose-Einstein correlation functions in the 0–20% most central Be + Be collisions at 150A GeV/c. We compared free parameter fits to fixed parameter fits, to reduce statistical uncertainty of the physical parameters ( α ,   λ   and   R ). We found that the results from the free parameter fits and the fixed parameter fits are similar, but the statistical uncertainty of each parameter is reduced by a large factor. The aim of this excercise was to show the trends of the parameters with fixing one parameter to a physically motivated value. Our results confirmed that in this collision system and at this collision energy, the Lévy exponent α is far from the Gaussian case, as well as from the conjectured value at the critical end point. We furthermore found that the R ( m T ) dependence is compatible with hydro predictions, and λ ( m T ) may show different patters at RHIC and SPS energies. These findings will be subsequently investigated in other collision systems as well.

Funding

Barnabás Pórfy was supported by the NKFIH grants FK123842 and FK123959.

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
CERNEuropean Organization for Nuclear Research
SPSSuper Proton Synchrotron
HBTHanbury–Brown and Twiss
CEPCritical end point
NA61/SHINESPS Heavy Ion and Neutrion Experiment

References

  1. Abgrall, N.; Andreeva, O.; Aduszkiewicz, A.; Ali, Y.; Anticic, T.; Antoniou, N.; Baatar, B.; Bay, F.; Blondel, A.; Blumer, J.; et al. NA61/SHINE facility at the CERN SPS: Beams and detector system. J. Instrum. 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–1048. [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–183. [Google Scholar] [CrossRef]
  4. 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. Phys. Rev. 2018, C97, 064911. [Google Scholar] [CrossRef]
  5. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  6. Csanad, M.; Csorgo, T.; Nagy, M. Anomalous diffusion of pions at RHIC. Braz. J. Phys. 2007, 37, 1002–1013. [Google Scholar] [CrossRef] [Green Version]
  7. Csorgo, T.; Hegyi, S.; Novak, 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]
  8. Csorgo, T.; Hegyi, S.; Novak, T.; Zajc, W.A. Bose-Einstein or HBT correlations and the anomalous dimension of QCD. Acta Phys. Polon. 2005, B36, 329–337. [Google Scholar]
  9. Csorgo, T. Particle interferometry from 40-MeV to 40-TeV. Acta Phys. Hung. 2002, A15, 1–80. [Google Scholar] [CrossRef]
  10. Csorgo, T.; Lorstad, B. Bose-Einstein correlations for three-dimensionally expanding, cylindrically symmetric, finite systems. Phys. Rev. 1996, C54, 1390–1403. [Google Scholar]
  11. 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]
  12. 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$$c-Minimization and Precise Critical Exponents. J. Stat. Phys. 2014, 157, 869–914. [Google Scholar] [CrossRef]
  13. 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] [PubMed] [Green Version]
  14. Sirunyan, A.M.; CMS Collaboration. Bose-Einstein correlations in pp, pPb, and PbPb collisions at s N N = 0.9–7 TeV. Phys. Rev. 2018, C97, 064912. [Google Scholar] [CrossRef]
  15. Sinyukov, Y.; Lednicky, R.; Akkelin, S.V.; Pluta, J.; Erazmus, B. Coulomb corrections for interferometry analysis of expanding hadron systems. Phys. Lett. 1998, B432, 248–257. [Google Scholar] [CrossRef]
  16. Bowler, M.G. Coulomb corrections to Bose-Einstein correlations have been greatly exaggerated. Phys. Lett. 1991, B270, 69–74. [Google Scholar] [CrossRef]
  17. Porfy, B. NA61/SHINE results on Bose-Einstein correlations. In Proceedings of the 12th International Workshop on Critical Point and Onset of Deconfinement (CPOD 2018), Corfu, Greece, 24–28 September 2018. [Google Scholar]
  18. Vertesi, R.; Csorgo, 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. 2011, C83, 054903. [Google Scholar]
  19. Beker, H.; Boggild, H.; Boissevain, J.; Cherney, M.; Dodd, J.; Esumi, S.; Fabjan, C.W.; Fields, D.E.; Franz, A.; Hansen, K.H.; et al. 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]
Figure 1. The setup of the NA61/SHINE detector system.
Figure 1. The setup of the NA61/SHINE detector system.
Universe 05 00154 g001
Figure 2. Lévy exponent α versus transverse mass: comparison between free parameter fit, fit with R fixed and fit with α fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Figure 2. Lévy exponent α versus transverse mass: comparison between free parameter fit, fit with R fixed and fit with α fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Universe 05 00154 g002
Figure 3. Correlation strength R versus transverse mass: comparison between free parameter fit, fit with α fixed and fit with R fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Figure 3. Correlation strength R versus transverse mass: comparison between free parameter fit, fit with α fixed and fit with R fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Universe 05 00154 g003
Figure 4. Lévy exponent λ versus transverse mass: comparison between free parameter fit, fit with α fixed and fit with R fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Figure 4. Lévy exponent λ versus transverse mass: comparison between free parameter fit, fit with α fixed and fit with R fixed. The boxes represent systematic uncertainties. For each bin, the results are slightly shifted to the right for visibility, but they are in the same bin.
Universe 05 00154 g004

Share and Cite

MDPI and ACS Style

Pórfy, B. Lévy HBT Results at NA61/SHINE. Universe 2019, 5, 154. https://doi.org/10.3390/universe5060154

AMA Style

Pórfy B. Lévy HBT Results at NA61/SHINE. Universe. 2019; 5(6):154. https://doi.org/10.3390/universe5060154

Chicago/Turabian Style

Pórfy, Barnabás. 2019. "Lévy HBT Results at NA61/SHINE" Universe 5, no. 6: 154. https://doi.org/10.3390/universe5060154

APA Style

Pórfy, B. (2019). Lévy HBT Results at NA61/SHINE. Universe, 5(6), 154. https://doi.org/10.3390/universe5060154

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