We explore the 511 keV emission associated with subGeV dark matter DM particles that can produce electronpositron pairs and form positronium after thermalizing. We use 16 yr of data from the Spectrometer on INTEGRAL to constrain DM properties including the full positron propagation and losses and the suppression of free electron density away from the Galactic plane. We show that the predicted longitude and latitude profiles vary significantly for different DM masses unlike previous assumptions and obtain the strongest limits on subGeV DM from the MeV to a few GeV so far excluding cross sections down to lt vgt  10SUP32SUP cmSUP3SUP sSUP1SUP for m SUB  SUB  1 MeV and lt vgt  10SUP26SUP cmSUP3SUP sSUP1SUP for m SUB  SUB  5 GeV and lifetimes up to   10SUP29SUP s for m SUB  SUB  1 MeV and   10SUP27SUP s for m SUB  SUB  5 GeV for the typical NavarroFrenkWhite DM profile. Our derived limits are robust within a factor of a few due to systematic uncertainties.
['Cosmic ray astronomy', 'Particle astrophysics', 'Milky Way dark matter halo', '324', '96', '1049', 'High Energy Physics - Phenomenology', 'Astrophysics - High Energy Astrophysical Phenomena']
New 511 keV Line Data Provide Strongest subGeV Dark Matter Constraints
{'New 511 keV line data provides strongest sub-GeV dark matter constraints': "Pedro De la Torre Luque, 1, 2, 3 Shyam Balaji, 4 and Joseph Silk 5, 6, 7  \n1 Departamento de F'ısica Te'orica, M-15, Universidad Aut'onoma de Madrid, E-28049 Madrid, Spain 2 Instituto de F'ısica Te'orica UAM-CSIC, Universidad Aut'onoma de Madrid, C/ Nicol'as Cabrera, 13-15, 28049 Madrid, Spain  \n- 3\n- The Oskar Klein Centre, Department of Physics, Stockholm University, Stockholm 106 91, Sweden\n- 4 Physics Department, King's College London, Strand, London, WC2R 2LS, United Kingdom\n- 5 Institut d'Astrophysique de Paris, UMR 7095 CNRS & Sorbonne Universit'e, 98 bis boulevard Arago, F-75014 Paris, France\n- 6 Department of Physics and Astronomy, The Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, U.S.A.  \n7 Beecroft Institute for Particle Astrophysics and Cosmology, University of Oxford, Keble Road, Oxford OX1 3RH, U.K.", 'ABSTRACT': 'We explore the 511 keV emission associated to sub-GeV dark matter (DM) particles that can produce electron-positron pairs and form positronium after thermalizing. We use ∼ 16 yr of SPI data from INTEGRAL to constrain DM properties, including the full positron propagation and losses, and the free electron density suppression away from the Galactic plane. We show that the predicted longitude and latitude profiles vary significantly for different DM masses, unlike previous assumptions, and obtain the strongest limits on sub-GeV DM (from the MeV to a few GeV) so far, excluding cross-sections down to ⟨ σv ⟩ ≲ 10 -32 cm 3 s -1 for m χ ∼ 1 MeV and ⟨ σv ⟩ ≲ 10 -26 cm 3 s -1 for m χ ∼ 5 GeV and lifetimes up to τ ≳ 10 29 s for m χ ∼ 1 MeV and τ ≳ 10 27 s for m χ ∼ 5 GeV for the typical Navarro-Frenk-White DM profile. Our derived limits are robust within a factor of a few due to systematic uncertainties.', '1. INTRODUCTION': 'The 511 keV γ -ray line from the galactic bulge indicates e + e -annihilation into γγ via positronium in the interstellar medium (ISM). This line has been measured by SPI on INTEGRAL (Siegert et al. 2016) and COSI (Kierans et al. 2020), and reviewed in Prantzos et al. (2011);  \npedro.delatorreluque@fysik.su.se shyam.balaji@kcl.ac.uk silk@iap.fr  \nKierans et al. (2019); Siegert (2023). The line has disk and bulge components, with the bulge component being narrow, bright and centered at the galactic center (GC), with a flux of ∼ 10 -3 cm -2 sec -1 (Siegert et al. 2016). This has been used to constrain possible positron sources, such as isotopes from stars (Prantzos et al. 2011; Bartels et al. 2018), low-mass X -ray binaries (Bartels et al. 2018) or neutron star mergers (Fuller et al. 2019). More recently stellar templates for the nuclear stellar bulge and a boxy bulge were  \nshown to have good matches with the observed signal Siegert et al. (2021).  \nThe bulge positrons could also come from dark matter (DM), which is abundant near the GC. This possibility has been explored, at least, since Boehm et al. (2004b). However, not all DM models work. For example, DM decay does not fit the signal shape Vincent et al. (2012); Ascasibar et al. (2006). Moreover, the positron injection energy must be around MeV scale Beacom et al. (2005); Boehm & Uwer (2006); Beacom & Yuksel (2006); Boehm (2009), which implies that the DM mass must be below some tens of MeV, unless there are intermediate steps in the annihilation. Such a low DM mass was thought to be inconsistent with cosmology, but recent studies have shown that it can be compatible with observations if there is also neutrino injection from DM annihilation in the early universe Escudero (2019); Sabti et al. (2020). We use these new findings to motivate study of light DM scenarios that may produce the 511 keV line.  \nVincent et al. (2012) showed that the DM annihilation rate needed to explain the 511 keV line emission in the bulge (from ∼ 8 years of SPI data in a 5 keV bin), was approximately  \n⟨ σv ⟩ 511 ≃ 5 × 10 -31 ( m χ 3 MeV ) 2 cm 3 sec ,  \nfor a Navarro-Frenk-White (NFW) profile DM density distribution for the best fit parameters in Vincent et al. (2012).  \nHowever, the previous studies have assumed that positrons follow the DM distribution and have ignored positron propagation in the Galaxy, that depends on their energy, E (with D ( E ) ∝ E δ , where D is the diffusion coefficient and δ > 0 means that higher energy positrons diffuse faster). This assumption was adopted since energy losses at MeV energies are dominant Strong & Moskalenko (1998), but diffusion is still important, and leads to a 511 keV profile that changes with the DM mass, as we show in Fig. 1, where we also compare with the profile expected from a  \nFigure 1. Diffuse 511 keV emission longitude profile in the latitude bin -10 . 75 · < b < 10 . 25 · from direct DM annihilation into electron-positron pairs. We show DM masses ( m χ ) ranging from 1 (cyan), 10 (green), 100 (blue) and 1000 MeV (black) respectively, and compare to the expected signal with no propagation (dashed). To facilitate the comparison with SPI data (shown in red), the emission for each mass is scaled by the factor indicated in the legend. We consider a typical NFW DM distribution.  \n<!-- image -->  \nNFWDMdensity distribution without any propagation. Moreover, the latest data from AMS-02 Aguilar et al. (2013) on secondary cosmic-rays (CRs) suggest a faster diffusion below the GeV scale (Weinrich, N. et al. 2020; Luque et al. 2021; Luque 2021; de la Torre Luque et al. 2022). Also, the interaction of charged particles with turbulence in the interstellar medium (ISM) causes energy exchange and diffusive reacceleration (Heinbach & Simon 1995; Strong & Moskalenko 1998; Cardillo 2019), which boosts low-energy particles to higher energies, making them diffuse faster. Therefore, we evaluate the 511 keV profile from the distribution of propagated positrons injected by DM, instead of assuming that it follows the DM distribution. In fact, some recent works (De la Torre Luque et al. 2023a,b,c), have shown the importance of adding propagation effects in constraining sub-GeV DM.  \nHere, we perform a detailed evaluation of the 511 keV line emission from DM and utilize the non-compatibility of these signals with the observed 511 keV line morphology, to constrain the DMannihilation cross-section and decay lifetime for the kinematically accessible channels. In doing so, we highlight three major novelties. i) We include the transport of positrons when calculating the 511 keV emission in the Galaxy. We show that the diffusion of positrons changes their spatial morphology significantly. ii) We include the fact that the free electron density drops quickly away from the Galactic Plane, which suppresses the latitude profile flux of the 511 keV emission. This is necessary for any positron source that could produce 511 keV radiation. iii) We derive strong limits on the lifetime of decaying DM. Remarkably, we find that the limits from the full SPI dataset are the most constraining ones below DM masses of ∼ 1 GeV, a mass range of interest given low sensitivity in DM direct detection experiments to light DM (Cirelli et al. 2023; De la Torre Luque et al. 2023c).', '2. EVALUATION OF THE 511 KEV SIGNAL': 'The literature assumes that the distribution of diffuse positrons from DM decay/annihilation follows the parent DM Galactic distribution. However, positrons can propagate and diffuse far from their source, changing their distribution and the profile of the observed diffuse 511 keV line in the Galaxy Jean, P. et al. (2009). Therefore, we evaluate the DM-induced 511 keV line flux considering the full propagation of positrons, including all sources of energy losses and other effects, such as reacceleration Seo & Ptuskin (1994); Osborne & Ptuskin (1987), triplet pair production Gaggero et al. (2013) and in-flight annihilation Heitler (1936). We calculate the steady-state solution for the distribution of positrons injected by DM (similar to Calore et al. (2022)) and integrate it over all possible energies.  \nWe compare our results with more than 16 years of data from INTEGRAL/SPI, which  \noffered the latitudinal and longitudinal profiles of the diffuse 511 keV line emission with unprecedented quality Siegert et al. (2019).', 'Positron propagation': 'We compute the distribution and energy spectra of the positrons produced by DM annihilation/decay in the Galaxy with a customized version De La Torre Luque (2023) of the DRAGON2 code Evoli et al. (2017, 2018), a CR propagation code that numerically solves the full diffusionadvection-loss equation for charged particles in the Galactic environment Ginzburg & Syrovatskii (1969). We use the same expression for the diffusion coefficient and propagation parameters as in De la Torre Luque et al. (2023c), which, at energies below 300 GeV, takes the form  \nD ( E ) = D 0 β η ( E 4 GeV ) δ , (1)  \nwhere δ = 0 . 49, η = -0 . 75, D 0 = 1 . 02 × 10 29 cm 2 /s and β is the ratio of the particle speed to the speed of light. This setup allows us to reproduce the current experimental CR data with very good accuracy, as broadly discussed and shown in Appendix A of De la Torre Luque et al. (2023c). These simulations include all the relevant sources of energy loss (namely, synchrotron, inverse Compton, bremsshtralung ionization and Coulomb losses), inelastic interactions, triplet pair production Gaggero et al. (2013) and inflight annihilation of positrons. We show in Appendix B a comparison of the relevant timescales of the simulation.  \nWe simulate the channels leading to the production of e ± final states through the processes: χχ → µ + µ -, π + π -, e + e -for DM annihilation and χ → µ + µ -, π + π -, e + e -for DM decay. The injection spectra dN e /dE e for each channel are calculated following Cirelli et al. (2021). We compute the signals from light DM from m χ = m i for annihilation and m χ = 2 m i for decays, up to  \nm χ = 5 GeV (where i = e, µ, π ), taking 3 masses per decade and a NFW DM distribution Navarro et al. (1996), as a baseline.  \nSimulations are performed with a spatial resolution of ∼ 150 pc, which is sufficiently fine for comparisons with SPI data and for the smooth NFW profile. We simulate electron-positron signals from DM annihilation and decay in the range of kinetic energies from 100 eV to 10 GeV, with an energy resolution of 5%. We need to consider such low energies for reasons explained below. Lowering the minimum energy or the spatial resolution does not impact our results significantly. We ensure convergence of the simulations using a variable time step (see Sect. 4.5 of Evoli et al. (2017)), with a minimum time step of 0 . 1 kyr and a maximum time step of 64 Gyr, that cover the relevant timescales involved.', '511 keV emission': 'Given the propagated and integrated flux of diffuse positrons in the Galaxy as a function of 3D position ( x, y, z ), ϕ e ( x, y, z ), we calculate the emission of 511 keV X -rays from a given direction by integrating over the line of sight, as  \ndϕ 511 γ d Ω = 2 k ps ∫ ds s 2 ϕ e ( x s,b,l , y s,b,l , z s,b,l ) 4 πs 2 , (2)  \nwhere k ps = 1 / 4 is the fraction of positronium decays corresponding to (singlet) para-positronium states) contributing to the 511 keV line signal, ϕ e = ∫ dϕ e dE dE is the energy-integrated flux of positrons, d Ω = dldb cos b is the solid angle element with l , b and s denoting the Galactic longitude, latitude and distance s (along the line of sight) from the Earth.  \nWe note that the probability of producing positronium away from the Galactic Plane depends on the positron density injected from DM and the ambient free electron density. Hence, we apply a scaling relation to the 511 keV profiles, following the vertical distribution of free electron density in the Galaxy. This is motivated by the fact that SPI observations of the  \nlatitude profile show a similar variation of the 511 keV emission as the one of the free electron density (see Appendix A and Fig. 3). This scaling leads to a more conservative but realistic evaluation of the predicted 511 keV emission. For the distribution of free electrons, we use the NE2001 model Cordes & Lazio (2003a,b), which contains the modelling of several Galactic sub-structures and is derived from a combined fit of pulsar dispersion measurements, temporal and angular broadening of radio pulses, and intensity measurements. We consider only the vertical dependence of the electron density, since the free electron density is always much higher than the diffuse positron density in the Galactic disk. This has usually been assumed true in all of the Galaxy (see, e.g. Keith & Hooper (2021); Calore et al. (2022)), however a different rate of positronium production is expected out of the Galactic plane, because of the steep reduction in free electron density when moving away from the Galactic plane. This scaling has an appreciable effect in the predicted latitude profiles, but has almost no effect on the predicted longitude profiles of the 511 keV line emission.  \nWe adopt the scaling predicted by the NE2001 model, applying the correction proposed in Gaensler et al. (2008), where the authors found that the height of the thick disk roughly doubles to ≃ 2 kpc: i) a thick disk with large Galactocentric scale height (whose scaling follows a exp[ -| z | /H 1 ] relation, where z is the Galactocentric height and the scale H 1 was found to be 1 kpc), ii) a thin, annular disk in the inner Galaxy (whose scaling decays as exp[ -| z | /H 2 ], with H 2 ≃ 140 pc), and iii) a GC component. We take this scaling to be maximal at the center of the Galaxy and exponentially decays like exp[ -| z | /H 1 ] + exp[ -| z | /H 2 ] at a height z , as shown in Fig. 3.', '3. SUB-GEV DARK MATTER CONSTRAINTS FROM SPI DATA': "In Fig. 1, we depict how predicted 511 keV line longitude profiles vary with different DM masses, specifically when directly annihilating into e + e -final states. Results for DM masses of 1, 10, 100, and 1000 MeV are compared to SPI data. A decreasing DM mass results in an increase of the signal normalization and a peakier profile around the GC. This effect is due to lower-energy injected positrons covering shorter distances before thermalization or annihilation, influenced by higher gas density toward the GC. Lower DM masses match better with SPI's longitude profile shape. However, the NFW profile fails to replicate the peakiness indicated by SPI measurements for DM masses above the electron mass. Only a highly peaked distribution like the Moore profile Moore et al. (1999) ( γ = 1 . 5) provides a satisfactory fit, as shown in Fig. 5. This emphasizes the importance of data away from central longitudes for deriving stringent DM bounds.  \nSimilar plots for other annihilation channels ( µ + µ -, π + π -, and direct e + e -) are presented in Fig. 6, including the predicted longitude and latitude profiles for various DM masses. For latitude profiles, scaling the free electron density yields a profile similar to SPI's, even at high ( ∼ GeV) DM masses. Considering positron annihilation in flight is crucial, impacting the distribution above or below the Galactic plane, particularly reducing emission in the plane with lower DM masses.  \nAdditionally, as confirmed by other studies Boehm & Fayet (2004); Vincent et al. (2012), decaying DM produces a flat profile inconsistent with observations, evident in Fig. 8.  \nWe show the bounds from SPI data and compare them to other DM constraints in Fig. 2, for DM annihilation (left panel) and decay (right panel). To derive these limits, we use a simple χ 2 fit to the SPI data, using the curve-fit Python package, and take the 2 σ uncertainty bands. We show the constraints from the SPI lat-  \nude profile as a dashed curve ('SPI Lat') and the SPI longitude profile as a solid curve ('SPI Long'). We use colour to distinguish the channel of annihilation or decay: e + e -(green), µ + µ -(red) and π + π -(blue). We restrain ourselves to computing the limits independently, since we expect that combining both datasets will not offer noticeable improvement and will be dominated by longitude data. In the left panel, we also compare with the limits from the e ± injection from DM annihilation on the CMB anisotropies, from Slatyer Slatyer (2016) and Lopez-Honorez et al. Lopez-Honorez et al. (2013), shown as dot-dashed lines ('CMB'). These limits are based on the fact that the injection of charged particles by DM during the cosmic Dark Ages must increase the residual ionization fraction, thus altering the anisotropies of the CMB. The CMB limits are weak for decaying light DM Liu et al. (2023a,b), unlike the case of annihilating DM. This is because DM clusters at redshifts z ≲ 100, which increases the DM annihilation rate and the e ± injection, while the decay rate stays the same. Since the CMB bounds are weak for decaying DM, we show the stronger bound from De la Torre Luque et al. (2023c) (labelled 'Voyager-1') as a dot-dashed curve in the right panel. This limit comes from the fact that decaying DM results in an e ± cascade that can be observed by Voyager 1 , without the solar modulation effect of the heliosphere, due to the satellite's distance from the solar system.  \nWe show the strong bounds from De la Torre Luque et al. (2023c), labelled as 'XMMNewton' in both panels, as dotted curves. These limits are derived from the fact that e ± products from an exotic DM injection of e ± may generate bremsstrahlung radiation and upscatter the low-energy Galactic photon fields via the inverse Compton process, producing a broad emission from X -ray to γ -ray energies observable by Xmm-Newton , which measured X -rays in galactocentric rings around the GC. These limits  \n<!-- image -->  \nFigure 2. Comparison of the 95% confidence bounds on annihilating (left panel) or decaying (right panel) DM derived in this work (thick dashed for SPI Lat and solid lines for SPI Long , respectively) with other existing constraints. We show the CMB bounds from Slayter Slatyer (2016) and Lopez-Honorez et al. Lopez-Honorez et al. (2013) (dot-dashed line), the bounds from XMM Newton De la Torre Luque et al. (2023c) (dotted line), the previous 511 keV line bounds from Vincent et. al Vincent et al. (2012) (red dashed line) and from Voyager 1 De la Torre Luque et al. (2023c) (dot-dashed line). In both panels, we show the π + π -(blue), µ + µ -(red) and e + e -(green) channels, respectively. Uncertainties from the propagation setup employed are shown as an orange and a blue bands, for the latitude and longitude profiles, respectively.  \n<!-- image -->  \ninclude best-fit CR propagation and diffusion parameters and represented the strongest astrophysical constraints for this mass range of DM of 1 MeV to a few GeV and surpass cosmological bounds across a wide range of masses.  \nWe show the bounds from Vincent et al. (2012), also using 511 keV data from SPI after a shorter ∼ 8 yr of operation. Our limits are stronger because their constraint used the spectrum of the emission at 511 keV, instead of the full profile. This, combined with the fact that we use the most constraining observations (from high longitude), explains the relative strength of our constraints.  \nIn fact, we see that the SPI longitude profile provides the strongest limit on annihilating DM, excluding cross-sections of ⟨ σv ⟩ ≃ 10 -32 cm 3 s -1 at masses of around an MeV, increasing to around the thermal relic values of ⟨ σv ⟩ ≃ 10 -26 cm 3 s -1 at masses of over several GeV. This SPI longitude profile only gets surpassed by Xmm-Newton X -ray bounds at masses ≳ 700 MeV. The CMB bounds only become more stringent than the SPI latitude profile  \nlimits at masses in excess of a few tens of GeV. Strikingly, limits derived from the SPI latitude profile are only marginally weaker (less than an order of magnitude across the entire DM mass range) than the SPI longitude profile limits. This was also found in Calore et al. (2022). Furthermore, this bound also surpasses the other DM limits across most of the mass range considered, except XMM-Newton at masses over a few hundreds of MeV and CMB bounds at slightly higher mass. Finally, we also compare our limits with those derived from Boehm et al. (2004a) from final state radiation emission.  \nFrom the right panel of Fig. 2, we see once again that the SPI longitude profile provides the strongest limit on decaying DM, excluding lifetimes of τ ≃ 10 29 s at masses of around an MeV decreasing to around 10 27 s at masses of over a GeV. This SPI longitude profile gets surpassed by Xmm-Newton X -ray bounds at masses ≳ 1 GeV. Limits derived from the SPI latitude profile are about an order of magnitude weaker than the SPI longitude profile limits across the entire mass range. The next strongest limit, coming  \nfrom Voyager 1 is significantly weaker (about three orders weaker than the SPI longitude profile limits). This illustrates the power of the 511 keV line in constraining sub-GeV DM. We remark that including the free electron density scaling changes our latitude bounds by a factor of 2-3, leading to a more conservative estimate.", '4. DISCUSSION AND CONCLUSION': "This study investigates sub-GeV dark matter (DM) particles annihilating or decaying into standard model (SM) particles, leading to an electron-positron cascade and positronium bound-state formation. About one-fourth of these states decay into 511 keV photons, detectable by the SPI spectrometer on the INTEGRAL satellite. SPI longitude and latitude profiles are used to establish new limits on DM properties, including mass, annihilation cross-section, and decay lifetime. We enhance previous analyses by considering two previously overlooked effects: positron diffusion and propagation in the interstellar medium, and the decrease in free electron density away from the Galactic plane. These considerations significantly impact the 511 keV photon line's shape and intensity, influencing derived DM constraints. Additionally, this study presents the first constraints on decaying DM using the SPI dataset. Assuming a NFW density distribution, our limits are the strongest across the sub-GeV DM mass range, ruling out thermally-averaged cross-sections in the range 10 -32 cm 3 s -1 ≲ ⟨ σv ⟩ ≲ 10 -26 cm 3 s -1 and decay lifetimes of 10 29 s ≲ τ ≲ 10 27 s, surpassing many cosmological and astrophysical limits. However, uncertainties persist: On the one hand, the total systematic uncertainties of the SPI data points are not totally under control, since they depend on some template models for the shape of the signal. The fact that these constraints are so sensitive to the high longitudes can be particularly affected by these systematic uncertainties, as pointed out in De la Torre Luque et al. (2023a), given that there are data points indicating no  \nflux detected at these longitudes. In that paper, we found that constraining ourselves to the inner 20 · instead (expected to be less affected by the detector's systematic uncertainties) leads to a weaker limit by a factor of 2. On the other hand, we expect that uncertainties in the CR propagation setup affect the predicted flux of positrons from DM. De la Torre Luque et al. (2023c) demonstrated that the main factors affecting the flux of positrons produced from sub-GeV DM are the halo height and the ' Alfven speed (dictating the level of reacceleration). Therefore, we have built a pessimistic and optimistic propagation setups to compute the uncertainty bands shown in Fig. 2. The pessimistic setup features the scenario with very slow diffusion, where the halo height is set to H = 3 kpc (given that the lower the halo height, the lower the positron flux (De la Torre Luque et al. 2023c), we take a very small value compared to the typical values obtained in CR analyses), keeping the ratio of the normalization of the diffusion coefficient to the halo height ( H / D 0 ) to the value that reproduces CR secondary ratios at Earth. In addition, we assume no reacceleration (i.e. V A = 0 km/s). In the optimistic case, we set V A = 40 km/s and H = 16 kpc. As we observe from Fig. 2, at low masses, uncertainties in reacceleration make our limits vary by up to an order of magnitude, while at high masses the effect of reacceleration is minor and uncertainties in the halo height make our limits uncertain by a factor of up to 2-3, compatible with what was also found in (De la Torre Luque et al. 2023c). Effects like strong convection or winds near the Galactic Center could also alter the positron spectra (Bartels et al. 2017), but accurate estimation of this with current data is challenging. On top of this, we have checked that these limits do not change by more than a factor of a few using a cuspier DM distribution than the NFW. An Einasto profile is expected to slightly strengthen these limits, given that it results in more DM den-  \nty at intermediate longitudes. Following the same reasoning, the use of a more cored profile, like the Isothermal one, could worsen the limits. However, even with these uncertainties included, our results show that the 511 keV line still provides the strongest astrophysical bounds so far, especially for DM decay across the full MeV range.  \nFuture data from the COSI (Siegert et al. 2020) with improved spatial resolution would provide insights into processes and structures producing the 511 keV line in the Milky Way. Additionally, testing these DM scenarios by characterizing diffuse MeV-scale γ -ray emission from the Milky Way halo, using proposed γ -ray telescopes like AMEGO (McEnery et al. 2019) or e-ASTROGAM (De Angelis et al. 2018), could  \nprovide useful background models and further strengthen bounds on sub-GeV DM.", 'Acknowledgements': 'We are grateful to Pierluca Carenza, Marco Cirelli, Daniele Gaggero, Dan Hooper, Jordan Koechler, Tim Linden, Pierrick Martin and Thomas Siegert for useful discussions. SB is supported by the STFC under grant ST/X000753/1. PDL is supported by the European Research Council under grant 742104 and the Swedish National Space Agency under contract 117/19. This project used computing resources from the Swedish National Infrastructure for Computing (SNIC) under project Nos. 2021/3-42, 2021/6326, 2021-1-24 and 2022/3-27 partially funded by the Swedish Research Council through grant no. 2018-05973.', 'REFERENCES': 'Boehm, C., & Uwer, P. 2006, Revisiting Bremsstrahlung emission associated with Light  \nDark Matter annihilations.  \nhttps://arxiv.org/abs/hep-ph/0606058  \nBoehm, C., et al. 2004b, Phys. Rev. Lett., 92, 101301, doi: 10.1103/PhysRevLett.92.101301  \nCalore, F., et al. 2022, Phys. Rev. D, 105, 063026, doi: 10.1103/PhysRevD.105.063026  \nCardillo, M. 2019, Galaxies, 7, 49,  \ndoi: 10.3390/galaxies7020049  \n- Cirelli, M., et al. 2021, Phys. Rev. D, 103, 063022, doi: 10.1103/PhysRevD.103.063022\n- -. 2023, JCAP, 07, 026,  \ndoi: 10.1088/1475-7516/2023/07/026  \nCordes, J. M., & Lazio, T. J. W. 2003a, NE2001. II. Using Radio Propagation Data to Construct a Model for the Galactic Distribution of Free Electrons.  \nhttps://arxiv.org/abs/astro-ph/0301598  \n- -. 2003b, NE2001.I. A New Model for the Galactic Distribution of Free Electrons and its Fluctuations.  \nhttps://arxiv.org/abs/astro-ph/0207156  \nDe Angelis, A., et al. 2018, Journal of High Energy Astrophysics, 19, 1-106,  \ndoi: 10.1016/j.jheap.2018.07.001  \n```\nStrong, A. W., & Moskalenko, I. V. 1998, The Astrophysical Journal, 509, 212, doi: 10.1086/306470 Vegetti, S., & Vogelsberger, M. 2014, Mon. Not. Roy. Astron. Soc., 442, 3598, doi: 10.1093/mnras/stu1284\n```', 'A. HEIGHT SCALE OF THE FREE ELECTRON DENSITY': "As explained above, in order to consider that different regions of the Galaxy will have a different rate (probability) of positronium formation, that also depends on the distribution of the ambient (free) electrons, we consider a scaling of the 511 keV line emission following the NE2001 model Cordes & Lazio (2003a,b) for the free electron density distribution. In particular, we consider that the disk of the Galaxy will always contain a much larger density of free electrons in the ISM gas than diffuse positrons (i.e. there will be always more electrons available than positrons to form positronium states, thus being 'saturated', and the probability of positronium formation will mainly depend on the diffuse positron distribution). However, this assumption must be broken well above and below the disk, since the density of free electrons falls very quickly as the galactic height, z , is increased, as shown by Fig. 3. In particular, in this figure we depict the electron density distribution, at the center of the Galaxy, as a function of the Galactic height, normalized to the electron density at the GC position. This is the scale factor that we apply to our evaluations of the 511 keV line. This scaling is particularly important for very low DM masses, when the positrons lose energy very fast in the disk (but not outside the disk, since the density of gas producing these losses is much lower) and the fraction of positrons annihilated in flight (directly proportional to the gas density where the positrons propagate) becomes relevant.  \nThis is also motivated by the fact that, as we observe from the SPI measurements of the latitude profile, the intensity of the 511 keV line decays by around one order of magnitude at 10 · above or below the disk. Assuming the the height of the halo is around 4-10 kpc, this means that with every ∼ 0 . 36-0 . 9 kpc, the intensity of the line decays by one order of magnitude, similar to the factor ≳ 6 that the free electron density decays every ∼ kpc above or below the disk. While the predicted latitude profiles from DM appear to be quite flat (compared to the SPI data) for all the masses explored here, using this scaling we find that the predicted latitude profiles become much more similar to the  \nFigure 3. Vertical profile of the free electron density normalized at the Galactic Center as a function of height above the Galactic Plane z . We take this distribution from the NE2001 model Cordes & Lazio (2003a,b), as indicated in the text.  \n<!-- image -->  \nshape of the SPI measurements. However, we remark that the longitude profiles predicted are roughly unaffected by this scaling.", 'B. RELEVANT TIMESCALES': 'Our calculations include all the relevant Galactic CR propagation processes in a broad energy range. For signals from sub-GeV DM, it is crucial to account for diffusion, reacceleration, and energy losses, that become dominant below a few tens of MeV. The most relevant cooling processes at these energies are Coulomb and ionization losses, although we include in our simulation synchrotron, bremsstrahlung, inverse Compton and catastrophic losses (annihilation) as well. A comparison of the relevant timescales is illustrated in Fig. 4, where we show a band around the benchmark diffusion time (i.e. diffusion time of particles with our benchmark diffusion setup) to indicate the difference between the optimistic and pessimistic scenarios. Moreover, cooling is much faster than diffusion. We note that the average distance travel by a positron injected at 10 MeV (which can be calculated as ⟨ l travel ⟩ ≈ √ 2 Dτ , being τ the dominant timescale) is larger than 100 pc in a 1 cm -3 gas, without accounting for wind advection or reacceleration (that can boost ⟨ l travel ⟩ to even more than 500 pc).  \nFigure 4. Main timescales involved in the propagation of electrons and positrons generated by sub-GeV DM particles, for a gas density of 1 cm -3 composed by H and a 10% of He. Gas densities. The magnetic field strength here is set to 5 µ G and the energy density of the interstellar radiation fields is taken from Evoli et al. (2020). The grey band indicates the difference in diffusion time in the Optimistic and Pessimistic setups, used to calculate uncertainties in the DM bounds.  \n<!-- image -->', 'C. 511 KEV EMISSION FROM COMMON GALACTIC DARK MATTER DENSITY DISTRIBUTIONS': 'In this appendix, we report a comparison of the longitude (left panels of Fig. 5) and latitude (right panels of Fig. 5) profiles predicted for DM annihilation, assuming different standard Galactic DM distributions, for DM masses of 10 and 100 MeV. In particular, we show the predictions obtained for the NFW profile (cyan lines), a contracted NFW profile (c-NFW), like the one that provides the best-fit to the Galactic Center Excess (GCE) Ackermann et al. (2017), which has a slope of γ = 1 . 25 (green lines), a Moore profile Moore et al. (1999) , which has a slope of γ = 1 . 5 (blue lines), and a cuspy profile, as an example of a distribution that accounts for spikes formed around the central black hole Sgr A*, taken from Refs. Gondolo & Silk (1999); Lacroix (2018); Balaji et al. (2023) (black lines).  \nWe defer the full study of the compatibility of the different DM distributions with the SPI observations to a future work. However, this comparison allows us to see how the profiles are sensitive to the  \ndifferent choice of the Galactic DM distribution. In particular, we see that profiles that are not very steep (with slope γ ≲ 1 . 25) are not able to produce a good simultaneous fit to the longitudinal and latitudinal SPI profiles for DM mass exceeding an MeV. However, a steeper profile would be able to provide relatively satisfactory fits for MeV masses. However, we remark that a high slope, like that adopted in the Moore profile, is already in tension with some astrophysical observations. This is because it can predict too much DM in the inner regions of galaxies and galaxy clusters to comply with rotation curve observations Spekkens & Giovanelli (2005), X -ray observations of galactic clusters Walker et al. (2012), gravitational lensing observations Foex et al. (2012); Vegetti & Vogelsberger (2014) and N-body simulations of ΛCDM halos Navarro et al. (2004). In general, this leads to the conclusion that DM cannot be the dominant source of positrons in the Galaxy, but, instead, it is expected that a combination of different sources is producing the peaked 511 keV emission profiles that SPI observes.  \nFigure 5. Comparison of the predicted diffuse 511 keV line emission from a low-mass DM particle of mass m χ = 10 MeV (top panels) and m χ = 100 MeV (bottom panels) annihilating into electron-positron pairs for different DM density profiles.The left panels are the predicted longitude profiles, for latitude bin -10 . 75 · < b < 10 . 25 · while the right panels are the latitude profile, for a longitude bin of -5 . 25 · < l < 3 . 75 · . We fix the annihilation rate to the values shown in the plots and scale the profiles for different masses to facilitate the comparison with the SPI data. We show the NFW profile (cyan), contracted NFW (green), Moore (blue) and cuspy DM profiles.  \n<!-- image -->', 'D. ANNIHILATION AND DECAY PROFILES FOR DIFFERENT CHANNELS': 'In this appendix, we report the predicted latitude and longitude profiles of the 511 keV emission from DM decay and annihilation. In particular, we show in Fig. 6 a comparison of profiles obtained assuming DM annihilation through the µ + µ -(top panels), π + π -(middle panelbos) and direct e + e -(bottom panels) channels for different DM masses (from twice the mass of the final state, around 100 MeV, to 5 GeV, in the case of the π + π -and µ + µ -channels, and from 1 MeV to 1 GeV for the direct e + e -channel). The left panels illustrate the predicted longitude profiles of the 511 keV emission and the right panels the predicted latitude profiles. These comparisons allow us to see that the data points corresponding to the higher longitude and latitude angles (in absolute value) are the most  \nFigure 6. Comparison of the predicted longitude (left panels) and latitude (right panels) of the 511 keV emission from DM annihilation for different DM masses, through the µ + µ -(top panels), π + π -(middle panels) and direct e + e -(bottom panels) channels. The lowest mass for the case of the π + π -annihilation channel is actually m χ = 2 m π ≃ 139 MeV and the minimum energy shown in the µ + µ -is m χ = 2 m µ ≃ 105 MeV. We fix the annihilation rate to the values shown in the plots and scale the profiles for different masses to facilitate the comparison with the SPI data. The DM profiles are shown in the same colours and the left and right panels are the same latitude bins as Fig. 5.  \n<!-- image -->  \nconstraining ones, which explains why the bounds derived from this dataset are much stronger than those derived from Ref. Vincent et al. (2012) (besides the fact that these measurements are reported after double the duration of event collection). As we see, the predicted longitude profiles show a clear trend: the lower the mass of the DM particle, the more concentrated the 511 keV emission around the GC. This is expected, since a lower DM mass means a lower injection energy for the positrons, and the lower the energy of the positrons, the slower they diffuse. Hence lower energy positrons have a distribution more similar to the DM distribution (NFW in these calculations) and lose their energy closer to their injection point. As we pointed out in the previous section, the predictions from a NFW DM distribution do not lead to a 511 keV emission longitude profile as peaky as that observed by the SPI data even for the low DM masses that are accessible in the direct e + e -channel. In turn, as we see from the predicted latitude profiles, masses of around 100 MeV reproduce the observed SPI latitude profile quite well. We remark that this is due to the scaling that accounts for the decay of the density of free electrons away from the GC, applied for first time in this work. Without accounting for this effect, we observe that the latitude profiles become much more diffuse (i.e. flatter profiles) than what is observed. In addition, for the latitude profiles, we see that the same trend wherein lower DM mass leads to a profile more concentrated around the Galactic Plane is observed above some tens of MeV. However, once the energy loss (mainly Coulomb and ionization energy losses) becomes the dominant process during the positron propagation, the 511 keV emission at the Galactic plane starts to decay quickly, given that these particles suffer from severe energy losses and annihilation in flight losses (since this effect is proportional to the gas density, it does not affect the positrons above and below the Galatic plane). We have investigated that the minimum energy employed in these simulations (100 eV) is low enough that our results do not change appreciably for a minimum energy of 1 eV.  \nWe show a direct comparison of the latitude profile including the density scaling (solid lines) and without including it (dashed lines) for different DM masses in Fig. 7, normalized to facilitate its visual comparison. This illustrates the importance of accounting for the density scaling to realistically evaluate the latitude profile of the 511 keV diffuse emission.  \nFigure 7. Comparison of the evaluated latitude profiles accounting for the scaling with the electron density (solid lines) and without considering it (dashed lines).  \n<!-- image -->  \nFinally, we show in Fig. 8 the predicted profiles from decay of DM for masses of 1 MeV to 1 GeV, in the direct e + e -channel. As already observed in the past, and discussed above, we observed that both the latitude and longitude profiles of the predicted 511 keV emission are much flatter (for every mass) than the ones measured by SPI, disfavouring DM decay as the main producer of the emission around the GC. This is simply due to the fact that the total positron emission from decay is proportional to the DM density, while in the case of annihilation the emission is ∝ ρ 2 DM .  \nFigure 8. Similar to Fig. 6 but for the case of DM decay into the direct e + e -channel.  \n<!-- image -->'}