Next Article in Journal
Traceable Scheme of Public Key Encryption with Equality Test
Next Article in Special Issue
Heavy-Ion Collisions toward High-Density Nuclear Matter
Previous Article in Journal
Simulating Solid-Liquid Phase-Change Heat Transfer in Metal Foams via a Cascaded Lattice Boltzmann Model
Previous Article in Special Issue
Study of Kinetic Freeze-Out Parameters as a Function of Rapidity in pp Collisions at CERN SPS Energies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Event-by-Event Investigation of the Two-Particle Source Function in Heavy-Ion Collisions with EPOS

1
Department of Atomic Physics, Eötvös Loránd University, H-1117 Budapest, Hungary
2
GSI Helmholtzzentrum für Schwerionenforschung, Planckstr. 1, 64291 Darmstadt, Germany
3
Faculty of Physics, Warsaw University of Technology, pl. Politechniki 1, 00-661 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(3), 308; https://doi.org/10.3390/e24030308
Submission received: 20 January 2022 / Revised: 18 February 2022 / Accepted: 18 February 2022 / Published: 22 February 2022

Abstract

:
Exploring the shape of the pair-source function for particles such as pions or kaons has been an important goal of heavy-ion physics, and substantial effort has been made in order to understand the underlying physics behind the experimental observations of non-Gaussian behavior. In experiments, since no direct measurement of the source function is possible, quantum-statistical momentum correlations are utilized to gain information about the space-time geometry of the particle emitting source. Event generators, such as EPOS, however, provide direct access to the freeze-out coordinates of final state particles, and thus the source function can be constructed and investigated. The EPOS model is a sophisticated hybrid model where the initial stage evolution of the system is governed by Parton-Based Gribov–Regge theory, and subsequently a hydrodynamic evolution is utilized, followed by hadronization and hadron dynamics. EPOS has already proven to be successful in describing several different experimental observations for systems characterized by baryon chemical potential close to zero, but so far the source shape has not been explored in detail. In this paper we discuss an event-by-event analysis of the two-particle source function in s N N = 200 GeV Au+Au collisions generated by the EPOS model. We find that when utilizing all stages of the model, Lévy-shaped distributions (unlike Gaussian distributions) provide a good description of the source shape in the individual events. Hence it is clear that it is not the event averaging that creates the non-Gaussian features in the pair distributions. Based on this observation, we determine Lévy-parameters of the source as a function of event centrality and particle momentum.

1. Introduction

A long-standing goal of high-energy nuclear physics has been to understand and explore the space-time geometry of the particle emitting source created in heavy-ion collisions [1]. One main observable that is of great interest is the two-particle source function, sometimes also called spatial correlation function or pair-separation distribution. Detailed investigation of this quantity is important for a multitude of reasons (connected to hydrodynamic expansion, critical behavior, etc.), however, it is not something that is easy to reconstruct experimentally. There is a whole sub-field of high-energy nuclear- and particle-physics called femtoscopy, which deals with such measurements of lengths and time intervals on the femtometer (fm) scale [2]. Since it was shown by G. Goldhaber et al. that intensity correlations of identical pions can be used to gain information about the pair-source function [3,4], femtoscopy has propelled to the forefront of investigations, and today it is still one of the most extensively studied field of high-energy physics. Besides the ample experimental studies, phenomenological investigations also placed great emphasis on describing the shape of the source function. Hydrodynamical model calculations suggest [5,6,7,8,9,10] that the source-shape is Gaussian, and this was adopted by many measurements as well [11,12].
Source imaging studies [13,14] on the other hand suggest that the two-particle source function of pions has a long-range component, obeying a power-law behavior. It was also shown recently by various experimental measurements, that a generalization of the Gaussian source shape, the Lévy distribution can provide a much more suitable description of the observed sources [15,16]. These kind of source shapes arise in many different scenarios [17] such as anomalous diffusion [18], jet fragmentation [19], critical behavior [20], or resonance decays. It has been shown that even averaging over many events with different sizes can contribute to the appearance of a power-law component [21,22]. In order to have a better understanding of the underlying processes behind the experimental results, more effort is needed from the phenomenology side. Among the important tools for such investigations are the event-generators that encompass different theoretical and phenomenological methods to model nuclear reactions. One of such event generators is the EPOS model [23]—the Energy conserving quantum mechanical multiple scattering approach, based on Partons (parton ladders), Off-shell remnants, and Splitting of parton ladders. In this paper we present a detailed event-by-event analysis of the two-pion source distribution in s N N = 200 GeV Au+Au collisions generated by EPOS. The event-by-event nature of our analysis helps in deciding if the role of event averaging is crucial in the apparent non-Gaussian but Lévy nature of the observed sources.
The paper is structured as follows. In Section 2 we introduce the EPOS model and its several stages of evolution. In Section 3 we discuss the basic definitions and properties of the pair source distribution. In Section 4 we discuss the details of the event-by-event analysis and the applied methods. In Section 5 we present the results of the analysis and discuss our findings. Finally, in Section 6 we summarize and conclude.

