Mildly boosted dark matter annihilation and reconciling indirect galactic signals
Steven J. Clark
sclark@hood.eduHood College, Frederick, MD 21701, USA
Abstract
The galactic center excess is a possible non-gravitational observation of dark matter; however, the canonical dark matter model (thermal freeze-out) is in conflict with other gamma-ray observations, in particular those made of the Milky Way’s satellite dwarf galaxies. Here we consider the effects of a two-component dark matter model which results in minimally boosted particles that must remain bound to their host galaxy in order to produce an observational signal. This leads to a signal that is heavily dependent on galactic scale and can help reconcile the differences in the galactic center and dwarf galaxy measurements under the dark matter paradigm.
I Introduction
The galactic center excess (GCE) [1, 2, 3] is a flux of gamma rays originating from the center of the Milky Way galaxy that is higher than predictions from astrophysical processes. One possible interpretation is that it is due to dark matter (DM) interactions with Standard Model particles (SM); if correct, this interpretation would be the first non-gravitational detection of DM [4, 5, 6]. However, the DM parameter space that best correlates with the DM interpretation is also in conflict with other measurements; in particular, it is in conflict with similar gamma-ray measurements of the Milky Way’s satellite spherical dwarf galaxies (dSph) [7, 8, 9, 10, 11, 12].
If the GCE does originate from DM interactions, then reconciling these two observations can shed light on DM properties. Multiple models have attempted to address this difference in light of the DM proposal. Approaches frequently revolve around modifying the SM spectra through different particle productions [13, 14] and adjusting the astrophysical interaction rates (commonly termed the J-factor). Some approaches for altering the J-factor include interactions with various velocity dependencies [15, 16, 17, 18, 19, 20] or signals originating from secondary highly boosted DM [21, 22, 23, 24]. For approaches that modify the J-factor, the central concept is that there is an inherent difference in the two environments (small and large galaxies) which leads to galactic dependencies not captured in the canonical value [24]. In order to reconcile the GCE and dSph signals, this would require either an enhancement in larger galaxies or a suppression in smaller ones.
Expanded dark sectors offer a possible approach at addressing these signals by introducing dynamics that play a crucial role in DM distributions. Note that galactic DM is non-relativistic; this implies that processes that impart small increases to a particle’s kinetic energy, compared with its rest mass, can lead to the particle achieving escape velocity, , from the host galaxy. Because is dependent on galactic size, the requisite energy is lower in smaller galaxies and thus easier to escape. If this “boosted” DM particle is SM active, larger fractions of escaping boosted DM correspond with lower observational galactic signals. This leads to an overall suppression in the observed rates that is more pronounced in smaller galaxies, thus providing a mechanism to reconcile the GCE and dSph results.
In this work, we investigate galactic signal rates from multi-component DM models where the SM active components are created with a mild boost from a dominant SM inert portion. These types of models lead to a strong galactic dependence in observational rates with a rapid transition where galaxies above a critical scale experience minimal alterations to canonical rates while those below can experience strong to total suppression.
II Model and Bounded Fraction
For illustrative purposes, we consider a basic two component dark matter toy model similar to [24] consisting of and with mass relationships and . annihilates to while annihilates to standard model particles.
(1)
(2)
where the “b” superscript indicates that the s are produced with extra kinetic energy. We assume annihilation is weak allowing for to serve as the dark matter candidate while the annihilation rate is comparatively much stronger. From the perspective of a galaxy, s comprise the majority of the dark matter, and s are produced through annihilation. These will either be produced with sufficient velocity to overcome the gravitational potential and escape the galaxy, or they will remain bound. We will assume that if they achieve escape velocity, the annihilation coupling is small enough that they escape without further interaction (for annihilation occurring while escaping the galaxy, see Ref. [24]). If s do not achieve escape velocity, they persist in the host galaxy until they annihilate with another . In this setup, the galactic population fluctuates until it reaches a steady state solution, balancing between injections from the first interaction (Eq. (1) adjusted by the fraction that escape the galaxy) and depletion from the second annihilation (Eq. (2)) producing a SM signal similar to canonical DM annihilation.111Because will only experience a mild boost, their annihilation spectra will be identical to similar final products as canonical DM annihilation. For this work, we assume that cross-sections are sufficient for all galaxies to reach this equilibrium state. The observable signal is directly proportional to the rate of annihilation, and when at equilibrium, it is also proportional to the injection rate.
Two quantities are required to determine the rate of injection into the host galaxy: the base rate of annihilation producing s and the fraction from any particular annihilation that does not achieve escape velocity. We first discuss the fraction that does not reach . Fig. (1) shows annihilation in the center of mass frame (COM). corresponds to the velocity of in the frame while is the velocity of the center of mass with respect to the galactic frame. is the angle between and .
Using conservation of energy, it is easy to show that
(3)
where in the non-relativistic limit. In the COM frame, both s have the same velocity; however, in the galactic reference frame, a difference develops depending on their orientation with . In the galactic frame, denoted by subscript “”, the daughter particle velocity is
(4)
where we have assumed that all velocities are non-relativistic. It should also be noted that in the non-relativistic limit where is the relative velocity between the two parent particles. This is true for both the COM and galactic reference frames.
For to remain bound to the galaxy, the total energy must be negative; this leads to the relationship where is the gravitational potential of the host galaxy. Combining this relationship with Eqs. (3) and (4) as well as , we arrive at the condition
(5)
where is the minimum angle that must make with in order for the to remain bound to the galaxy. If the right hand side of Eq. (5) is greater than 1, then the condition is always satisfied, all products are bound, and ; if it is less than , and all products escape. Note that Eq. (5) is valid only for . For , the binding condition is ; otherwise, all products escape.
For simplicity, we assume that annihilation is isotropic in the COM frame. In this isotropic example, the bound on translates to the fraction of bounded annihilation products () through the fractional area of a 2-sphere between .
(6)
III Initial Annihilation Rate
To determine the initial annihilation rate from , we follow the approach from Ref. [18] with the addition of as discussed above to restrict the effective rate to include just the fraction of daughter particles which remain bound to the galaxy. For this work, the relevant quantity is . In canonical velocity independent annihilation rates, is analogous to the square of the DM density () and evaluated from the DM velocity distribution, capturing the relative rate of annihilation occurring at a particular position in the galaxy.
(7)
Subscripts correspond to the two individual parent particles involved in the annihilation. (Note that this differs from the preceding section where subscripts indicated different particle species.) is the phase-space distribution for in the galaxy assumed here to be isotropic. is the relative velocity between the two parent particles with being a model parameter. is the fraction of daughters bound to the host galaxy in an annihilation. Also note the relationship .
For convenience in analyzing multiple galaxies later in this work, we have introduced the scaled distances, densities, and velocities in Eq. (7)
(8)
where and are the galactic scale and density parameters, and is the gravitational constant. Furthermore, the scaled phase-space distribution and gravitational potential are
(9)
(10)
where .
As stated before, Eq. (7) encapsulates the relative rate of DM annihilation, and captures the velocity dependence. For this work, we consider only the velocity independent interaction and leave for future studies. For and , Eq. (7) reduces to as expected.
J-factors (a measure of the expected flux) adjust the annihilation rate by accounting for the distance to the object and the observation window through a line of sight (l.o.s.) and region of interest (ROI) integration over .
(11)
III.1 Potentials and Phase-Space Distributions
From a known DM distribution, the potential can be calculated using Newtonian gravity. If the distribution is spherically symmetric, then the scaled potential is [25, 26, 18]
(12)
In addition, if the DM velocity distribution is assumed to be isotropic and the halo is in equilibrium, the velocity distribution can also be determined from the density and the potential in terms of , the energy per unit mass through [18]
(13)
(14)
where we assume and go to zero at . In this manner, we are able to define a fully self-consistent velocity distribution for the DM halo. The density can be found from the velocity distribution through
(15)
(16)
where . For this work, we assume the NFW profile [27] for all galaxies. This leads to the gravitation potential where is the natural logarithm. From and , we solve for numerically.
IV Adjusted Annihilation Rates
Due to the boost received during the first annihilation, only a fraction of remains bound to the galaxy. Fig. (2) shows for various kick velocities. The ratio between and is effectively the number of daughters that remain bound to the galaxy and the number that are produced. This ratio captures the fraction of the first annihilation products that are bound to the galaxy and which can participate in future interactions.
As would be expected from Eqs. (3 - 5), for , all products remain in the host galaxy. This is due to the minimal changes to the energy distribution in the system. As increases, escape occurs at large galactic radii due to the lower gravitational potential where even a small velocity change can provide the required energy. With increasing , the outer edges continue to become suppressed with the development of a critical radius above which there is total suppression.
This total suppression is due to the kick velocity providing the required energy for all valid velocity combinations to escape at the particular radius. From Eq. (4), , and , it is easy to find this maximum kick velocity where all annihilation products are unbound from the galaxy,
(17)
For each curve shown in Fig. (2), this value is observed by the location of the sharp right cutoff.
This pattern continues until where increased suppression begins at small radii. This suppression is due to a lack of to allow for back-scattering events () that sufficiently reduce their energy. At , the increase in kinetic energy is equal to the deepest portion of the NFW potential. This requires all bound products to be back-scattered with respect to in order to reduce their speed. The amount of energy reduction in back-scattering is more pronounced for larger . At small galactic radii in the NFW distribution, a large proportion of the velocity distribution has low velocities compared with distributions at larger radii. This results in parent particles having lower average at small radii, leading to an incapability to sufficiently back-scatter to keep products bound to the galaxy even with the larger gravitational potential. Instead, these back-scattering events still have sufficient energy to escape. Rates at small radii thus experience a more dramatic suppression when compared to larger distances, and a peak in the rates is introduced for large as observed in Fig. (2).
For , all products escape and the galaxy is completely suppressed to 0. For the NFW profile, this kick velocity corresponds to ().
IV.1 Total Rates
The total scaled annihilation rate for a galaxy can be found through . In Fig. (3), we show the ratio between the galactic rate and the canonical result for the NFW distribution. Due to the assumption that and are in equilibrium, this ratio is equal to the change in expected signal from the second annihilation and the canonical result. For , there is minimal variation from the canonical result as expected. At , the rate begins to dramatically decrease such that it is negligible by . At , the rate reaches zero as no annihilation products are bound to the galaxy. As discussed earlier, for the NFW distribution, this occurs at .
For a single model, the kick velocity () will be constant for all galaxies; however, scaled velocities () are galaxy dependent. Due to the sharp transition from full canonical expectations and a complete reduction to zero, some galaxies may experience almost no variation while others could experience dramatic departures from the canonical value. The main contributor to how different galaxies behave is their physical dimensions. In Fig. (3), we have included for the Milky Way (0.451) and the dwarf galaxy Draco (4.232) assuming .222For the Milky Way (Draco), we used and . [28] Other dSph galaxies commonly used in dark matter searches have similar values to Draco (). For this choice of , we would expect minimal alterations from the canonical signal for the Milky Way while expecting a complete suppression from Draco. Similar results will occur for other dSph galaxies. This suppression could account for the discrepancy between GC and Fermi dSph measurements.
Observing a signal from the Milky Way while experiencing a departure from the canonical result for dSph necessitates , which corresponds to a mass splitting . This range is at the extreme of expected differences; a more conservative range would be to allow for a sizable Milky Way signal while experiencing large variations from the canonical result for dSphs. Interestingly, this is the same mass splitting range as the 2cDM model needed to explain the missing galaxy, core-cusp, and too-big-to-fail problems of -body simulations [29, 30, 31]; this is not completely unexpected due to the similarity in the models.
Care should be taken when interpreting these results, however, as they are a measure of the total annihilation rate of the galaxy. When converting to J-factors, they will be most accurate for distant galaxies, which the Milky Way is not. Further work is needed to understand how the secondary annihilation distribution in the galaxy will differ from canonical annihilation to identify if there are additional features when considering the angular J-factor along with the effect of galactic DM over-densities. In addition, more work is needed to better understand the equilibrium requirements necessary to reach the steady state situation assumed here and how it can differ in galactic environments.
Overall, this model introduces a mechanism that can explain the GCE as well as the lack of an observation in dSph under the DM interpretation. A characteristic feature to distinguish it from other models is a dramatic drop in excess SM products at a critical galactic scale. If measured, this scale can be used to measure the mass splitting between the two DM species.
Acknowledgements
We thank Bhaskar Dutta and Savvas Koushiappas for helpful discussion and suggestions in the preparation of this work.
References
[1]
L. Goodenough and D. Hooper, “Possible Evidence For Dark Matter Annihilation In The Inner Milky Way From The Fermi Gamma Ray Space Telescope,” arXiv:0910.2998 [hep-ph].