Collision energy dependence of source sizes for primary and secondary pions at NICA energies

We study the evolution of the two-pion correlation function parameters with collision energy in the context of relativistic heavy-ion collisions within the NICA energy range. To this end, we perform UrQMD simulations in the cascade mode to produce samples of pions from $5\times 10^6$ Bi+Bi collisions for each of the studied energies. The effects of the quantum-statistical correlations are introduced using the correlation afterburner code CRAB. We fit the correlation function using Gaussian, exponential and symmetric L\'evy shapes and show that for all collision energies the latter provides the best fit. We separate the sample into pions coming from primary processes and pions originating from the decay of long-lived resonances, and show that the source size for the latter is significantly larger than for the former. The source size for the secondaries, is similar but in general larger than the size for the whole pion sample. To further characterize the pion source, we also simulate the effects of a non-ideal detector introducing a momentum smearing parameter, representing the minimum pair momentum and thus a maximum source size that can be resolved. The values of the correlation function intercept parameter are therefore modified from the values they attain for the perfect detector case. Using the core-halo picture of the source, we show that the values of the intercept parameter are influenced by the presence of a significant fraction of core pions coming from the decay of long-lived but slow-moving resonances. These findings serve as a benchmark to compare with future Monte Carlo studies that consider an Equation of State and thus allow for a phase transition within the studied energy domain.