2. The EPOS Model

EPOS, Energy conserving quantum mechanical multiple scattering approach, based on Partons (parton ladders), Off-shell remnants, and Saturation of parton ladders, is a phenomenological model based on Monte Carlo techniques. It opens up the possibility of investigating various phenomena and observables such as particle production, momentum distributions or flow correlations, providing a better understanding of the evolution of the system created in elementary (proton-proton) collisions and also during complex reactions involving heavy-ions. The theoretical framework included in the model provides a coherent description of the space-time expansion of matter based on a precise spectrum of studies of both elementary processes such as electron-positron annihilation or lepton-nucleon scattering and more compound collisions of protons or nuclei. The model was designed to describe processes appearing in collisions at μ B 0 , at very high (top RHIC or LHC) energies and for various systems, such as Au+Au, Pb+Pb or p+p.
The EPOS model consists of several phases of evolution:
  • initial stage (based on the Parton Gribov–Regge theory);
  • core/corona division;
  • hydrodynamical evolution;
  • hadronization;
  • hadron rescattering;
  • resonance decays.
In this Section, all of these stages are described.

2.1. Initial Stage of the Evolution

In the theoretical framework of the model the crucial element is the sophisticated treatment of both the hadron-hadron scattering and the initial stage of the collisions at ultra-relativistic energies. It is highly relevant in the understanding of possible parton–hadron phase transitions. In EPOS, a merged approach of the Gribov–Regge Theory (GRT) and the eikonalised parton model is utilised to provide proper treatment of the first interactions happening just after a collision. This approach satisfies conservation laws, and treats the subsequent Pomerons (interactions) equally (as opposed to other multiple interaction approaches, for example, Pythia, where the first interaction is not treated exactly the same way as the others) [24].
The formalism used for the calculation of the cross-sections is the same as the one used for calculating particle production. It is based on the Feynman diagrams of the QCD-inspired effective field theory and provides energy conservation. The nucleons are divided into a certain number of “constituents” carrying the incident momentum fraction. The fractions sum to unity in order to ensure momentum conservation. A nucleon is called a spectator if it is not part of the interaction region of the colliding nuclei. If a nucleon is not a spectator, then its constituents can either be participants taking part in the elementary interactions with constituents from the opposite side, or a remnant, which although part of the interaction region, does not take part in the elementary interactions. This is illustrated in Figure 1.
The particle production is based on the String Model approach [26,27]. The parton ladders are recognized as a quasi-longitudinal color field (elementary flux tubes) and are treated as classical strings [23]. The intermediate gluons introduce the transverse motion into the kinky string evolution. The schematic picture of the flux tube with the transverse kink is shown on Figure 2.

2.2. Core-Corona Approach

If the density of the strings is very high, they cannot decay independently. This situation is characteristic for heavy-ion collisions and high-multiplicity p p collisions. Henceforth, in EPOS a dynamical process of the division of the strings segments into core and corona is introduced [23,28,29].
The core-corona division is based on the abilities of a given string segment to leave the “bulk matter”. The transverse momentum of the element and the local string density are considered as criteria for the division. If the string segment belongs to a very dense area, it will not escape but will contribute to the core, which will be governed in the next step by the hydrodynamical evolution. When the segment originates from the part of the string close to the kink, it is characterized by high transverse momenta. It escapes the bulk matter and will contribute to the corona. It consequently will show up as a hadron in a jet. There is also a possibility that the string segment is close to the surface of the dense part of the medium, and its momentum is high enough to leave it; then it also becomes a corona particle.

2.3. Viscous Hydrodynamical Evolution, Event-by-Event Treatment and EoS

In EPOS a 3D+1 viscous hydrodynamics approach is applied, called vHLLE (viscous relativistic Harten-Lax-van Leer-Einfeld Riemann solver-based algorithm). In the simulations, the separate treatment of individual events is highly important—smooth initial conditions for all events are not applied. The event-by-event approach in hydrodynamical evolution is based on the random flux tube initial conditions [23]. It has a relevant impact on the final observables such as spectra or various harmonics of flow. The viscous hydrodynamics uses Equation of State X3F (“cross-over” and “3 flavor conservation”) which is compatible with lattice QCD data from Ref. [30]. It corresponds to μ B = 0 MeV, hence this feature limits the applicability of the model to describe the region of the QCD phase diagram characterized by finite baryon density [23].

2.4. Hadronization and Hadronic Cascades

The expanding medium in the processes of hydrodynamical evolution reaching the given freeze-out condition is transformed into the particle spectra. In EPOS 3, the criterion characterizing the hadronization hypersurface is that the energy density equals 0.57 GeV/fm 3 . EPOS 3 furthermore utilizes the Cooper–Frye formula [31] when determining distributions. The final part of the simulation uses a so-called hadronic afterburner, based on UrQMD [32,33]. The hadronic scattering has a significant impact on the final observables [34].

