Michael A. Fedderke \faOrcidmfedderke@perimeterinstitute.caPerimeter Institute for Theoretical Physics, Waterloo, Ontario, N2L 2Y5, Canada
The William H. Miller III Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA
David E. Kaplan \faOrciddavid.kaplan@jhu.eduThe William H. Miller III Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA
Anubhav Mathur \faOrcida.mathur@jhu.eduThe William H. Miller III Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA
Surjeet Rajendran \faOrcidsrajend4@jhu.eduThe William H. Miller III Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA
Erwin H. Tanin \faOrcidehtanin@stanford.eduStanford Institute for Theoretical Physics, Stanford University, Stanford, California, 94305, USA
The William H. Miller III Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA
Abstract
The tentative identification of approximately ten relativistic anti-helium () cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy anti-nuclei.
We propose a novel scenario for the enhanced production of such anti-nuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) anti-quarks in our galaxy by physics beyond the Standard Model (BSM).
We demonstrate that SM anti-nucleosynthetic processes that occur in the resulting rapidly expanding, thermalized fireballs of SM plasma can, for a reasonable range of parameters, produce the reported tentative ratio of to events at AMS-02, as well as their relativistic boosts.
Moreover, we show that this can be achieved without violating anti-deuterium or anti-proton flux constraints for the appropriate anti-helium fluxes.
A plausible BSM paradigm for the catastrophic injections is the collision of macroscopic composite dark-matter objects carrying large net anti-baryon number.
Such a scenario would require these objects to be cosmologically stable, but to destabilize upon collision, promptly releasing a fraction of their mass energy into SM anti-particles within a tiny volume.
We show that, in principle, the injection rate needed to attain the necessary anti-helium fluxes and the energetic conditions required to seed the fireballs appear possible to obtain in such a paradigm.
We leave open the question of constructing a BSM particle physics model to realize this, but we suggest two concrete scenarios as promising targets for further investigation.
I Introduction
The AMS-02 Collaboration111The Alpha Magnetic Spectrometer (AMS-02) is a particle-physics detector located on the International Space Station [1]. has unofficially reported [2, 3, 4, 5, 6, 7] highly relativistic cosmic-ray events detected in years of data that are consistent with tentative identification as anti-helium [6, 7].222These unofficial reports [2, 3, 4, 5, 6, 7] have taken the form of public oral presentations on behalf of the Collaboration in the context of scientific conferences or major colloquia, as well as the associated publicly available presentation slides.
We stress however that these data have not to date been published, have never been officially claimed by the AMS-02 Collaboration to present a formal detection of anti-helium cosmic rays, and have always been accompanied by disclaimers and caveats that the origin of these candidate events requires more study.
Additionally, only partial data is available for these events. Although publicly available mass determinations are uncertain [6], the data are reported to be consistent with tentative identification of both and candidate events, with an event ratio of roughly (albeit with small statistics and large uncertainties) [6, 7].
Additionally, tentative identification of 7 anti-deuterium candidate events has recently been reported [7].
While the tentative identifications of these candidate events require more work to confirm [7], the anti-helium candidates in particular have been the subject of extensive recent interest in the literature [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18] because, taken at face-value, they are surprising: the formation of complex anti-nuclei in astrophysical environments is challenging and the rates of production implied by these candidate events are hard to reconcile with known Standard Model (SM) physics.
Within the SM, a known source of anti-nucleus cosmic rays is spallation induced by primary cosmic rays (hydrogen or helium) in the interstellar medium (ISM) [12].
Spallation is, however, inefficient at producing anti-nuclei; due to kinematics, this is particularly true for those anti-nuclei with higher atomic mass number , such as .
Coalescence of two anti-nucleons or anti-nuclei (or of an anti-nucleon and an anti-nucleus) produced in spallation into a higher- anti-nucleus is probable only if their relative kinetic energy is comparable to or below the difference in the nuclear binding energies of the initial and final states; roughly, .
On the other hand, the threshold energy for production of nuclear anti-particles via spallation of a primary cosmic ray against the ISM is .
Because primary cosmic-ray fluxes tend to be power laws as a function of energy (see, e.g., Ref. [1]), it follows that the kinetic energy of nuclear anti-particles produced by spallation of sufficiently energetic primary cosmic rays is also of typically , except very near threshold.
The formation rates of heavier anti-nuclei via coalescence of such products therefore tend to suffer significant phase-space suppressions, leading to strong hierarchies between the expected numbers of anti-nucleus events with subsequently higher that would be observed at AMS-02: i.e., .
These conventional astrophysics predictions for anti-nucleus cosmic ray fluxes from spallation were quantified recently in Ref. [12], where the expected numbers of anti-nucleus events at AMS-02 were found to scale roughly as , in line with the phase-space argument above.
For model parameters that reproduce the anti-proton flux observed by AMS-02, Ref. [12] thus found that the predicted anti-helium fluxes are always orders of magnitude below AMS-02 sensitivity.
Conversely, in order to reproduce, e.g., the tentatively identified flux, one would have to overproduce, e.g., the observed anti-proton flux by many orders of magnitude.
Similar challenges in producing anti-nuclei also occur in decaying or annihilating particle dark-matter (DM) models commonly studied in the context of indirect searches [12, 10, 11, 8, 9].
The anti-nucleus production rates in these models suffer from similar phase-space suppressions as for spallation of primary cosmic rays.333The phase-space suppression argument made above, and the resulting hierarchy of anti-nucleus fluxes with higher , should apply for any anti-nucleus production scenario that starts with high-energy () processes. One should keep in mind however that there are considerable variations in the predicted anti-nucleus fluxes stemming from uncertainties in the parameters of the nuclear-coalescence model [8, 9] and different choices of the cosmic-ray propagation model [14].
With optimistic assumptions [14, 11], and possibly
with enhancements from decays [15], it has been found that it might be possible for annihilating particle dark matter to be the origin of the tentatively identified flux.
However, even a single confirmed event at AMS-02 would be challenging to explain.
See also Ref. [19] for a recent detailed investigation of these points.
If the candidate events are confirmed, explaining the presence of both the and events at AMS-02, with their comparable rates, seems to require the absence of the severe phase-space suppressions that inexorably lead to a strong hierarchical relation of the anti-nucleus fluxes.
These suppressions can be ameliorated if the anti-nucleons that combined to form the anti-nuclei have low relative momenta, .
Various beyond-the-Standard-Model (BSM) scenarios have been proposed to achieve this.
Refs. [12, 16, 18] considered anti-nucleus production occurring in anti-matter–dominated regions of primordial origin that have cooled down significantly by the time of the anti-nucleus production.
These scenarios require a BSM mechanism in the early Universe to cause the requisite matter–anti-matter segregation.
Anti-nucleus production scenarios via the decay of a new particle carrying anti-baryon number (which may or may not be the dark matter) have been also considered.
In Ref. [13], the mass of the decaying particle was tuned to be very close to the mass of the desired anti-nucleus in order to restrict the final-state phase space, such that the produced particles are non-relativistic.
Such a scenario however requires several new decaying particles, each with its own mass tuning to separately enhance the production of , , and perhaps .
In Ref. [17], a strongly coupled dark sector was considered, where dark hadron showers triggered by the decay of a new particle simultaneously increase the multiplicity of the decay products and decrease their relative momenta .
The final decay products (i.e., the lightest dark bound states) then decay to SM anti-quarks which subsequently form anti-nuclei.
A challenging aspect of this scenario is the need to model strong-coupling phenomena such as dark hadron showers.
Another challenging aspect of the tentatively identified events at AMS-02 is their relativistic nature (i.e., large Lorentz boosts).
Overcoming phase-space suppressions of nuclear coalescence rates by considering scenarios where the colliding particles have low relative momenta usually also results in anti-helium products that are non-relativistic.
There are however ways around this: the anti-helium products in scenarios that start with the decay of a particle, such as those considered in Refs. [13, 17], can be made relativistic if the decaying particle is already boosted in the galactic rest frame.
This can be achieved if, e.g., the decaying particle is produced in turn from the earlier decay of another heavier particle.
Other works [12, 16, 18] have either suggested acceleration mechanisms based on supernova shockwaves (similar to the Fermi mechanism), or did not address this question.
In this paper, we propose a scenario where the anti-helium nuclei observed at AMS-02 originate from sudden and localized “injections” (of BSM origin) of anti-baryon number in our Galaxy in the form of energetic SM anti-quarks.
These particles subsequently thermalize into relativistically expanding, optically thick fireballs with a net anti-baryon number.
Anti-nuclei are produced in each of these fireballs thermally through a nuclear reaction chain similar to that operating during Big Bang Nucleosynthesis (BBN), albeit with remarkable qualitative and quantitative differences due to the very different (anti-)baryon–to–entropy ratios, timescales, and expansion dynamics involved.
The evolution of the fireball plasma after the initial injection of SM particles is purely dictated by SM physics and, moreover, since this process involves thermalization, depends only on certain bulk properties being achieved by the injection, and for the most part not on the exact details of the latter.
However, owing to the inefficiency of weak interactions, there are regions of parameter space where the anti-neutron–to–anti-proton ratio in the plasma at the onset of the anti-nucleosynthetic processes may depend on the details of the initial injection of SM particles.
With that one exception, the predicted anti-particle cosmic ray fluxes produced in our proposed scenario are therefore both predictive and largely agnostic as to the microphysical origin of the injections that seed these fireballs.
Assuming that the requisite injections can occur, we show that this scenario could explain not only the tentatively identified anti-nucleus fluxes at AMS-02, but also the relativistic Lorentz boosts of the detected particles.
The seemingly paradoxical requirements of a low-energy environment to foster production of higher- anti-nuclei and the high energies required for relativistic anti-nucleus products are reconciled naturally in our scenario by the expansion dynamics of the fireball plasma: as the plasma expands, its temperature falls adiabatically while its thermal energy is converted to bulk kinetic energy by the work of its internal pressure.
This allows a low-energy environment for anti-nucleosynthesis to proceed, while at the same time accelerating its products to relativistic speeds.
Of course, injections of the requisite amounts of SM anti-particles with the correct properties to generate these fireballs cannot occur spontaneously: a BSM mechanism is required.
Suggestively, we show that collisions of certain very heavy, macroscopic, composite dark objects (possibly a sub-fraction of all of the dark matter) that carry SM anti-baryon number could at least achieve a high-enough rate of injections with appropriate parameters to explain the tentatively identified anti-helium events at AMS-02, provided that a substantial fraction of the dark objects’ mass energy can be converted into SM anti-quarks as a result of dark-sector dynamics triggered by the collision (in some parameter regions, there may be a requirement to inject also a comparably sized asymmetry of charged leptons).
While this is encouraging, we have not yet developed a detailed microphysical model that achieves the necessary dark-sector dynamics in a way that is amenable to robust and controlled understanding; however, we offer some speculative initial thoughts on certain model constructions that we believe are promising avenues to explore toward that goal.
For the purposes of this paper, we ultimately leave this question open; as it is crucial to providing a concrete realization of the scenario we advance, however, we both intend to return to it in our own future work and we encourage other work on it.
The remainder of this paper is structured as follows.
In Sec. II, we review the candidate anti-nuclei events observed by AMS-02.
In Sec. III, we discuss synthesis of anti-nuclei in an expanding relativistic fireball, discussing first questions of thermalization after energy injection [Sec. III.1], then turning to the fireball expansion dynamics [Sec. III.2] and its termination at the point of photon decoupling [Sec. III.3], and the nuclear reactions at play during the expansion [Sec. III.4], before summarizing [Sec. III.5].
In Sec. IV, we then discuss how the anti-nuclei thus produced would propagate to AMS-02 and give our scenario’s projections for the anti-nuclei spectra and event rates [Sec. IV.1]; we also discuss other potential observables [Sec. IV.2].
In Sec. V, we discuss a potential origin for the fireballs in collisions of composite dark-matter states, showing first that the rates could work [Sec. V.1], that the fireballs could be appropriately seeded if certain benchmarks can be met [Sec. V.2], and then offering some speculative thoughts toward particle physics models that may be worth further investigation to see it they are able to achieve the necessary conditions [Sec. V.3].
We conclude in Sec. VI.
A number of appendices add relevant detail.
Appendix A gives derivations of various scaling laws for fireball expansion that we rely on in the main text.
Appendix B discusses the ratio of anti-nucleons from which the anti-nucleosynthesis is initiated.
Appendices C and D, respectively, give details of the nuclear reaction networks and cross-sections we have used.
Appendix E discusses whether dynamical changes in the number of degrees of freedom during fireball expansion are relevant.
Appendix F discusses prompt versus slow injections.
Finally, Appendix G reviews an estimate of the AMS-02 rigidity-dependent sensitivity to anti-helium events.
II The AMS-02 candidate anti-helium events
In this section, we summarize the data that is publicly available [2, 3, 4, 5, 6, 7] regarding the AMS-02 candidate anti-helium and anti-deuterium events.
Table 1: Parameters for individual candidate anti-helium events displayed in the identified references: is momentum, is mass, is charge ( is the elementary charge), is the Lorentz factor, and is rigidity.
For identification purposes across the references, the reported event date and corresponding day of year are given (the date of event 2 is not given in the references; we give the date of its first public presentation in Ref. [2]).
Unavailable data are denoted by “—”.
We have set .
Data for , , and are given in the references [speed data, , is additionally given for event 2]; we derived [or for event 2] and and propagated uncertainties naïvely.
For comparison, and .
Event # 3 is clearly tentatively identified in Refs. [3, 7] as a candidate event.
Event # 4, which is similarly identified as an candidate in the oral presentation of Ref. [7], is not included in our analysis; see discussion in text.
To our knowledge, the most recent scientific presentations on these candidate events are Refs. [6, 7].
Ref. [6] provides graphical mass and rigidity histogram data for 9 candidate events collected in of AMS-02 observations.
Ref. [7] also appears to show an additional anti-helium candidate event dated after Ref. [6], with a mass measurement consistent with either or , but closer to the former.
Table 1 shows the detailed basic parameter values for 3 of the 9 candidate events that have been shown publicly [2, 3, 4, 5, 7], as well as the data available for the additional candidate event [7]; such detailed data has not been presented for the other candidate events.
We base our analysis on the data for the 9 candidate events presented in Ref. [6], which are consistent with the identification of 6 candidate events, and 3 candidate events in of AMS-02 data (see discussion about below).555This identification is also consistent with earlier presentations based on a smaller dataset with a total of 8 candidate anti-helium events [3, 5], of which 2 were tentatively identified as candidates [3, 4], giving a 3:1 ratio for .,666We do not explicitly consider the “additional” event shown in Ref. [7] (event # 4 in Table 1).
This event is labeled in the slide deck for Ref. [7] to have occurred on September 20, 2022, which is after the date of presentation of the histograms in Ref. [6] on February 28, 2022; it is thus highly likely to be a new event not previously discussed in past references before Ref. [7].
However, at the level of uncertainty regarding these events that we work in this paper, whether or not we include this event has almost no relevant impact on our discussion of the overall or relative event rates for the and events.
Assuming either that or for the data in Ref. [6] (see discussion in main text), it is also entirely consistent within statistical errors for one additional event to occur in the additional integration time of relevance for the of data that appears to be discussed in Ref. [7]. The mass histogram shown in Ref. [6] is however more than broad enough to support an inference that some of the candidate events thus classified could in fact be events, and vice versa (the additional event shown in Ref. [7] could also be identified as either isotope within the uncertainties).
Nevertheless, throughout this paper, we adopt a event ratio for vs. as fiducial.
Of the candidate events in Ref. [6], 3 have rigidities (where is the particle momentum and its charge) in the approximate range , while the other 6 candidate events have rigidities in the range .
The data in Table 1 however make clear that care should be taken not to associate the 3 candidate events with larger shown in the histogram in Ref. [6] with the 3 candidate events and the remainder with the 6 candidate events: indeed, event # 3 in Table 1 has a rigidity in the low range, but a mass most consistent with .
One robust inference however is that all of the candidate events for which the necessary data are available to make this determination are highly relativistic, with being a typical fiducial value (an extreme range of possible Lorentz factors for the other candidate events based on the rigidity and mass data shown in Ref. [6] is roughly , assuming ).
Given the confirmed helium events in the AMS-02 data as of Ref. [6], the ratio of the 9 candidate anti-helium events in that reference to the confirmed helium events is approximately .
Likewise, Ref. [7] reports confirmed helium events; with 10 total anti-helium events, this yields the same ratio within uncertainties.
Finally, for the purposes of converting total candidate anti-helium event numbers into rate estimates, we assume that the relevant AMS-02 data-taking period over which all 9 events discussed in Ref. [6] were observed is .
There is some uncertainty on this given available information: the data-taking period of relevance to the 9 events reported as of Ref. [6] may actually be slightly shorter, (this point is not made unambiguously clear in Ref. [6]).
There is thus an uncertainty on the required rates of arising from the incompleteness of the publicly available information on this point; this is however significantly smaller than the uncertainty on the rates owing to the small statistics of the relevant event samples.
Additionally, Ref. [7] makes a new report of 7 candidate anti-deuterium () events in of AMS-02 data.
However, only scant information is available regarding these events: a single histogram showing that the charge-sign–mass product for these events lies in the range .
III Fireball anti-nucleosynthesis
In this section, we show how an abrupt localized injection in a region of space of a large amount of energy and anti-baryon number in the form of SM anti-quarks can lead to formation of a locally thermalized fireball comprised of a plasma mixture of free anti-nucleons, pions, leptons, and photons.
As this fireball expands hydrodynamically [20, 21, 22, 23, 24, 25], it cools adiabatically, eventually permitting Standard Model anti-nucleosynthetic processes to produce bound anti-nuclei, including anti-helium, in the hot and dense environment.
Owing to the dynamics of the expansion in the interesting region of parameter space, the radial bulk expansion velocity of the fireball constituents also becomes relativistic by the time of anti-nucleosynthesis, resulting in the anti-helium thus produced being released into the interstellar medium relativistically (assuming the fireball is located within our galaxy).
We also demonstrate that the fireball expansion shuts off anti-nucleosynthetic processes prior to the attainment of complete nuclear statistical equilibrium, thereby allowing the amount of produced, after decays of unstable products, to be larger than the amount of .
The thermalized nature of the initial fireball caused by the requisite injection of energy and anti-baryon number conveniently erases most of the history of how such a thermal state comes to exist.
As such, the conclusions we reach in this section as to the anti-nuclear outputs of the fireball expansion are largely independent of the model details of how such a fireball comes to exist; instead, they depend only on bulk physical properties of the fireball state, such as its temperature, anti-baryon–to–entropy ratio, and radius at certain critical points in its evolution.
The exception to this is that, due to incomplete thermalization via inefficient weak interactions, the results we obtain can depend on the net charge on the hadronic sector that is injected primarily via the BSM process that seeds the fireball; this essentially becomes another parameter we must consider.
Of course, this history-erasure does not alleviate the requirement for a concrete BSM mechanism by which the requisite energy and anti-baryon number injection could occur; we comment on this aspect of the problem in Sec. V but ultimately defer this to future work.
In the remainder of this section, we discuss: the thermalization following energy and anti-baryon number injection [Sec. III.1], the dynamics of the fireball expansion [Sec. III.2], and the anti-nucleosynthetic and other outputs of the fireball expansion [Sec. III.4].
We summarize in Sec. III.5.
III.1 Thermalization
We will be mainly concerned with SM fireballs whose initial size is of order , and which double in size on characteristic timescales of order . The injected particles will interact through various Standard Model processes which overall tend to bring themselves toward local thermal equilibrium.
For simplicity, we limit ourselves to the parameter space where the injected SM energy density is such that the would-be temperature of the SM plasma is below the QCD scale,777At higher energy densities, the plasma would lie above the QCD phase transition, a quark-gluon plasma (QGP) would form, and one would need more careful treatment of the evolution of the fireball back through the phase transition as it cools.
.
The path toward thermalization in this regime will involve a process of hadronization where the injected anti-quarks confine and fragment into anti-nuclei and copious pions.
The rate of hadronization is set by the QCD scale in the rest frame of the confining or fragmenting particles and will be time-dilated at some level in the center of mass frame of the injected gas of particles. Unless the initial Lorentz factors of the injected anti-quarks are extremely high, (corresponding to energies ), a case we do not consider, hadronization occurs essentially instantaneously compared to the initial expansion timescale of the fireball, .
The subsequent thermalization should proceed similarly to that in analogous setups involving hadrons found in the contexts of heavy-ion colliders [26, 27, 28] and in the early Universe before [29, 30] and during [31, 32] BBN.
A proper description of the thermalization process would require solving a complex set of Boltzmann equations dictating the time-evolution of the energy spectra of relevant particles and resonances.
We instead provide rough estimates of the typical rates of the processes involved.
The typical rates of strong interactions involving pions (e.g., ), ; electromagnetic (EM) interactions for charged particles that are relativistic at a given temperature (e.g., ), ; and weak interactions for relativistic particles (e.g., ), , are respectively given by
(1)
(2)
(3)
with and [32, 33].
These rates suggest that within the timescale of expansion of the initial fireball , the SM particles can interact efficiently through strong and electromagnetic processes, but generically not through weak processes, unless the temperature is significantly higher than .
In the cases that we consider, strong and electromagnetic processes are sufficient for thermally populating all SM particle species with masses below a given temperature.
For instance, anti-protons can be created from charged pions and anti-neutrons through strong interactions, photons can be produced via bremsstrahlung from charged pions and the decay of neutral pions, and leptons pairs can be produced in photon annihilations.
However, if the weak interactions are indeed always inefficient after SM particle injection, charge would need to be separately conserved in the leptonic and hadronic sectors, because reactions such as (and similar crossed or charge conjugated reactions, as well as similar reactions with the anti-baryons) that would allow charge to be exchanged between those two sectors would not be efficient.
Therefore, were the initial injection of particles to be such that the the net charge is zero in both sectors (e.g., only net charge- and color-neutral combinations of anti-quarks are injected), the charged pions would have a chemical potential that fixes their population asymmetry to the number of anti-protons.
If the plasma temperature falls, this may have what we will see to be undesirable consequences for our scenario, such as a dramatic depletion of the anti-proton abundance when the symmetric, thermal pion abundance becomes Boltzmann suppressed while strong interactions between the asymmetric abundance and the population remain efficient for some time.888One possible injection-model-dependent solution to this issue is an initial injection of SM particles that is still net color- and charge-neutral, but which in addition to having a large net negative baryon number, also has a large charged lepton asymmetry and therefore baked-in opposite-sign EM charge asymmetries in the hadronic and leptonic sectors that cannot be removed by the inefficient weak interactions. We discuss this issue in more detail in Sec. III.4.3 and Appendix B; see also comments in Sec. V.3.
Immediately upon the completion of (partial) thermalization, the fireball energy is dominated by radiation comprised of a subset of photons and relativistic , , and (depending on the temperature); its mass is dominated by anti-baryons in the form of free and .
The thermal pressure of the trapped radiation drives an adiabatic expansion of the fireball, which can be described hydrodynamically owing to the short mean free path of the constituent particles, .
This is the topic of the next subsection.
III.2 Relativistic fireball expansion
We now turn to considering the dynamics of the expanding thermalized plasma.
We treat the plasma as a spherically symmetric perfect fluid and consider only its radial expansion.
In what follows, quantities defined in the comoving rest-frame of the radially moving fluid will be marked with a prime ′, while those defined in the fireball center-of-mass (c.m.) frame are unprimed.
If is the CoM time co-ordinate and is the CoM-frame radial coordinate centered on the fireball, then the radial expansion of the fireball is described by two relativistic fluid equations that arise from the covariant conservation of the energy-momentum tensor for a perfect fluid assuming spherical symmetry: , where with for radial fluid flow () and is the metric for spherical co-ordinates .
Respectively, the equation encodes energy conservation and the equation encodes radial momentum conservation of the fireball fluid [34, 35]:
(4)
(5)
where is the Lorentz factor associated with the CoM-frame bulk radial velocity of the fluid at a given radius , and and are the comoving energy density and pressure of the fluid.
While Eqs. (4) and (5) are not analytically solvable, it is understood that the fireball will undergo an initial phase of rapid acceleration to relativistic speeds [34, 36], whereupon analytically tractable evolution takes over.
The acceleration is initially limited to a thin shell near the surface of the fireball where the pressure gradient is strong.
This surface expansion then generates inward-traveling rarefaction waves which accelerate the bulk of the fireball.
Numerical simulations [37, 38, 39] (see also Refs. [40, 41]) suggest that the radial layers comprising an initially static thermal fireball will accelerate to relativistic radial velocities within the time it takes for sound waves (whose speed in a radiation-dominated fluid is ) to cover the initial radius of the fireball.
Upon arrival at the center of the of the fireball, the waves get reflected outward, creating a strong underdensity at that point.
This underdensity creates a hollow structure at the center of the fireball that becomes more pronounced over time, essentially turning the thermal plasma into a radially moving shell of density concentration with an initial thickness .
The radial bulk velocities of all the radial layers of the plasma shell soon approach the speed of light and the resulting nearly flat velocity profile keeps the shell thickness constant until much later times.
Once the radial layers of the plasma shell are moving relativistically with Lorentz factors , its subsequent evolution follows simple scaling laws, which can be parameterized by the initial temperature and the initial radius (and thickness) when , as well as the (conserved) anti-baryon–to–entropy ratio of the plasma shell [34, 36].999Because the plasma temperature will be related to the shell radius through simple scaling laws, the precise definitions of and beyond what we have described are not important.
One can more generally define as the temperature when the radius of the shell is (or vice versa), where the reference (or ) is arbitrarily chosen.
So long as the reference point is sufficiently early in the shell evolution that the energy of the plasma is still dominated by radiation, the specific choice of reference point is unimportant.
For simplicity, in this paper we choose the reference point to be when . Note that and naturally have radial variation over the shell; however, thanks to the smoothing effect of sound waves, this variation is typically mild [34].
For simplicity, our parametric estimates will assume a single (average) value of and for the entire shell.
We will however take into account the radial variation of the Lorentz factor between the inner and outer radius of the shell because, as we will see, it has an important consequence at the late stages of the shell evolution; where relevant, we assume that varies monotonically across the shell, being larger on the outer edge of the shell.
The expansion in the relativistic regime proceeds in three stages (see Appendix A for derivations):
1.
Acceleration. As long as the shell energy density is radiation-dominated, the internal pressure of the radiation continually accelerates the radial velocity of the shell, converting radiation energy into anti-baryon kinetic energy in the process.
In this phase, the shell thickness in the center-of-mass frame of the fireball remains approximately constant, , and the bulk Lorentz factor scales with the radius of the shell as .
See also Refs. [42, 43].
2.
Coasting. Once the radiation energy drops below the anti-baryon kinetic energy, the radiation can no longer accelerate the shell appreciably and from then on the shell simply coasts at its terminal Lorentz factor
(6)
where is the proton mass.
3.
Spreading. The (assumed monotonic) variation of the Lorentz factor between the inner and outer radii of the shell translates to velocity variation . This leads to increasing shell thickness which becomes important when and is well captured by
(7)
To sum up, the Lorentz factor of the shell scales as
(8)
Neglecting changes in the number of degrees of freedom as various species (pions, muons, electrons, positrons) fall out of thermal equilibrium and get Boltzmann suppressed,101010This amounts to setting (from photons) throughout.
We discuss in Appendix E how accounting for changes in yields only mild, quantitative changes to the picture presented here, but leaves the qualitative evolution unchanged.
the average comoving temperature of the shell scales as
(9)
These scaling laws are in agreement with the numerical simulations in Refs. [44, 34, 36].
Note that the co-moving dynamical expansion timescale (i.e., -folding timescale) for the fireball, as measured by the change in its comoving temperature,111111This turns out to be a convenient measure of fireball expansion for our later purposes; see also Appendix C at Eq. (97). can then be approximated as
(10)
where the three cases here map to the three different expansion regimes for the fireballs discussed above.
We display the evolution of this timescale for a set of benchmark parameters that will be of interest (see Sec. III.4) in Fig. 1.
The evolution discussed here is valid so long as photons remain tightly coupled to the plasma; we turn to this topic in Sec. III.3.
Figure 1: The comoving expansion timescale (solid black line) as a function of the comoving temperature [cf. Eq. (10)] of a fireball with , , and .
The three vertical dashed lines mark, from left to right, the transition at from the acceleration phase to the coasting phase, the transition at from the coasting phase to the spreading phase, and the point at which photon decoupling occurs [corresponding to ; cf. Eq. (20) with ].
The orange shaded region is thus the region in which the plasma expands as a single tightly coupled fluid through its various expansion phases (as annotated), while the unshaded region where is the region where the plasma has decoupled; see Sec. III.3.
III.3 Photon decoupling
The fluid evolution described in Sec. III.2 applies as long as the photons remain tightly coupled with the charged particles in the plasma.
Eventually, this assumption is violated and the photons decouple.
In this section, we estimate when this occurs.
Photons in the plasma scatter dominantly with the electrons and positrons.
In our parameter space of interest, we will show that it is the case that photon decoupling always occurs at a temperature well below the electron mass, , after any symmetric thermal population of charged leptons have annihilated away.
The opacity of the plasma is then due to the remaining charged leptons that must exist to guarantee net neutrality of the plasma; while their identity depends on some details of the BSM injection process, residual positrons are most important to this estimate.121212As we discuss in the next section, typical expansion timescales are too short for any muon asymmetry that is present to be guaranteed to decay by the time of this decoupling.
Because of the scaling however, residual muons are not relevant to our estimates unless . Let us therefore parametrize the population of residual positrons available for the photons to scatter from as having a CoM-frame density , where can be a BSM-injection dependent parameter, but which we typically expect to be .
Note that our estimates will therefore conservatively underestimate the population of light charged particles available to scatter if they are naïvely extended to .
While all of the light charged particles are of course initially highly relativistic upon thermalization in the comoving fireball frame (i.e., long before photon decoupling), because we will show that the fireball becomes optically thin only when , it will be appropriate in this section (when working in the comoving fireball frame) for us to self-consistently use the Thompson cross section as the relevant scattering cross-section for photons from the residual positrons around the time of photon decoupling.
In fact, we will also employ as the relevant cross section even when naïvely extending these results to ; the optical depth we derive should therefore not be understood to be accurate for except insofar as it indicates an optical depth (i.e., an opaque fireball), which qualitative conclusion we do not expect would be modified were we to instead use the full Klein–Nishina cross-section.
Now, consider a photon moving in the CoM of the fireball at an angle relative to the radial direction.
If , its mean free path in the comoving frame of the fireball plasma is .
This can be related via a Lorentz transformation to the mean-free path in the CoM frame as .
Solving for and using , we find .
The optical depth for a photon emitted from a radius inside the shell and escaping to infinity is thus given by
(11)
where is the displacement of the photon in the CoM frame.
Note that is in principle radius-dependent; however, we can assume for the following reason.
Due to relativistic beaming, most of the thermal photons are concentrated within in the CoM frame, yielding and .
Therefore, we can take and, up to a numerical factor that is not important for this parametric estimate, approximate .
The optical depth for a typical photon in the shell emitted from (when ) is then [45]
(12)
The antiproton density appearing in this expression needs to understood with some caution: it is the density experienced by the photon when it is located at CoM radial co-ordinate , which differs from the density of the antiprotons as a function of the CoM radial co-ordinate when evaluated at a fixed instant of time; we denote that latter density here as , using the fireball radius in the CoM frame as a proxy for time.
For the purposes of this estimate, suppose that is approximately constant over the interval , where is the -dependent thickness of the fireball shell when the inner edge of the fireball has radius .
It will turn out (and we will show this a posteriori) that the moment at which the photons decouple from the fireball plasma will be deep in the spreading phase of its expansion (consistent with ).
As a result, (for ) is approximately constant, and .
Now, as the photon travels outward at speed from initial radius , the fireball will also be expanding, so the anti-proton density it will experience is , where .
This density is non-zero on the interval from where is defined implicitly by .
For the purposes of this estimate, let us take .131313Note that this estimate actually fails if we insert ; we assess that this is due to various factors that we have neglected here.
For instance, it would work for .
While we wish to be transparent about this issue, it is not a serious problem for this estimate, which is intended to be parametric only. In the limit , this yields .
Note that this yields a range of support for the integral in Eq. (12) over a range .
This is really the important point that this parametric estimate leads to: the integral in Eq. (12) will run over a range .
Physically, the relativistic expansion of the fireball at implies that a photon spends much longer inside the fireball plasma than it would were the same (CoM-frame) thickness of plasma stationary (in the CoM frame).
The constant value that we assume that takes over its range of support is approximately , where (without an argument) is the CoM fireball antiproton density when it is approximately spherical and has radius (this estimate assumes that anti-protons are not significantly depleted by nuclear reactions during expansion, which is true in the parameter region interest to us; see Sec. III.4).
Now, because will increase by more than an order of magnitude as increases to , it follows that will fall rapidly over the range of its support in the integral in Eq. (12).
As such, up to numerical factors that we neglect in this estimate, we will approximate the integral in Eq. (12) using the value of the integrand evaluated at , multiplied by the range over which the integrand takes approximately that value, which will be according to the estimates above (i.e., the true integrand falls by over this characteristic change in the integration variable, and the upper limit of the integration is parametrically large enough as so to not cut it off before then).
That is, we will approximate Eq. (12) as
(13)
where we used , (spreading phase), and we have defined to be
(14)
(15)
(16)
(17)
(18)
where we inserted , which is understood here to be the baryon number density when the fireball is at radius and temperature ; applied the definition of the anti-baryon–to–photon ratio; and dropped a numerical constant for (see Appendix E), because this estimate should be understood parametrically only.
If we define to be the approximate moment in time where the plasma becomes optically thin to photons emitted in its bulk, and the corresponding radius of the fireball at that time to be , we can write , leading to
(19)
(20)
At the fiducial parameter point, , which is less than the rest-frame muon lifetime, as noted above.
This is also deep in the spreading regime for these parameters and we also have (cf. Fig. 1), validating our assumptions above.
When the fireball radius hits , the bulk of the plasma becomes optically thin and a burst of photons is released.
This is also the approximate moment at which the anti-baryons and other particles that were coupled to the fireball plasma (see Sec. III.4) are released into the interstellar medium (see Sec. IV), assuming that the fireball was located in our Galaxy.
III.4 Nuclear physics
The monotonically decreasing co-moving temperature of the fireball plasma given at Eq. (9) implies that it eventually becomes thermodynamically favorable for bound anti-nuclei to form.
Were thermodynamic equilibrium among the lightest few anti-nuclei species to be achieved (analogous to the situation in Big Bang Nucleosynthesis [BBN]), almost the entirety of the available anti-neutron abundance would be converted into , leading to a final configuration dominated by and , with only trace amounts of other complex light anti-nuclei.
Given that AMS-02 has tentatively identified similar numbers of and candidates (up to a factor of a few; statistics are small), such an outcome would not be phenomenologically viable.
To understand how to avoid this outcome, consider a key feature of how the analogous process of ordinary BBN proceeds.
The most efficient nuclear-reaction pathway to has as an initial step free-neutron capture on hydrogen to form deuterium D: [46], with the D then being processed by further nuclear burning to other light elements.
However, the relatively low deuterium binding energy and the low baryon-to-entropy ratio during BBN make deuterium prone to photodissociation back to free neutrons and protons: this is the famous “deuterium bottleneck” [47, 48, 49, 50].
Consequently, production in BBN was delayed until the temperature of the primordial plasma cooled down significantly below , whereupon the abundance of photons capable of dissociating deuterium was considerably Boltzmann suppressed, enabling the deuterium abundance to rise.
Crucially, in the BBN realized in our Universe, the deuterium bottleneck was overcome while thermodynamic equilibrium was still being maintained among the light species: nuclear reaction rates were still sufficiently fast compared to Hubble expansion that, once deuterium was capable of being created without being photodissociated, it was rapidly burned to tritium and , and then further to .
But standard BBN successfully overcame the deuterium bottleneck while maintaining thermodynamic equilibrium only marginally.
If the expansion rate of the Universe were to have been sufficiently larger, neutron capture on hydrogen would have decoupled (i.e., frozen out) before the bottleneck could have been overcome.
In that case, the output of the nucleosynthesis would not have been dictated by thermodynamic equilibrium among the light nuclei.
Instead, the immediate nucleosynthesis products in this scenario would have been mostly free protons and neutrons, with smaller abundances of deuterium, tritium, , and being produced in amounts controlled by the relative rates of nuclear reactions that produce them, which can be comparable to one another.
Accounting for the fact that unstable neutrons later decay to protons and that tritium later decays to , the final output in this counterfactual case could easily have been such that .
In what follows, we show that parameter space exists for which the analogous anti-nucleosynthesis occurring in our expanding fireball remains in this “stuck in the bottleneck” regime, yielding phenomenologically viable amounts of as compared to .
III.4.1 Preliminaries
In order to separate changes in the number density of an element due to nuclear reactions from that due to the fireball expansion, in this section we describe the evolution of a nuclear species in terms of its fractional abundance , defined as
(21)
where is the number density of element and is the anti-baryon number density.
We consider only .
We obtain the evolution of the abundances of nuclear elements by numerically solving the Boltzmann equations detailed in Appendix C describing the simplified network of nuclear reactions among these species,141414We solve a simplified, partial reaction network accounting only for light species, which is acceptably accurate for our purposes.
In principle, more accurate results could be obtained by using a modified version of BBN nucleosynthesis codes such as PRyMordial [51], PRIMAT [52], or AlterBBN [53].
using the nuclear cross sections shown in Appendix D.
The results of these numerical computations are summarized in
Figs. 2–5.
While those numerical results are of course more accurate, we also wish to gain an understanding of, and intuition for, the most important nuclear processes at work, and determine the dependencies of the final anti-helium isotope abundances on the fireball parameters .
In what follows, we therefore develop an analytical understanding that reproduces the gross features of the numerical results.
The main approximation we employ in our analytical arguments is as follows.
In general, nuclear species with higher mass numbers are produced from those with lower in a sequence of successive two-body nuclear reactions.
Our numerical analysis shows that in the fireball parameter space that yields , a hierarchy is maintained between the nuclear abundances with successive mass numbers, namely , as depicted in Figs. 2 and 3.
Furthermore, that analysis shows that anti-nucleosynthesis occurs mainly at temperatures low enough that all of the endothermic reverse nuclear reactions to the exothermic ones considered in this subsection are negligible, except for the photodissociation of anti-deuterium responsible for the bottleneck in the first place,151515This endothermic process is an exception because (a) it involves photons, which are highly abundant compared to anti-baryons due to the low anti-baryon–to–entropy ratio we consider; and (b) the binding energy of anti-deuterium is unusually small, . which we thus take into account in our analysis.
These observations suggest that, instead of solving the whole nuclear reaction network at once, we can treat the nuclear reactions sequentially; that is, we can consider elements produced earlier in the chain of successive two-body reactions as fixed sources for reactions later in the chain of successive reactions, neglecting back-reaction on those sources arising from those later reactions.
In what follows, we first discuss relevant expansion timescales, and then proceed to discuss in turn heavier and heavier anti-nuclei synthesized in this approximate sequential paradigm.
Finally, we summarize and discuss other, non–anti-nucleosynthetic outputs.
III.4.2 Expansion timescales
Anti-nucleosynthesis in the fluid rest frame of the expanding fireball is qualitatively similar to BBN in that there are various nuclear reactions occurring in an adiabatically expanding background [47, 48, 49, 50].
However, it differs from the BBN in important ways.
The plasma in our scenario is spatially finite and expands relativistically into vacuum.161616Our anti-nucleosynthesis process resembles in this respect that occurring in the context of gamma-ray burst [54, 55, 56, 57] or heavy-ion collision [58], but is otherwise very different. This of course leads to a non-trivial dependence of the dynamical expansion timescale on the comoving temperature , as shown at Eq. (10).
The parameter space that is viable for our model is roughly (i.e., ), ,
and (i.e., ).
We displayed the scalings of in Fig. 1 for benchmark parameters in these ranges.
The results are quantitatively and qualitatively very different compared to the Hubble time as a function of temperature during BBN; as we will see, this leads to important differences between fireball anti-nucleosynthesis and BBN.
As we will show, fireball anti-nucleosynthesis in this parameter space commences at temperature , which satisfies and thereby always falls in the “spreading” phase of the fireball expansion.
Figure 2: The abundances of nuclear species (solid and dashed colored lines, as identified in the legend) as a function of the comoving fireball temperature , computed by numerically solving the Boltzmann equations for a simplified nuclear-reaction network detailed in Appendix C, for a fireball with , , ().
Also shown is (dotted black line), the abundance of photons with energies above the anti-deuteron photodissociation threshold .
Initially, nuclear bound states are essentially non-existent, apart from whose abundance is kept at an exponentially small value due to the high abundance of photodissociating photons, .
As the fireball cools down, decreases while increases. Eventually at [cf. Eq. (32)], drops sufficiently low that photodissociation becomes inefficient, thus marking the onset of anti-nucleosynthesis: the absence of photodissociation allows to rise significantly, which enables nuclear reactions involving to produce whose presence, in turn, enables the production of .
The rise in also enables production but the produced quickly converts into .
Since the nuclear reaction rates are proportional to the antibaryon number density , most of the anti-nucleus production occurs within the first few -foldings of expansion after drops below .
Nuclear species depicted in dashed lines are stable on the timescale of fireball evolution (until it becomes optically thin) but are expected to decay to the species shown by the solid line of the same color during their journey to the solar system.
Figure 3: The abundances of nuclear species (with as annotated in the legend) released when a fireball with an initial radius becomes optically thin, computed by numerically solving the Boltzmann equations for a simplified nuclear-reaction network, as detailed in Appendix C.
Here, we set and ().
While the hierarchy is maintained, the abundances of , , released by the fireball scale in reasonable agreement with our analytical predictions, when the appropriate comparisons are made: (a) numerically, , which is only slightly shallower than the analytical scaling from Eq. (35); (b) the relative scaling holds reasonably well [cf. the form of expressed at Eq. (43)] leading to when combined with the numerical result , in reasonably good agreement with these numerics [actually, is a slightly better numerical fit for the relative scaling, leading to , which is also a slightly better numerical fit]; and (c) the relative scaling also holds very well [cf. Eq. (46)], leading to when combined with the numerical result , again in very good agreement with these numerics. Note however, in connection with (b) and (c), that the naïve analytical expectations based on our discussion in the main text would be and , respectively, if one took from Eq. (35); see also further discussion in Sec. III.4.8.
Besides that, is produced promptly in negligible amount, [cf. Eq. (38) and surrounding discussion].
Nuclear species depicted in dashed lines are stable on the timescale of the fireball evolution (until it becomes optically thin) but are expected to subsequently decay to the species shown in with solid lines of the same color on timescales short compared to their propagation time in the Milky Way (e.g., later decays to with a half-life of years in the rest frame).
III.4.3 Anti-nucleon abundance
At temperatures , nuclear bound states have not formed and all the anti-baryons reside in unbound anti-neutrons and anti-protons.
Anti-neutrons can, in principle, convert to and from anti-protons through both weak (e.g., ) and strong (e.g., ) processes [59, 32, 33].
If at least one of these processes is efficient, the relative abundance of these anti-nucleons is initially kept at its chemical equilibrium value, .
In the early Universe, the matter analogs of both processes were efficient at some point after hadronization.
Then, strong processes decoupled first as pions rapidly decayed and annihilated away, and hence the pre-BBN freeze-out abundances of neutrons and protons were determined by the later-occurring decoupling of weak interactions.
By contrast, in the fireball anti-nucleosynthesis scenario we consider, the typically short timescales of the fireball expansion render weak interactions inefficient at all times.
Consequently, the anti-nucleons freeze out as soon as the pion-mediated strong interconversion processes become inefficient.
After the fireball has thermalized, the following strong-mediated charge exchange reactions (SMCER) are initially in equilibrium [31] (see discussion about below)
Additionally,171717The reactions , , , and are all in equilibrium.
, , and .
These imply
(22)
The chemical equilibrium anti-neutron–to–anti-proton ratio for is thus given by
(23)
where .
Note that the chemical potential of the charged pions, , depends on the physics before the fireball has thermalized via efficient strong and EM interactions.
Therefore, it is model-dependent.
For instance, it depends on whether electroweak interactions were ever efficient in this pre-thermalization stage, and on some details of the BSM particle injection process that seeds the fireball.
For simplicity, we neglect the chemical potential of charged pion in our analysis here by assuming181818As we discuss Appendix B, this specific assumption is equivalent to assuming that there is a net negative charge in the hadronic sector of the plasma that has a certain very specific value: defining as at Eq. (76), we would have ; cf. Eq. (29).
This charge is compensated by opposite charge in the leptonic sector so that the plasma as a whole is net EM-neutral as expected from fireballs seeded by EM-neutral dark states.
However, our results as stated in the main text are unchanged qualitatively, and change quantitatively by only factors, so long as it is approximately true that by the time that the SMCER become inefficient.
We also show in Appendix B that we may even be able to tolerate values as small as , although that changes some conclusions stated in the main text in a qualitative fashion.
(24)
We discuss the model-dependence of and the case when is non-negligible in Appendix B.
The pion-mediated strong interactions decouple at a temperature where the pion abundance becomes sufficiently Boltzmann suppressed that the found in Eq. (1) goes below the fireball expansion rate found in Eq. (10).
We find that, numerically, invariably in the whole parameter space that is viable for our scenario.
The freeze-out value of the anti-neutron–to–anti-proton ratio can be approximated by its chemical-equilibrium value at that time
(25)
If , unlike what we have assumed, then the freeze-out value of and our subsequent results would change; however, as long as , these changes are only and most of our conclusions remain valid.
See Appendix B for further discussion.
Moreover, while neutron decay is an important phenomenon in the BBN that was realized in the early Universe, in our scenario anti-neutrons do not decay until well after anti-nucleosynthesis finishes.
Assuming that only a small fraction of the and are burned to higher nuclei (true throughout our parameter space of interest), we will thus have, for all times relevant for the anti-nucleosynthesis in the expanding fireball, the following:
(26)
(27)
where
(28)
implying that
(29)
III.4.4 Anti-deuterium production
Anti-deuterium is produced primarily through the reaction191919The anti-deuterium formation releases some amount of energy density to the plasma, given by the total binding energy of the anti-deuterium formed: , where . This amounts to a tiny fraction of the radiation energy density in the parameter space of our interest, where , , and when the anti-deuterium forms. .
Initially, however, the reverse reaction (photodissociation) is in equilibrium and the high abundance of photons with energies above the anti-deuterium binding energy suppresses the (quasi-equilibrium) anti-deuterium abundance, which is given by the Saha equation:
(30)
This continues until the abundance of photons with sufficient energy to photodissociate anti-deuterium,
(31)
starts to fall below the anti-deuterium abundance; i.e., .
The temperature at that point can be estimated as
(32)
where the displayed range of values corresponds to the range of viable anti-baryon–to–entropy ratios .
As mentioned earlier, in our the parameter space of interest falls in the spreading phase of the fireball expansion (i.e., it satisfies ) and, neglecting the mild logarithmic dependence on , the fireball expansion timescale at decoupling of the photodissociation reactions is
(33)
where we have set (corresponding to ).
After the photodissociation of decouples at , at which point the fireball expansion timescale is , anti-deuterium production through is no longer thwarted, and so the anti-deuterium abundance rises monotonically.
At around the same time, heavier elements that rely on anti-deuterium burning as an initial step begin to be populated sequentially.
Since the product of the fireball expansion timescale and the nuclear reaction rates that form any of the light elements scales as during the spreading phase (), these nuclear reactions are most efficient in populating the light elements in the first fireball-expansion -fold or so after the decoupling of photodissociation, before the anti-baryon density is significantly diluted by the expansion.
Moreover, as we operate in the regime where anti-deuterium is not efficiently burned to more complex nuclei (see the next sub-subsection), a simple estimate for the final anti-deuterium abundance can be obtained by assuming the anti-deuterium abundance is that generated by neutron–proton fusion reactions operating in a single dynamical expansion timescale at the point of anti-deuterium photodissociation freeze-out.
We estimate that abundance to be
(34)
(35)
where we took the value of cross-section to be that at (the lower end of the range of values for ); see Appendix D.
In writing the above results, we have used and we manually inserted an prefactor in Eq. (34), such that the final result is in better agreement with what we obtained by numerically solving the Boltzmann equations at this benchmark point.
Note also that, in our numerical results, we find the scaling of with is actually closer to (cf. Fig. 3); this is important in the context of later results that will raise this result to large powers [cf. Eqs. (43) and (46)].
We are interested in the regime where : i.e., the anti-deuterium production decouples before its abundance rises to .
III.4.5 Anti-deuterium burning
During and slightly after the production of anti-deuterium, a small fraction of it also burns through the following dominant channels:
with essentially equal branching fractions (both are strong-mediated nuclear reactions).
Because we are analyzing production in the regime where anti-deuterium is not efficiently burned to more complex anti-nuclei, we may treat the abundance of anti-deuterium as a fixed source which acts to populate the more complex nuclei over roughly a single dynamical expansion timescale after the anti-deuterium are produced.
As such, the prompt production of and can be estimated as
(36)
where we have evaluated all the quantities at , assumed around the time of this production, and manually included an numerical prefactor [cf. the factor introduced in Eq. (34)].
Alternative production channels for and are and , respectively; however, we verified numerically that the former production channel is negligible as long as and the latter is negligible as long as , which are always satisfied in the parameter space we consider.
It is understood that these process are inefficient because they both suffer from photon-emission suppression (i.e., they are electromagnetic-mediated, rather than strong-mediated, nuclear reactions).
The reaction can also proceed with a branching ratio of (it is electromagnetically mediated).
Because of this small branching fraction, prompt production through this channel, , is negligible compared to other, deuterium–tritium-burning channels that we discuss below so long as around the time of production.
III.4.6 Anti-helium-3 burning
The strong-mediated reaction is also present.
It is also extremely efficient in part because it is not Coulomb suppressed.
It thus gives the one counterexample to our earlier statement that we can ignore back-reaction on sequentially produced species: because it is exothermic ( [60, 61], assuming completely ionized nuclei as appropriate at the relevant temperatures), this reaction burns essentially all the anti-helium-3 that are produced primarily (via ) to anti-tritium.
As such, the anti-helium-3 abundance is maintained at a very low, quasi-equilibrium level: we numerically found that the residual never exceeds .
At the same time, the abundance is roughly doubled because our estimate at Eq. (36) indicated roughly equal production abundances for the two anti-nuclei before this depletion reaction was accounted for.
Our estimates of the anti-tritium and anti-helium-3 abundances after this burning should therefore be revised to
(37)
(38)
III.4.7 Anti-tritium burning
Anti-tritium burns efficiently to anti-helium-4 through the following dominant process:
assuming that throughout the burning, we find that the prompt production is
(39)
where we again manually introduced an prefactor to better match our numerical results.
III.4.8 Final anti-nucleosynthesis products
Figure 4: The fireball parameter space.
The initial temperature and radius of the fireball are defined at the point right after the fireball plasma has reorganized itself into a shell moving with an average Lorentz factor .
The various styles of blue lines show contours of constant values of ratio of the injected anti-helium-4 and anti-helium-3 abundances, .
These contours were obtained by numerically solving the Boltzmann equations for a simplified nuclear-reaction network detailed in Appendix C, while fixing the anti-baryon–to–entropy ratio such that the terminal Lorentz factor of the shell is .
The orange lines show contours of constant , and provide estimates, up to numerical factors, of the total injection energy in the form of anti-quarks required to create the fireball corresponding to a parameter space point of interest.
The black star indicates the benchmark parameter point defined at Eq. (49).
The purple region is the parameter space where the fireball would thermalize at a temperature above the QCD phase transition (making it a quark–gluon plasma [QGP]), a regime we avoid to keep our analysis tractable.
In the the gray region, the would-be temperature of the fireball is so low that pions are too Boltzmann suppressed to facilitate the thermalization of the injected anti-quarks; i.e., .
Figure 5: Output space.
The injected anti-helium isotope ratio and the typical injected Lorentz factors of the anti-particles are two key observables we seek to explain with the fireball anti-nucleosynthesis scenario we propose in this paper.
Publicly available AMS-02 tentative data currently favor observed values of and (this figure does not include propagation effects; see Sec. IV).
Dashed and dotted lines show the injected anti-helium isotope ratio obtained by numerically solving the Boltzmann equations for a simplified nuclear-reaction network detailed in Appendix C, while dot-dashed lines show the analytical approximation, Eq. (47).
The black star indicates the benchmark parameter point defined at Eq. (49).
As far the numerical results are concerned, this figure is constructed as follows: for each indicated value of (as annotated in the legend), we vary with so as to keep [cf. Eq. (48)] fixed to the annotated constant value.
The analytical prediction on the anti-helium isotope ratio has a precise quadratic dependence
, and we do not therefore need to make further specific numerical assumptions about the values of and ; although it is only strictly speaking valid when [see discussion below Eq. (48)], we also show it here for , where it is still reasonably accurate.
The numerical computation of confirms the expected strong dependence on , but it also shows mild sensitivity to variation orthogonal to .
Anti-neutrons have a mean rest-frame lifetime of minutes (decaying to an anti-proton), while anti-tritium decays to anti-helium-3 via the beta decay with a rest-frame half-life of 12.3 years.
Even accounting for Lorentz factors , the and will decay on timescales (at most) as viewed in the fireball center-of-mass frame.
After these decays, the remaining light anti-nuclei outputs of our scenario are, at late time, given by
(40)
(41)
(42)
(43)
(44)
(45)
(46)
while all other products are negligible.
Note that, in Eqs. (42) and (44), we made use of Eqs. (34) and (41) to rewrite the factors of in Eqs. (36) and (39) in terms of , and .
We also inserted the numerical values of the coefficients here: , , and .
It is important to note that the analytically predicted relative scalings of and are found to be reasonably accurate in numerical results obtained from solving the Boltzmann equations (cf. Appendix C), at least for the hierarchy (cf. Fig. 3): numerically, we actually find a scaling somewhat closer to , but is found numerically to be very accurate.
However, the naïve absolute scaling with of these results that is implied by combining them with the analytically predicted scaling from Eq. (35) should be understood with some caution owing to the high powers to which is raised and the fact that we find numerically that is a more accurate scaling result, at least the same hierarchy of the (again, cf. Fig. 3).
That is, the naïve predictions would be and , whereas we observe numerical scalings in Fig. 3 more consistent with the and results that would follow from combining the quite accurate relative scalings predicted analytically with the more accurate numerical result [actually, would be a slightly more accurate result, reflecting the numerically obtained scaling ; cf. Fig. 3].
If we take , for which , , and , the anti-helium isotope ratio injected into the interstellar medium for a given is given by
(47)
where in the last line we have substituted from Eq. (35).
Hence, to obtain comparable anti-helium-3 and anti-helium-4 abundances, the fireball parameters must be such that the combination of parameters
(48)
is an number; see Fig. 4.
According to the simple analytical estimates of this section, the tentative AMS-02–observed anti-helium isotope ratio corresponds to .
Note that this estimate is made here for the injected ratio, without regard to the impact of propagation of the anti-nuclei in the Galaxy or the isotope-dependent AMS-02 sensitivity, which we discuss and account for in Sec. IV.
There, we will show that is required, corresponding to .
Note that the sequential production approximation we used to derive the analytical predictions for [Eqs. (40)–(46)] is justified if .
This translates to , or .
We have also numerically computed the anti-helium isotope ratios for different fixed values of using the set of Boltzmann equations described in Appendix C; see Fig. 5.
While these numerical results show reasonable overall agreement with the analytical scaling derived at Eq. (47), we find numerically that still varies mildly with if is held fixed.
Nevertheless, because we have tuned the numerical constants , , and in the analytical results to the numerical computations, we find that this mild violation of the scaling makes the analytical results for that isotope ratio at worst an factor discrepant from the numerical results throughout the viable parameter space.
We also find that there is still reasonable agreement of the analytical and numerical results for as large as , notwithstanding the limitation that we noted previously.
III.4.9 Non-nuclear outputs
In addition to anti-nuclei, the fireball evolution described above will result in the injection of light SM particles throughout the galaxy.
When the fireball becomes optically thin, a burst of photons is released, with an average energy .
We assume (and this is almost certainly the case) that this dominates the integrated emission from the photosphere throughout the previous expansion.
Since the pions present in the fireball remain in thermal equilibrium until at least , they are Boltzmann-suppressed prior to decoupling and therefore no significant gamma-ray signal is expected from their decay.
Anti-neutrinos are continually produced through weak interactions, via both inelastic weak interactions and decays.
The dominant scattering production occurs when the fireball first thermalizes and its temperature is the highest. At this point, all the produced anti-neutrinos (which possess an average kinetic energy ) will escape the fireball since their mean free path is , where is taken to be the (thermal) electron or positron number density since .
Additional neutrinos are produced as the pions within the fireball decay.
These have an average kinetic energy .
There may also be a neutrino contribution immediately on injection of the SM species that seed the fireball, but this is model dependent.
Charged leptons present in the fireball may also be injected into the interstellar medium after the plasma becomes optically thin. The number of such leptons remaining upon annihilation is model dependent; it is set by the initial charge asymmetry in the lepton sector. However, if the fireball as a whole is electrically neutral, this injection should be dominated by positrons with an average energy , whose number cannot exceed the total number of anti-baryons injected.
The observability of X-ray and lepton bursts is discussed in Sec. IV.2.
We find that for the benchmark parameters sufficient to explain the AMS-02 candidate anti-helium events, detection of these additional particles is infeasible due to the low count of particles arriving at the Earth and, in some cases, their low energies.
III.5 Summary
We have identified a parameter space (see Fig. 4) where a sudden and spatially concentrated BSM injection of energetic anti-quarks in our Galaxy triggers (subject to certain properties of the injection) a series of events, dictated purely by Standard Model physics, that lead to relativistic anti-helium anti-nuclei being released with number ratios and Lorentz boosts roughly consistent with AMS-02 observations (we discuss in Sec. IV how propagation effects modify the observed number-ratios from the injected values we have thus far discussed).
Here, we summarize this predicted series of events and provide benchmark values for key quantities at various points of the process; we denote these benchmark quantities with a tilde and an appropriate subscript.
Following the anti-quark injection, the anti-quarks rapidly hadronize and thermalize mainly via strong and electromagnetic processes into an optically thick, adiabatically expanding fireball with conserved anti-baryon–to–entropy ratio .
This fireball then undergoes a period of rapid acceleration which turns it into a plasma shell moving with relativistic radial speed, with an average Lorentz factor .
Right at the onset of this phase of its evolution, the temperature, outer radius, and thickness of the shell are in the ballpark of
(49)
The subsequent evolution of the plasma shell proceeds in three sequential stages (see Fig. 1):
1.
Acceleration: the shell continues to radially accelerate under its own thermal pressure with its average Lorentz factor increasing linearly with the shell’s outer radius , and keeping its thickness approximately constant, .
2.
Coasting: as is approaching close to the terminal radial bulk Lorentz boost , the shell enters the second stage of expansion where it simply coasts with an approximately constant Lorentz factor , again keeping its thickness approximately constant, .
3.
Spreading: the expansion timescale becomes long enough that the radial velocity difference between the innermost and outermost layers of the shell causes the shell’s thickness to increase significantly over time.
We now describe how the fireball’s particle content evolves as it expands and cools down.
Initially, while anti-nuclei are still absent due to rapid photodissociation of , anti-neutrons and anti-protons are kept in detailed balance by pion-mediated interconversion processes such as .
This continues until they finally decouple at a comoving temperature , at which temperature their relative abundance freezes out at202020
For the purposes of this summary discussion, we are assuming the appropriate charge asymmetry on the hadronic sector is achieved at injection (i.e., that at ); qualitatively similar results are however obtained so long as anti-neutron–to–anti-proton ratio remains .
See discussion in Sec. III.4.3 and Appendix B.
.
Rapid photodissociation of any fusion-produced ceases only deep in the final (spreading) expansion stage, when the comoving temperature, outer radius, and thickness of the plasma shell are about
(50)
The following nuclear reactions then proceed to produce light anti-nuclei (see Fig. 2):
•
production through .
•
production either (1) directly through , or (2) indirectly through , followed by the highly efficient .
The latter process depletes and keeps its abundance low.
•
production through .
In a way somewhat analogous to how dark matter is produced in freeze-in scenarios, these processes sequentially produce nuclear anti-particles with their final (frozen) numbers satisfying , and essentially no other elements (see Fig. 3).
As the shell further expands and decreases in density, the bulk of the plasma eventually becomes transparent to photons when its comoving temperature, outer radius, and thickness are around
(51)
At that point, the relativistic anti-nucleosynthetic products and a burst of X-ray photons are released from the plasma shell.
While traversing the interstellar medium, the decay to and the decay to (these decay timescales are very short compared to the galactic dwell-time).
Each fireball seeded with a total anti-baryon number therefore contributes to the nuclear anti-particle population in the interstellar medium as follows:
(52)
Note that these specific numerical results depend on the the benchmark values of the parameters that were chosen at Eq. (49) such that the resulting injected ratio of anti-helium isotopes reproduces the current AMS-02 candidate-event observed value of , and the terminal Lorentz boost of the plasma shell is at injection (see Fig. 5).
These injected values are however somewhat modified by Galactic propagation effects that we discuss and account for in the next section.
IV Propagation and detection
In the previous section, we showed how anti-nucleosynthesis occurring in an expanding thermal fireball state characterized by a certain temperature, radius, anti-baryon content, and net hadronic charge asymmetry could generate, after unstable elements have decayed, both and in an isotopic ratio broadly consistent with the candidate AMS-02 events.
In this section, we discuss how the properties of the anti-helium (and other species) injected by such fireballs at locations within the Milky Way (MW) are processed by propagation from the source to the AMS-02 detector, as well as the necessary parameters to generate event rates consistent with the candidate AMS-02 observations.
Our analysis is predicated on the following basic assumptions: (1) all seeded fireballs have similar , , and parameters (and hadronic charge asymmetries); and (2) a large enough number of anti-helium producing fireballs have been, and continue to be, seeded at random times up to the present day that we can neglect both spatial and temporal clumpiness in the injection and instead model it as a temporally constant and spatially smooth source.
Additionally, motivated by having dark-matter collisions seed the fireballs (see Sec. V), we assume that (3) the spatial distribution of the injections is , where is a Navarro–Frenk–White (NFW) profile [62],212121Following Ref. [14], we take the NFW profile to be normalized such that the local average DM density is [63, 64] at [65], and use a scale radius [14]. so that the anti-helium injection is occurring dominantly within the Milky Way itself and is peaked toward its center.
IV.1 Cosmic rays
The fireball injection model discussed in Sec. III is such that all cosmic-ray species are initially injected in the vicinity of the fireball with a narrow range of velocities centered around the bulk Lorentz factor222222Although many anti-helium nuclei are produced through decay processes (and never thermalize with the fireball) [e.g., the bulk of production is from decay], the associated nuclear decay values are sufficiently small that the product nuclei are always non-relativistic in the decay rest frame. As a result, their Lorentz factor in the galactic rest frame at the time of injection does not differ appreciably from . .
The subsequent motion of these injected cosmic rays to Earth (and hence the AMS-02 detector in low-Earth orbit) is of course diffusive in both position and momentum space [66].
In principle, we should thus pass the fireball-injected species to galprop [67, 68] to solve the necessary transport equations and account for various propagation effects; see also Ref. [12].
Instead of using a modified version of galprop to study the propagation of injected anti-particles, we argue as follows.
Galprop natively solves the transport equation for positively charged nuclei.
Of course, the opposite sign of the charge for the anti-nuclei does not impact the diffusive nature of the transport [66, 14].
To be sure, there are additional annihilation reactions that can occur for anti-particles interacting with the (dominantly ordinary matter) interstellar medium (ISM), and the inelastic cross-sections scattering with the ISM also differ somewhat for particles vs. anti-particles.
Nevertheless, at the level of precision at which we work, the annihilation cross-sections can however reasonably be ignored for our purposes, as they constitute a negligible correction to the total inelastic cross-section of the anti-particle species at energies [69, 70, 71] and can therefore be absorbed into the uncertainties associated with the propagation [14].
Ignoring also any other differences in the anti-particle vs. particle inelastic cross-sections for interaction with the ISM, we employ galprop results for the corresponding positive charged nuclei as an approximation to the desired results for the negatively charged anti-nuclei (e.g., we inject primary instead of , and read off results accordingly, etc.).
This approximation could of course be revisited; however, as we shall see, the results of this approximate treatment indicate a ratio of observed anti-helium fluxes that differs by only an numerical factor as compared to the ratio injected by the fireballs; it is therefore unclear whether a modification of propagation code as in Ref. [12] to more correctly treat the anti-particle propagation is justified given other, larger uncertainties in our scenario.232323We note that such modification was however important for Ref. [12] to achieve accurate results, as one of the main issues addressed in that work was to refine predictions for the fully propagated anti-particle secondary fluxes produced by the primary ordinary matter cosmic-ray spectra.
Specifically, we model the injection of anti-cosmic rays of species by specifying galprop source terms for the corresponding positively-charged, ordinary-matter species, which we denote here as :
(53)
where is taken to be the isotopic abundance for the anti-particle species from Sec. III.4, is taken to be a narrow top-hat function centered at the rigidity corresponding242424While the fireballs inject each species at a single rigidity, the width of the top-hat (chosen here to be of the central value for numerical reasons) is inconsequential as long as it is subdominant to the momentum-space diffusion occurring during propagation, which we verify a posteriori.
In order to specify a fixed , and hence a different dependence for each species, galprop had to be run multiple times: in each run, a single species was injected at the required rigidity ; the output spectra from each such run were then summed with weights . to Lorentz factor for species as injected by the fireball, and is the NFW profile.
We of course then also read off local flux results for the species , and impute those to species .
The source terms are normalized such that the imputed total injection rate of anti-baryon number, summed over all species and integrated over the whole galprop simulation volume, is .
Furthermore, we approximate the decay of to and to as occurring instantaneously at the fireball location, and we thus consider only of the anti-particle species when running galprop via the above procedure, in the ratios specified at252525As discussed in Sec. III.4, while the numerical solutions to the Boltzmann equations are more accurate than the analytical results at Eqs. (40)–(46), the tuning of those analytical results to the numerics via the constants , , and , makes the analytical results sufficiently accurate for our purposes here, notwithstanding the minor violation of the scaling of discussed in Sec. III.4.8; see also Fig. 5.
Eqs. (40)–(46) as a function of .
For the galprop diffusion model, we adopt the propagation parameters for “ISM Model I” in Ref. [72]; within the range of parameters consistent with existing cosmic-ray observations, our results are largely insensitive to the choice of transport model.
We also neglect the modulation of the cosmic-ray fluxes at Earth due to the heliospheric magnetic field: in the force-field approximation, solar modulation is governed by a single parameter known as the Fisk potential which, at high energies, corrects the flux at most by a factor [73].
Note also that the typical mean-free path for an anti-nucleus traveling with is , meaning that it stays within the galaxy for a duration where is the thickness of the MW disk.
A set of example post-propagation spectra at the position of the Earth is shown in Fig. 6, alongside AMS-02 sensitivity curves or observations.
In each case, the propagated flux of injected anti-particles peaks at kinetic energies corresponding to the Lorentz parameter [i.e., ].
Figure 6: Predicted spectra of anti-nuclei upon propagation to the Earth as a function of kinetic energy per nucleon , for different choices of the initial Lorentz factor , as annotated by the different colored lines.
The best available sensitivity of AMS to anti-helium (top panel) and anti-deuterium (middle panel) events [74, 75, 76, 77, 78] are presented for comparison in gray (see text).
In the top panel, the solid and dashed lines correspond to the isotopes and , respectively.
The anti-proton flux observed by AMS is shown (bottom panel) by a solid black line [79], while the additional flux required to exceed the uncertainties on the measured flux is shown as the dashed black line in that same panel; fireball production is a negligible source of galactic at these parameter values.
Throughout this figure, we fix , while the isotopic abundances at the fireball are set by ; see Eqs. (42) and (44).
The expected number of anti-helium events accumulated at AMS over 10 years (at the 95% CL) can be computed using the output spectra from galprop and the published AMS sensitivity to the flux ratio between anti-helium and helium [74, 75].
This is done by recasting the anti-helium acceptance of each energy bin in terms of this sensitivity in conjunction with published helium data [76], following the procedure outlined in Appendix B of Ref. [15]; see our Appendix G for a brief review.
Assuming the number of events follows Poissonian statistics, and allowing for the joint probability of the predicted numbers of and to deviate from their AMS-02 tentative observed values within the 68% confidence interval, the parameter space consistent with the AMS-02 candidate anti-helium events is shown in Fig. 7.
Explicitly, we require [80]
(54)
The events for each isotope are assumed to be independent, satisfying
(55)
where is the number of events for species predicted by galprop given a set of model parameters , whereas is the fiducial number of candidate events reported by AMS-02.
For the benchmark shown by the black dot in Fig. 7, the isotope ratio changes from at injection to an observed value of – (depending on the choice of ).
Part of this change in the isotope ratio from injection to observation has to do with physical propagation effects.
For instance, spallation of onto the interstellar medium; however, only % of the abundance arises from this effect, which moves the isotopic ratio by only an factor.
Moreover, that effect would tend to drive the ratio in the other direction (i.e., it reduces and increases ).
We also expect there to be some mild differences in diffusion for different-mass isotopes at the same kinetic energy per nucleon owing to different rigidities.
However, we identify the larger part of the change to arise not from propagation effects at all, but rather from the fact that the AMS-02 anti-helium sensitivity, which we take to be given as the same function of rigidity for all anti-helium isotopes (see Appendix G), is not flat as a function of ; on the other hand, the post-propagation energy-per-nucleon spectra of and are almost the same (since they have the same energy-per-nucleon at injection), resulting in rigidity distributions for the two species that peak at different values of (i.e., when stated in terms of the kinetic energy per nucleon, the AMS-02 sensitivity that we assume differs for species with different charge-to-mass ratios; see the top panel of Fig. 6).
Ultimately, however, this -factor change could easily be absorbed into a slightly different parameter point if any of the assumptions leading to this effect are found to be inaccurate.
For the case of anti-deuterium, a robust detection relies on the rejection of backgrounds (particularly and He) that are substantially more abundant.
As described in Ref. [78], the latest published AMS-02 anti-deuterium sensitivity curve (shown in Fig. 6) is based on an earlier superconducting-magnet configuration for AMS-02, rather than the permanent-magnet configuration actually in use, and cuts off at .
Nevertheless, taking that sensitivity, only injections taking place with would result in expected events in the parameter space of interest, but it would be challenging on the basis of that sensitivity to simultaneously account for the 7 candidate anti-deuterium events reported in Ref. [7] in addition to the candidate anti-helium events.
That said, an updated study of the AMS-02 anti-deuteron sensitivity would be required to accurately determine whether a single choice of parameters could achieve this.
Finally, anti-proton events arising from the fireball are also expected to be observed at AMS-02; indeed, as seen in Fig. 6, this flux dominates those of the other species in the parameter regime of interest.
However, conventional astrophysical sources are responsible for an anti-proton flux that is greater than that created by the fireballs by a factor of [79].
Producing sufficiently many events to exceed the uncertainty on this measurement (tantalizing in light of the anti-proton excess of at similar energies [72]) would require a greater anti-nucleus injection rate than is favored by the candidate anti-helium observations; see however Fig. 9 in Appendix B (and related discussion) for an alternative parameter point that may be more interesting from this perspective.
In contrast to the heavier species, secondary production contributes significantly to the flux at low kinetic energy per nucleon, resulting in an approximately flat spectrum that falls off only for (i.e., below the range shown in the lower panel of Fig. 6).
However, this plateau is observationally unimportant as compared to the peak in the spectrum arising from the primary fireball-injected , as the former lies further below the detectable flux level than the latter.
This feature is absent from the spectra for heavier antinuclei in the ranges plotted owing to the inefficiency of secondary production of heavy nuclei (as discussed in, e.g., Sec. I).
Figure 7: The region of parameter space for which AMS-02 would be expected to observe 3 events of and 6 events of in years, with the injection taking place at .
(The allowed region is approximately identical over the range .)
The anti-nucleus injection rate normalizes the total rate of fireball injection of anti-cosmic rays into the MW, whereas the individual source isotopic ratios are set by the fireball anti-deuterium abundance : via Eqs. (40)–(46).
The predicted number of anti-helium events is allowed to independently vary from the fiducial quantity within the 68% confidence interval of the joint probability distribution [see Eq. (54)].
The black dot shows the fiducial parameter values referred to in Fig. 6 and Eq. (61).
IV.2 Other signatures (indirect detection)
In this section, we examine the detectability of X-ray, anti-neutrino, and positron bursts which arise from the fireball model (discussed in Sec. III.4.9).
While fireball injection of anti-nuclei is envisaged as a continuous process compared to the galactic dwell-time for these (diffusively transported) particles, the rapid expansion timescale and low galactically-integrated injection event rates involved (see Sec. V.1 for this estimate) mean that the non-nuclear outputs from independent fireball injections are temporally well-separated at the Earth and so are best treated independently (with the exception of positrons, whose transport is also diffusive).
At the fluxes required to explain the AMS-02 candidate anti-helium events, there should be events occurring within the galaxy over .
Suppose that we made an observation to look for their non-nuclear products that has a duration and that observes a fraction of the whole sky.
In this time, we estimate that the closest observable injection event would occur at a distance from Earth, if we assume that the injections arise from collisions of the sub-component composite DM states whose benchmark parameters we discuss below in Sec. V.1.
Photons are rapidly released in a burst when the fireball arising from the injection becomes optically thin, at size ; see Eq. (51) for this benchmark.
We estimate the energy released in this burst to be [using the benchmark at Eq. (51), and ], which results in the emission of photons with average energy .
X-ray telescopes sensitive to such events have an effective area [81, 82], a field of view of at most square degrees (i.e., ), and mission lifetimes of .
The total expected number of photons arriving at the detector from the closest release that would occur would be .
Even if we conservatively ignore the finite field of view (i.e., set ) and take , this is too faint to observe.
There may additionally be a diffuse background of X-ray photons due to the cumulative injections taking both intra- and extra-galactically over extended periods of time, but we have not estimated this here.
Anti-neutrinos may be injected promptly when the fireball thermalizes, or indirectly through the decay of injected pions.
The fraction of the total energy [see Eq. (49) for this benchmark; we set ] that is carried by the prompt anti-neutrinos is (for as at early times; see Fig. 1), corresponding to particles with an average energy of .
As far as (anti-)neutrinos from decay are concerned, thermal pions have a number density within the fireball, so the number of neutrinos injected due to their decay is , where is the charged pion lifetime.
Their average energy is also .
The decay contribution is thus sub-dominant.
Given a neutrino detector in this energy range, with an area [83], roughly neutrinos would pass through the detector as a result of the nearest injection in the observation time.262626We took here to reflect that Earth is quite transparent to neutrinos with energies below a few TeV, so the detector has full-sky coverage; see, e.g., Ref. [84]. The probability of a single neutrino interacting in the detector is where is the target (e.g., water [83]) number density, and is the target thickness.
Taking and as relevant for our assumptions, we find , so the expected number of detectable events272727By way of comparison, the Kamionkande-II detector detected only 11 neutrino interactions of energies of from SN1987A [85], which released neutrinos of all flavors [85] at a distance of [86], leading to electron anti-neutrinos passing through the detector area (varies slightly depending on the orientation relative to the source) [85].
Taking and gives , which is broadly consistent with these numbers, demonstrating the consistency of our estimate in the text.
is for any conceivable observation duration.
Positrons injected by the fireball are potentially observable in two ways: directly as cosmic rays, and indirectly through photons produced when they annihilate with electrons in the ISM.
However, the injected positron energies are , which falls below the AMS-02 sensitivity threshold (even without accounting for further energy loss during propagation), likely making them unobservable.282828Unsurprisingly given that solar-modulation effects severely impact sub-GeV positrons [87], this energy also falls below the lowest-energy positron measurements reported by PAMELA [88], HEAT [89], CAPRICE94 [90], AMS-01 [91], or FERMI [92, 93].
Relevant Voyager 1 data for the sum of electrons and positrons down to energies of from periods after it crossed the heliopause are available [94] (see also Ref. [95]); we do not however pursue these constraints further in this work.
For the 511 keV emission, the total integrated positron injection rate from fireballs cannot exceed the total integrated anti-baryon injection rate (cf. the discussion about in Appendix B), which from the numbers shown in Fig. 7 is maximally ; this is much smaller than the integrated positron injection rate of that explains 511 keV emission in the MW [96].
The fireball-injected positrons annihilating therefore likely only contributes a sub-dominant 511 keV flux.
On the basis of these estimates, these non-nuclear products thus do not appear to be observable; however, a more detailed investigation of these (indirect-detection) signatures of this class of models may be worthwhile in future work.
V A dark-matter origin for the fireballs?
Thus far, we have operated under the assumption that a rapid, localized injection of energetic Standard Model anti-baryons can be achieved in order to seed the fireballs.
A BSM mechanism is required to explain these injections.
In this section, we discuss whether the collisions of large, composite dark-matter states that carry anti-baryon number may be able to provide such a mechanism.
Assuming that a substantial fraction of the mass energy of such colliding DM states can be promptly converted to SM anti-quarks as a result of dynamics triggered by the collision, we demonstrate in Sec. V.1 that the requisite injection rate of anti-baryon number could be achieved for DM states with certain bulk physical properties (i.e., total mass, number of constituents, and physical size).
In Sec. V.2 we then show that, provided certain benchmarks can be realized in the conversion of the dark-state mass energy to the SM, fireballs with appropriate bulk parameters could also be seeded.
As this is suggestive, in Sec. V.3 we then advance some speculations toward specific microphysical DM models in which we suspect this could possibly occur.
We emphasize, however, that we have not settled the question of whether one or more of these models actually do in fact realize the required dynamics.
While we intend to return to this open question in future work, we also encourage other work on this point.
V.1 Collision rates
Suppose that a fraction of the DM energy density is comprised of large, cosmologically stable composite objects.
The existence, formation, and signatures (or lack thereof) of a variety of objects of this type have been subjects of extensive study in the literature [97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117] (see also Refs. [118, 119]).
For the purposes of the estimate we give here, we assume292929Slightly different assumptions could be made (and indeed are made in the following subsection), but would mostly just lead to parameter-space remappings in the discussion that follows here: for instance, we could take the baryon number per constituent to be something other than , but still of ; e.g., per constituent.
Likewise, the total mass of the composite object could get corrections from binding energy or relativistic motion of its constituents.
that these objects have a mass dominated by constituents, each with mass and baryon number of , giving them a total mass and total baryon number ; we take their radius303030Note that the initial fireball radius discussed in Sec. III can be very different from the dark-blob radius .
to be .
Taking a collision cross-section
(56)
where is a Sommerfeld-like enhancement to the geometrical cross-section, the collision rate of two such blobs integrated over the whole simulation volume used in galprop (see Sec. IV.1), can be estimated as313131The blob collision rate can be significantly modified if the blobs formed binaries in the early Universe [120, 121, 122].
(57)
(58)
(59)
At Eq. (57), is the DM energy density at galactocentric radius , as specified in the NFW model [123] (see Sec. IV.1), and is the typical velocity of blobs in the galaxy (we take this to be constant over ).
At Eq. (58), we have have defined
In the above numerical estimate, we considered benchmark blobs with radius , number of constituents , and constituent mass , giving a total blob mass of (roughly the mass of a typical Main Belt asteroid with a diameter of a few-hundred meters) and an average density (while much more dense than ordinary matter, this is still less dense than a non-extremal SM white dwarf).323232Such an object, while very dense, is still many orders of magnitude larger than its own Schwarzschild radius: . We have also assumed that they constitute only about % of the DM, and have conservatively neglected any Sommerfeld-like enhancement to the collisions rate.
We assume further that, upon such pair-wise collisions, the DM states can be destabilized in such a way that a fraction of their combined number of constituents is converted to Standard Model anti-quarks, seeding into the SM a total baryon number per collision.
The injection rate of anti-baryon number (cf. Fig. 7) can then be estimated as
(60)
(61)
Importantly, for these benchmark bulk DM parameters, we find that this could achieve the injection luminosity required to obtain the correct fiducial number of candidate anti-helium events at AMS-02, (cf. Fig. 7), for .
That is, the requisite injection rate could be achieved even with large composite DM that is a sub-component of the total DM density and that is not very efficient in converting constituent anti-baryon number to SM anti-quarks upon collision, all without any Sommerfeld-like enhancement to the collision rate.
Of course, these conclusions change for different DM bulk parameter ranges, per Eq. (61).
Note also that the typical inter-collision time here, , is orders of magnitude shorter than the diffusion time for the charged cosmic rays to escape the MW [72], which also justifies a posteriori treating the injection as a roughly constant-in-time source term in galprop (see Sec. IV.1).
V.2 Mapping to fireball parameters
While we cannot robustly estimate the fireball parameters that would obtain from such collisions without a model, we can show that these dark composite DM states would be able, at the order of magnitude level, to generate approximately the benchmark fireball parameters discussed at Eq. (49).
If we assume that is also the fraction of energy injected into the SM anti-quarks, then the energy injected is .
If we take (verified a posteriori), this energy must go to the mass-energy of anti-nucleons (they dominate the fireball anti-baryonic output) and the remainder to relativistic species.
That is, .
Therefore,
(62)
(63)
where we have assumed that the fireball injection happens in a region of radius333333We limit our analysis in this paper to the “burst” injection regime in which the timescale for injection is smaller than the region of size over which the injection occurs: ; we have also taken for simplicity, although we could also have in this regime.
However, the fireball anti-nucleosynthesis should proceed similarly in the opposite, wind regime (), but with some differences in how the injection properties are related to the properties and outputs of the resulting thermalized fireball.
Parametrically, one can think of the anti-quark injection in the wind regime as a sequence of anti-quark bursts occurring in a region of size continuously one after another, each with an injection energy and anti-baryon number smaller by a factor of relative to their total values and .
Consequently, in the wind regime the resulting initial fireball temperature will be lower by a factor of and at the same time the overall anti-nuclei output for the same will be times higher compared to that of the burst regime.
We refer the reader to Appendix F for more discussion on this point.
that is much smaller than (i.e., ), consistent with that injection occurring as a result of some catastrophic collapse dynamics in the post-collisional evolution (see Sec. V.3).
Meanwhile, in the final post-fireball evolution, we have roughly that , so , corresponding to .
Moreover, the constant defined at Eq. (48) is , implying an injected anti-helium isotope ratio of , roughly consistent with an observed isotope ratio of , as discussed in Sec. IV.
In order to achieve the requisite net negative charge on the hadronic sector of the fireball in order to obtain a value around the onset of anti-nucleosynthesis (see the discussion in Sec. III.4.3 and Appendix B) while also maintaining net fireball EM neutrality (as required for injections from net-neutral dark objects), it would be sufficient for the injection to take place in such a way that a fraction of the anti-baryons are anti-protons, with a corresponding number of positively charged leptons being injected.
On energetics grounds, because and , this can easily be accomplished with minimal change to the above discussion.
It would however require the initial dark object to carry net-negative lepton number in addition to net-negative baryon number.
While the precise values shown in both this subsection and the previous one should not be taken too literally given the roughness of the estimates, what this nevertheless indicates is that collisions of large, macroscopic dark objects can in principle attain both the requisite injection rates and required fireball conditions to make our anti-helium anti-nucleosynthesis mechanism operate in a phenomenologically viable fashion, provided that the benchmarks and broad model features we have given can be attained in a concrete model.
V.3 Towards a particle physics model
The development a detailed microphysical model for the composite dark states, and their evolution in and after collisions is a question we ultimately leave open in this paper.
Nevertheless, in this section, we sketch the outlines of two particle physics models that we believe are promising targets for future investigation, and indicate where the difficulties in understanding their evolution lie.
We intend to return to this point in future work, and also encourage other work on it.
We also note at the outset of this discussion that obtaining a model with a consistent and observationally allowed cosmological evolution in light of the greatly increased density of dark matter in the early Universe343434We remind the reader that the local DM density in the MW, [63, 64], is approximately 5.5 orders of magnitude larger that the present-day average cosmological abundance, [124], which is in turn roughly 10.5 orders of magnitude down from the abundance at matter-radiation equality (MRE), .
The DM density at MRE is thus approximately 5 orders of magnitude larger than locally in the MW today. is an additional model-building constraint that would need to be carefully evaluated in the context of a full model, taking into account also the redshift dependence of the average speed of the collisions.
V.3.1 Imploding fermion+Yukawa model
The requisite injection may result from rapid conversions of fundamental DM particles , each carrying baryon number ( [and possibly lepton number], into SM anti-particles via a higher-dimensional effective operator; e.g., something of the form if only anti-baryon number needs to be injected, or something like if both anti-baryon and anti-lepton number need to be injected primarily (here, and represent an anti-quark and anti-lepton, respectively); see discussions in Secs. III.1 and III.4.3, and Appendix B.
Simultaneous and spatially concentrated conversions of a large number of particles can be naturally achieved if the particles participating in the conversion existed in the form of blobs; i.e., large composite bound states.
While these blobs, presumably formed in the early Universe, must survive over timescale until the present epoch, they must also rapidly convert to anti-quarks in timescale when the time is ripe (within the last timescale for anti-nuclei to escape the Galaxy).
To bridge these extremely long and extremely short timescales required for the DMSM process, we can imagine the following scenario.
These blobs may have existed in a metastable state in the sense that they are individually stable, with the rate of the DMSM conversion process within each blob satisfying , but the merger of a pair of blobs would trigger a runaway collapse of the merger product accompanied with many orders of magnitude increase in , eventually causing most of the particles to convert to SM anti-particles when .
The stable blob can be kept in equilibrium, for instance, by the balance between the repulsive effect of the Fermi pressure of particles and the compressing effect due to a higher vacuum pressure outside compared to inside the blob.
Such a vacuum pressure difference arises naturally in scenarios where the blobs are formed through a cosmological first-order phase transition from a higher-energy false vacuum to a lower-energy true vacuum in the early Universe [108, 97, 106, 109, 119, 118].
In those scenarios, the present-epoch setup is that the inside of the blobs remains in the false vacuum while the rest of the Universe is in the true vacuum.
Correspondingly, the pressure associated with the vacuum energy is more negative inside than outside the blob.
Runaway collapse can be achieved by incorporating additional Yukawa forces between particles with an intermediate-range mediator: long compared to the typical spacing between but short compared to the radius of the blob.
We additionally assume that the are non-relativistic to avoid suppressions in their Yukawa forces (cf. Refs. [100, 102] and the discussion in the next subsection).
The total Fermi kinetic energy of the blob, the vacuum energy difference between inside and outside the blob, and the attractive Yukawa potential for an intermediate-range mediator then scale with the blob’s number of constituents and radius as , , and , respectively.
For a given sufficiently small , the blob’s total energy as a function of has a metastable minimum at a radius that is set by the balance between the non-relativistic Fermi pressure and vacuum-pressure difference.
However, when exceeds a certain threshold the metastable minimum ceases to exist and it becomes energetically favorable for to decrease indefinitely, until new effects that deplete the or reverse the tendency to collapse turns on.353535We note also that this short-range Yukawa regime has different qualitative behavior as compared to the cognate “saturated” regime considered in Ref. [102] because the fermions we consider here are still non-relativistic, whereas those considered in Ref. [102] are relativistic by the time the Yukawa becomes short range compared to the size of the blob.
This modifies both the scaling of the degeneracy pressure with blob radius and the behavior of the Yukawa force, as compared to the saturated case in Ref. [102].
Such a runaway collapse can be triggered by the merger of two near-critical blobs.
Aside from the complexity of accurately modeling the envisaged scenario, one can already see competing effects or constraints that make it challenging to present a concrete realization of this scenario.
First, the requirement to reduce the blob radius by many orders of magnitude during the blob collapse in order to increase from to can increase the Fermi momentum of the particles to the point that they become relativistic.
Once the particles are relativistic, the Yukawa forces become progressively suppressed [100, 102] and at the same time the Fermi pressure increases faster as the blob radius decreases [125].
Both these effects tend to stop the collapse, potentially preventing from reaching .
Moreover, consistency with collider data (and to obtain ) requires raising the cutoff of the higher-dimensional DMSM operator to at least above .
The latter implies a strong suppression on unless the blob becomes very dense at the time of conversion.
Preliminary estimates suggest that in order to obtain , the energy density of the blob at the time of conversion must be , which is much higher than the typical energy density of the fireball we seek to create.
To dilute such a high energy density one may need to introduce another operator that converts the particles to a lighter dark particles much more efficiently than the DMSM conversion, such that most of the convert into with a rate , and only a small fraction of with energy density converts to SM anti-particles during the timescale.
V.3.2 Dark-dwarf model
The second model that we believe may be interesting to consider is based on the dark states behaving somewhat like white dwarfs in the SM.
The particle content of this model would be a heavy “dark proton” , a slightly heavier “dark neutron” , a much lighter “dark electron” , and a massless “dark neutrino” .
The and carry baryon number , and both couple with equal strength to a Yukawa force mediated by a scalar .
Additionally, the and are taken to both couple to a dark electromagnetism, with opposite charges.
We also assume the existence of a two higher dimensional 4-fermion operators allowing (and crossed) reactions, and permitting , where is a SM anti-quark.
With such particle content, it is possible to form dense macroscopic composite objects (“dark dwarfs” [109]) with very large constituent number, consisting of held together by the Yukawa force, and supported against collapse by the (relativistic) degeneracy pressure of the .
The mass hierarchies of the particles can be arranged such that, at zero temperature, these objects would initially contain no or constituents.
Similar objects, albeit with just a single uncharged constituent fermion supported against collapse due to attractive Yukawa forces by its own degeneracy pressure, have been considered in, e.g., Refs. [99, 100, 102, 105]; see also Ref. [109].
A straightforward extension of these works to the two-component dark dwarf we have in mind indicates that these objects, while not possessing a singular collapse instability (i.e., there is no true Chandrasekhar-like limit),363636The avoidance of the singular collapse instability arises from the suppression of the coupling of scalar-mediated forces to relativistic fermions, as was pointed out in Ref. [100] and examined further in Ref. [102].
possess of threshold constituent number for which the following is true: for , their radii are large, say ; on the other hand, for , their radii are much smaller.373737One way to understand this is that the relationship for behaves similarly to the mass–radius relationship for a SM white dwarf (as computed in purely Newtonian gravity, using the full equation of state for degenerate electrons) [125]: a slow function of for , but showing a vary rapid decrease toward for .
However, very near , the become relativistic, leading to the softening of the Yukawa coupling to the discussed in Ref. [100, 102], which avoids a singular crunch but leads to a minimum size near and a slow (power law) increase in away from above .
Roughly, the hierarchy of these radii is in ratio of the dark electron and proton masses: .
If the cosmologically stable dark objects have with , then collisions of two such objects could in principle trigger a large change in the bulk properties of the pre- and post-collision states because (assuming no ejection of constituents), the post-collision object would than have constituents.
If the collision were to play out in such a way that both the pre- and post-collision states were at exactly zero (dark) temperature, the idea would heuristically be the following:
before collision, the Fermi energy of the (relativistic) is much smaller than the mass difference of the and .
However, when these objects collide and then dynamically relax down to their much smaller post-collision equilibrium size, the density increases greatly, leading to their Fermi energy becoming much larger.
Once the Fermi energy exceeds the – mass difference, a large fraction of the can suddenly “neutronise” to (we adopt the nomenclature of the cognate process that occurs in the collapse toward a SM neutron star [125]).
The in turn are unstable to decay to SM anti-quarks.
As a result, the sudden dynamical “collapse” of the incoming dark states to a much smaller and more dense state upon collisions could trigger a catastrophic decay of a large fraction of the mass energy of the initial colliding dark objects to SM anti-quarks.383838Of course, the existence of decay channels for the dark dwarfs through annihilation reactions that proceed through an off-shell must also be considered both for blob stability and for understanding the rapidity of decay to SM quarks during this collapse.
Unfortunately, the dynamics of this system are actually more complicated than this naïve picture.
Initial estimates indicate that for scalar forces with a range larger than the blob sizes (the regime we assume), there is a speed up of the dark states as they collide, leading to the collision taking place in a supersonic regime with respect to the estimated speed of sound in the blobs.
As a result, two things are true: (1) the collision is a violent process that is not amenable to simple analytical treatment, and could lead to the ejection of constituents, or the formation of a large bound object, or (more likely) something more complicated; and (2) as the collision occurs, it is likely that the blob constituents are heated significantly (via, e.g., shock heating), possibly to initial post-collision dark temperatures that are a significant fraction of the dark electron mass .
As a result, the zero-temperature dynamics of the object are likely significantly modified by resulting thermal pressures.
It appears that the resulting state may still lose energy and shrink via surface dark-photon emission, volume relativistic emission of the particle, and later potentially either volume or surface emission of dark neutrinos.
The resulting state may thus still evolve toward the neutronization regime sufficiently rapidly for this model to be viable.
However, the entirety of this evolution is complicated and requires further study to establish in detail whether it works.
While this model may potentially be able to achieve the requisite destabilization dynamics, obtaining a net-negative EM charge on the anti-hadrons for fireballs seeded by this model (see discussions in Secs. III.4.3 and V.2, and Appendix B) may be challenging in the exact formulation advanced here.
A modification to this picture, or further model building, may be required.
The scope of work required to investigate this model in full and place the speculative statements in this sub-section on a firm footing is such that we defer its consideration to a future paper.
VI Discussion and conclusion
In this paper we have studied the anti-nucleosynthesis of elements up to anti-helium-4 in rapidly expanding fireballs of SM plasma that carry net anti-baryon number and that we assume to be seeded within the Milky Way by a BSM process.
For appropriate initial conditions set by the initial radius, temperature, and anti-baryon–to-entropy ratio of the plasma, the evolution of these fireballs is such that their thermal pressure drives the system toward a regime where there is relativistic bulk radial motion of a thin shell of plasma, in which the temperature of plasma falls as the expansion proceeds.
This permits purely SM thermal anti-nucleosynthesis of elements (similar to BBN) to occur in the expanding, cooling thin shell, while the products obtain relativistic boosts with respect to the rest frame of the fireball.
Eventually, for appropriate parameters, the expansion rate shuts off the anti-nucleosynthesis in a regime where the anti-nucleosynthetic products have not reached their thermodynamic equilibrium values, which allows the abundances of and anti-tritium (which later decays to ) produced to not be too dissimilar, and not highly suppressed with respect to anti-deuterium or anti-protons; see Figs. 3–5.
Once this expanding fireball becomes optically thin, its products cease to be driven by thermal pressures, and are launched at relativistic speeds into interstellar space, all with the same bulk speed (i.e., Lorentz boost).
The anti-tritium subsequently decays to on a short timescale, leading to injected amounts of and that are not too dissimilar.
Assuming that these injections of high-energy anti-nucleons and anti-nuclei are spatially distributed as the square of an NFW profile, and making use of galprop to approximate the galactic transport of these injected products, we showed that the relative fluxes received at Earth are such that one could obtain a ratio of roughly 2:1 for the isotopes to , without being excluded by AMS-02 anti-deuterium or anti-proton constraints; see Figs. 6 and 7.
We showed that other products of this scenario (e.g., photons) do not appear to supply additional constraints, but we did not undertake an exhaustive indirect detection study.
We also computed the overall required injection luminosity of anti-baryons into the fireballs we considered in order to obtain anti-helium events in the AMS-02 10-year exposure.
Our conclusion that the sizes of the anti-proton (larger) and anti-helium (smaller) fluxes are not too dissimilar so as to violate observational constraints on the anti-proton flux, while also explaining the anti-helium isotope ratio, does depend to some extent what the anti-proton–to–anti-neutron ratio is in the partially thermalized fireball plasma prior to the onset of anti-nucleosynthesis.
Because weak interactions are inefficient in the plasma and weak decay timescales for the charged pions are long compared to the dynamical expansion timescale of the fluid, this value is however dependent on the details of the BSM injection process that seeds the fireball.
In reaching our conclusions as stated above in the main text, we assumed that the charged pions have a small chemical potential (i.e., ) until SMCER reactions such as became inefficient, locking in the –to– ratio.
This was equivalent to assuming that some combination of the BSM injection and possible SM processes following it were such that they established a certain very specific overall net negative electromagnetic charge (i.e., ; see Appendix B) on the (anti-)hadronic constituents of the fireball no later than the time at which the SMCER begin to be inefficient, with a compensating overall positive charge held by charged leptons in order to allow overall neutrality.
But we also showed in Appendix B that we could relax this assumption on the small pion chemical potential somewhat, which allows to vary in some range around 0.5 with only changes to our results.
However, if is too small around the time the SMCER reactions become inefficient, a larger anti-proton abundance results for the same anti-helium flux, which can be observationally challenging.
While we argue that SM processes impose , we found that is likely required for phenomenological viability (see Fig. 9).
Realizing this thus likely requires that the BSM injection process needs to be able to primarily inject both anti-quarks and anti-leptons of compensating charge.
Finally, we showed that it would be plausible (at least in principle) for the required injection luminosity to be obtained via the collisions of supermassive, composite dark states (possibly a subcomponent of the dark matter), provided that these otherwise individually cosmologically stable states can become destabilized in the collision in such a way that activates a decay channel that converts a non-negligible fraction of their mass energy to SM anti-quarks (and possibly some fraction of positively charged leptons), with this taking place both rapidly and in a localized region of space.
We also showed that, assuming certain benchmark bulk behavior in the post-collisional evolution can be realized, this same scenario could seed fireballs with the correct bulk physical values to realize anti-nucleosynthesis scenario we advance in a phenomenologically viable fashion.
We did not supply a full particle physics model that could computably realize the required behavior, but we offered two potential model paradigms for these dark objects that we believe may be interesting to investigate further in future work.
The first paradigm is based on a model with a single species of fermions held together by a combination of vacuum pressure differences inside and outside a large agglomeration of such objects, and an intermediate-range Yukawa force, but supported against collapse by their nonrelativistic degeneracy pressure.
Collisions of pairs of such objects may in certain regimes provoke a catastrophic collapse of the resulting combined state if it exceeds a critical threshold of constituents.
We can imagine operators that allow some number of these fermions to then mutually annihilate to a collection of SM anti-quarks; modified versions of this behavior could likely realize the opposite-sign lepton injection too.
The dramatic density difference in the pre- and post-collision states of these objects may allow them to be cosmologically stable before collision, but for the annihilation rates to spike high enough during the collapse to permit the requisite injection in a sufficiently rapid timescale.
The second paradigm is a based on a collision of dark analogs of degenerate white dwarfs, that may undergo rapid collapse after collision, provided the collision produces a merger product above a certain threshold number of constituents.
During this collapse, constituents of the dark dwarfs may be converted to other dark species that are unstable to prompt decay to SM anti-quarks.
In this case, additional model building or a modification of the picture we presented are likely required to achieve the opposite-charge lepton injection.
Neither of these scenarios is however amenable to full analytical control, and important aspects of the evolution of each are thus still unclear.
To be completely clear: it is thus still an open question whether, and how, the requisite BSM injection of SM anti-quarks (and opposite-sign charged leptons) can be achieved via dark-state collisions in order to seed the fireballs whose post-injection SM behavior we have studied.
The most obvious and pressing follow-up is therefore a fuller investigation of these (or other) model paradigms, so as to place this BSM aspect of our scenario on a rigorous footing.
We anticipate returning to this point in future work.
On the observational side, further AMS-02 data will continue to be taken until 2030 [7]; combined with the existing data, this may yield further clarity about the status of the AMS-02 candidate anti-helium events.
Additionally, the balloon-borne GAPS experiment [126] is approved to fly in late-2024 [127], and is anticipated to provide data on lower-energy anti-deuterium and anti-helium fluxes that will be complementary to those from AMS-02.
These data will continue to inform model-building.
The fireball anti-nucleosynthesis scenario we have advanced in this paper provides an interesting and novel alternative SM pathway to anti-helium formation, provided that the BSM seeding issue can be appropriately addressed.
We therefore view this work a partial step toward understanding the origin of the tentatively identified AMS-02 anti-helium events.
Acknowledgements.
We thank Anne-Katherine Burns, Junwu Huang, Tim Linden, Vivian Poulin, and Stefano Profumo for useful discussions and correspondence.
E.H.T. thanks Xuheng Luo and Ngan H. Nguyen for useful discussions on another project.
We also thank an anonymous referee for constructive comments that improved the clarity of the presentation of aspects of this work.
This work was supported by the U.S. Department of Energy (DOE), Office of Science, National Quantum Information Science Research Centers, Superconducting Quantum Materials and Systems Center (SQMS) under Contract No. DE-AC02-07CH11359.
E.H.T. acknowledges support by NSF Grant PHY-2310429 and Gordon and Betty Moore Foundation Grant No. GBMF7946.
D.E.K. and S.R. are supported in part by the U.S. National Science Foundation (NSF) under Grant No. PHY-1818899.
S.R. is also supported by the Simons Investigator Grant No. 827042, and by the DOE under a QuantISED grant for MAGIS.
D.E.K. is also supported by the Simons Investigator Grant No. 144924.
Research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science, and Economic Development, and by the Province of Ontario through the Ministry of Colleges and Universities.
M.A.F. gratefully acknowledges the hospitality of the Simons Center for Geometry and Physics at Stony Brook University.
M.A.F. also thanks the Aspen Center for Physics, where parts of this work were undertaken, supported by National Science Foundation (NSF) Grant Nos. PHY-1607611 and PHY-2210452.
Appendix A Derivation of fireball scaling laws
As described in Sec. III.2, following a transient and localized injection of energy and anti-baryon number, a thermalized fireball with anti-baryon–to–entropy ratio is formed.
In our parameter space, the fireball immediately accelerates to semirelativistic bulk velocities corresponding to Lorentz factors , at which point its radius and temperature are and respectively, turning into a dense shell with a central underdensity in the process.
At the point where the radius of the fireball is , this shell structure is just beginning to form, and the shell thickness can be taken to be .
The subsequent evolution of the shell is dictated by the following relativistic fluid equations which, respectively, encode conservation of anti-baryon number, energy, and momentum in the comoving fluid frame:
(64)
(65)
(66)
where , , (), , , and are, respectively, the radial position, time, bulk velocity (and its associated Lorentz factor), anti-baryon number density, energy density, and pressure of a fluid element in the fireball.
Primed quantities are defined in the comoving rest frame of the fireball fluid; unprimed quantities are defined in the fireball CoM frame.
In this Appendix we show how, in the ultra-relativistic bulk velocity limit , the fireball expansion follows the simple scaling laws used Sec. III.2.
See also Refs. [34, 35] for alternative explicit derivations, and Refs. [43, 42] for related studies.
In the fireball CoM frame, the energy and anti-baryon number of a thin radial slice of the shell with radius and thickness are given by
(67)
(68)
The time derivatives of and can be written as
(69)
(70)
where we have used along the radial trajectory of a shell, assumed that the shell evolves adiabatically with a constant equation of state , and made use of the fluid equations, Eqs. (64)–(66).
For small , we can perform the expansion .
Hence, the content of the round brackets in the above equations is and thus negligible for sufficiently small .
Moreover, for sufficiently large , the term with the square brackets is also negligible.
Therefore, in these thin-shell and ultrarelativistic limits we can assume that and are approximately constant.
We now focus on a differential radial layer whose initial radius is and neglect for the moment the time evolution of its thickness .
The change in the thickness of the shell will be important but (as we will show) only at later times.
In the interest of explaining the AMS-02 relativistic anti-helium candidate events, we are particularly interested in scenarios with , where each shell is initially radiation dominated (RD) and thus obeys .
That latter scaling, taken together with and the conservation of and imply that and during RD.
This ‘accelerating phase’ lasts until , whereupon the shell becomes matter dominated (MD).
Following that, , which implies that and during MD.
We refer to this as the ‘coasting phase’.
To summarize, as long as the shell thickness is approximately constant, the following scaling laws hold:
(71)
(72)
Due to the differences in the velocities of neighboring radial layers, the constant-thickness approximation employed above breaks down at some point.
Assuming the velocity gradient is monotonic and approximately uniform over the full thickness of the shell, , and that , we find .
This implies that the thickness of a differential shell will behave as
(73)
where we have used ; the second term in the -bracket encodes the shell-thickness spreading effect.
Since during the accelerating phase, the spreading effect in that phase decreases with increasing time/radius, leading to .
On the other hand, in the coasting phase , the spreading term grows linearly with .
The spreading term becomes dominant when , implying the existence of a second expansion phase during MD where the shell thickness grows as :
(74)
This ‘spreading phase’ continues until the plasma becomes optically thin to the photons, at which point the assumption of strong coupling of the photon and anti-baryon fluids breaks down. As described in Sec. III.2, this occurs when the radius of the fireball is , found in Eq. (20).
Appendix B Freeze-out ratio of anti-nucleons
The discussion in this Appendix supplements that in Sec. III.4.3 and concerns the time period after the fireball has thermalized via strong and EM interactions () and before the anti-neutron–to–anti-proton ratio has frozen out () at SMCER decoupling.
While strong and electromagnetic interactions are efficient given the fireball expansion timescales, weak interactions typically operate at timescales that are significantly longer.393939For the benchmark fireball with and , we initially have . This means that the fireball we are considering is only partially thermalized.
Some assumptions must be made about the chemical potentials of certain particle species, such as leptons and other particles that are in (relative) chemical equilibrium with them, in order to fully characterize the thermal state of the fireball.
In the main analysis, we made such an assumption by neglecting the chemical potential (i.e., setting ) in calculating the chemical-equilibrium anti-neutron–to–anti-proton ratio at Eq. (25).
In this appendix, we depart from that assumption and instead use charge asymmetry in the lepton (or hadron) sector to parameterize the non-thermalized part of the fireball.
We then describe the requirements to achieve the assumed negligible in terms of the charge asymmetry parameter and comment on the impact of relaxing this assumption on the outputs of fireball anti-nucleosynthesis.
Let the comoving number densities of , , and charge residing in the lepton sector, be denoted respectively as , , and . Assuming the fireball is locally charge neutral, we must then have404040Charged kaons are stable over the timescales of our interest and so should also be present, albeit with much lower abundances compared to since they are significantly heavier.
We expect their influences to be qualitatively similar to that , but that they will not change our conclusions due to their relatively low abundances compared to the pions.
(75)
The inefficiency of weak interactions in the fireball expansion timescales implies that the comoving net charge density in the hadronic sector and in the lepton sector are separately approximately conserved.
Therefore, we can treat414141Note that can in principle be negative, however that implies a negative , which leads to exponentially suppressed anti-proton abundance in most cases.
We therefore do not consider this case.
(76)
as a fourth parameter, in addition to the three parameters , characterizing the evolution of the incompletely thermalized fireball.
The anti-neutron–to–anti-proton ratio is a key input for our anti-nucleosynthesis analysis. Its freeze out value is well approximated by its chemical equilibrium value at the moment of SMCER decoupling.
The impact of the charge asymmetry enters the through dependence of and on it, as dictated by equilibrium thermodynamics.
To quantify this dependence, we start by assuming
(77)
(78)
(79)
(80)
where the first and second lines are justified by efficient (see discussion below) strong and EM interactions such as , , and being in equilibrium; the third and fourth lines are justified as long as the charge densities carried by and are .
Hence, we can approximate the chemical-equilibrium abundances and as [128, 129]
(81)
(82)
The charge neutrality relation Eq. (75) thus reduces to
(83)
For simplicity, we will set .
Eq. (83) can be solved to give .
For and , we found that varies monotonically from to as is varied from to , providing an a posteriori justification toward the lower range of for the assumption at Eq. (80) that was used to justify the approximation at Eq. (82).
Figure 8: Chemical-equilibrium abundances for different charged-pion chemical-potential assumptions.
The benchmark parameters , , are assumed in these plots.
Panel (a): The chemical-equilibrium abundance of anti-protons as a function of the comoving temperature for different values of [defined at Eq. (76)]; the for is also shown for comparison.
Panels (b)–(d): The chemical-equilibrium abundances of anti-protons , symmetric population of charged pions , and anti-symmetric population of charged pions as a function of the comoving temperature for the selected parameter values , , and , respectively.
For comparison, we also show and for (as assumed in Sec. III.4.3).
The gray-shaded regions are where we expect and to be given by their freeze-out values, instead of the chemical-equilibrium values, and , that are shown in (or inferable from) these plots.
Knowing allows us to compute the chemical-equilibrium abundances of anti-protons and charged pions; see Fig. 8.
At high temperatures, chemical potentials are insignificant and consequently particles are thermally populated in a democratic way, yielding and (to keep charge neutrality) .
On the other hand, at low temperatures, thermodynamics favors storing charge asymmetry in rather than . Very approximately, this can be understood heuristically in the following way.
We have assumed a net anti-baryon number for the plasma.
Because baryon number is conserved in the plasma, a fixed number of anti-nucleons (forming the net anti-baryon asymmetry) therefore automatically exist in the plasma at any temperature, either as or ; specifically, the mass of these particles need not be extracted from the thermal bath in order for them to be in existence.
We have also assumed a non-zero negative charge asymmetry on the hadronic sector, and charge is also conserved in the plasma.
An anti-proton kills two birds with one proverbial stone: it can carry both baryon number and charge.
By contrast, pions do not carry baryon number and therefore are not automatically in existence in the plasma in any abundance: they must be created thermally; i.e., their rest-mass energy must be extracted from the thermal bath.
Effectively, it therefore costs much more energy to create a massive thermally from the bath in order to store negative charge on it, than it does to simply store the charge on .
As a result, unless , most of the charge asymmetry typically ends up in anti-protons, such that well before the SMCER interactions decouple.
Figure 9: The abundances of nuclear species (solid and dashed colored lines, as identified in the legend) as a function of the comoving fireball temperature , computed by numerically solving the Boltzmann equations for a simplified nuclear-reaction network detailed in Appendix C, for a fireball with , , ().
Note that this plot is similar to Fig. 2; however, instead of neglecting the charged pion chemical potential, here we set such that , , and , all evaluated at .
Note that here is an order of magnitude larger than in Fig. 2.
Starting at a sufficiently high (), the actual comoving number densities of , , and closely track their chemical-equilibrium values (obtained from equilibrium thermodynamics as we have computed above) until they eventually freeze out once become too Boltzmann suppressed to keep the SMCER interactions efficient.
We found that the decoupling temperature424242We computed using the criterion , where we took which generalizes Eq. (1).
Furthermore, we compare the rate of the pair production process and find that .
Therefore, the strong and electromagnetic interactions freeze out at similar temperatures. varies slightly in the range when is varied in the range .
Importantly, for , we find that already for ; see Fig. 8.
This means that we can indeed rely on our equilibrium thermodynamics arguments above to establish that the charge asymmetry is maintained on the anti-protons in this regime.
In general, we find that having results in at SMCER decoupling; see Fig. 8.
Because there are only quantitative changes to the anti-nucleosynthetic abundance results in the main text for more general values for than the value of that we used in Sec. III.4.3 (under the more restrictive assumption of ), this implies that the same qualitative conclusions that we reached in the main text will continue to hold for any values of .
This is turn is significant because, in order to achieve exactly at (as we assumed in Sec. III.4.3) we would need to choose such that it almost exactly matches the freeze-out anti-proton abundance in the absence of : i.e., .
Were that exact value a crucial input to the anti-nucleosynthesis analysis, this would constitute a fine tuning because is supposed to be set by pre-thermalization physics which has no connection to the post-thermalization physics (although see the discussion in the last paragraph of this Appendix below).
Our analysis in this Appendix thus establishes that our main-text analysis is not finely tuned.
Next, we discuss how the fireball outputs are affected if .
In this case, we generically have at SMCER decoupling.
However, even then, a anti-helium isotope ratio can still be obtained without overproducing anti-protons; i.e., when the fireball becomes optically thin, in order to avoid violating the uncertainties on the AMS-02 anti-proton flux [79]; see Fig. 6.
To demonstrate this, we show in Fig. 9 the evolution of nuclear abundances in this case with and other parameters chosen to explain the AMS-02 candidate anti-helium events.
The freeze-out abundances shown in Fig. 9 are marginally consistent with the anti-proton flux observed at AMS-02.
To explain the anti-helium isotope ratio with would however lead to an overproduction of anti-protons. Therefore, compatibility with the AMS-02 data requires .
Note, however, that the marginal set parameters we picked for Fig. 9 to complement could be interesting from the perspective of the anti-proton excess [72].
We have so far treated the charge asymmetry as a free parameter.
In actuality, it depends both on the properties of the initial particle injection, which are model-dependent but could naturally give , and post-injection SM processes.
In case the initial injection does not produce a considerable charge asymmetry, some degree of can arise spontaneously from non-equilibrium Standard Model processes occurring before or after thermalization.
Electroweak process such as and its variants can transfer charges between the hadron and lepton sectors; however, they have a preferred charge-transfer direction only at sufficiently low temperatures, , at which point electroweak interaction rates are extremely suppressed (given the dynamical timescales involved in the fireball expansion).
A higher contribution to arises from the decay of (a small fraction of) charged pions.
As shown in Fig. 8(b), even if we start with being virtually zero, a charge asymmetry of is automatically present in at high temperatures (SMCER reactions are efficient, so if , we have ).
Within the initial fireball expansion timescale (during which the bulk plasma motion is non-relativistic), a fraction (for ) of the charged pions decay to charged muons [ with [80]], thereby generating in the range .
Unfortunately, at least for the benchmark parameters we used in the main text (or the benchmark used in Fig. 9), this spontaneously generated is too small to be phenomenologically viable in light of the AMS anti-proton results (see discussion above).
Appendix C Simplified nuclear reaction network
Anti-neutrons and anti-protons in the fireball plasma are initially held in chemical equilibrium by charged-pion–mediated anti-neutron–anti-proton interconversions (i.e., SMCER) at .
As found in Sec. III.4.3, such interconversion process decouples relatively early, at , when there are virtually no nuclear bound states present.
On the other hand, photodissociation stalls the synthesis of anti-deuterium (and the whole nuclear chain) until the temperature is considerably below the anti-deuterium binding energy , at , by which point the anti-neutron–to–anti-proton ratio has completely frozen out.
Hence, we can treat the anti-neutron–anti-proton decoupling separately from the whole nuclear reaction chain.
In our numerical procedure, we thus first solve for the freeze-out ratio of ignoring nuclear reactions, finding per Eq. (25) [assuming ; similar results could be derived for more general as discussed in Appendix B], and use this as an input when solving the Boltzmann equations governing the subsequent nuclear reactions.
We now describe our simplified nuclear reaction network governing the evolution of relative abundances of elements, .
We consider only elements with , namely , , , , , and .
Heavier species with such as and isotopes are virtually absent and have negligible impacts on the elements that we consider.
Conservation of baryon number ensures , with being the atomic mass number of the nuclear species .
This can be used to solve for the abundance of one element, which we take to be , given the other abundances
(84)
We set as the initial conditions
(85)
per Eq. (25) at an initial temperature434343This value of the initial temperature is arbitrarily chosen. As long as it lies in the range , the exact value of the temperature at which the initial conditions are set is not important. , and evolve the abundances of elements other than with the following simplified nuclear reaction network [130]:
(88)
(91)
(92)
(95)
(96)
where we have defined the dimensionless quantity
(97)
which quantifies the efficiency of nuclear reactions on the timescale of fireball expansion.
Apart from anti-deuterium photodissociation, we include in the above set of equations only thresholdless (forward) reactions.
The reverse processes to these reactions are endothermic, and are suppressed by factors of the form (where is the threshold energy of the reverse reaction; i.e., the energy liberated in the forward reaction) in their thermally averaged cross section as well as the smallness (verified a posteriori) of the factors relative to those of the forward reactions.
It is nevertheless important to take into account anti-deuterium photodissociation due to its uniquely low threshold energy and the potentially large photon abundance.
The abundance of photons energetic enough to photodissociate the nuclear species is
(98)
where is the threshold energy of that process.
The photodissociation cross-sections as functions of energy typically rise sharply at and fall of smoothly at higher .
We approximate the thermal averaged photodissociation cross-sections with their near-threshold values.
We consider only the photodissociation processes , with and [131]. The other ones have significantly higher threshold energies and are consequently much less efficient due to the exponentially suppressed photon abundance (as considered at the temperature when anti-nucleosynthesis is efficient).
Appendix D Nuclear cross sections
Figure 10: The assumed values of thermal-average cross sections for (anti-)nuclear reactions as a function of temperature .
We assume that the nuclear cross-sections for purely anti-matter processes are equal to those of the analogous reactions involving matter.
For completeness, we list here the thermal averaged nuclear cross sections from [46, 48, 132, 133] and express them in terms of the temperature ratio :
•
()
•
() and (). For each,
•
()
•
()
•
()
•
()
•
()
We show these cross sections in Fig. 10, rescaled as appropriate to fit them all on the same plot.
Most important to our analysis are cross-sections at temperatures – where anti-nucleosynthesis mainly occurs.
We neglect the mild temperature dependences of the cross-sections for processes that do not suffer from Coulomb-barrier suppression.
Appendix E Changes in the number of degrees of freedom
Here we account for the effect of changes in the number of degrees of freedom on the fireball expansion dynamics (see Sec. III.2), neglecting the small differences between the for entropy and energy density [134]; i.e., we assume .
We first write the energy and entropy conservation equations, respectively, for each of the shells (assumed relativistic, with ) in more general forms:
(99)
(100)
In cases of our interest, the fireball can be approximated as consisting of only radiation (photons and relativistic massive species) and (anti-)matter (dominated by anti-baryons), allowing us to write
(101)
where and are the contributions from radiation and (anti-)matter, respectively, and we used .
Eliminating from the energy and entropy conservation equations, and using , we find
(102)
Substituting , Eq. (102), and Eq. (73) into Eq. (100), we obtain an equation that determines the evolution of with :
(103)
which shows how the temperature-radius relation is modified in the presence of changes in compared to the for constant assumed in the main text.
During RD (), we thus have
(105)
(106)
while during MD () we have
(107)
(108)
Notice that the impact of a change in is weak if it occurs during matter domination.
In our scenario, the effective number of relativistic degrees of freedom of the fireball is given by
(109)
reflecting, respectively, the contributions from the sets of particles , , and in each temperature range.
That is, we have two possible “jumps” in , happening at and (the transitions of course are smooth, but abrupt).
Since we require and , only the jump at has a chance to spoil the accuracy of the constant- approximation adopted in the main analysis.
This jump at amounts to an error in the analytical estimate for by a factor of and that for by a factor of .
These are within the level precision we are aiming for, given that our analysis does not fully capture the effects of the inhomogeneity and initial evolution of the fireball, which would also introduce uncertainties in our estimates for .
Note that the Lorentz factor remains to be related to the temperature as given in Eq. (102).
Although appears at various places in the Boltzmann Eqs. (88)–(96) [it plays the role of time, appears in the cross sections, and controls ], its relation to the shell radius , namely [which is affected by changes in ], only enters via , as defined in Eq. (97), through the comoving expansion timescale .
The latter is given by Eq. (10) in the absence of changes in and would change by a mild factor when the effects of changes in are included.
This amounts to a slight shift in the value of in relation to that yields a given .
Since controls the output nuclear abundances, the lines of constant anti-helium isotope ratio in, e.g., Fig. 4, would be displaced, though only slightly.
Appendix F Prompt (‘burst’) or slow (‘wind’) injection
The duration of the anti-quark injection can be shorter (burst regime) or longer (wind regime444444See Ref. [135] for a discussion of thermalization and hydrodynamics in the wind regime in an analogous setup (but completely different context).) than the spatial size of the injection .
Following the injection, the anti-quarks will have the following number density profile if they do not interact with one another and simply free stream away
(110)
where is the radius from the center of a spherical injection region of radius , is the total energy of the injected SM particles, and the factor
(111)
accounts for the burst vs. wind branching of cases: in the wind case the density of the injected anti-quarks within is diluted by a factor of .
Regardless of the relative sizes of and , the probability that an anti-quark undergoes a process of any sort before arriving at a radius , , decreases with for , where we have assumed454545When geometrical effects are taken into account, one would find that goes down with as the motion of the injection particles become increasingly radial [136], which means the probability would actually reduce even faster with , further strengthening the argument we are making.
that the relative velocity is relativistic and the cross section is independent of .
It follows that in order for these particles to thermalize, the thermalization rate needs to be efficient inside the injection region; i.e., at .
Once the fireball can achieve thermalization and has attained a semirelativistic bulk radial velocity , its subsequent evolution is completely described by the fireball temperature and radius when it first thermalizes (as well as ), and these are given by
(112)
(113)
Note that the anti-baryon–to–entropy ratio does not depend on the relative size between and .
As explained in Appendix A, the plasma outflow can be treated as a series independent differential radial slices which separately go through nearly the same thermal and hydrodynamical evolution.
In each of these slices, anti-nucleosynthesis proceeds as described in Sec. III.
The only difference is that and are now more generally given by Eqs. (112) and (113).
Since the abundances of nuclear species released by each radial slice are completely determined by the set of parameters which are essentially the same for all the shells, the numbers of anti-nuclei or anti-nucleons released into the ISM are still given by regardless of whether the injection is in the burst or wind regime.
We review here the procedure referred to in the main text of Sec. IV.1 that we used to obtain an estimate of the AMS-02 anti-helium sensitivity as a function of rigidity .
This procedure was specified in Appendix B of Ref. [15]; our exposition here merely reproduces the argument in that reference (adding some additional detail per Ref. [137]) and is given only to make our presentation self-contained.
For a particle species with rigidity-dependent flux at AMS-02 given by , the number of events in the rigidity range observed at AMS-02 in a time can be written as [15]
(114)
where the species-dependent “acceptance” is defined here [as in Ref. [15], where it was called ] to fold in all relevant effects such as integration time , detector effective area, trigger efficiency, etc. [i.e., it subsumes the factors , , and in Eq. (1) of Ref. [76] (or the factors of , , and the isotope-dependent integration time implicit in the rate in Eq. (1) of Ref. [138])].
Following Ref. [15], we consider first ordinary helium and look to the Supplemental Material of Ref. [76] (hereinafter, “Ref. [76, Suppl.]”) (see also Ref. [1]), which gives values for and its statistical uncertainty over narrow rigidity bins in the range .
These data can be used to extract in the following manner [15, 137].
Suppose that one of the aforementioned rigidity bins is centered at , and has width ; it then follows from Eq. (114) that we can write the number of events in that rigidity bin as
(115)
correspondingly, the Poisson statistical uncertainty on that number is
(116)
(117)
But we also have from Eq. (115) that the statistical uncertainty on the flux is
All the quantities on the RHS of the last line are known for from Ref. [76, Suppl.], allowing us to extract as
(121)
We also then have that
(122)
It remains to translate this to an anti-helium sensitivity.
We follow the argument given in Ref. [15] to do this.
In Ref. [74], a projection for the 95% confidence upper limit on the ratio of the anti-helium flux to the measured helium flux, , assuming no observed events in an integration time in the range is given as
(123)
Because the Poisson distribution function has , it follows that the number of events in that rigidity range corresponding to the flux of anti-helium that saturates that 95% CL upper limit in an integration time would be
(124)
We then assume that the acceptances for anti-helium and helium are proportional to each other, with a rigidity-independent proportionality constant :
(125)
We thus have that
(126)
(127)
(128)
(129)
where we used at the last line that and that for steady-state helium fluxes.
Note that Eq. (129) no longer makes reference to the timescale .
Solving Eq. (129) for and substituting into Eq. (125), we have
(130)
(131)
again using at the last line that for steady-state helium fluxes.
Note that as expected because .
Although we should take from the above argument, we shift this to , the lower limit of available data from Ref. [76, Suppl.]; likewise, we shift , the nearest upper bin edge to in Ref. [76, Suppl.].
We can then (a) use Eq. (122) to construct the number of helium events in each rigidity bin given in Ref. [76, Suppl.] and sum them up to find ; (b) use Eq. (121) to construct , again using the data in Ref. [76, Suppl.]; and (c) construct using Eq. (131).
Once we have the acceptance for anti-helium, we can compute event numbers using Eq. (114).
Note that we have assumed here throughout that the same acceptance as a function of rigidity applies for all the species with [i.e., all the (anti-)helium isotopes].
The analysis in Ref. [139, §4.7] supports that assumption at the level of a few tens of percent at low rigidity, and better than 15% at high rigidity ().
We have also verified that the reconstruction procedure for the helium acceptance based on Eq. (121) and the data in Ref. [76, Suppl.] reproduces the acceptance that can be constructed from Fig. 4.29 in Ref. [139, §4.7] (after accounting for differing integration times) to within at worst %, and usually within 30% or better (with agreement generally becoming better for higher rigidity) in most of the rigidity bins from Ref. [76, Suppl.].
This is acceptably accurate for our purposes.
We also tried to apply this acceptance reconstruction technique to the separated and isotope data presented for a much smaller range of low rigidities () in the Supplemental Material of Ref. [138] (see also Ref. [1]).
However, a naïve application of Eq. (121) to those data for each isotope separately gives results for acceptances for that are generally a factor of smaller than those for .
However, we believe that to be a spurious result for two reasons, and therefore disregard it: (1) the raw number of events selected for analysis in Ref. [138] were events for and events for , while the (somewhat rigidity-dependent) flux ratio for to was also reported to be around 10–15%; because the event ratio is similar to the flux ratio, it seems impossible for the isotope acceptances to differ by as much as a factor of 10; and (2) the result conflicts with the effective acceptance curves464646The effective acceptance in Ref. [139], , is defined as ; it folds in both effective area and trigger-efficiency effects as defined in Ref. [138], but not integration time.
shown in Ref. [139], which that reference used to closely reproduce official (combined) AMS-02 helium flux results of Refs. [1, 76] to within .
We suspect that the acceptance reconstruction procedure of Ref. [15] that we have reviewed here is simply inaccurate as applied to the isotope-separated data in the Supplemental Material of Ref. [138] because those data are a small, , sub-component of the total helium flux data that must be unfolded by AMS-02 to obtain the individual isotope results.
Moreover, as applied to the data in the Supplemental Material of Ref. [138], it gives results with much larger variations from rigidity bin to rigidity bin than those obtained from applying it to the earlier (combined) helium data in Ref. [76, Suppl.] (with integration-time differences accounted for).
Because of these issues with these reconstructions and because they also only cover a lower rigidity range than where we need the results, we do not use as reconstructed from the isotope-separated data in Ref. [138].
A comment is also in order on the assumption of steady-state helium fluxes.
AMS-02 has reported data on the time variation of helium fluxes in Refs. [138, 140, 1].
Between the years 2011 and 2019, the fluxes at very low rigidities, have increased by a factor of [140], while those for rigidities have varied by [140]; see also Ref. [138] for alternative presentation showing changes of similar magnitude between 2011 and 2017.
The variations of the fluxes over the timescale , which covers the years 2011–3, are much smaller: they change by only a few tens of percent around their average values, even at low rigidity [140].
In the high-rigidity regime of most interest to the AMS-02 anti-helium candidate events, the assumption of steady-state helium flux in the above derivation is thus well justified.
And while there may be a mild violation of the scaling of arising from the changing flux at the low end of the rigidity range, we estimate that this effect has only -factor overall impact on our analysis.
\justify
References
[1]AMS Collaboration, M. Aguilar et al., The Alpha Magnetic
Spectrometer (AMS) on the International Space Station: Part II —
Results from the first seven years,
Phys. Rep.894 (2021) 1–116.
[2]
S. Ting, The First Five Years of the Alpha Magnetic Spectrometer on the
International Space Station: Unlocking the Secrets of the Cosmos, CERN
colloquium presentation, December 8,
2016 \faExternalLink*.
[3]
S. Ting, Latest Results from the AMS Experiment on the International
Space Station, CERN colloquium presentation, May 24,
2018 \faExternalLink*.
[4]
S. Schael, Latest Results from the AMS Experiment, conference
presentation at ‘COSPAR 2018’, July
2018 \faExternalLink*.
[5]
A. Oliva, Observations of cosmic-rays and search for anti-nuclei with
AMS-02, conference presentation at ‘Light Anti-Nuclei as a Probe for New
Physics (LAN2019)’, October 14, 2019
\faExternalLink*.
[6]
P. Zuccon, AMS-02 results & upgrade, conference presentation at
MIAPbP conference ‘Antinuclei in the Universe?’, February 28,
2022 \faExternalLink*.
[7]
S. Ting, Latest Results from the AMS Experiment on the International
Space Station, CERN colloquium presentation, June 8,
2023 \faExternalLink*.
[9]
M. Cirelli, N. Fornengo, M. Taoso and A. Vittino, Anti-helium from Dark
Matter annihilations,
JHEP08 (2014)
009 [arXiv:1401.4017].
[10]
M. Korsmeier, F. Donato and N. Fornengo, Prospects to verify a possible
dark matter hint in cosmic antiprotons with antideuterons and antihelium,
Phys. Rev. D97 (2018) 103011 [arXiv:1711.08465].
[18]
A. Bykov, K. Postnov, A. Bondar, S. Blinnikov and A. Dolgov, Antistars
as possible sources of antihelium cosmic rays,
JCAP2023 (2023) 027 [arXiv:2304.04623].
[19]
P. De La Torre Luque, M. W. Winkler and T. Linden, Cosmic-Ray
Propagation Models Elucidate the Prospects for Antinuclei Detection,
arXiv:2404.13114.
[26]
J. Letessier and J. Rafelski, Hadrons and Quark–Gluon
Plasma.
Oxford University Press, 2002,
10.1017/9781009290753.
[27]
J. Noronha-Hostler, M. Beitel, C. Greiner and I. Shovkovy, Dynamics of
Chemical Equilibrium of Hadronic Matter Close to T(c),
Phys. Rev. C81 (2010) 054909 [arXiv:0909.2908].
[33]
A. Boyarsky, M. Ovchynnikov, O. Ruchayskiy and V. Syvolap, Improved big
bang nucleosynthesis constraints on heavy neutral leptons,
Phys. Rev. D104 (2021) 023517 [arXiv:2008.00749].
[39]
R. Yaresko, M. G. Mustafa and B. Kampfer, Relativistic Expansion of
Electron-Positron-Photon Plasma Droplets and Photon Emission,
Phys. Plasmas17
(2010) 103302 [arXiv:1008.3495].
[40]
D. F. G. Fiorillo, G. G. Raffelt and E. Vitagliano, Supernova emission
of secretly interacting neutrino fluid: Theoretical foundations,
Phys. Rev. D109 (2024) 023017 [arXiv:2307.15122].
[42]
M. Diamond, D. F. G. Fiorillo, G. Marques-Tavares, I. Tamborra and
E. Vitagliano, Multimessenger Constraints on Radiatively Decaying
Axions from GW170817,
Phys. Rev.
Lett.132 (2024) 101004 [arXiv:2305.10327].
[44]
P. Meszaros, P. Laguna and M. J. Rees, Gas dynamics of relativistically
expanding gamma-ray burst sources: Kinematics, energetics, magnetic fields
and efficiency, Astrophys.
J.415 (1993) 181–190
[arXiv:astro-ph/9301007].
[48]
P. D. Serpico, S. Esposito, F. Iocco, G. Mangano, G. Miele and O. Pisanti,
Nuclear reaction network for primordial nucleosynthesis: A Detailed
analysis of rates, uncertainties and light nuclei yields,
JCAP12
(2004) 010 [arXiv:astro-ph/0408076].
[49]
F. Iocco, G. Mangano, G. Miele, O. Pisanti and P. D. Serpico, Primordial
Nucleosynthesis: from precision cosmology to fundamental physics,
Phys. Rep.472 (2009) 1–76 [arXiv:0809.0631].
[60]
W. Huang, M. Wang, F. Kondev, G. Audi and S. Naimi, The ame 2020 atomic
mass evaluation (i). evaluation of input data, and adjustment procedures*,
Chinese Physics C45 (mar 2021) 030002.
[61]
M. Wang, W. Huang, F. Kondev, G. Audi and S. Naimi, The ame 2020 atomic
mass evaluation (ii). tables, graphs and references*,
Chinese Physics C45 (mar 2021) 030003.
[65]GRAVITY Collaboration, R. Abuter, A. Amorim, M. Bauböck,
J. P. Berger, H. Bonnet, W. Brandner et al., A geometric
distance measurement to the Galactic center black hole with 0.3%
uncertainty,
Astron.
Astrophys.625 (2019) L10
[arXiv:1904.05721].
[68]
A. Vladimirov, S. Digel, G. Jóhannesson, P. Michelson, I. Moskalenko, P. Nolan
et al., Galprop webrun: An internet-based service for calculating
galactic cosmic ray propagation and associated photon emissions,
Comput.
Phys. Commun.182 (2011) 1156–1161.
[69]
I. V. Moskalenko, A. W. Strong, J. F. Ormes and M. S. Potgieter,
Secondary anti-protons and propagation of cosmic rays in the galaxy
and heliosphere, Astrophys.
J.565 (2002) 280–296
[arXiv:astro-ph/0106567].
[70]
A. Moiseev and J. Ormes, Inelastic cross section for antihelium on
nuclei: an empirical formula for use in the experiments to search for cosmic
antimatter,
Astropart.
Phys.6 (1997) 379–386.
[72]
I. Cholis, T. Linden and D. Hooper, A Robust Excess in the Cosmic-Ray
Antiproton Spectrum: Implications for Annihilating Dark Matter,
Phys. Rev. D99 (2019) 103026 [arXiv:1903.02549].
[73]
A. Reinert and M. W. Winkler, A precision search for WIMPs with charged
cosmic rays,
JCAP2018 (2018) 055.
[74]
A. Kounine, Status of the AMS Experiment,
arXiv:1009.5349.
[76]AMS Collaboration, M. Aguilar, D. Aisa, B. Alpat, A. Alvino,
G. Ambrosi, K. Andeen et al., Precision Measurement of the Helium Flux
in Primary Cosmic Rays of Rigidities 1.9 GV to 3 TV with the Alpha Magnetic
Spectrometer on the International Space Station,
Phys. Rev.
Lett.115 (2015) 211101.
[78]
T. Aramaki et al., Review of the theoretical and experimental status of
dark matter identification with cosmic-ray antideuterons,
Phys. Rep.618 (2016) 1–37 [arXiv:1505.07785].
[79]AMS Collaboration, M. Aguilar et al., Antiproton Flux,
Antiproton-to-Proton Flux Ratio, and Properties of Elementary Particle Fluxes
in Primary Cosmic Rays Measured with the Alpha Magnetic Spectrometer on the
International Space Station,
Phys. Rev.
Lett.117 (2016) 091103.
[81]
N. Lund, C. Budtz-Jørgensen, N. J. Westergaard, S. Brandt, I. L.
Rasmussen, A. Hornstrup et al., JEM-X: The X-ray monitor aboard
INTEGRAL,
Astron.
Astrophys.411 (2003) L231–L238.
[85]
K. Hirata, T. Kajita, M. Koshiba, M. Nakahata, Y. Oyama, N. Sato et al.,
Observation of a neutrino burst from the supernova SN1987A,
Phys. Rev. Lett.58 (1987) 1490–1493.
[86]
G. Pietrzyński, D. Graczyk, A. Gallenne, W. Gieren, I. B. Thompson,
B. Pilecki et al., A distance to the Large Magellanic Cloud that is
precise to one per cent,
Nature567
(2019) 200–203 [arXiv:1903.08096].
[88]
O. Adriani, G. C. Barbarino, G. A. Bazilevskaya, R. Bellotti, A. Bianco,
M. Boezio et al., Cosmic-Ray Positron Energy Spectrum Measured by
PAMELA,
Phys. Rev.
Lett.111 (2013) 081102 [arXiv:1308.0133].
[89]
M. A. DuVernois, S. W. Barwick, J. J. Beatty, A. Bhattacharyya, C. R. Bower,
C. J. Chaput et al., Cosmic-Ray Electrons and Positrons from 1 to 100
GeV: Measurements with HEAT and Their Interpretation,
Astrophys. J.559
(2001) 296.
[90]
M. Boezio, P. Carlson, T. Francke, N. Weber, M. Suffert, M. Hof et al.,
The Cosmic-Ray Electron and Positron Spectra Measured at 1 AU during
Solar Minimum Activity,
Astrophys. J.532
(2000) 653.
[91]AMS Collaboration, J. Alcaraz, B. Alpat, G. Ambrosi, H. Anderhub,
L. Ao, A. Arefiev et al., Leptons in near Earth orbit,
Phys.
Lett. B484 (2000) 10–22.
[92]Fermi-LAT Collaboration, M. Ackermann, M. Ajello, A. Allafort, W. B.
Atwood, L. Baldini, G. Barbiellini et al., Measurement of Separate
Cosmic-Ray Electron and Positron Spectra with the Fermi Large Area
Telescope,
Phys. Rev.
Lett.108 (2012) 011103 [arXiv:1109.0521].
[93]Fermi-LAT Collaboration, S. Abdollahi et al., Cosmic-ray
electron-positron spectrum from 7 GeV to 2 TeV with the Fermi Large Area
Telescope, Phys.
Rev. D95 (2017) 082007 [arXiv:1704.07195].
[94]
E. C. Stone, A. C. Cummings, F. B. McDonald, B. C. Heikkila, N. Lal and W. R.
Webber, Voyager 1 Observes Low-Energy Galactic Cosmic Rays in a Region
Depleted of Heliospheric Ions,
Science341
(2013) 150–153.
[109]
C. Gross, G. Landini, A. Strumia and D. Teresi, Dark Matter as dark
dwarfs and other macroscopic objects: multiverse relics?,
JHEP09 (2021)
033 [arXiv:2105.02840].
[112]
J. F. Acevedo, J. Bramante and A. Goodman, Accelerating composite dark
matter discovery with nuclear recoils and the Migdal effect,
Phys. Rev. D105 (2022) 023012 [arXiv:2108.10889].
[132]
V. A. Rubakov and D. S. Gorbunov, Introduction to the Theory of the
Early Universe: Hot big bang theory.
World Scientific, Singapore, 2017,
10.1142/10447.
[133]
R. Esmailzadeh, G. D. Starkman and S. Dimopoulos, Primordial
nucleosynthesis without a computer,
Astrophys. J.378
(1991) 504–518.
[138]AMS Collaboration, M. Aguilar, L. Ali Cavasonza, G. Ambrosi,
L. Arruda, N. Attig, A. Bachlechner et al., Properties of cosmic helium
isotopes measured by the alpha magnetic spectrometer,
Phys. Rev.
Lett.123 (2019) 181102.
[139]
R. Sonnabend, Search for antihelium nuclei in cosmic rays with the
AMS-02 experiment on the International Space Station.
PhD thesis, RWTH Aachen U., 2023.
10.18154/RWTH-2023-10864.
[140]AMS Collaboration, M. Aguilar, L. A. Cavasonza, G. Ambrosi,
L. Arruda, N. Attig, F. Barao et al., Properties of Daily Helium
Fluxes,
Phys. Rev.
Lett.128 (2022) 231102.