Introduction
Two-particle correlation studies have become a prime tool to determine the size and lifetime, as well as to infer the global properties, of the strongly interacting systems produced in relativistic heavy-ion collisions [1][2][3][4][5][6][7][8][9][10].Since the two-particle correlation function is related to the Fourier transform of the spatio-temporal component of the phase-space density of the emitting source, measurements of the correlation function provide access in particular to the space-time features of this source.Pions are by far the most abundant particles produced in these collisions.Therefore, it is common to perform correlation studies by experimentally measuring the two-pion correlation function.The technique is closely related to the photon intensity interferometry measurements introduced by Hanbury-Brown and Twiss to determine the size of stellar objects [11,12].In the context of collisions of hadron systems, Goldhaber, Goldhaber, Lee and Pais [13] developed a similar technique to study the interaction region formed in these systems.Lednicky [14] coined the term femtoscopy to emphasize that the technique is in this case applied to measurements at the fermi or femtometer scale.Femtoscopic studies have nowadays become a sophisticated and ever more precise tool, both experimentally as well as theoretically.
From the phenomenological point of view, the twopion correlation function can be described in terms of a set of parameters that contain not only information about the space-time size of the pion emitting source but also about the kind of processes that drive pion production.These parameters include the average source size  and the intercept  at vanishing relative momentum.Since particle production processes depend on the collision energy, it is expected that such parameters likewise evolve with energy.Of particular importance is to find whether the variation of these parameters show signals of criticality in the range of collision energies, where strongly interacting matter undergoes a transition from hadron to quark and gluon-dominated degrees of freedom.In this sense, femtoscopic studies provide also a tool for the exploration of the QCD phase diagram [15,16].Experimentally, a promising energy domain where signals of criticality can be found is planned to be scanned by the Multipurpose Detector (MPD) [17] at the Nuclotron-based Ion Collider fAcility (NICA), currently under construction at the Joint Institute for Nuclear Research (JINR).
When the correlation function is described as a function of the relative invariant pair momentum for a fixed average pair momentum, it is common to assume that the shape of the correlation function can be parametrized in terms of a Gaussian.Since the Fourier transform of a Gaussian is also a Gaussian, this assumption provides a simple description of the space-time source from the information obtained in momentum space.However, in recent times, it has become clear that this simple parametrization is not adequate and that a better description is achieved if the correlation function is parametrized based on a source described by a symmetric Lévy distribution [18][19][20].This is based on the realization that particle sources may show a large tail in configuration space and thus a description in terms of a distribution containing only one characteristic length may not be appropriate.The symmetric Lévy distribution is a generalization of the Gaussian distribution where the exponent  is called the Lévy exponent.The case  = 2 corresponds to the Gaussian distribution and  = 1 corresponds to a Cauchy distribution.
Two-particle correlation studies based on Lévy shape fits have been performed to describe experimental data for a single, or simultaneously at most for a couple of collision energies, mainly in the large energy domain.These studies correspond to SPS [21], RHIC [22][23][24][25][26] and LHC [27] energies for nucleusnucleus (A-A) collisions and also to LHC [28] energies for proton-proton collisions.A recent compilation in this energy range for the A-A case has been reported in Ref. [29].However, these studies have not yet been performed, neither fitting experimental data nor at Monte Carlo (MC) level, for energies where the putative Critical End Point (CEP) is thought to exist.In this work we aim to set up these studies within the NICA energy range where the CEP could be found.We report on correlation studies using MC generated data for central Bi+Bi collisions, which is one of the beam species planed to be used at the startup of the MPD data taking.The collision energies considered are within the range that is planned to be explored by the MPD experiment: √    = 4.0, 5.8, 7.7 and 9.2 GeV.For our present purposes, the MC generated data does not include an Equation of State (EoS) and thus it does not consider a possible phase transition at a given energy.In this sense, this work represents a benchmark study that we aim to use to compare with future studies that will consider an EoS in the MC generated data.To contrast the commonly used shapes that describe the correlation functions, these are fitted with Gaussian, exponential and symmetric Lévy shapes and the fits are compared.We show that, as it happens for the large energy domain, in the considered range and for all the studied cases, the symmetric Lévy shape provides the best fits.From the fits we extract the source parameters  and , and in the case of the Lévy fit, also  and find their dependence with the collision energy.We also separate the pion sample into pions coming from the decay of long-lived resonances, (secondary) and pions coming from direct processes (primary) and analyze the parameters describing their Lévy shape fits.We show that, as expected, primary pions come from a source with a smaller average size than secondary pions.However, the source size for secondary pions, although similar, is in general larger than the source size for the whole pion sample.To further characterize the source, we then introduce in the MC simulation the effects of a smeared momentum determination, that can be translated into a minimum pair momentum resolution, to extract the fraction of pions that come from the core using the core-halo picture of the source [30][31][32][33].This momentum smearing hampers resolving particles coming from the halo, and thus affects the extracted values of the intercept  in the correlation function.In this case, which simulates a real detector with a finite momentum resolution, the core-halo picture identifies the intercept  with the square of the fraction of pions coming from the core [30,34].We show that, for the energy range considered, the core still contains a significant fraction of pions that we interpret as coming from long-lived but slow-moving resonances and thus is not made up mainly from primary particles.We point out that future comparisons of MC studies that include an EoS with our results, can help to identify the emergence of particle production processes that change the relative abundance of primary and secondary pions in the core that in turn could signal critical behavior, which can be identified analyzing the properties of the Lévy stability index  as the collision energy changes.
The work is organized as follows: In Sec. 2 we describe the generalities of the two-pion correlation function and the implementation of the MC simulation.In Sec. 3 we fit the MC generated correlation functions with Gaussian, exponential and symmetric Lévy shapes and extract the correlation parameters.In Sec. 4, we separate the source into two components; one corresponding to primary and another corresponding to secondary pions and analyze the two samples finding the average source sizes for each of these sets.We also show an example of the source image for a collision producing a pion source at a time where 90% of the pions have attained freeze-out and compare the average size of the freeze-out hyper-surface with the Lévy size, finding an excellent agreement.We then introduce a smeared momentum resolution to simulate the effects of a real detector.We find the values of the intercept parameter and use the core-halo picture to link this parameter with the square of the fraction of core pions.We find that the core contains a significant component of secondary pions that we interpret as coming from long-lived but slow-moving resonances.We finally summarize and conclude in Sec. 5.

