Bremsstrahlung as a probe of baryon stopping in heavy-ion collisions

In collisions between heavy ions at ultra-relativistic energies the participating protons lose energy, which is converted into new particles. As the protons slow down, they emit bremsstrahlung radiation. The yield and angular distribution of the emitted radiation are sensitive probes of how much energy the incoming protons have lost. In this paper, the spectrum of bremsstrahlung radiation is calculated for different stopping scenarios, and the results are compared with the expected yield of photons from hadronic interactions.


Introduction
In collisions between heavy ions at relativistiv energies there is convincing evidence that a new state of matter, a quark-gluon plasma, is formed [1].In this state, the quarks and gluons are no longer confined to nucleons but can move freely over distances large compared with the size of a single nucleon.The energy density in the quark gluon plasma formed in Pb+Pb collisions at the LHC has been estimated from measurements of the total transverse energy at midrapidity to be on the order of 12-14 GeV/fm 3 at a time of 1 fm/c after the collision [2,3].This is far above the densitites of 0.2-0.5 GeV/fm 3 lattice QCD calculations find are required for deconfinement [4].
The energy deposited in the quark-gluon plasma comes from the energy lost by the incoming nuclei, and one of the most fundamental question one can address in the study of high energy heavy-ion collisions is therefore how much energy the incoming baryons lose.This is usually referred to as the amount of baryon stopping.One can have scenarios ranging from complete stopping, where the incoming baryons lose all their energy, to full transparency, where the baryons lose no or very little energy.Full stopping would imply that all baryons end up close to midrapidity, whereas full transparency would leave the baryons near beam rapidity.Since baryon number is conserved, the fate of the baryons in the colliding nuclei can be determined from the rapidity distribution of net baryons, that is dn B /dy−dn B /dy.For experimental reasons, one is often restricted to study the net proton rather than the net baryon distributions.
Results from the Relativistic Heavy-Ion Collider (RHIC) [5][6][7] and fixed target experiments at the CERN SPS [8] show that the amount of stopping decreases with increasing collision energy in the range √ s N N = 7 − 200 GeV.There have been attempts to explain the energy loss in this energy range from hadron transport models [9] and models based on the Color Glass Condensate [10,11].
At the CERN SPS and RHIC, identified protons and anti-protons could be measured down to low transverse momenta p T ∼ 0 over a wide rapidity range, and the net-proton rapidity distributions could thus be determined.Such measurements were performed by the NA49 [8] and BRAHMS [6,7] experiments.At the LHC, the situation is different.The only experipment which has measured identified protons and anti-protons at low p T is ALICE.The results have shown that in the central rapidity region |y| ≤ 0.5 there are no net protons [12][13][14].But beyond that, there are no experimental constraints on how the net-protons are distributed.
To improve this situation, we propose to use the bremsstrahlung photons emitted when the nuclei slow down.This idea was first suggested before the start of the relativistic heavy-ion programs at CERN and RHIC [15,16].Before the start-up of RHIC, several studies were made where this process was considered [17][18][19].These also included a proposal to build a dedicated detector to study this radiation [18], but those plans were never realized.Recently, the idea was brought up again in the context of the LHC [20].
The previous studies mentioned above all use a similar, semi-classical approach to calculate the bremsstrahlung spectrum, based on the description in [21].Our calculations will follow the same path.We will, however, implement improved stopping scenarios which are either based on model calculations or phenomenological and consistent with existing data.The stopping scenarios considered in [20] are simplified and do not take into account the fact that the central region (|y| ≤ 0.5) at the LHC is baryon free.This was, however, recently followed up by a study of phenomenological stopping scenarios where the central region is almost baryon free [22].
We will also, for the first time, make a detailed estimate of the background from hadronically produced photons, primarily from the decay of π 0 mesons.This background is obtained from simulations with PYTHIA 8.3 [23].The goal is to determine in what regions of phase space one can expect bremsstrahlung photons to provide a realistic measure of the nuclear stopping.One can, in addition to photons from hadronic interactions, also expect a large background from secondary photons produced in the detector material.The latter is specific to a certain experiment and its material budget and is thus beyond the scope of this paper.