3. the Two-Particle Source Function

In this section we discuss the basics definitions and properties of the two-particle source function. The pair source distribution is defined as the auto-correlation of the single particle phase-space density S ( x , p ) :
D ( r , K ) = S ( x 1 , K ) S ( x 2 , K ) d 4 ρ = S ( ρ + r / 2 , K ) S ( ρ r / 2 , K ) d 4 ρ ,
where instead of the single particle variables x and p the pair-variables appear – the pair center of mass four-vector ρ = ( x 1 + x 2 ) / 2 , the pair separation four-vector r = x 1 x 2 , and the average momentum K = ( p 1 + p 2 ) / 2 . The D ( r , K ) distribution is the quantity that can be reconstructed indirectly from femtoscopic momentum correlation measurements, and experiments usually investigate the source-parameters that describe the shape of this distribution, see details, for example, in [15]. It was recently shown by different experiments that for pions this pair-source exhibits a power-law behavior, and can be described with a Lévy-stable distribution [15,16]. In case of spherical symmetry, the symmetric, centered stable distribution is defined as
L ( r ; α , R ) = 1 ( 2 π ) 3 d 3 q e i q r e 1 2 | q R | α ,
where the temporal dimension is removed from the dependence, made possible by the mass-shell condition, as detailed in Ref. [15]. The two important parameters that describe such a distribution are the Lévy-scale parameter R and the Lévy-exponent α . One can think of the latter as the parameter that is responsible for “how far” the distribution is from the Gaussian. In the α = 2 case L ( r ; α , R ) is identical to a Gaussian distribution, while in case of α < 2 it exhibits a power-law behavior. An illustration of the shape of such distributions can be seen on Figure 3. Since this distribution retains the same α exponent under convolution of random variables, if the single-particle source densities have a Lévy-shape then it follows that the two-particle source will also have such a shape, only the scale-parameter will be different. This can be summarized via the three-dimensional single-particle source S ( x ) and the pair-source D ( r ) (where we now suppress the momentum dependence, which is in turn contained in the parameters of the distribution):
S ( x ) = L ( x ; α , R ) D ( r ) = L ( r ; α , 2 1 / α R ) .
There are already multiple experimental measurements for the Lévy source parameters. Most notably the PHENIX experiment published results [15] for 0–30% centrality, s N N = 200 GeV Au+Au collisions, in the m T region of 0.23 GeV/ c 2 to 0.87 GeV/ c 2 . They found the α parameter to be very slightly dependent on m T , with an average value around 1.2. The NA61/SHINE experiment also has measured the Lévy-exponent in 150 A GeV Be+Be collisions, and obtained average α values of around 1.2 [16].

4. Analysis Method