Monte Carlo simulations for the two-pion correlation function
The two-pion correlation function is defined as where  1 and  2 are the single-and two-pion momentum distributions, respectively, and  1 and  2 are the four-momenta of each of the pions.Usually, the two-pion correlation function is analysed as a function of the relative four-momentum  =  1 −  2 and for fixed values of the average pair momentum  = 1 2 (  1 +  2 ).The two-pion correlation function can also be related to the pion emitting source function in phase space, (, ), by noting that if dynamical correlations, such as Coulomb or strong final state interactions are neglected, the single-and two-pion momentum distributions can be written as where Ψ  and Ψ  1 ,  2 are the single-and two-pion symmetrized wave functions.
With the purpose of identifying possible non-Gaussian structures in the correlation function, which becomes difficult in a three-dimensional study, we perform a one-dimensional analysis in terms of the variable  inv = √︃  2 0 − | ì | 2 [34].Then, the two-pion correlation function can be written as where S is the Fourier transform of .If the relative momentum between the particles is much smaller than the average pair momentum, then the two-pion correlation function can be written as [35,36] which is a function of  inv , for a fixed average pair momentum .Since the simulations require a large amount of events to produce a statistically significant sample, hereby we consider all possible values of  that contribute to a given  inv .In this sense, our results have to be regarded as describing an effective source that contains all possible sizes corresponding to all values of the pair momenta.The resulting correlation function can thus be parametrized by different shapes, for instance Gaussian, exponential or Lévy shapes [20].
From the previous equation, it can be seen that the twopion correlation function can reach a maximum value of 2 at zero relative momentum, where the correlation function intercept,  2 ( inv → 0), is usually, denoted by  2 ( inv → 0) = 1 + , with  also known as the chaoticity or intercept parameter.However, different effects such as final state interactions and a finite experimental resolution, can prevent the intercept parameter from reaching the value 1 and can be understood in terms of the core-halo picture, whereby particles that come from the decays of long-lived resonances create a component of the source with a size that may not be resolved when the corresponding width of the pair momentum difference becomes smaller than the detector resolution.Accounting for this possibility the phase space emitting source can be modeled as consisting of two components  =  core +  halo , where each component has a Fourier transform and the core is composed of pions that come from primary processes.Notice that and hence, S(0) =  core +  halo .Thus, for experimentally resolvable values of the relative momentum, it can be assumed that S( inv ) ≃ Score ( inv ).Therefore, the  1: Parameters resulting from fits of the two-pion correlation function to exponential, Gaussian and Lévy forms for the complete pion set obtained at √    = 5.8 GeV.
two-pion correlation function can be expressed as As a consequence, in the core-halo picture, we can identify for any given pair invariant average momentum [20,30] On the other hand, the previous definition of the two-pion correlation function, Eq. ( 1), is barely used when dealing either with experimental or MC produced data.Instead of Eq. ( 1), one usually defines the measured two-pion correlation function as [36,37] where () is the relative momentum distribution of pions created in the same event which contains the quantum statistical effects, and () is the relative momentum distribution of pions created in different events, which does not contain Bose-Einstein correlations.N  and N  are normalization factors for  and .
q inv [MeV] In this work, we obtain the relative momentum distributions by means of MC simulations of relativistic heavy-ion collisions, from the Ultra-relativistic Quantum Molecular Dynamics Model (UrQMD) [38,39] in the cascade mode, which is used to simulate five million Bi+Bi central collisions at different collision energies within the NICA range ( √    = 4.0, 5.8, 7.7 and 9.2 GeV) with the simulation stopped at a time 200  fm.Since UrQMD does not include quantum statistical correlations [40,41], these are added by the formalism included in the "correlation after-burner" (CRAB) analyzing program [42], which uses the phase space distributions at their freeze-out positions to implement correlation weights.CRAB performs a boost for each of the pion pair momenta to the pair center of mass frame, where the pair wave function is easier to symmetrize, to then produce the quantum correlation.This is expressed as a function of the relative momentum which can then in turn be written in terms of  inv .

Correlation function fits and parameters
In general, it is assumed that the phase space source distribution, , can be factorized into a space-time distribution and momentum distribution.The Fourier transform of the space-time part, which is often referred to as the characteristic function, is assumed to be an analytic function around zero relative momentum and its second order Taylor expansion characterizes its behavior, even for large values of  inv [20].Thus, the two-pion correlation function can be approximately written as where  inv is the characteristic size of the source.Eq. ( 11) will be referred to as the Gaussian form.The previously mentioned assumptions can be translated into the stochastic nature of the several (independent) pion emission process.If one assumes that there are many independent processes that shift the emission position and that the final production point is a sum of many, similarly distributed, random shifts whose variance is finite, then according to the Central Limit Theorem, the distribution tends to a Gaussian.However, if the processes are characterized by large fluctuations originating power-like tails and a nonanalytic behaviour of the characteristic function, then the limiting distribution is not a Gaussian but instead a Lévy distribution.A special case of this distribution is the one named symmetric stable Lévy distribution.Utilizing such a distribution for the source function, the correlation function is written as where  is called the index of stability and can be related to the correlation critical exponent of QCD [43][44][45].Equation ( 12) will be referred to as the Lévy form of the correlation function.As a last example of the different functions that can describe the two-pion correlation function, we consider the exponential shape Equation ( 13) will be referred to as the exponential form of the correlation function.