The bremsstrahlung spectrum
The energy radiated per solid angle from a current J(r, t) is given by [21] The vector n is a unit vector in the direction of the photon, and here and throughout we use units where = c = 1.We choose coordinates where the incoming beams move along the z-axis and the photon is emitted in the xz-plane with angle θ, giving n = (sin(θ), 0, cos(θ)).
In the center of mass, the incoming nuclei are Lorentz contracted in the longitudinal direction to a size ∼ R/γ, where R is the nuclear radius and γ the Lorentz factor of the beam.For a lead nucleus at the LHC, this corresponds to a longitudinal size of about 0.003 fm.It is thus justified to ignore the longitudinal relative to the transverse extension and write the currents for the incoming nuclei as Here, σ is the nuclear electric charge density in the transverse plane, v 0 the velocity, v 0 = tanh(y b ), where y b is the beam rapidity, and θ(t) the Heaviside step function.The charge density is normalized to σd 2 r ⊥ = Z, where Z is the number of protons in the beam nucleus.The outgoing protons will have a distribution in the transverse plane, which we assume is the same as for the incoming particles, and a distribution in velocities, which may or may not depend on the position in the transverse plane.Writing the velocity in terms of the rapidity, v(y) = tanh(y), and the corresponding density as ρ(y, r ⊥ ) one gets for the outgoing current The density of the outgoing protons is normalized to since it includes the contribution from both incoming nuclei.With the total current J = J + + J − + J f , the radiated energy becomes Integrating by parts in the integral over time one obtains The time, ∆t, and longtudinal distance over which the protons are slowed down can be expected to be small compared with the transverse size, ∆t << R. For low energy photons (ω << 1/∆t), it is therefore justified to neglect the time and longitudinal components in the phase factor and make the assupmtion With these assumptions the integrals over time and z in (6) can be performed, leading to Here, α = e 2 /4πǫ 0 is the fine structure constant.This result is in agreement with [18].To do the integral over the transverse dimensions, one has to know the function ρ(y, r ⊥ ).It is conceivable that there is a dependence of the rapidity loss of the protons on the transverse coordinate; protons close to the center of the nuclei can be expected to lose more energy than those on the periphery.Previous studies have, however, found that the dependence of ρ(y, r ⊥ ) on the transverse position has only a minor effect on the spectrum of bremsstrahlung photons [18,19].Moreover, the models we will use for the rapidity loss of the protons provide the average loss, independent of position in the transverse plane.We will therefore ignore the dependence on r ⊥ here and assume ρ(y, r ⊥ ) = ρ(y).The integral over r ⊥ can then be performed and the spectrum of emitted photons can be written Here, F (Q) is the nuclear form factor obtained from a Fourier transform of the nuclear charge distribution, ρ A : We use a form factor where R A = 6.62 fm, ρ 0 = 0.161 fm −3 , and a = 0.70 fm for a Pb nucleus.This parameterization has been shown to reproduce the Fourier transform of a Woods-Saxon distribution in configuration space very well [24].One can note that for low energy photons and small emission angles, ω sin(θ) << 1/R, the Form Factor in Eq. 9 is approximately 1.The energy and angular dependencies then factorize, with the energy dependence given by 1/ω and the angular dependence by To proceed with the calculations one has to define the function ρ(y), the rapidity distribution of the net-protons in the final state.We consider 3 scenarios for central Pb+Pb collisions at the LHC.
1.The net proton distribution as given by PYTHIA 8.3 [23].Heavy-ion collisions have been implemented in PYTHIA through the Angantyr model [25].
The scaling from proton-proton collisions is based on the Glauber model and the original PYTHIA framework is used to describe the individual nucleon-nucleon sub-collisions.It thus extrapolates the dynamics of pp collisions to heavy-ion collisions, without introducing any collective effects between the nucleon-nucleon collisions.It does reproduce the measured charged particle pseudorapidity distributions in Pb+Pb collisions at the LHC. 2. The net proton distribution as given by the hadron transport model SMASH-2.2 [9,26].The model combines a string model, where the colliding hadrons are excited to strings which fragment, with elastic and inelastic interaction between hadrons in the later stages of the collisions.This approach leads to a considerably larger amount of stopping compared with PYTHIA.In fact, the model predicts a non-zero number of net-protons at midrapidity.We use this result for the current calculation anyway, since it represents a valid result of the model at LHC energies [H.Elfner, private communication].3. A phenomenological model where the central region, |y| ≤ 0.5, contains no net protons, but where the protons have a considerable shift away from beam rapidity.This is modelled by the sum of two skewed Gaussians.
The net-proton rapidity distributions, ρ(y), for these 3 scenarios are shown in Fig. 1.Of the three models, Pythia clearly shows the least amount of stopping.While SMASH-2.2 exhibits a non-zero yield of net-protons at midrapidity, the protons are on average not shifted so much toward y = 0 as in scenario 3.