For the analysis presented below we used s N N = 200 GeV Au+Au events generated by EPOS359. Using like-sign pion pairs, we measured the one-dimensional pair-source distribution in the longitudinal co-moving system (LCMS). The LCMS pair-separation vector can be expressed in lab-frame single-particle coordinates as
r L C M S = x 1 x 2 , y 1 y 2 , z 1 z 2 β ( t 1 t 2 ) 1 β 2 , where β = p z , 1 + p z , 2 E 1 + E 2 .
Using this variable, one can construct the spatio-temporal distance distribution D ( r L C M S , t ) . After angle- and time-integration we obtain the one-dimensional distance distribution as:
D ( r L C M S ) = d Ω L C M S d t D ( r L C M S , t ) ,
where we now suppress the K dependence indicated in Equation (1). Note that the dependence on the lab-frame time-coordinate disappears after the time integral of Equation (5), since we only keep the dependence on r L C M S , the final variable. In fact when analyzing the EPOS output, we only calculate the number of pairs in a given r L C M S bin, hence dependence on all other coordinates is naturally integrated out. When selecting pions we used the single-particle rapidity and transverse momentum requirements of | η | < 1 and 0.2   G e V / c < p T < 1.0 G e V / c . For each individual event we constructed the D ( r L C M S ) distribution for 5 different average transverse momentum k T classes in equal bins ranging from 0.2 to 0.4 GeV / c . Note that k T = 0.5 K x 2 + K y 2 is the transverse component of K used in Equation (1). We chose this k T region to be around the peak of the pair k T distribution to have adequate statistics (number of pairs) in the individual k T bins. To investigate centrality dependence as well, we separated the measurements to the centrality classes of 0–5%, 5–10%, 10–20%, 20–30%. In total we used 63,000 EPOS events, 10,500 for the first two centrality classes and 21,000 for the rest.
As mentioned before, EPOS has different stages of evolution including hydrodynamic expansion and hadronic rescattering. In order to identify the effect of the different stages on the shape of the pair-source distribution, as well as the contribution from the resonance decay products, we separated our investigation to four different cases as follows:
(a)
CORE with only primordial pions;
(b)
CORE with primordial + decay pions;
(c)
CORE+CORONA+UrQMD with only primordial pions;
(d)
CORE+CORONA+UrQMD with primordial + decay pions,
where primordial pions include pions coming from the thermal medium, that is, primordial pions are those that are not decay products. For each single event we fitted a Lévy distribution to the constructed D ( r L C M S ) distribution within the range of 2 fm to 100 fm. The fit was considered good if the confidence level calculated from the χ 2 and N D F values was greater than 0.1%. An example of such fits for the four different cases can be seen on Figure 4.
In case (a) we found that the events exhibit sharp cutoff features and mostly can be fitted well with a Gaussian. In case (b) the inclusion of the decay product pions results in power-law like structures, appearing at different regions in r L C M S . The shape of the events can be very different depending on the number (and origin) of decay pions in the sample. Due to the increased fluctuations and different event shapes the event-by-event extraction of the source parameters could not be done in the previous two cases. The fit settings would have to be fine-tuned for each event separately for this to work, which makes it impossible to do for thousands of events. In case (c) and (d) however, distinct non-Gaussian structures (power-law tails) are present in all events, shapes can be described by Lévy distributions in a statistically acceptable manner, and the extraction of the event-by-event source parameters is feasible. Furthermore, their distribution in the event sample can also be determined.
An example for such a distribution can be seen on Figure 5. In this example the distribution was reconstructed from fitting 21,000 events, out of which 18,460 fits were successful for case (c) and 18,768 for case (d) according to the confidence level criteria. Note that the distribution of source parameters was quite the same for non-acceptable fits as well; however, those do not necessarily represent the acquired source distributions, hence they were omitted from the further calculations. As mentioned before, we measured these two-dimensional R vs. α distributions for four different centrality classes, and five different k T regions. From these we can extract the mean and standard deviation values, and investigate their centrality and k T dependence. There are multiple ways to determine the mean and standard deviation parameters—on one hand, one could do normal distribution fits to the obtained 2D histograms, and on the other hand, one can simply calculate the first and second momenta of the distributions. We checked both and since the results were quite similar, for the sake of simplicity we chose the latter one.
Let us reiterate the point here that we analyzed individual EPOS events, and determined the R and α parameters of the pair source distribution in those individual events. We saw Lévy-shaped distributions when we included hadronic scattering, with slightly different α parameter values when decay pions were also included (besides primordial pions).

5. Results and Discussion

As described above, we performed fits to individual EPOS events from the final stage of EPOS (CORE+CORONA+UrQMD), and investigated averages of the resulting R and α parameters. We repeated this exercise for various centralities (0–5%, 5–10%, 10–20% and 20–30%) and k T regions (five equal bins from 0.2 to 0.4 GeV/c). We analyzed two cases separately: first the case of using only primordial pions, and then a case where both primordial and decay pions were included in the sample. Results for the mean R and α values versus m T = m 2 + k T 2 are shown in Figure 6.
One can observe a clear decreasing trend in R with both m T and centrality, for both the case of primordial pions, as well as the case where decays were also included. This trend with m T is similar to the observed R Gauss 2 m T trend observed universally across collision centrality, particle type, colliding energy, and colliding system size [11,35], even though it is based on Gaussian source radii. The decrease of R with increasing centrality shows the relation of the Lévy-scale to the initial fireball size. One can also observe that R is only weakly affected by the inclusion of decay pions; the values are slightly higher in the latter case.
The Lévy-stability index α shows less prominent centrality dependence, although a small decrease for more peripheral events is visible in case we include decay pions as well. This feature, that is, the centrality dependence of α , was not yet investigated in experimental publications (with the exception of preliminary data in Ref. [36]). One may also observe a weak decrease with m T similarly to what was observed in Ref. [15]. Furthermore it is clearly visible that when decay pions are also included, the α parameter decreases. This is expected as decay pions produce an even stronger tail, creating a smaller α value.
Concrete values of R can be compared to measured values from Ref. [15]. There R values of 7–8 fm were measured for the 0–30% centrality class and the m T = 0.25–0.45 GeV/ c 2 window. Our calculation yields similar values in the 10–20% and 20–30% centrality class. In Ref. [15], however, α values around 1.2 with a weak m T dependent decrease were found. The trend in m T is similar in our analysis, but the magnitude of the α values is somewhat different. The reasons for this could be multifold, they can range from the unavoidable event averaging present in the experiment to initial or final state effects not present in our simulations. The exploration of this difference is beyond the scope of present paper.
Note that the filled bands on both plots indicate the standard deviation of the R and α distributions over the investigated event sample (let us remind the reader that we investigate and fit pair source distributions in individual events). The statistical uncertainty of these data points is basically negligible, due to the large number of events this average was performed over. This means that the trends observed in Figure 6 and discussed above are true features of the EPOS event sample investigated in this paper.