Source features
If the source can be thought of as consisting of two components, one made of pions coming from primary and another one from secondary processes, on average the distance from the center of the fireball, where these particles are produced, is different.This is the core-halo picture of particle production.In this section we test whether the simulated two-particle correlation functions can be analyzed separating the pion sample into primary and secondary pions, with the former �sNN [GeV]

Complete set Primaries Secondaries
Fig. 4: Source size  inv as a function of the collision energy, √    , for the complete set of pions (black circles), primary pions (red rectangles) and secondary pions (blue diamonds), obtained from Bi+Bi collisions with impact parameter  = 0 − 1 fm. inv is obtained from fits to a Lévy form.coming from the core and the latter from the halo.As we proceed to show, this does not happen.To see this, we first study the case of a perfect resolution detector, and show that for the collision energies studied, although primary pions do indeed come from a small size source, the source size for secondary pions is similar albeit in general larger than the overall source size, which may seem puzzling.The picture is further elucidated when introducing a finite momentum resolution.
As we also show, the fraction of pions coming from the core contains a significant portion of secondary pions which we attribute to the decay product of long-lived but slow-moving resonances.In this work we include all the resonances available in UrQMD.However the resonances that give rise to most of the secondary pions are , Δ, ,  * ,  (1440), Δ(1600), Δ(1700), Δ(1950), Σ,  1 and (1700) [31].Most of these resonances have not such a long life-time.However they are produced all over the life-time of the system.It is in this sense that we consider them to be "long-lived".The picture that emerges is that the overall source size is the average between a small size core, containing a large population of secondaries, and a larger size halo.We extract the corresponding fireball parameters from fits to the central value of the correlation function in each bin.The rationale is to simulate the case of a very large statistical sample.We study the evolution of these parameters as a function of the collision energy �sNN [GeV]

Complete set Primaries Secondaries
Fig. 5: Lévy index of stability  as a function of the collision energy, √    , for the complete set of pions (black circles), primary pions (red rectangles) and secondary pions (blue diamonds), obtained from Bi+Bi collisions with impact parameter  = 0 − 1 fm. is obtained from fits to a Lévy form.

Ideal resolution case
Figure 1 shows as an example the case for √    = 5.8 GeV.The resulting parameters from fits to Gaussian, exponential and Lévy shapes are shown in Table 1.Notice that the fit that better describes the correlation function is obtained with the Lévy shape.For all the studied shapes, from the fits we should obtain  = 1.Nevertheless, the obtained values slightly deviate from 1 due to the limited statistics.The deviation from 1 within errors is an indication of how good the sample is for statistical purposes.By increasing the statistics  should tend to 1 and the error bar become smaller.
The results of the source size obtained from the fit to the Lévy form for the complete set of pions can be compared to the source image obtained from the freeze-out coordinates ( FO , ì  FO ) given by UrQMD.An example of the three-dimensional source image is shown in Fig. 2. The image is obtained in the center-ofmass system for a collision energy with √    = 9.2 GeV from the spatial freeze-out coordinates given by UrQMD at time  = 35 fm when 90% of the charged pions have frozen out.The source is nearly spherical in coordinate space.A long tail with a small density of q inv [MeV]
frozen out pions can also be inferred from this image.The average space-time interval of the pion freeze-out coordinates Δ = √︁ ⟨(FO 2 )  − (|ì  | 2 FO )  ⟩ coincides, within a few percent, with the Lévy fit value.We interpret this as the average width of the freeze-out hyper-surface, divided by a correction factor that accounts for the three-dimensional nature of this hyper-surface needed to compare with the one-dimensional Lévy analysis.
To explore the origin of pions that populate the core and the halo, we first perform the separation of the sample into primary and secondary pions.We accomplish this using the UrQMD parent process identification, whereby secondary pions are those produced as a result of decays of long-lived resonances.Figure 3 shows the separate contribution from primary and secondary pions to the the total correlation function shown in Fig. 1.Since the best fit is obtained using a Lévy form, hereafter we consider only the parameters obtained from this kind of fit, which are shown in Table 2.The results are consistent with the picture whereby primary pions come mainly from the core, since their source has a size significantly smaller than that of secondary pions.
Figures 4 and 5 show the evolution of the source radii and of the Lévy index of stability, respectively, with the collision energy for the complete, primary and secondary sets of pions.Notice that for most of the studied energies the following hierarchy of source size parameters holds:  inv, prim <  inv, all <  inv, second .except for the largest energy considered where we have instead  inv, prim <  inv, second <  inv, all .This may seem somewhat confusing if one thinks that the larger portion of secondary pions should come from the halo, since they originate from long-lived resonances.We also notice that in general the following hierarchies for the intercept and Lévy index:  prim <  all <  second and  second ≲  all <  prim , and that the correlation function for primary pions has the closest behaviour to a Gaussian.For the whole sample of pions, there is a general tendency for  inv to increase as the collision energy increases.For the separate samples of primary and secondary pions, although  inv grows with the collision energy from the lowest to the largest energy considered, this growth is non-monotonic.The index  slightly decreases with the collision energy from the lowest to the largest energy considered but, except for the case of the primary pion sample, the overall decrease is non-monotonic.To summarize, for the perfect detector case, the general tendency for the source radii is to grow with the collision energy, while the Lévy index of stability does not have a clear general tendency.The values of  inv for secondary pions are very similar but in general larger than those of the complete set of pions for the whole energy range; the opposite happens in general for the index .