Results
The angular distributions of bremsstrahlung photons in the forward direction, integrated over the energy range 0.1 ≤ ω ≤ 0.5 GeV, from the 3 scenarios are shown in Fig. 2 a).The distributions are peaked close to 1/γ, as expected, and the yields increase with an increasing amount of stopping.The curves also exhibit different angular dependencies, the scenarios with more stopping having more persistent tails toward larger angles.Thus, the models are differentiated in both total photon yield and angular dependence.
To facilitate a comparison with experimental acceptances, which are usually defined in terms of pseudorapidity, η = − ln(tan(θ/2)), the angular distribution in Eq. 9 can be rewritten as The pseudorapidity distribution, integrated over azimuthal angle and energy interval 0.1 ≤ ω ≤ 0.5 GeV, is shown in Fig. 2 b).The difference between the spectra is more pronounced in the pseudorapidity distribution than in the angluar distribution.This is due to the factor sin 2 θ = 1/ cosh 2 (η) which varies rapidly for large |η|.Covering the entire range of emission angles, the figure illustrates the difference in total number of radiated photons between the stopping scenarios.Furthermore, while the peaks of all three spectra lie at large pseudorapidities, the scenarios give significantly differing photon yields at lower η as well, making such scenarios potentially discernable within experimental acceptances.To put these numbers in context, we compare them with the photon yield from the 5% most central Pb+Pb collisions from PYTHIA 8.3 in Fig. 3.These hadronically produced photons, most of which come from the deacy π 0 → γ+γ, constitute a background to the bremsstrahlung photons we are considering here.The background yield is shown by the black histograms in the figure.The sum of the yield of the background and bremsstrahlung photons is shown by the solid, blue histograms for scenario 1-3 in Fig. 3 a), b), and c), respectively.The number of bremsstrahlung photons is calculated from Eq. 13, integrated over azimuthal angle and photon energy 0.1 ≤ ω ≤ 0.5 GeV.
The bremsstrahlung calculations above assume that all protons participate in the collision.This will not be the case for collisions within a finite impact parameter range.We therefore calculate a correction factor ( N p, part /2Z) 2 , where N p, part = 151 is the average number of participating protons in the 5% most central collisions in Pythia.The result of applying this correction factor to the yield of bremsstrahlung photons is shown by the dashed, blue histograms in the figure .From the figure one can see that at low pseudorapidities, the background completely dominates.In the very forward direction, however, the background falls off quickly while the bremsstrahlung peak emerges.The yield of bremsstrahlung photons differ widely between the different scenarios, which emphasizes that this is indeed a very sensitive probe of the amount of nuclear stopping.The shapes of the pseudorapidity distributions are also quite different between the scenarios.This means that the limit in pseudorapidity where one can expect a significant signal over background is lower the larger amount of stopping one has.
As mentioned, a detailed discussion of in which experiments one might extract a bremsstrahlung signal is beyond the scope of this paper.We nevertheless indicate in Fig. 3 the experimental acceptances of the current and future experiments where it might be possible.The existing LHCb experiment has an electromagnetic calorimeter coverage between 2.0 ≤ η ≤ 4.5 [27].We include it here, although it might not be able to reach low enough photon energies at large pseudorapidities [28], [R.McNulty, private communication].During the Next Long Shutdown at the LHC (2026 -2029) it is foreseen to install a forward calorimeter (FoCal) in the ALICE experiment [29].It will consist of a high resolution electromagnetic and hadronic calorimeter covering 3.4 ≤ η ≤ 5.8.Finally, beyond LHC Run 4 there are plans to upgrade the ALICE experiment to ALICE-3 [30].The current design of ALICE-3 includes a Forward Conversion Tracker, which should have the possibility to measure photons with energies down to or below 100 MeV in the pseudorapidity range 3.0 ≤ η ≤ 5.0.The pseudorapidity coverages of these detectors are shown by the red lines in Fig. 3.
From the figure one can see that for scenario 3, "No central charge", there is a visible excess over the hadronic background within the pseudorapidity coverage of all the detectors mentioned above.For scenario 1, "PYTHIA 8.3",  the situation is less favorable and one would have to go to the most forward regions of ALICE-3 and FoCal to find a good signal to background ratio.
Since the Photon Conversion Tracker in ALICE-3 aims at measuring photons with energies below 100 MeV, we also include a plot of the photon pseudorapidity distributions for the photon energy range 0.01 ≤ ω ≤ 0.1 GeV in Fig. 4. In this energy range, there is a significant excess inside the ALICE-3 acceptance for all stopping scenarios.In addition to being peaked in the forward direction, the bremsstrahlung spectrum increases rapidly with decreasing photon energy, approximately as 1/ω, as was mentioned above.This is contrary to the background from hadronically produced photons, which decrease with decreasing ω in the energy range considered here.To illustrate this, we plot the energy spectrum integrated over azimuthal angle and pseudorapidity range 4 ≤ η ≤ 5 in Fig. 5.As in Figs. 3  and 4, the background is given by the black histograms, and the uncorrected and corrected signal plus background by the solid and dashed blue histograms, respectively.As for the pseudorapidity distributions, the energy below which one can expect a significant signal over background is highly dependent on the stopping scenario.Also the yield is strongly dependent on the stopping scenario.
The inset in Fig. 5 a) shows the low energy region, and it emphasizes that if one can go to low enough photon energies, a signal will be visible also in scenarios with a small amount of stopping.One should keep in mind that it might be possible to extract a bremsstrahlung signal, even with a rather low signal to background ratio, by subtracting the hadronic background.The hadronic background should be well constrained from measurements of charged and neutral particle spectra.