6. Summary and Conclusions

We investigated a sample of EPOS events individually, and using identical pion pairs we reconstructed the pair source function in every individual event. In the case when only primordial pions were analyzed before hadronic scattering, a Gaussian shape was observed. However, when decay pions were also included, already power-law like structures appeared, and after the inclusion of hadronic scatterings (via UrQMD), Lévy-shaped pair distributions arose in the individual events. It is hence clear that it is not the event averaging that creates the non-Gaussian features in the pair distributions (and the arising correlation functions in femtoscopical measurements).
Subsequently, using the final stage of EPOS events (CORE+CORONA+UrQMD), we analyzed the event sample mean of the event-by-event Lévy-scale R and Lévy-index α values for the case of using only primordial pions and the case of including both primordial and decay pions. We observed clear trends as a function of m T and centrality for both cases. These observations show that in a realistic hydrodynamics-based simulation deviations from the Gaussian source shape appear when one includes hadronic scattering and decays. The values and trends of R are compatible with experimentally measured values, although we did not perform a detailed data comparison here. The weak decrease of α with m T is similar to what was observed experimentally, but the values of α are somewhat larger than measured values.
In the future we plan to utilize similar techniques to explore the dependence of these results on particle species as well. Further investigations might also include expanding the analysis to multiple dimensions, different collision energies, as well as reconstructing femtoscopical correlation functions.
Finally let us note that one of the important conclusions of our analysis is that the Lévy-shaped source assumption provides an acceptable description of the pion pair-source in EPOS3.

Author Contributions

Conceptualization, M.C.; formal analysis, D.K.; resources, M.S.; software, D.K., M.S.; writing—original draft preparation, D.K., M.C., M.S.; writing—review and editing, M.C., D.K.; visualization, D.K.; supervision, M.C. All authors have read and agreed to the published version of the manuscript.

Funding

Our research has been partially supported by the Hungarian NKFIH Grants No. FK-123842, K-138136, and 2019-2.1.11-TÉT-2019-00080. D.K. was also supported by the ÚNKP-21-4 New National Excellence Program of the Hungarian Ministry for Innovation and Technology from the source of the National Research, Development and Innovation Fund. We also acknowledge the usage of computer cluster DWARF at Warsaw University of Technology supported by the Polish National Science Center (NCN) under Contracts No. UMO-2017/26/E/ST3/00428 and UMO-2017/27/B/ST2/02792. These studies were also funded by IDUB-POB-FWEiTE-1 and IDUB-POB-FWEiTE-3 project granted by Warsaw University of Technology under the program Excellence Initiative: Research University (ID-UB). M.S. also wishes to acknowledge the financial support of the German Academic Exchange Service (DAAD).

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank H. Zbroszczyk, K. Werner and M.I. Nagy for useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The data from Figure 6 is listed in Table A1.
Table A1. Mean and standard deviation values extracted from the source parameter distributions measured in EPOS3 Au+Au collisions at s N N = 200 GeV. For the average transverse momentum k T the range is indicated, while for the average transverse mass m T the central value is given. The correlation coefficients and the covariance values are also indicated.
Table A1. Mean and standard deviation values extracted from the source parameter distributions measured in EPOS3 Au+Au collisions at s N N = 200 GeV. For the average transverse momentum k T the range is indicated, while for the average transverse mass m T the central value is given. The correlation coefficients and the covariance values are also indicated.
Primordial Pions
Centr. k T [GeV/c] m T [GeV/ c 2 ] R [fm] σ R [fm] α σ α cor( R , α )cov( R , α ) [fm]
0–5%0.20–0.240.2619.3810.3221.6700.056−0.225−0.004
0.24–0.280.2959.2220.3531.6410.056−0.273−0.005
0.28–0.320.3319.0210.3701.6210.056−0.305−0.006
0.32–0.360.3688.8110.3681.6100.056−0.319−0.007
0.36–0.400.4058.6190.3611.6030.056−0.324−0.007
5–10%0.20–0.240.2618.8920.4251.6690.064−0.328−0.009
0.24–0.280.2958.6630.4191.6450.063−0.326−0.009
0.28–0.320.3318.4310.4011.6300.062−0.329−0.008
0.32–0.360.3688.2230.3851.6210.062−0.332−0.008
0.36–0.400.4058.0490.3721.6140.061−0.316−0.007
10–20%0.20-0.240.2618.1210.5211.6710.071−0.229−0.008
0.24–0.280.2957.8880.5041.6500.070−0.235−0.008
0.28–0.320.3317.6830.4801.6360.068−0.230−0.008
0.32–0.360.3687.5100.4631.6270.067−0.224−0.007
0.36–0.400.4057.3700.4441.6200.066−0.208−0.006
20–30%0.20-0.240.2617.3500.5021.6550.078−0.142−0.006
0.24–0.280.2957.1350.4711.6390.076−0.128−0.005
0.28–0.320.3316.9710.4381.6270.073−0.137−0.004
0.32–0.360.3686.8440.4071.6190.072−0.119−0.003
0.36–0.400.4056.7490.3901.6100.071−0.114−0.003
Primordial + Decay Pions
Centr. k T [GeV/ c ] m T [GeV/ c 2 ] R [fm] σ R [fm] α σ α cor( R , α )cov( R , α ) [fm]
0–5%0.20–0.240.2619.3850.3011.5860.049−0.287−0.004
0.24–0.280.2959.2690.3181.5620.045−0.224−0.003
0.28–0.320.3319.1010.3291.5420.045−0.265−0.004
0.32–0.360.3688.9110.3261.5290.045−0.277−0.004
0.36–0.400.4058.7320.3201.5200.045−0.309−0.004
5–10%0.20–0.240.2618.9950.3751.5730.054−0.311−0.006
0.24–0.280.2958.8050.3681.5500.051−0.292−0.005
0.28–0.320.3318.6050.3551.5340.051−0.305−0.005
0.32–0.360.3688.4130.3441.5220.050−0.309−0.005
0.36–0.400.4058.2400.3361.5140.050−0.313−0.005
10–20%0.20–0.240.2618.2840.4721.5570.063−0.236−0.007
0.24–0.280.2958.1090.4601.5370.061−0.237−0.007
0.28–0.320.3317.9220.4461.5230.059−0.250−0.007
0.32–0.360.3687.7440.4331.5140.059−0.253−0.006
0.36–0.400.4057.5860.4211.5080.059−0.263−0.007
20–30%0.20–0.240.2617.5000.5171.5440.074−0.252−0.010
0.24–0.280.2957.3050.4891.5270.071−0.256−0.009
0.28–0.320.3317.1320.4651.5160.070−0.261−0.008
0.32–0.360.3686.9820.4371.5090.068−0.256−0.008
0.36–0.400.4056.8560.4111.5050.067−0.242−0.007