Finite resolution case
To account for finite resolution effects, recall that it has been reported that the NICA-MPD will have a minimum momentum resolution of about 1.5 % for particles with total momentum around 0.2 GeV [46].Therefore, the relative momentum resolution of MPD will be of about Δ  = 10 MeV.This effect can be included in our studies by fixing the smearing parameter of CRAB to 10 MeV.This means that the momentum  of a given particle will be assigned the same value, provided it lies in the range  ± Δ  /2, which in turn can be translated into a limit for the minimum pair momentum resolution.Our results show that for the studied energy range and for a non-perfect detector with a smearing, or equivalently, a track resolution Δ  = 10 MeV, the resolvable size, in invariant variables, is of order 10 fm and the sample contains a non-resolvable halo which, as opposed to the case at higher collision energies, is basically a smooth continuation of the core, that is to say, it is not made of a very largely separated component.The existence of this non-resolvable 2 / ndf Primary 3.426 ± 0.012 0.905 ± 0.003 1.709 ± 0.016 1.615 Secondary 6.591 ± 0.054 0.647 ± 0.005 1.447 ± 0.022 0.259 Table 4: Results of the fit to the two-pion correlation function with a Lévy form accounting for a finite resolution of the detector, setting a smearing of 10 MeV for primary and secondary pions at √    = 5.8 GeV.
component is signaled by a value  < 1.By means of the Heisenberg uncertainty relation, this finite resolution therefore translates into the fact that source components whose characteristic size is larger than  ≳ 1/Δ  , cannot be resolved [32] and we hereby refer to these as the halo.
Figure 6 shows the effect of this finite resolution on the two-pion correlation function, shown in Figure 1, together with a fit to a Lévy form.The results of the fit are shown in Table 3.Notice that the finite resolution has the effect of significantly diminish the source size (by about 23%) and the value of the intercept (by about 39%), while increasing the value of the Lévy stability index (by about 14%).
Figure 7 shows the effect of a finite resolution (smearing of 10 MeV) on the separation of primary and secondary pions shown in Figure 3, together with a fit to a Lévy form.The results of the fits are shown in Table 4. Notice that, once again, the finite resolution has the effect of diminishing the source size (by about 2% for primary pions and about 21%, for secondary pions) and the value of the intercept (by about 8% for primary pions and about 39% for secondary pions), while the value of the Lévy index decreases for primary pions (by about 8%) and increases for secondary pions (by about 10%).This means that the effect of the finite resolution is of the same order for the set containing all the pions than for the set containing only the �sNN [GeV]