Summary
To conclude, we have shown that even with realistic stopping scenarios the bremsstrahlung spectra show a strong sensitivity to the amount of nuclear stopping.Comparisons with Pythia show that a significant signal over the hadronic background is obtained in the range η 4−5 and ω 300−500 MeV.Again, the exact limits depend on the amount of stopping one has.
Considering the importance of determining the amount of stopping in heavy-ion collisions at the LHC, and given that no alternative methods are available, we believe the possibility to use bremsstrahlung photons should be considered seriously.Hopefully this paper can help in the design of future detectors to accomplish such a measurement.

2 Fig. 1
Fig.1Net-proton rapidity distributions for the 3 scenarios described in the text.

2 Fig. 2 a
Fig.2a) The angular distribution of bremsstrahlung photons with energies between 0.1 ≤ ω ≤ 0.5 GeV for the 3 scenarios.b) The pseudorapidity distributions, dNγ /dη, of bremsstrahlung photons within the same energy range.

Fig. 3
Fig.3The pseudorapidity distributions for photons with 0.1 ≤ ω ≤ 0.5 GeV integrated over the azimuthal angle.The black histogram shows the background from hadronically produced photons.The solid blue histogram shows the sum of the photons from bremsstrahlung radiation and hadronic production.The correction applied to the bremsstrahlung spectrum to obtain the dashed blue histogram is described in the text.

Fig. 5
Fig.5The energy distributions for photons with 4.0 ≤ η ≤ 5.0 integrated over the azimuthal angle.The black histogram shows the background from hadronically produced photons.The solid blue histogram shows the sum of the photons from bremsstrahlung radiation and hadronic production.The correction applied to the bremsstrahlung spectrum to obtain the dashed blue histogram is described in the text.The inset in a) shows the low energy region.