References

  1. 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]
  2. Lednicky, R. Femtoscopy with Unlike Particles; International Workshop on the Physics of the Quark Gluon Plasma, 2001; Available online: https://arxiv.org/abs/nucl-th/0112011 (accessed on 17 February 2022).
  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. 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–312. [Google Scholar] [CrossRef]
  5. Akkelin, S.V.; Sinyukov, Y.M. The HBT-interferometry of expanding inhomogeneous sources. Z. Phys. C 1996, 72, 501–507. [Google Scholar] [CrossRef]
  6. Akkelin, S.V.; Sinyukov, Y.M. The HBT interferometry of expanding sources. Phys. Lett. B 1995, 356, 525–530. [Google Scholar] [CrossRef]
  7. Csörgo, T.; Lorstad, B.; Zimányi, J. Quantum statistical correlations for slowly expanding systems. Phys. Lett. B 1994, 338, 134–140. [Google Scholar] [CrossRef]
  8. Csizmadia, P.; Csörgo, T.; Lukacs, B. New analytic solutions of the nonrelativistic hydrodynamical equations. Phys. Lett. B 1998, 443, 21–25. [Google Scholar] [CrossRef]
  9. Csörgo, T.; Lorstad, B. Bose-Einstein correlations for three-dimensionally expanding, cylindrically symmetric, finite systems. Phys. Rev. C 1996, 54, 1390–1403. [Google Scholar] [CrossRef]
  10. 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]
  11. PHENIX Collaboration; Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Amirikas, R.; Aphecetche, L.; Aronson, S.H.; et al. 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]
  12. STAR Collaboration; Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Badyal, S.K.; Bai, Y.; et al. Pion interferometry in Au+Au collisions at S(NN)**(1/2) = 200-GeV. Phys. Rev. C 2005, 71, 044906. [Google Scholar] [CrossRef]
  13. PHENIX Collaboration; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Al-Jamel, A.; Aoki, K.; Aphe-cetche, L.; Armendariz, R.; et al. Source breakup dynamics in Au+Au Collisions at s(NN)**(1/2) = 200-GeV via three-dimensional two-pion source imaging. Phys. Rev. Lett. 2008, 100, 232301. [Google Scholar] [CrossRef] [PubMed]
  14. PHENIX Collaboration; Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Amirikas, R.; Aphecetche, L.; Aronson, S.H.; et al. Evidence for a long-range component in the pion emission source in Au + Au collisions at s(NN)**(1/2) = 200-GeV. Phys. Rev. Lett. 2007, 98, 132301. [Google Scholar] [CrossRef]
  15. PHENIX Collaboration; Adare, A.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Alexander, J.; Alfred, M.; Al-Ta’ani, H.; Angerami, A.; et al. Lévy-stable two-pion Bose-Einstein correlations in sNN=200 GeV Au+Au collisions. Phys. Rev. C 2018, 97, 064911. [Google Scholar] [CrossRef]
  16. Pórfy, B. Lévy HBT Results at NA61/SHINE. Universe 2019, 5, 154. [Google Scholar] [CrossRef]
  17. Csörgo, 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]
  18. Csanád, M.; Csörgo, T.; Nagy, M. Anomalous diffusion of pions at RHIC. Braz. J. Phys. 2007, 37, 1002–1013. [Google Scholar] [CrossRef]
  19. Csörgo, 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]
  20. Csörgo, 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]
  21. Cimermaň, J.; Plumberg, C.; Tomášik, B. The Shape of the Correlation Function. Phys. Part. Nucl. 2020, 51, 282–287. [Google Scholar] [CrossRef]
  22. Cimerman, J.; Plumberg, C.; Tomášik, B. The Shape of the Correlation Function. POS 2021, ICHEP2020, 538. [Google Scholar] [CrossRef]
  23. Werner, K.; Karpenko, I.; Pierog, T.; Bleicher, M.; Mikhailov, K. Event-by-Event Simulation of the Three-Dimensional Hydrodynamic Evolution from Flux Tube Initial Conditions in Ultrarelativistic Heavy Ion Collisions. Phys. Rev. C 2010, 82, 044904. [Google Scholar] [CrossRef]
  24. Drescher, H.J.; Hladik, M.; Ostapchenko, S.; Pierog, T.; Werner, K. Parton based Gribov-Regge theory. Phys. Rep. 2001, 350, 93–289. [Google Scholar] [CrossRef]
  25. Stefaniak, M. Beam Energy Scan Dependence of Elliptic and Triangular Flow of Identified Hadrons in the STAR Experiment and the EPOS Model. Ph.D. Thesis, Warsaw University of Technology, Warsaw, Poland, 2021. [Google Scholar]
  26. Andersson, B. The lund string model. In Proceedings of the 7th European Symposium on Antiproton Interactions: From LEAR to the Collider and Beyond, Durham, UK, 9–13 July 1984. [Google Scholar]
  27. Ferreres-Solé, S.; Sjöstrand, T. The space–time structure of hadronization in the Lund model. Eur. Phys. J. C 2018, 78, 983. [Google Scholar] [CrossRef] [PubMed]
  28. Werner, K.; Guiot, B.; Karpenko, I.; Pierog, T. Analysing radial flow features in p-Pb and p-p collisions at several TeV by studying identified particle production in EPOS3. Phys. Rev. C 2014, 89, 064903. [Google Scholar] [CrossRef]
  29. Werner, K. Core-corona separation in ultra-relativistic heavy ion collisions. Phys. Rev. Lett. 2007, 98, 152301. [Google Scholar] [CrossRef]
  30. Borsanyi, S.; Endrodi, G.; Fodor, Z.; Jakovac, A.; Katz, S.D.; Krieg, S.; Ratti, C.; Szabo, K.K. The QCD equation of state with dynamical quarks. JHEP 2010, 11, 077. [Google Scholar] [CrossRef]
  31. Cooper, F.; Frye, G. Comment on the Single Particle Distribution in the Hydrodynamic and Statistical Thermodynamic Models of Multiparticle Production. Phys. Rev. D 1974, 10, 186. [Google Scholar] [CrossRef]
  32. Bleicher, M.; Zabrodin, E.; Spieles, C.; Bass, S.A.; Ernst, C.; Soff, S.; Bravina, L.; Belkacem, M.; Weber, H.; Stöcker, H.; et al. Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model. J. Phys. G 1999, 25, 1859–1896. [Google Scholar] [CrossRef]
  33. Bass, S.A.; Belkacem, M.; Bleicher, M.; Brandstetter, M.; Bravina, L.; Ernst, C.; Gerland, L.; Hofmann, M.; Hofmann, S.; Konopka, J.; et al. Microscopic models for ultrarelativistic heavy ion collisions. Prog. Part. Nucl. Phys. 1998, 41, 255–369. [Google Scholar] [CrossRef]
  34. Stefaniak, M. Examination of Heavy-ion Collisions Using EPOS Model in the Frame of BES Program. Acta Phys. Polon. Supp. 2018, 11, 695. [Google Scholar] [CrossRef]
  35. PHENIX Collaboration; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Al-Jamel, A.; Aoki, K.; Aphecetche, L.; Armendariz, R.; et al. Kaon interferometric probes of space-time evolution in Au+Au collisions at s(NN)**(1/2) = 200-GeV. Phys. Rev. Lett. 2009, 103, 142301. [Google Scholar] [CrossRef] [PubMed]
  36. Lökös, S. Centrality dependent Lévy-stable two-pion Bose-Einstein correlations in s N N = 200 GeV Au+Au collisions at the PHENIX experiment. Universe 2018, 4, 31. [Google Scholar] [CrossRef]