Complete set Primaries Secondaries
Fig. 8: Source size  inv as a function of the collision energy, √    , for the complete set of pions (black circles), primary pions (red rectangles) and secondary pions (blue diamonds), obtained from Bi+Bi collisions with impact parameter  = 0 − 1 fm and a smearing of 10 MeV. inv is obtained from fits to a Lévy form.secondary pions, while the effect of a finite resolution barely affects the set of primary pions.
Figures 8 and 9 show the evolution of the source radii and the Lévy index, respectively, with the collision energy for the complete, primary and secondary sets of pions.The hierarchy of the fit parameters for the invariant radii is  inv, prim, FR <  inv, all, FR <  inv, second, FR , except for le lowest energy considered where we have instead  inv, prim, FR <  inv, second, FR ≲  inv, all, FR .The hierarchy for the intercept parameter becomes now  all, FR <  second, FR <  prim, FR , while it remains the same for the Lévy index  second, FR <  all, FR <  prim, FR .Notice that for the set of primary pions,  inv is basically constant for the energy range considered, whereas  shows a moderate decrease with energy.This means that a finite momentum resolution has a small effect on the pions of primary origin.The values for  inv and  for the set of secondary pions are closer to the corresponding parameters when the whole pion sample is considered.This means that when no separation of the sample between primary and secondary pions is made, the full sample is dominated by the secondary pions.For the whole sample of pions, as well as for the set of secondary pions,  inv shows an overall tendency to increase with the collision energy, whereas for both of these samples,  is basically constant around the same value  ≃ 1.5.
�sNN [GeV]     6: Average charged pion multiplicity evolution with the collision energy for the complete, primary and secondary sets of pions.
Table 5 shows the evolution of the intercept parameter  with respect to the collision energy for the complete set of pions, as well as for the primary and secondary pions.Notice that for all the cases, the value of  decreases as the energy increases.Nonetheless, for the complete set and the set of secondary pions, the decrease is around 10%, whereas for the set of primary pions the decrease is only marginal and around 2%.In the core-halo picture, the intercept parameter can be related to the square of the fraction of pions coming from the core.This implies, according to Eq. (9), that between 77% and 82% of the pions should come from the core and only a smaller fraction should come from the halo.According to our definition, pions of secondary origin are those coming from the decay or long-lived resonances.The results show that simulating a finite resolution detector with a smearing Δ ∼ 10 MeV, this fraction of pions of secondary origin together with the pions of primary origin come from a spacetime region within  inv , such that Δ inv ≲ 1 implying that the space region that can be accessed is restricted to  inv ≲ 20 fm.We have verified from the simulations that indeed, about 20% of the pions are produced in a region with a distance from the collision center larger than 20 fm.Therefore, the spatial region where about 20% of secondary pions come from, corresponds to an average size  halo ≳ 20 fm.This explains why, when considering a finite resolution detector, the core pions (as defined by the core-halo picture) do not contain this fraction of secondary pions and thus, the average size for the source of secondaries decreases from  inv ∼ 8 fm to  inv ∼ 6.5 fm.Since the intercept parameter for secondaries decreases more than for the primaries as the energy increases, this means that for larger energies, long-lived resonances decay further away from the center, as expected, and thus contribute less to the population of core pions.From the UrQMD simulation, it is possible to directly identify that the average fraction of primary pions, which always populate the core, increases marginally with energy to be between 6% and 12%.Therefore, we can conclude that for lower energies, the core is mainly populated by pions from resonance decays and that this population slowly decreases as the collision energy increases.For completeness, and in order to fully characterize the pion sample, Table 6 shows the average multiplicity per event for the complete, the primary and the secondary sets of pions.

