Abstract
In nuclear collisions the incident protons generate a Coulomb field which acts on produced charged particles. The impact of these interactions on charged-pion transverse-mass and rapidity spectra, as well as on pion–pion momentum correlations is investigated in Au + Au collisions at \(\sqrt{s_\text {NN}}\) = 2.4 GeV. We show that the low-\(m_t\) region (\(m_t < 0.2\) GeV/\(c^2\)) can be well described with a Coulomb-modified Boltzmann distribution that also takes changes of the Coulomb field during the expansion of the fireball into account. The observed centrality dependence of the fitted mean Coulomb potential energy deviates strongly from a \(A_{\text {part}}^{2/3}\) scaling, indicating that, next to the fireball, the non-interacting charged spectators have to be taken into account. For the most central collisions, the Coulomb modifications of the HBT source radii are found to be consistent with the potential extracted from the single-pion transverse-mass distributions. This finding suggests that the region of homogeneity obtained from two-pion correlations coincides with the region in which the pions freeze-out. Using the inferred mean-square radius of the charge distribution at freeze-out, we have deduced a baryon density, in fair agreement with values obtained from statistical hadronization model fits to the particle yields.
Similar content being viewed by others
1 Introduction
The fireball produced in a relativistic heavy-ion collision has a net positive electric charge due to the protons from both target and projectile nuclei. Charged particles emitted from the expanding fireball are subjected to the long-range Coulomb interaction caused by this electric charge resulting in distortions of their emission spectra. Such modifications can be observed in the spectra of charged hadrons, ideally by comparing the effects on oppositely charged states. Positively and negatively charged pions are consequently ideal probes of the Coulomb effects: they have the same mass but opposite charge, and are produced copiously with comparable yields already in few-GeV collisions. The force exerted by the electric field accelerates the \(\pi ^+\) and decelerates the \(\pi ^-\) resulting in a change of their kinetic energies. The corresponding average energy kick is determined by the charge distribution of the incident protons, which create the Coulomb field, as well as by the distribution of pion emission points, i.e. the pion source. Investigating the Coulomb effects has the potential to reveal information on the characteristics of the pion source at freeze-out. An increase of the \(\pi ^-/\pi ^+\) yield ratio at low momenta was first observed in early Bevalac experiments [1,2,3] and was quickly interpreted as being related to the Coulomb field [4, 5]. Indeed, these effects turned out to be ubiquitous in heavy-ion collisions, ranging from SIS beam energies [6,7,8], to AGS [9, 10], SPS [11], and even up to RHIC energies, as discussed in an extensive report [10]. The early theoretical work done to interpret the Bevalac data [4, 5] was taken up again as the SIS, AGS, and SPS results became available [12,13,14,15,16,17,18]. In particular at low bombarding energies, the \(\pi ^{-}/\pi ^{+}\) ratio is expected to provide information on the symmetry term of the nuclear equation of state [19], thus requiring the Coulomb effect to be under full control [20]. Finally, the Coulomb field also influences the distributions of relative momenta used in two-pion interferometry [21,22,23,24] and this must be taken into account to properly interpret the measured HBT radii of charged-particle sources [13, 25,26,27,28,29,30].
The High Acceptance Di-electron Spectrometer (HADES) experiment [31] at the SIS18 accelerator has extensively studied particle production in 1.23 \(A\) GeV Au + Au collisions, equivalent to a center-of-mass energy of \(\sqrt{s_\text {NN}}\) = 2.4 GeV. In particular, charged-pion spectra were measured with very high statistics and the pion phase-space distribution could be reconstructed with good accuracy [32], limited mostly by the systematic uncertainties of the detector efficiency and acceptance of order 5–10%. Especially, the extrapolation of the measured differential particle yields into the low-\(m_t\) region, not fully covered by HADES, is a potential source of systematic uncertainties. In the past, Boltzmann distributions were often used for this extrapolation, neglecting the spectral distortions caused by the Coulomb interaction. One motivation for the analysis presented here was to provide a more accurate parameterization of pion distributions at low \(m_t\) (or low \(p_t\)), allowing to reduce systematic uncertainties due to the extrapolation of measured yields into this region. This is particularly important owing to the fact that the geometry of the HADES detector and its toroidal magnetic field lead to a somewhat different low-momentum cut-off of the \(\pi ^-\), bent to large polar angles, and the \(\pi ^+\), bent towards the beam axis [31, 32].
This paper is organized as follows. In Sect. 2 we introduce the formalism used to describe the Coulomb effect on charged-particle spectra and we validate our fit procedure by simulations. In Sect. 3, we present an improved analysis of the charged-pion spectra measured in central and semi-peripheral Au+Au collisions at \(\sqrt{s_\text {NN}}\) = 2.4 GeV. The extracted values of the average Coulomb potential energy \(V_C\) as well as the resulting source radius are discussed as a function of the collision centrality. We then compare these results with estimates of the Coulomb potential energy derived from the HBT radii of the measured like-sign pion-pion momentum correlations. From the pion source volume, we extract a baryon density at freeze-out which is compared with the Statistical Hadronization Model (SHM) fits to the particle yields in the same event sample. Using published data from other experiments, we then present the excitation function of \(V_C\) from 1 to 10 \(A\) GeV kinetic beam energy. In Sect. 4, finally, we summarize our findings and finish with a brief outlook.
2 Methodology
2.1 Coulomb field acting on charged pions
In a central collision at \(\sqrt{s_\text {NN}}\) = 2.4 GeV, the colliding Au nuclei are stopped in a fireball which afterwards expands radially with a velocity of about \(\langle \beta \rangle \approx \) 0.3–0.4 [33, 34], with the produced particles moving away in the Coulomb field of the positive net charge of the fireball. The influence of the Coulomb field is most noticeable in the velocity or momentum distributions of low mass charged particles, in particular pions: \(\pi ^{+}\) are sped up and \(\pi ^{-}\) are slowed down, causing substantial modification of their differential yields. Due to the resulting reshuffling, at low center-of-mass momenta, the \(\pi ^{+}\) yield is reduced whereas the \(\pi ^{-}\) yield is increased relative to the uncharged pions. In collisions in which the nuclei are not fully stopped their longitudinal motion as well as possible charged spectators need to be taken into account.
The nucleus–nucleus center-of-mass energy in the final state \(E_{f}\) of pions can be expressed in terms of their initial energy \(E_{i}\), i.e. the energy corrected for the presence of the electrostatic field, and the Coulomb potential energy \(V_{C}\) as
where \(p_f\) and \(p_i\) are the corresponding final and initial pion center-of-mass momenta, and the ± sign corresponds to the different pion charges. As the total charge of the fireball is due to the incoming proton charges, \(V_C\) is a positive quantity. In this simple static picture, the Coulomb potential energy leads to a shift of the pion total energy by an amount \(+V_C (-V_C)\) for positively (negatively) charged pions. However, as discussed in [14], the expansion of the charged fireball causes an attenuation of the Coulomb effect. Indeed, the central electric field felt by a given pion is only produced by those charges that are slower than this pion, i.e. the charged matter shell that has overtaken the pion in the outward expansion does not contribute. For a complete description of the expansion, a fully dynamic calculation would have to be used, e.g. within a hadronic transport model allowing for different freeze-out times and loci depending on the charge and momentum of the pions. Here, we follow instead a more data-driven ansatz similar to the one applied in [10, 14], where all pions are assumed to freeze-out at the same time, and parameterize the attenuation of the electrostatic force by replacing \(V_{C}\) with an effective potential \(V_{\text {eff}}\) expressed as a function of the pion kinetic energy \(E_{\pi }- m_{\pi }\).Footnote 1 Defining \(x = \sqrt{(E_{\pi }/m_{\pi } - 1) \, m_p/T_p}\), which is a measure for the relative velocities of pions and protons, and integrating over the proton velocity distribution, the authors of [14] found
where \(m_p\) is the proton mass, \(T_{p}\) is the inverse slope parameter from a fit to the proton \(m_{t}\) distribution, and \({{\,\mathrm{erf}\,}}(x)\) is the error function. The 2D case corresponds to a cylindrical geometry with a transverse expansion (i.e. for boost-invariant systems, relevant at high collision energies) and the 3D case stands for a spherical expansion. The latter one is more appropriate for the low beam energies at which HADES operates. The rationale behind Eq. (2) is the following: the protons provide the dominant part of the fireball net charge and the attenuation term is given by the fraction of protons that are slower than a given outgoing pion. Note that the proton energy distribution can be described alternatively by a blast-wave fit, using e.g. Siemens-Rasmussen [35], although in that case no closed formula is available for \(V_{\text {eff}}\) in Eq. (2). The integration over the proton velocity distribution would have to be carried out numerically, increasing massively the cost for repeatedly computing \(V_{\text {eff}}\) in the iterative fitting procedure. Ultimately, the attenuation of the Coulomb field as a function of pion energy is realized in our fit function by replacing the parameter \(V_{C}\) in Eq. (1) with \(V_{\text {eff}}\) of Eq. (2), and \(V_{\text {eff}}\) taking the sign of \(V_C\).
To extrapolate pion spectra into phase-space regions not covered by experiment, it is customary to use a relativistic Boltzmann distribution adjusted to the data. In Ref. [32], we had shown that a proper description of the transverse-mass spectrum requires in fact the sum of two Boltzmann distributions, expressed as
where E is the total energy of the pion in the center of mass,Footnote 2A is the normalization, \(T_{1}\) and \(T_{2}\) are the slopes of the two spectral components, and f and \(1-f\) are their corresponding fractional amplitudes (\(0 < f \le 1\)). Following Ref. [14], we insert the initial total pion energy \(E_{i}\) into Eq. (3) and, with the help of Eq. (1), express it as a function of the observed final energy \(E_{f}\) to obtain
This distribution contains the full Jacobian \(J \times J_{\mathrm {eff}}\) of the transformation from “initial” to “final” kinematic variables. It consists of two parts: the first factor J corresponds to the Jacobian proposed in Ref. [13] for a constant Coulomb potential energy, namely
The second factor \(J_{\mathrm {eff}}\) results from the explicit dependence of \(V_{\text {eff}}\) on the pion kinetic energy introduced via Eq. (2). As shown in Appendix A, it is expressed as
with x as defined above for Eq. (2).
Note that this second factor has been omitted in previous work, e.g. in [10, 14] and is here applied for the first time in the analysis. As demonstrated on simulated pion spectra in the next subsection, we find that the extracted values for the Coulomb potential energy would come out almost twice larger if this factor was not included.
2.2 Effect on the kinematic distribution of charged pions
In order to illustrate the effect of a Coulomb field on the kinematic distributions of charged pions, we have used the event generator Pluto [36] to simulate thermal spectra of all three pion species, \(\pi ^+\) \(\pi ^0\), and \(\pi ^-\). To do that, the Pluto code was modified to include the Coulomb potential energy, basically by implementing Eqs. (4)–(6) in the energy sampling routine. In the simulation, we have set the parameters to \(T_1 = 50\) MeV, \(T_2 = 90\) MeV, \(f = 0.95\), \(V_c = 15\) MeV, \(T_p = 130\) MeV, and used as well angular distribution coefficients set to \(A_2\) = 0.54/0.60/0.66 for \(\pi ^{+}/\pi ^0/\pi ^{-}\), which are typical values observed for the energy regime where HADES operates [32, 34]. The left panel of Fig. 1 depicts the simulated distributions around mid-rapidity (\(0.74 \pm 0.05\)) for \(\pi ^-\) and \(\pi ^+\) as a function of \(m_t-m_0\), compared to the pure \(\pi ^{0}\) case, i.e. without the effect of the Coulomb field, naturally realized by neutral pions. In order to focus on the Coulomb effect, isospin related differences were ignored, i.e. the three pion flavors were simulated with equal multiplicities. Fluctuations at the higher \(m_t\) values are caused by the limited event statistics. The lower panel depicts the ratios of charged pions with respect to the \(\pi ^{0}\) distribution. The effect of the Coulomb potential is most prominent at very low \(m_t-m_0\) leading to an enhancement of the \(\pi ^-\) and a depletion of the \(\pi ^+\) by almost 50%. At high \(m_t-m_0\), the effect is opposite for both \(\pi ^-\) and \(\pi ^+\), amounting to about 15%. Note that, for the determination of the total pion yield from measured data, the low part of the transverse-mass spectrum is of particular interest: it holds a substantial fraction of the yield, while often requiring extrapolation to correct for incomplete detector acceptance.
The right panel of Fig. 1 displays the simulated \(\pi ^-\) (blue) and \(\pi ^+\) (red) center-of-mass rapidity distributions dN/dy, again compared to the \(\pi ^0\) case (black); the lower panel depicts the ratios relative to the \(\pi ^{0}\) distribution. At mid-rapidity, the effect of the Coulomb field on the total pion yields amounts to about 5%. In the tails, however, the relative differences are larger, with the Coulomb force leading to a narrowing of the \(\pi ^-\) distribution and a broadening of the \(\pi ^+\) one relative to \(\pi ^{0}\).
2.3 Test of the fit procedure
From data, the parameter \(V_{C}\) can be obtained most robustly by adjusting Eq. (4) simultaneously to both the measured \(\pi ^{+}\) and \(\pi ^{-}\) transverse-mass (transverse-momentum) spectra. Before turning to measured data, however, we have cross checked the performance of our fitting procedure using a Monte Carlo simulation. Employing the event generator Pluto [36], we have set up a pion source according to Eq. (4) with the realistic values of the parameters \(T_{1}\), \(T_{2}\), f, and \(V_{C}\), as given before. In this simulation we have assumed a spherical expansion of the fireball, as described by Eq. (2) for the 3D case.
Figure 2 shows the simulated charged-pion \(dN/dm_t\) spectra at mid-rapidity, together with fits using the Coulomb-modified two-component Boltzmann distribution of Eq. (4). As already pointed out, the influence of the Coulomb field on the charged pions manifests itself mostly at low \(m_t - m_0\). From the figure one can see that the values of the fitted parameters agree within statistical errors with the input, validating our procedure. In order to illustrate the impact of using an attenuated potential, we present two cases: first, using a \(V_{\text {eff}}\) to account for 3D expansion of the fireball (Eq. (4)) and, second, ignoring the expansion, i.e. setting \(V_{\text {eff}}= V_{C}\) and \(J_{\text {eff}}= 1\). The two fits are displayed in Fig. 2 as solid and dashed colored curves, respectively. As the inserted close-up shows, the two fit curves differ substantially at very low \(m_t\), giving also largely different values of the extracted Coulomb potential energy \(V_C = 15.2\) MeV (\(\chi ^2/\mathrm {ndf} = 1.07\)) compared to \(V_C = 20.9\) MeV (\(\chi ^2/\mathrm {ndf} = 3.08\)). We investigated as well the fitting scheme proposed in Ref [10], realized by Eq. (4) with \(J_\mathrm {eff} = 1\), resulting in \(V_C = 32.8\) MeV (\(\chi ^2/\mathrm {ndf} = 1.38\)), i.e. more than double the value used to simulate the events.
2.4 Role of the spectator protons
In the picture of an expanding charged fireball, only the total charge of the participant nucleons contributes to the Coulomb effect. The implicit assumption is that the Coulomb field of the spectator protons is small, and that the latter move away fast enough to be neglected. As a consequence, \(V_C\) would scale with \(A_{\mathrm {part}}\) like the ratio of volume over radius of the charge distribution, i.e. \(V_C \propto Z_{\mathrm {part}}/R \propto A_{\mathrm {part}}^{2/3}\). However, in order to estimate the influence of spectator charges, we have implemented in our simulation the three-source ansatz of Gyulassy and Kauffmann [5] by defining in Pluto a charged fireball at rest in the center of mass and two receding charged spectators. In this scheme, the pions are emitted from the participant zone only, but are affected by the combined Coulomb field of all three charges. By fitting Eq. (4) to the generated pion spectra at mid-rapidity, the evolution with \(A_{\mathrm {part}}\) of the resulting effective Coulomb potential energy can be determined. The fitted parameter \(V_C\) is shown in Fig. 3 for calculations performed with and without the spectators included. While still not being dynamic, this simulation illustrates the trend of \(V_C\) with \(A_{\mathrm {part}}\), suggesting that, while spectator contributions are indeed small in central collisions, already in semi-peripheral and more so in peripheral collisions, the Coulomb potential energy deviates from the simple, central-source \(A_{\mathrm {part}}^{2/3}\) scaling. Evidently, the interpretation of Coulomb effects in peripheral events or off mid-rapidity will be more demanding than for central events. The influence of spectator charges on pion spectra has also been investigated in Refs. [17, 18], but for SPS energies only.
3 Results and discussion
3.1 Effect of the Coulomb field in Au + Au collisions
We have applied the fit procedure presented in the previous section to the charged-pion spectra measured with HADES in Au+Au collisions at \(\sqrt{s_\text {NN}}\) = 2.4 GeV [32]. In that experiment, events were selected online with a hardware trigger based on the multiplicity of hits registered in the HADES time-of-flight detectors. In total, \(10^{8}\) events in the \(0-40\%\) most central collisions were selected for the present analysis. By tracking the charged particles through the HADES magnetic field, their momentum was reconstructed and by recording their time of flight, a velocity was obtained. Particle identification was finally achieved by cutting on the characteristic momentum-velocity relation. More details of the event reconstruction, centrality selection, and pion identification are given in [32]. The rapidity acceptance of the detector is about 0.65 for charged pions with slight differences between low-momentum \(\pi ^{+}\) and \(\pi ^{-}\) caused by the toroidal geometry of the magnetic field. To obtain the total pion yields, extrapolation to zero \(m_t\) is required as well as (model-dependent) extrapolation outside of the acceptance. The uncertainties on the yields are dominated by systematics of the efficiency correction (3%) and the extrapolation to full solid angle (5–7%). The statistical uncertainties on the yields are negligible (see Ref. [32] for details).
We start by discussing the mid-rapidity charged pion transverse-mass spectra (\(|y_{\text {c.m.}}|< \) 0.05) obtained for the 10% most central events, where the centrality of the Au+Au collisions has been selected by cutting on the number of hits in the HADES time-of-flight detectors. The spectra for both charges have been fitted simultaneously in the range of 0.025 \(< m_{t}<\) 0.65 GeV/\(c^2\) assuming an attenuated potential \(V_{\text {eff}}\) to account for 3D expansion of the fireball (fixing \(T_p\) at 142, 125, 115, and 106 MeV for 0–10%, 10–20%, 20–30%, and 30–40% centrality classes, respectively). The data points and the combined fit are displayed in Fig 4 as well as the ratios of data to fit functions. Except for the lowest \(m_t\) bin of the negative pion distribution, the ratios differ by less than 5% from unity, demonstrating that the fit function based on Eq. (4) is adequate. The adjusted parameters \(T_1\), \(T_2\), and \(V_C\) are listed in Table 1, together with the average number of participants \(\langle A_\text {part}\rangle \) obtained from a Glauber calculation [37]. In particular, the mean Coulomb potential energy is found to be \(V_C = 13.6 \pm 0.6\) MeV. By integration of the adjusted functions, the pion yields and their fit errors are obtained, leading to a \(\pi ^-/\pi ^+\) yield ratio of \(2.05 \pm 0.10\) at mid-rapidity.
The influence of the Coulomb field on the charged pion spectra has been investigated with a BUU transport model in Ref. [38]. A similar theoretical study, done with the QMD model, can be found in [39]. The authors calculated a mean Coulomb potential energy in the range 20 – 30 MeV for 1\(A\) GeV Au+Au collisions, i.e. substantially larger than what we observe in our data. However, in the BUU calculation, the potential was determined by averaging over the full space-time of the heavy-ion collision as well as over the full rapidity range, whereas our result is obtained at mid-rapidity. The fits of the pion spectra used in the present work are most sensitive to low-momentum pions most likely emitted from \(\Delta \) resonances late in the expansion of the fireball and thus emerging from a dilute charge distribution. In contrast, the authors of Ref. [38] discuss only pions moving faster than the expanding charge cloud. Furthermore, in their model calculation, fast pions are emitted from higher-lying resonances and freeze out earlier, at a density of \(\approx 0.8 \rho _0\), seeing thus a much more compact charge distribution which translates into a correspondingly larger mean Coulomb energy.
3.2 Centrality dependence of the Coulomb potential
The mid-rapidity charged pion spectra from other centrality classes (10–20%, 20–30%, and 30–40% ) are fitted in a similar fashion. They are shown in Fig. 4 together with the corresponding data-over-fit ratios. The latter demonstrate again that, except for the lowest \(m_t-m_0\) bin, the data is very well described by Eq. (4). The extracted fit parameters are listed in Table 1 and their evolution with centrality is summarized in Fig. 5 where \(V_C\), the \(\pi ^{-}/\pi ^{+}\) ratio, and the inverse slope parameters \(T_1\) and \(T_2\) are plotted as a function of mean \(\langle A_{\mathrm {part}} \rangle \). One observes that the Coulomb potential energy \(V_C\) decreases smoothly from central to peripheral collisions. This decline goes in hand with the decreasing overlap volume, resulting in a smaller fireball and less net charge contributing to the Coulomb effect. As can be seen in the figure, the \(A_{\mathrm {part}}^{2/3}\) scaling of \(V_C\) is not fully realized, suggesting that spectator contributions are indeed present in the most peripheral events. Note, however, that the observed effect appears to be weaker than in the schematic model calculations [5] underlying Fig. 3. In Fig. 5, one furthermore sees that both inverse slope parameters increase with increasing centrality, whereas the pion ratio remains basically constant. To avoid complications due to spectator effects, we focus our further analysis on central collisions.
3.3 Pion rapidity densities and total yields
The Coulomb field affects not only the transverse momentum but also the rapidity density distributions. However, here the situation is more complicated due to the observed forward-backward peaked polar anisotropies of the emitted pions [32]. As discussed, the extrapolation and integration of the pion \(m_t\) spectra with the use of the Coulomb-modified fit function, the dN/dy distributions of \(\pi ^{+}\) and \(\pi ^{-}\) can be obtained inside the rapidity range covered by HADES. The left panel of Fig. 6 shows the integrated pion rapidity density of central events as a function of the center-of-mass rapidity. For comparison, the yields based on the Boltzmann fits done in a previous analysis [32] are shown as well. It is apparent that the inclusion of the Coulomb field affects the pion yields, decreasing \(\pi ^{+}\) and increasing \(\pi ^{-}\). To obtain the total pion multiplicities per event, an extrapolation of the data points outside of the HADES rapidity coverage is required. Total pion multiplicities are hence obtained by extrapolating the measured yields with the help of various transport calculations (see Ref. [32]), resulting in \(M(\pi ^{-}) = 17.5 \pm 1.0\) and \(M(\pi ^{+}) = 8.6 \pm 0.5\), superseding the \(17.1 \pm 1.2\) and \(9.3 \pm 0.7\) values obtained in Ref. [32] without considering the impact of Coulomb effects. The error bars given are systematic and are in both cases dominated by the extrapolation to full solid angle (for details see Ref. [32]). From these numbers, the ratio of the total pion multiplicities is found to be \(2.03 \pm 0.14\). This ratio can be compared to the predictions of the isobar model [40], namely 1.95 when pion production is mediated solely by \(\Delta \) resonance excitation, and 1.70 for production solely via \(N^{*}\) resonances. The observed ratio favors clearly a scenario where the \(\Delta \) dominates pion production, also in agreement with the findings of a detailed study of resonance excitation and decay in the Au+Au collision system [41].
The rapidity dependence of the measured pion yield ratio, symmetrized around mid-rapidity, is displayed in the middle panel of Fig. 6, together with a Pluto simulation using source parameters set to the values obtained from the \(m_t\) fit, and the polar angular anisotropies \(A_2(\pi ^{-})\) = 0.66 and \(A_2(\pi ^{+})\) = 0.54 obtained in Ref. [32]. The main features of the observed ratio, in particular the shape of the marked dip around mid-rapidity caused by an interplay of the pion intrinsic angular distribution and the Coulomb field are reproduced.
The right panel of Fig. 6 shows the dependence of the fitted Coulomb potential energy on the center-of-mass rapidity, symmetrized around mid-rapidity. It is apparent that \(V_C\) has a minimum at mid-rapidity, with the strong rise at backward and forward rapidities being caused by the presence of the target and beam spectator charges, respectively.
3.4 Geometry of the charge distribution
The Coulomb potential energy obtained by fitting Eq. (4) to the \(m_t\) spectra can be used to estimate the hadronic density of the pion emitting system at freeze-out. To do this, we have to make assumptions about the geometry of the charge distribution, on the one side, and of the distribution of pion emission points, i.e. shape of the pion source, on the other side. To get a handle on this problem, we assume that the pions are produced in a spherically symmetric fireball. Furthermore, we consider two extreme cases: (1) a constant charge density and (2) a Gaussian charge density.
Following Ref. [27], the electrostatic potential energy of a uniformly charged sphere of radius \(R_{\circ }\) and total charge Ze is obtained as
For volume emission, i.e. assuming that the pions decouple instantaneously from the fireball, the average potential energy \(\langle V_C \rangle \) is obtained by integrating the potential over a constant source distribution \(f({\mathbf {r}})\propto \Theta (R_{\circ }-r)\), where \(\Theta \) is the Heaviside function, such that
with \(Z = Z_{\mathrm {part}}\) being the number of participating protons, \(Z_{\mathrm {part}} = A_{\text {part}} \times Z_{\text {Au}}/A_{\text {Au}}\). For the 10\(\%\) most central Au+Au events, \(\langle A_{\text {part}} \rangle \) = 303 ± 12 and \(\langle Z_{\text {part}}\rangle =122 \pm 5\) [37]. Note that, doing so, we neglect the spectator charges, which is a valid approximation for the most central collisions (see Fig. 3 and discussion above). If, however, only surface emission is assumed, we must replace in Eq. (8) the Heaviside function with the Dirac delta \(f({\mathbf {r}})\propto \delta (r-R_{\circ })\), recovering the relation \(\langle V_C \rangle = e^2 \, Z_{\text {part}} / R_{\circ }\), as used e.g. in Ref. [7]. Thus, in the 10% most central collisions, from the fitted potential \(\langle V_C\rangle = 13.6~\pm ~0.6\) MeV, we find a source radius \(R_{\circ } = 15.4 \pm 0.8 ~\text {fm}\) (\(12.9 \pm 0.7\) fm) for volume (surface) emission, respectively.
In two-particle momentum correlation analyses, usually a different picture of the fireball at freeze-out is assumed, namely a 3D Gaussian density profile. The widths, which in the most general case can be different along the three spatial directions, are denoted by \(R_{\mathrm {out}}\), \(R_{\mathrm {long}}\) and \(R_{\mathrm {side}}\). Restricting ourselves to a spherically symmetric shape, the Coulomb potential energy of a pion emitted at a distance r from the center of a Gaussian charge distribution of width \(\sigma \) is given by [26]
Assuming that the pions are emitted from a source following the same density profile, their average potential energy is given by
With this expression, we obtain for the 10% most central collisions a Gaussian radius of \(\sigma = 7.2 \pm 0.4\) fm. As discussed in the next section, this value is also consistent with the radii obtained from a two-pion correlation analysis using HADES data [29, 30].
By equating the mean potentials obtained by both relations, (8) and (10), we can define an equivalent hard-sphere radius [42] of the Gaussian distribution as
In order to make a meaningful comparison between sizes of the fireball calculated in the hard-sphere and Gaussian cases, one can use instead the root mean square radii \(R_{\text {rms}}\) of the two density distributions. In the (volume emission) hard-sphere case one has
and in the Gaussian case
Consequently, expressing \(R_{\circ }\) in Eq. (8) and \(\sigma \) in Eq. (10) in terms of \(R_{\text {rms}}\), we find that the two density profiles lead to very similar (within \(\le 7\)%) r.m.s. radii for a given fitted potential \(V_C\). It shows that the Coulomb effect is not very sensitive to details of the charge and pion source distributions, which is due to the long-range nature of the Coulomb force.
3.5 Coulomb potential from HBT radii
The Coulomb field of the expanding fireball also acts on the relative momentum of like-sign pion pairs, affecting in a characteristic way the \(R_{\pi ^{-}\pi ^{-}}\) and \(R_{\pi ^{+}\pi ^{+}}\) radii extracted from HBT intensity-interferometry [13, 26, 27, 43, 44], namely
where \(q_\mathrm {i}\) (\(p_\mathrm {i}\)) is the initial relative (absolute) momentum of the pair, \(q_\mathrm {f}\) (\(p_\mathrm {f}\)) is its final, i.e. with Coulomb push, relative (absolute) momentum, and \(V_{\mathrm {eff}}\) is the effective Coulomb potential energy. Relative and absolute pair momenta are formed from the momenta of the two individual pions, \(p_\mathrm {1,(i,f)}\) and \(p_\mathrm {2,(i,f)}\) as \(p_\mathrm {i,f} = (p_\mathrm {1,(i,f)} + p_\mathrm {2,(i,f)})/2\) and \(q_\mathrm {i,f} = (p_\mathrm {1,(i,f)} - p_\mathrm {2,(i,f)})/2\).
Typically, one has \(V_{\mathrm {eff}} / k_\mathrm {t} \ll 1\), where \(k_t\) is the average transverse momentum of the two pions forming the pair, so that the second-order term \(V_\mathrm {eff}^2/p_\mathrm {f}^2\) is very small.Footnote 3 Neglecting this term, the source radii of constructed neutral-pion pairs (denoted here by \(\tilde{\pi }^{0} \tilde{\pi }^{0}\)) had been obtained in Refs. [29, 30] as the average of the squared charged radii
Likewise, by taking the difference of the squared radii, the effective potential \(V_{\text {eff}}\) can be determined as a function of \(m_t\), where \(m_t\) is the average pion transverse mass around mid-rapidity
Fitting the resulting \(V_{\text {eff}}(m_t)\) with Eq. (2), one obtains from the HBT radius another estimate of the Coulomb potential energy \(V_C\). The situation is, however, more complex in HBT interferometry because the full 3D analysis of the emitting source gives various radii, namely \(R_{\mathrm {long}}\), \(R_{\mathrm {side}}\), \(R_{\mathrm {out}}\), and, finally, \(R_{\mathrm {int}}\) which can be considered as an average over the three axes (see Ref. [30] for a definition). Moreover, the interpretation of these radii in terms of a volume of homogeneity rather than a geometric charge volume has to be kept in mind when comparing with the \(R_{\text {rms}}\) radius of the single-pion analysis.
Applying the procedure to the various charged-pion HBT radii published in Ref. [30], the \(V_\mathrm {eff}^\mathrm {HBT}\) values depicted in Fig. 7 were obtained as a function of the mean pion \(m_t\). Results are shown for the 0–10% centrality, a pion rapidity coverage of \(-0.35 \le y_{\text {c.m.}}\le 0.35\), and for \(m_\mathrm {t} < 250\) MeV\(/c^{2}\) only, as at larger transverse mass the sensitivity of the HBT approach to \(V_\mathrm {C}\) dwindles. All effective potentials, \(V^\mathrm {inv}_\mathrm {eff}\), \(V^\mathrm {side}_\mathrm {eff}\), and \(V^\mathrm {long}_\mathrm {eff}\), were found to range between 5 and 15 MeV. Just as the various radii are not equal, the corresponding \(V_\mathrm {eff}\) also differ. We have excluded \(V^\mathrm {out}_\mathrm {eff}\) from the discussion because of the more complex influence of the Coulomb field on \(R_{\text {out}}\) (discussed in Ref. [27]) which renders Eq. (16) inapplicable.
To study the centrality dependence of the potential, the fits were performed for all four measured centrality classes. The values of \(V^\mathrm {inv}_\mathrm {C}\), \(V^\mathrm {side}_\mathrm {C}\) and \(V^\mathrm {long}_\mathrm {C}\) are similar compared to each other with variations of less than 25% within one centrality class. Again, we do not show the ’out’ direction and focus here on ’side’ and ’long’ noting that \(V^\mathrm {side}_\mathrm {C}\) varies over all centrality classes by less than 20% while \(V^\mathrm {long}_\mathrm {C}\) increases from peripheral to central collisions by more than 50%. Both are plotted in Fig. 8 as a function of \(\langle A_\mathrm {part}\rangle ^{2/3}\) (green and blue triangles) together with the values obtained from the single-pion spectra (shown also in Fig. 5). The Coulomb potential energy resulting from both methods, single-particle spectra and HBT, are in reasonable agreement. In particular, the HBT ’long’ direction values are very close to the single-pion spectra result. The similarity of the Coulomb potential energies obtained from the transverse-mass spectra and from the HBT analysis supports the assumption that the region in which the pions freeze-out and the region of homogeneity of two-pion correlations overlap. In Fig. 8, a deviation between the expected \(A_\mathrm {part}^{2/3}\) scaling and the measured values of \(V_{C}\) is again visible towards the more peripheral collisions. In the figure this scaling is shown by the dashed red line normalized to the most central data point. We can only speculate that this behavior may be caused by the contribution of the spectator protons to the Coulomb field, as was already pointed out when discussing Fig. 3.
3.6 Net-baryon density and chemical potential
Having confirmed that the source radius extracted from single-particle spectra is in fair agreement with the HBT radii, we use this radius as an approximation of the fireball size and estimate the corresponding baryon density \(\rho _B\) at freeze-out. Setting the baryon number equal to \(A_{\mathrm {part}}=303\) and using the hard-sphere radius for volume emission, \(R_{\circ } = 15.4\pm 0.8\) fm, the average baryon density is found to be \(\langle \rho _{B} \rangle = A_{\mathrm {part}}/\frac{4}{3}\pi R_{\circ }^3=0.12 \pm 0.02~\rho _{0}\), where \(\rho _{0}\) is the nuclear saturation density. Alternatively, we can assume a Gaussian source distribution with a maximal density at the center (\(r = 0\)) of \(\rho _{B,max} = 0.32\pm 0.05~\rho _0\), corresponding to an average density \(\langle \rho _B \rangle = \rho _{B,max}/2^{3/2} = 0.11\pm 0.02~\rho _0\), very close to the hard-sphere value. Note that using a Gaussian source also implies volume emission of the pions. This is warranted because the Coulomb effects manifest themselves mostly at the very low momenta where pion re-absorption in the medium is known to be weak [45, 46]. With the \(V_C\) values extracted for the four centrality classes, we obtain the evolution of \(\rho /\rho _0\) with the number of participating nucleons \(A_{\text {part}}\) shown in Fig. 9(left). The power law \(\rho /\rho _0 = n_0 \, A_{\text {part}}^{\alpha }\) fitted to the extracted densities, shown as red solid line, yields \(\alpha \simeq -1/3\) and \(n_0 \simeq 1\).
The freeze-out density can be converted to the baryon chemical potential \(\mu _B\) of an ideal hadron resonance gas (HRG) of temperature T by the following relation (see e.g. [47]):
where, \(m_N\) is the nucleon mass and \(g=4\) is the degeneracy factor for nucleons. In order to apply this relation, the temperature T needs to be fixed. However, an unambiguous determination of T from the pion spectra is challenging due to distortions by resonance decays and collective expansion (see Ref. [48] for a recent study). Therefore, we calculate \(\mu _B\) with the temperature extracted from a systematic study of freeze-out points [49] where T and \(\mu _B\) have been parameterized as a function of the reaction center-of-mass energy. The resulting relation is displayed in Fig. 9(right), together with freeze-out points extracted in various SHM analyses of measured hadron yields in heavy-ion collisions over a range of collision energies spanning from 2 GeV up to 2.76 TeV [49,50,51,52,53,54]. The temperature corresponding to our collision energy is found to be \(T_{\text {HRG}} = 54.5 \pm 0.5\) MeV and, inserted into Eq. (17), it leads to a baryochemical potential of \(\mu _B = 772 \pm 10\) MeV. The resulting freeze-out point is indicated by the red filled symbol in the phase diagram. Note finally that the source radius \(R_{\circ } = 15.4\pm 0.8\) found in our analysis is also very close to the freeze-out radius \(R_{\text {f.o.}}\) typically extracted from SHM analyses of the measured hadron yields, e.g. Ref. [48] found \(R_{\text {f.o.}}\) = 16 fm.
3.7 Collision energy dependence of the Coulomb effect
To investigate the beam-energy dependence of the Coulomb effects, we also analyzed published mid-rapidity charged-pion spectra from the KaoS experiment at the SIS18 accelerator, from the E895 experiment at the AGS, and from the NA49 experiment at the SPS. The KaoS data [7] were measured in 1\(A\) GeV Au+Au collisions for the 14 \(\%\) most central events at a fixed laboratory polar angle of \(\theta _{\mathrm {lab}} = 44^{\circ } \pm 4^{\circ }\), corresponding to a narrow, but non-uniform center-of-mass acceptance bin located around \(y_{\text {c.m.}}\simeq 0.15\). Figure 10 shows a fit of the pion center-of-mass kinetic-energy distributions \(d\sigma /dE^{\mathrm {kin} }_{\mathrm {c.m.}}\) with an adaptedFootnote 4 version of Eq. (4), setting \(T_p = 120\) MeV. The resulting Coulomb potential energy is \(V_C = 15.6 \pm 2.2\) MeV and the energy-averaged \(\pi ^{-}/\pi ^{+}\) ratio is \(2.20\pm 0.10\). In case of E895, the published data correspond to the 5\(\%\) most central Au+Au collisions measured at beam energies of 2, 4, 6, and 8\(A\) GeV, respectively [55]. Finally, the analyzed NA49 data were collected for the 7% most central Pb+Pb collisions at 20 and 30\(A\) GeV [56]. When extending our analysis to higher bombarding energies (i.e. up to 30\(A\) GeV), the assumption of a spherical fireball may not anymore be applicable. At high energies, the shape of the particle source close to mid-rapidity is usually assumed to be cylindrical, implying a 2D expansion rather than the 3D one expressed in Eq. (2). To account for this effect, we consider both cases, 2D and 3D, in our analysis. Figure 11 shows the fits done to the E895 and NA49 mid-rapidity charged-pion \(m_t\) spectra using Eq. (4), with \(V_{\mathrm {eff}}\) corresponding to a 3D fireball up to 8\(A\) GeV and 2D above; for 6 and 8\(A\) GeV, both results are presented. All fit parameters and their statistical uncertainties are listed in Table 2. For the E895 data, the \(T_2\) slopes were kept fixed to their published values as the limited \(m_t\) acceptance for positive pions of this experiment did not allow stable fits with both slopes varying freely.
One can see from Figs. 10 and 11 that the Coulomb-modified fit function gives overall a good description of the available data sets. Furthermore, when comparing all presented pion \(m_t\) distributions, it appears that, going from low to high beam energy, the spacing between negative and positive pion yields decreases. This can be attributed to the gradual decrease of both the \(\pi ^{-}/\pi ^{+}\) ratio and \(V_\mathrm {eff}\) with increasing bombarding energy.
Combining the results from all fits, the collision-energy dependence of the extracted Coulomb potential energy can be investigated for central events. Although the centrality selections of the various experiments are not identical, ranging between 0–5% and 0–14%, Fig. 12 shows that \(V_C\) decreases overall with increasing \(\sqrt{s_\text {NN}}\). Such a trend is expected as, with rising energy, the center-of-mass momenta of the incident protons will be increasingly focused into the longitudinal direction, diluting the electric field felt by the pions at mid-rapidity. On the other hand, the limited accuracy of the AGS data precludes to firmly decide whether the 3D or 2D expansion scenario matches better the general trend of the fitted Coulomb energy with \(\sqrt{s_\text {NN}}\). Finally, note that an upper limit on \(V_C\) can be estimated as given by the Coulomb potential energy of a charged sphere corresponding to two fully overlapping gold nuclei, i.e. \(V_C < 1.44 \, Z/R_{\circ } \simeq 33\) MeV, where \(Z=2 \times 79\) and \(R_{\circ } \simeq 7\) fm.
4 Summary
In this publication we have addressed the Coulomb interaction exerted on charged pions by the cumulated electric charge of the expanding fireball produced in relativistic heavy-ion collisions. We have illustrated the observable effects with charged-pion spectra measured by the HADES spectrometer in \(\sqrt{s_\text {NN}}\) = 2.4 GeV Au + Au reactions.
Fitting the pion \(m_t\) spectra with a Coulomb-modified two-slope Boltzmann distribution, we could reduce the systematic uncertainties arising in the determination of the total pion yields. Besides the \(m_t\) inverse-slope parameters \(T_1\) and \(T_2\), those fits in addition deliver information on the average Coulomb potential energy experienced by the produced particles. Comparisons with results obtained by the former KaoS, E895, and NA49 experiments show that the latter quantity follows a smooth trend with center-of-mass energy.
Relating the Coulomb potential energy to the spatial geometry of the expanding fireball, a freeze-out density can be calculated and, within the statistical hadronization model, constraints on the baryochemical potential can be obtained. Assuming volume emission of the pions, from either a hard-sphere or a Gaussian source, we deduced a range of baryon densities of 0.11–0.12 \(\rho _{0}\) and a baryochemical potential \(\mu _{B}\) of 772 ± 10 MeV, consistent with values found in various SHM analyses of measured hadron yields.
The cumulative electric charge of the participant nucleons also affects the two-pion momentum correlations, as already pointed out in the HBT analysis of our Au+Au data presented in [29, 30]. We have shown that these second-order Coulomb manifestations are compatible with the effects seen and quantified in the transverse-mass distributions of the charged pions.
The present analysis of high-statistics charged-pion spectra has revealed the importance of low-momentum pions for a better understanding of the expansion phase of a heavy-ion collision. These are pions moving slower than or with the expanding charged matter. This finding may encourage new fully dynamical calculations, e.g. transport simulations going beyond the assumptions of the present analysis, to elucidate further the role of low-momentum pions for the dynamics governing the freeze-out of the fireball produced in relativistic heavy-ion collisions.
Data Availability Statement
The manuscript has associated data in a data repository. [Authors’ comment: The datasets generated during and/or analyzed during the current study are available from the corresponding author on request.]
Notes
In a classical approximation, the average kinetic energy of thermal protons is \(\langle E^p_{\text {kin}} \rangle = 3/2 \, k \, T_p\), e.g. for an effective proton temperature \(T_p = 130\) MeV, \(\langle E^p_{kin} \rangle = 195\) MeV. The pions of the same velocity, i.e. with a kinetic energy \(E^{\pi }_{kin} = m_{\pi }/m_p \times 195 = 29\) MeV, experience a Coulomb field attenuated by a factor \(1 - e^{-1.5} \simeq 3/4 (3/5)\) in the 2D (3D) case.
Ignoring the Coulomb force, \(E = E_f = E_i\).
In our analysis, the smallest used average \(k_\mathrm {t}\) value is about 80 Mev/c, leading to \(V_{\mathrm {eff}} / k_\mathrm {t} < 1/5\). Furthermore, one has \(p_\mathrm {f} = \sqrt{k_\mathrm {t}^2 + k_\mathrm {l}^2} > k_\mathrm {t}\).
For this fit, the KaoS cross-section data was transformed from its kinetic-energy representation at fixed laboratory angle into a transverse-mass representation.
References
W. Benenson, G. Bertsch, G.M. Crawley, E. Kashy, J.A. Nolen, H. Bowman, J.G. Ingersoll, J.O. Rasmussen, J. Sullivan, M. Koike et al., Phys. Rev. Lett. 43, 683–686 (1979). https://doi.org/10.1103/PhysRevLett.43.683
S. Nagamiya, M.C. Lemaire, E. Moller, S. Schnetzer, G. Shapiro, H. Steiner, I. Tanihata, Phys. Rev. C 24, 971–1009 (1981). https://doi.org/10.1103/PhysRevC.24.971
J.P. Sullivan, J.A. Bistirlich, H.R. Bowman, R. Bossingham, T. Buttke, K.M. Crowe, K.A. Frankel, C.J. Martoff, J. Miller, D.L. Murphy et al., Phys. Rev. C 25, 1499–1517 (1982). https://doi.org/10.1103/PhysRevC.25.1499
K.G. Libbrecht, S.E. Koonin, Phys. Rev. Lett. 43, 1581–1584 (1979). https://doi.org/10.1103/PhysRevLett.43.1581
M. Gyulassy, S.K. Kauffmann, Nucl. Phys. A 362, 503–533 (1981). https://doi.org/10.1016/0375-9474(81)90507-8
C. Müntz [E802, E866 and KaoS], Acta Phys. Polon. B 29 (1998), 3253-3257 [arXiv:nucl-ex/9807001 [nucl-ex]]
A. Wagner, C. Müntz, H. Oeschler, C. T. Sturm, R. Barth, M. Cieslak, M. Debowski, E. Grosse, P. Koczon, M. Mang, et al. Phys. Lett. B 420, 20–24 (1998). https://doi.org/10.1016/S0370-2693(97)01549-9. arXiv:nucl-ex/9712004 [nucl-ex]
B. Hong et al., [FOPI]. Phys. Rev. C 71, 034902 (2005). https://doi.org/10.1103/PhysRevC.71.034902. arXiv:nucl-ex/0501021 [nucl-ex]
J. Barrette et al., [E877]. Phys. Rev. C 62, 024901 (2000). https://doi.org/10.1103/PhysRevC.62.024901. arXiv:nucl-ex/9910004 [nucl-ex]
D. Cebra, S. G. Brovko, C. E. Flores, B. A. Haag, J. L. Klay, arXiv:1408.1369 [nucl-ex]
H. Bøggild et al., [NA44]. Phys. Lett. B 372, 339–342 (1996). https://doi.org/10.1016/0370-2693(96)00114-1
B. A. Li, Phys. Lett. B 346, 5–8 (1995). https://doi.org/10.1016/0370-2693(95)00010-I. arXiv:nucl-th/9501001 [nucl-th]
G. Baym, P. Braun-Munzinger, Nucl. Phys. A 610, 286C–296C (1996). https://doi.org/10.1016/S0375-9474(96)00363-6. arXiv:nucl-th/9606055 [nucl-th]
H. W. Barz, J. P. Bondorf, J. J. Gaardhoje, H. Heiselberg, Phys. Rev. C 57, 2536–2546 (1998). https://doi.org/10.1103/PhysRevC.57.2536. arXiv:nucl-th/9711064 [nucl-th]
A. Ayala, J. I. Kapusta, Phys. Rev. C 56, 407–411 (1997). https://doi.org/10.1103/PhysRevC.56.407. arXiv:hep-ph/9701401 [hep-ph]
A. Ayala, S. Jeon, J. I. Kapusta, Phys. Rev. C 59, 3324–3328 (1999). https://doi.org/10.1103/PhysRevC.59.3324. arXiv:nucl-th/9807062 [nucl-th]
A. Rybicki, A. Szczurek, Phys. Rev. C 75, 054903 (2007). https://doi.org/10.1103/PhysRevC.75.054903. [arXiv:nucl-th/0610036 [nucl-th]]
V. Ozvenchuk, A. Rybicki, A. Szczurek, A. Marcinek, M. Kiełbowicz, Phys. Rev. C 102, no.1, 014901 (2020) https://doi.org/10.1103/PhysRevC.102.014901. arXiv:1910.04544 [nucl-th]
J. Estee et al. [S\(\pi \)RIT], Phys. Rev. Lett. 126(16), 162701 (2021). https://doi.org/10.1103/PhysRevLett.126.162701. arXiv:2103.06861 [nucl-ex]
J.R. Stone, P. Danielewicz, Y. Iwata, Phys. Lett. B 826, 136915 (2022). https://doi.org/10.1016/j.physletb.2022.136915
W.A. Zajc, J.A. Bistirlich, R.R. Bossingham, H.R. Bowman, C.W. Clawson, K.M. Crowe, K.A. Frankel, J.G. Ingersoll, J.M. Kurck, C.J. Martoff et al., Phys. Rev. C 29, 2173–2187 (1984). https://doi.org/10.1103/PhysRevC.29.2173
Lanczos, C.: The Variational Principles of Mechanics. Unabridged Dover (1986) republication of the fourth edition published by University of Toronto Press, Toronto (1970)
A.D. Chacon, J.A. Bistirlich, R.R. Bossingham, H. Bossy, H.R. Bowman, C.W. Clawson, K.M. Crowe, T.J. Humanic, M. Justice, P. Kammel et al., Phys. Rev. C 43, 2670–2688 (1991). https://doi.org/10.1103/PhysRevC.43.2670
W.B. Christie, W.F.J. Mueller, D.L. Olson, T.J.M. Symons, H.H. Wieman, D. Beavis, F.P. Brady, J.L. Romero, C.E. Tull, T. Abbott et al., Phys. Rev. C 45, 2836–2853 (1992). https://doi.org/10.1103/PhysRevC.45.2836
S. Pratt, Phys. Rev. D 33, 72–79 (1986). https://doi.org/10.1103/PhysRevD.33.72
Lounesto, P.: Clifford Algebras and Spinors. Cambridge University Press, Cambridge (2001)
H. W. Barz, Phys. Rev. C 59 (1999), 2214-2220 https://doi.org/10.1103/PhysRevC.59.2214[arXiv:nucl-th/9808027 [nucl-th]]
R. Maj, S. Mrowczynski, Phys. Rev. C 80, 034907 (2009). https://doi.org/10.1103/PhysRevC.80.034907. arXiv:0903.0111 [nucl-th]
J. Adamczewski-Musch et al. [HADES], Phys. Lett. B 795, 446–451 (2019). https://doi.org/10.1016/j.physletb.2019.06.047. arXiv:1811.06213 [nucl-ex]
J. Adamczewski-Musch et al. [HADES], Eur. Phys. J. A 56 (2020) no.5, 140 https://doi.org/10.1140/epja/s10050-020-00116-w[arXiv:1910.07885 [nucl-ex]]
G. Agakishiev et al. [HADES], Eur. Phys. J. A 41 (2009), 243-277 https://doi.org/10.1140/epja/i2009-10807-5[arXiv:0902.3478 [nucl-ex]]
J. Adamczewski-Musch et al. [HADES], Eur. Phys. J. A 56 (2020) no.10, 259 https://doi.org/10.1140/epja/s10050-020-00237-2[arXiv:2005.08774 [nucl-ex]]
W. Reisdorf et al. [FOPI], Nucl. Phys. A 848 (2010), 366-427 https://doi.org/10.1016/j.nuclphysa.2010.09.008[arXiv:1005.3418 [nucl-ex]]
H. Schuldes, Charged kaon and\(\phi \)reconstruction in Au+Au collisions at 1.23\(A\)GeV, PhD thesis, G.W. Goethe University, Frankfurt am Main, Germany (2016)
P.J. Siemens, J.O. Rasmussen, Phys. Rev. Lett. 42, 880–887 (1979). https://doi.org/10.1103/PhysRevLett.42.880
I. Fröhlich, L. Cazon, T. Galatyuk, V. Hejny, R. Holzmann, M. Kagarlis, W. Kühn, J. G. Messchendorp, V. Metag, M. A. Pleier, et al. PoS ACAT (2007), 076 https://doi.org/10.22323/1.050.0076[arXiv:0708.2382 [nucl-ex]]
J. Adamczewski-Musch et al. [HADES], Eur. Phys. J. A 54 (2018) no.5, 85 https://doi.org/10.1140/epja/i2018-12513-7[arXiv:1712.07993 [nucl-ex]]
S. Teis, W. Cassing, M. Effenberger, A. Hombach, U. Mosel, G. Wolf, Z. Phys. A 359 (1997), 297-304 https://doi.org/10.1007/s002180050405[arXiv:nucl-th/9701057 [nucl-th]]
V. S. Uma Maheswari, C. Fuchs, A. Faessler, L. Sehn, D. S. Kosov, Z. Wang, Nucl. Phys. A 628 (1998), 669-685 https://doi.org/10.1016/S0375-9474(97)00646-5[arXiv:nucl-th/9706004 [nucl-th]]
R. Stock, Phys. Rept. 135, 259–315 (1986). https://doi.org/10.1016/0370-1573(86)90134-1
J. Adamczewski-Musch, O. Arnold, C. Behnke, A. Belounnas, A. Belyaev, J.C. Berger-Chen, A. Blanco, C. Blume, M. Böhmer, P. Bordalo et al., Phys. Lett. B 819, 136421 (2021). https://doi.org/10.1016/j.physletb.2021.136421. [arXiv:2012.01351 [nucl-ex]]
G. Goldhaber, S. Goldhaber, W.Y. Lee, A. Pais, Phys. Rev. 120, 300–312 (1960). https://doi.org/10.1103/PhysRev.120.300
G. Baym, Acta Phys. Polon. B 29 (1998), 1839-1884 [arXiv:nucl-th/9804026 [nucl-th]]
D. Hardtke, T. J. Humanic, Phys. Rev. C 57 (1998), 3314-3318 https://doi.org/10.1103/PhysRevC.57.3314 [arXiv:nucl-th/9710010 [nucl-th]]
O. Buss, L. Alvarez-Ruso, P. Mühlich, U. Mosel, Eur. Phys. J. A 29 (2006), 189-207 https://doi.org/10.1140/epja/i2006-10075-y[arXiv:nucl-th/0603003 [nucl-th]]
V. Metag, M. Nanova, E. Y. Paryev, Prog. Part. Nucl. Phys. 97 (2017), 199-260 https://doi.org/10.1016/j.ppnp.2017.08.002[arXiv:1706.09654 [nucl-ex]]
L. D. Landau, E. M. Lifshitz Statistical Physics 2nd ed. Course of Theoretical Physics, ed2, vol5, 1969,
S. Harabasz, W. Florkowski, T. Galatyuk, M. Gumberidze, R. Ryblewski, P. Salabura, J. Stroth, Phys. Rev. C 102 (2020) no.5, 054903 https://doi.org/10.1103/PhysRevC.102.054903[arXiv:2003.12992 [nucl-th]]
A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, Nature 561 (2018) no.7723, 321-330 https://doi.org/10.1038/s41586-018-0491-6[arXiv:1710.09425 [nucl-th]]
F. Becattini, J. Steinheimer, R. Stock, M. Bleicher, Phys. Lett. B 764 (2017), 241-246 https://doi.org/10.1016/j.physletb.2016.11.033[arXiv:1605.09694 [nucl-th]]
J. Cleymans, H. Oeschler, K. Redlich, S. Wheaton, Phys. Rev. C 73, 034905 (2006). https://doi.org/10.1103/PhysRevC.73.034905. [arXiv:hep-ph/0511094 [hep-ph]]
L. Adamczyk et al., STAR. Phys. Rev. C 96, 044904 (2017). https://doi.org/10.1103/PhysRevC.96.044904. [arXiv:1701.07065 [nucl-ex]]
X. Lopez et al., FOPI. Phys. Rev. C 76, 052203 (2007). https://doi.org/10.1103/PhysRevC.76.052203. [arXiv:0710.5007 [nucl-ex]]
G. Agakishiev et al. [HADES], Eur. Phys. J. A 52 (2016) no.6, 178 https://doi.org/10.1140/epja/i2016-16178-x[arXiv:1512.07070 [nucl-ex]]
J.L. Klay et al., E895. Phys. Rev. C 68, 054905 (2003). https://doi.org/10.1103/PhysRevC.68.054905. [arXiv:nucl-ex/0306033 [nucl-ex]]
C. Alt et al., [NA49]. Phys. Rev. C 77, 024903 (2008). https://doi.org/10.1103/PhysRevC.77.024903. [arXiv:0710.0118 [nucl-ex]]
C. Alt et al., [NA49]. Phys. Rev. C 73, 044910 (2006). https://doi.org/10.1103/PhysRevC.73.044910. [arXiv:nucl-ex/0512033 [nucl-ex]]
Acknowledgements
We thank Andreas Wagner for clarifying discussions and for providing the KaoS data points. We are furthermore grateful to Kai Gallmeister and Ulrich Mosel for elucidating discussions of the Coulomb effects in their GiBUU transport calculations. This work is in memory of our colleague Robert Greifenhagen. Financial support by the following agencies and Grants is acknowledged: SIP JUC Cracow, Cracow (Poland), 2017/26/M/ST2/00600; TU Darmstadt, Darmstadt (Germany), VH-NG-823, DFG GRK 2128, DFG CRC-TR 211, BMBF:05P18RDFC1; Goethe-University, Frankfurt (Germany), BMBF:06FY9100I, BMBF:05P19RFFCA, GSI F &E, HIC for FAIR (LOEWE); Goethe-University, Frankfurt (Germany) and TU Darmstadt, Darmstadt (Germany), ExtreMe Matter Institute EMMI at GSI Darmstadt; NRNU MEPhI Moscow, Moscow (Russia), RFBR funding within the research project no. 18-02-00086, Fundamental properties of elementary particles and cosmology No 0723-2020-0041; JLU Giessen, Giessen (Germany), BMBF:05P12RGGHM; IPN Orsay, Orsay Cedex (France), CNRS/IN2P3; NPI CAS, Rez, Rez (Czech Republic), MSMT LTT17003, MSMT LM2018112, MSMT OP VVV CZ.02.1.01/0.0/0.0/18-046/0016066.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Consortia
Corresponding author
Additional information
Communicated by Silvia Masciocchi.
The preparation of this manuscript was completed before February 24, 2022.
Appendix A: Derivation of the Jacobian \(J_{\mathrm {eff}}\)
Appendix A: Derivation of the Jacobian \(J_{\mathrm {eff}}\)
In this appendix we derive the Jacobian \(J_{\mathrm {eff}}\), expressed in Eq. (6), introduced by the explicit momentum dependence of the effective Coulomb potential energy \(V_{\mathrm {eff}}\) given by Eq. (2). To start, we recall that the Jacobian of the total energy transformation \(E = E_0 \pm V_C\) caused by a constant Coulomb potential energy \(V_C\) is given by
The first equality follows from the spherical symmetry of the potential, i.e. by integration over polar and azimuthal angles, and the last one from the constancy of \(V_C\) leading to \(dE_0/dE = 1\).
If the constant \(V_C\) is replaced by an energy-dependent effective potential \(V_{\mathrm {eff}}(E)\), \(E = E_0 + V_{\mathrm {eff}}(E)\) and the \(dE_0/dE\) term changes into
As discussed in Sect. 2.1, the expansion of the charged fireball leads to an attenuation of the average Coulomb potential acting on a charged pion of given energy \(E = E(v)\), which we model, following Ref. [27], by the fraction of protons that have a lower velocity than the pion, \(v_p < v_{\pi }\). This results in the effective potential \(V_{\mathrm {eff}}(E)\) expressed by Eq. (2) for 2D and 3D expansions, respectively. Setting the variable \(x = \sqrt{(E_{\pi }/m_{\pi } - 1) \, m_p/T_p}\), we obtain by derivation
from which the second Jacobian term in Eq. (6) follows.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Adamczewski-Musch, J., Arnold, O., Behnke, C. et al. Impact of the Coulomb field on charged-pion spectra in few-GeV heavy-ion collisions. Eur. Phys. J. A 58, 166 (2022). https://doi.org/10.1140/epja/s10050-022-00796-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epja/s10050-022-00796-6