skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: An INTEGRAL/SPI view of reticulum II: particle dark matter and primordial black holes limits in the MeV range
ABSTRACT Reticulum II (Ret II) is a satellite galaxy of the Milky Way (MW) and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain γ-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the γ-ray band 25–8000 keV, and we derive a flux limit of ${\lesssim}10^{-8}\, \mathrm{erg\, cm^{-2}\, s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of ${\lesssim}1.7 \times 10^{-4}\, \mathrm{ph\, cm^{-2}\, s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of e+s produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses ${\lesssim}8 \times 10^{15}\, \mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into e+e− are $\langle \sigma v \rangle \lesssim 5 \times 10^{-28} \left(m_{\rm DM} / \mathrm{MeV} \right)^{2.5}\, \mathrm{cm^3\, s^{-1}}$. We conclude that analysing isolated targets in the MeV γ-ray band can set strong bounds on DM properties without multi-year data sets of the entire MW, and encourage follow-up observations of Ret II and other dwarf galaxies.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
914 to 924
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The nucleosynthesis in classical novae, in particular that of radioactive isotopes, is directly measurable by its γ-ray signature. Despite decades of observations, MeV γ-rays from novae have never been detected – neither individually at the time of the explosion, nor as a result of radioactive decay, nor the diffuse Galactic emission from the nova population. Thanks to recent developments in modelling of instrumental background for MeV telescopes such as INTEGRAL/SPI and Fermi/GBM, the prospects to finally detect these elusive transients are greatly enhanced. This demands for updated and refined models of γ-ray spectra and light curves of classical novae. In this work, we develop numerical models of nova explosions using sub- and near-Chandrasekhar CO white dwarfs as the progenitor. We study the parameter dependence of the explosions, their thermodynamics and energetics, as well as their chemical abundance patterns. We use a Monte Carlo radiative transfer code to compute γ-ray light curves and spectra, with a focus on the early time evolution. We compare our results to previous studies and find that the expected 511-keV-line flash at the time of the explosion is heavily suppressed, showing a maximum flux of only $10^{-9}\, \mathrm{ph\, cm^{-2}\, s^{-1}}$ and thus making it at least one million times fainter than estimated before. This finding would render it impossible for current MeV instruments to detect novae within the first day after the outburst. Nevertheless, our time-resolved spectra can be used for retrospective analyses of archival data, thereby improving the sensitivity of the instruments. 
    more » « less
  2. ABSTRACT A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a ‘velocity kick’ parameter. The simplest 2cDM model with cross-section $0.01\lesssim \sigma /m\lt 1\, \textrm {cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{\mathrm{ k}}\simeq 100\, \rm {km\, s}^{-1}$ have been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in N-body cosmological simulations tested on Milky Way (MW)-like haloes of a virial mass ${\sim}5 \times 10^{11}\, {\rm M_{\odot }}$ (Papers I & II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of ∼107−108 and ${\sim}10^{13} - 10^{14}\, {\rm M_{\odot }}$, respectively. We find that σ0/m ≳ 0.1 cm2g−1 is preferentially disfavoured for both dwarfs and galaxy cluster haloes in comparison with observations, while σ0/m = 0.001 cm2g−1 causes little perceptible difference from that of the CDM counterpart for most of the cross-section’s velocity dependence studied in this work. Our main result is that within the reasonable set of parameters, the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes, and the model leaves us an open window for other possible alternative DM models. 
    more » « less
  3. null (Ed.)
    Abstract Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 $$\hbox {C}_3\hbox {F}_8$$ C 3 F 8 superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (JCAP 1509: 052, 2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $$\sigma _{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm {cm}^2$$ σ SD ≲ 3 × 10 - 39 cm 2 ( $$6 \times 10^{-38} \mathrm {cm}^2$$ 6 × 10 - 38 cm 2 ) at $$\gtrsim 90\%$$ ≳ 90 % C.L. for a DM particle of mass 1 TeV annihilating into $$\tau ^+ \tau ^-$$ τ + τ - ( $$b\bar{b}$$ b b ¯ ) with a local density of $$\rho _{\mathrm{DM}} = 0.3~\mathrm {GeV/cm}^3$$ ρ DM = 0.3 GeV / cm 3 . The constraints scale inversely with $$\rho _{\mathrm{DM}}$$ ρ DM and are independent of the DM velocity distribution. 
    more » « less
  4. Abstract

    We report on a search for electron antineutrinos (ν¯e) from astrophysical sources in the neutrino energy range 8.3–30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we find no significant excess over background model predictions. Assuming several supernova relic neutrino spectra, we give upper flux limits of 60–110 cm−2s−1(90% confidence level, CL) in the analysis range and present a model-independent flux. We also set limits on the annihilation rates for light dark matter pairs to neutrino pairs. These data improve on the upper probability limit of8B solar neutrinos converting intoν¯e,Pνeν¯e<3.5×105(90% CL) assuming an undistortedν¯eshape. This corresponds to a solarν¯eflux of 60 cm−2s−1(90% CL) in the analysis energy range.

    more » « less
  5. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant. 
    more » « less