Summary and Conclusions
In the context of relativistic heavy-ion collisions within the NICA energy range, we have studied the evolution with collision energy of the parameters that describe the two-pion correlation function.We have performed MC simulations using the UrQMD event generator in the cascade mode, to produce 5×10 6 events for each considered energy.In each case, the quantum correlations are included using the CRAB analyzing code.No other source of correlations but the quantum ones have been included.We studied the correlation function as a function of the invariant relative pair momentum for a fixed value of the average pair momentum.To find the parameters that describe the correlation function we performed fits using Gaussian, exponential and symmetric Lévy shapes.We have shown that, as is the case when considering larger energies, the Lévy shape provides that best description of the correlations for the different settings and across the considered energy range.The most likely origin of the Lévy behavior is that the particle source has a large tail in configuration space and thus a description in terms of a distribution containing only one characteristic length is not appropriate.This is confirmed by our study of the three-dimensional source image.The pion sample is separated into its primary and secondary components.The latter is defined as the set of pions coming from the decay of long-lived resonances, mostly , Δ, ,  * ,  (1440), Δ(1600), Δ(1700), Δ(1950), Σ,  1 and (1700).We have shown that the source size for the sample of secondaries is similar but larger than that for the whole sample and significantly larger than the source size of the primaries.The intercept parameter exhibits the same hierarchy whereas the Lévy index exhibits the opposite one.For the case of the primary pion sample, the Lévy index shows an overall slight decrease with collision energy, however the secondary and full set of pions do not show a clear tendency with collision energy.We also notice that when restricting only to the set of primary pions, the correlation function already requires  < 2. Adding to the sample of primary pions the pions coming from resonance decays produces in general a further decrease of  [47].
In order to obtain a more accurate picture of the space-time characteristics of the pion producing sources, we have simulated the case of a non-ideal detector introducing a smearing parameter in the CRAB code to mimic a minimum resolution for the determination of the relative pair momentum.From the uncertainty relation between momentum and position, this translates into a maximum source size from where pion pairs can be identified.This smearing produces that the intercept of the correlation function becomes smaller than 1.Within the core-halo picture, the impossibility to determine source sizes larger than the inverse of the smearing momentum can be turned into an advantage since the size of the intercept can be directly identified with the square of the fraction of pions coming from the core.Our results indicate that the core pion sample has a large component that comes from the decay of long-lived but slow-moving resonances, as well as a small component of pions coming from primary processes.The former decreases whereas the latter increases with collision energy.In this sense, the analysis of the relative abundance of pions in the core coming from resonance decays and from primary processes, as the collision energy changes, becomes more important as a tool to study signals of criticality within the NICA energy range, when future comparisons are made with results from similar analysis but using an event generator that includes a phase transition within the same energy range.Indeed, when additional particle producing processes within the core introduce extra sources of correlations with lengths of order of the size of the system, such as when an EoS is considered in the MC generator, it is expected that these leave an imprint that can show up in particular as a non monotonic evolution of the Lévy index with collision energy [47].In this sense, or results represent a benchmark to use for comparison with future studies where the event generator includes and EoS.We are performing this kind of analysis based on the findings of this work and the results will soon be reported elsewhere.

Fig. 1 :
Fig. 1: Two-pion correlation function for Bi+Bi collisions at √    = 5.8 GeV, with impact parameter  = 0 − 1 fm.The dots represent the output of CRAB which is compared with exponential, Gaussian and Lévy fits.

Fig. 2 :
Fig. 2: (a) Example of the three-dimensional source image in the center-of-mass system for an event with √    = 9.2 GeV obtained from the spatial freezeout coordinates given by UrQMD at time  = 35 fm when 90% of the charged pions have attained their freeze-out.The source is nearly spherical in coordinate space.(b) One-dimensional distribution of Δ  = √︁ ( 2 FO )  − (|ì  | 2 FO )  for each particle for a sample event obtained from the freeze-out coordinates given by UrQMD at time  = 200 fm.The average value is obtained as Δ = 9.19 fm.

Fig. 3 :
Fig. 3: Two-pion correlation function of primary and secondary pions, produced in Bi+Bi collisions at √    = 5.8 GeV, with impact parameter  = 0 − 1 fm.The dots represent the output of CRAB, while the solid lines represent the fit of Lévy forms.

Fig. 6 :
Fig. 6: Two-pion correlation function of Bi+Bi collisions at √    = 5.8 GeV, with impact parameter  = 0 − 1 fm.The blue dots represent the same output of CRAB of Figure 1, while the red dots include the finite resolution effect of MPD, with a smearing of 10 MeV.Solid lines represent the fit to a Lévy form.

Fig. 7 :
Fig. 7: Two-pion correlation function of primary and secondary pions, produced in Bi+Bi collisions at √    = 5.8 GeV, with impact parameter  = 0 − 1 fm.The dots represent the output of CRAB with a smearing of 10 MeV, while the solid lines represent the fit of Lévy forms.

Table 2 :
Parameters resulting from fits of the two-pion correlation function to a Lévy form for primary and secondary pions obtained at √    = 5.8 GeV.

Table 5 :
Evolution of the intercept parameter  with the collision energy for the complete, primary and secondary sets of pions, obtained from fits to a Lévy form accounting for a finite resolution of the detector with a smearing of 10 MeV.