Figure 1. An illustration of the nucleon-nucleon rescattering with two projectiles A and two targets B. The splitting between participants and remnants shows the momentum sharing between constituents ensuring conservation of the given variable [25].
Figure 1. An illustration of the nucleon-nucleon rescattering with two projectiles A and two targets B. The splitting between participants and remnants shows the momentum sharing between constituents ensuring conservation of the given variable [25].
Entropy 24 00308 g001
Figure 2. The flux tube with the kink. The flux tube is mainly longitudinal but the kink part of the string moves transversely, here in the y-direction [25].
Figure 2. The flux tube with the kink. The flux tube is mainly longitudinal but the kink part of the string moves transversely, here in the y-direction [25].
Entropy 24 00308 g002
Figure 3. Lévy-stable source distributions with S ( r ) = L ( | r | ; α , R ) for α = 2, 1.5, and 1. The dependence on R is scaled out.
Figure 3. Lévy-stable source distributions with S ( r ) = L ( | r | ; α , R ) for α = 2, 1.5, and 1. The dependence on R is scaled out.
Entropy 24 00308 g003
Figure 4. EPOS3 example single event fits for 10–20% Au+Au collisions at s N N = 200 GeV. The average transverse momentum k T of the pion pairs is within the range of 0.28 GeV/c to 0.32 GeV/c. The measured D( r L C M S ) pion pair source distributions are denoted with downward-pointing triangles, and the fitted Lévy distributions are plotted with a red line (continuous in the fit region, dashed outside). The four panels from left to right, top to bottom are (a) CORE with only primordial pions, (b) CORE with primordial + decay pions, (c) CORE+CORONA+UrQMD with only primordial pions, and (d) CORE+CORONA+UrQMD with primordial + decay pions.
Figure 4. EPOS3 example single event fits for 10–20% Au+Au collisions at s N N = 200 GeV. The average transverse momentum k T of the pion pairs is within the range of 0.28 GeV/c to 0.32 GeV/c. The measured D( r L C M S ) pion pair source distributions are denoted with downward-pointing triangles, and the fitted Lévy distributions are plotted with a red line (continuous in the fit region, dashed outside). The four panels from left to right, top to bottom are (a) CORE with only primordial pions, (b) CORE with primordial + decay pions, (c) CORE+CORONA+UrQMD with only primordial pions, and (d) CORE+CORONA+UrQMD with primordial + decay pions.
Entropy 24 00308 g004
Figure 5. Example source parameter distributions for EPOS3 CORE+CORONA+UrQMD with only primordial pions (left) and with primordial and decay pions (right). The distribution was reconstructed from 10–20% centrality s N N = 200 GeV Au+Au events for pion pairs with an average transverse momentum k T between 0.28 GeV/c and 0.32 GeV/c.
Figure 5. Example source parameter distributions for EPOS3 CORE+CORONA+UrQMD with only primordial pions (left) and with primordial and decay pions (right). The distribution was reconstructed from 10–20% centrality s N N = 200 GeV Au+Au events for pion pairs with an average transverse momentum k T between 0.28 GeV/c and 0.32 GeV/c.
Entropy 24 00308 g005
Figure 6. Mean R and α values vs. average transverse mass m T for four different centrality classes (0–5%, 5–10%, 10–20%, 20–30%), in case of s N N = 200 GeV Au+Au collisions generated by EPOS3. The left column corresponds to the case of using primordial pions only, while the right column corresponds to the case of including both primordial and decay pions in the sample. The mean values are plotted with different filled markers. The corresponding colored boxes are representing the standard deviation values. The data from this Figure is listed in Table A1 for both cases.
Figure 6. Mean R and α values vs. average transverse mass m T for four different centrality classes (0–5%, 5–10%, 10–20%, 20–30%), in case of s N N = 200 GeV Au+Au collisions generated by EPOS3. The left column corresponds to the case of using primordial pions only, while the right column corresponds to the case of including both primordial and decay pions in the sample. The mean values are plotted with different filled markers. The corresponding colored boxes are representing the standard deviation values. The data from this Figure is listed in Table A1 for both cases.
Entropy 24 00308 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/e24030308

AMA Style

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(3):308. https://doi.org/10.3390/e24030308

Chicago/Turabian Style

Kincses, Dániel, Maria Stefaniak, and Máté Csanád. 2022. "Event-by-Event Investigation of the Two-Particle Source Function in Heavy-Ion Collisions with EPOS" Entropy 24, no. 3: 308. https://doi.org/10.3390/e24030308

APA Style

Kincses, D., Stefaniak, M., & Csanád, M. (2022). Event-by-Event Investigation of the Two-Particle Source Function in Heavy-Ion Collisions with EPOS. Entropy, 24(3), 308. https://doi.org/10.3390/e24030308

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