NuSTAR tests of sterile-neutrino dark matter: New Galactic bulge observations and combined impact Brandon M. Roach,1, ∗ Kenny C. Y. Ng,2, y Kerstin Perez,1, z John F. Beacom,3, 4, 5, x Shunsaku Horiuchi,6, { Roman Krivonos,7, ∗∗ and Daniel R. Wik8, yy 1Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 2Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot, Israel 3Center for Cosmology and AstroParticle Physics (CCAPP), Ohio State University, Columbus, OH 43210, USA 4Department of Physics, Ohio State University, Columbus, OH 43210, USA 5Department of Astronomy, Ohio State University, Columbus, OH 43210, USA 6Center for Neutrino Physics, Department of Physics, Virginia Tech, Blacksburg, VA 24061, USA 7Space Research Institute of the Russian Academy of Sciences (IKI), Moscow 117997, Russia 8Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA (Dated: Received 4 October 2019; accepted 2 April 2020; published 8 May 2020) We analyze two dedicated NuSTAR observations with exposure ∼190 ks located ∼10◦ from the Galactic plane, one above and the other below, to search for x-ray lines from the radiative decay of sterile-neutrino dark matter. These fields were chosen to minimize astrophysical x-ray backgrounds while remaining near the densest region of the dark matter halo. We find no evidence of anomalous x-ray lines in the energy range 5{20 keV, corresponding to sterile neutrino masses 10{40 keV. Interpreted in the context of sterile neutrinos produced via neutrino mixing, these observations provide the leading constraints in the mass range 10{12 keV, improving upon previous constraints in this range by a factor ∼2. We also compare our results to Monte Carlo simulations, showing that the fluctuations in our derived limit are not dominated by systematic effects. An updated model of the instrumental background, which is currently under develop- ment, will improve NuSTAR's sensitivity to anomalous x-ray lines, particularly for energies 3{5 keV. doi: https://doi.org/10.1103/PhysRevD.101.103011 I. INTRODUCTION the detector, allowing for a rejection of photons from known astrophysical sources. Final states with mono- Multiple lines of cosmological evidence indicate that energetic photons are particularly valuable for indirect 80% of the matter density of the Universe, and 25% DM searches, as they result in line-like signals atop a of∼ its energy density, is nonbaryonic and nonluminous,∼ (usually) smooth continuum background. hence its name, dark matter (DM) [1]. At present, the A popular DM candidate with mχ keV is the ster- effects of DM are only measurable via its gravitational ile neutrino, with models such as the∼ νMSM provid- effects on astronomical scales, ranging from the motions ing explanations for the particle nature of DM, neutrino of galaxies and galaxy clusters to the power spectrum masses, and baryogenesis [9{12]. The radiative decay of of the Cosmic Microwave Background [2{7]. The lack sterile neutrinos via χ ν + γ would produce a mono- of a viable Standard Model candidate for particle DM energetic x-ray photon! and an active neutrino, each with (hereafter symbolized χ) has led to a plethora of theo- E = mχ = 2 [13{20]. retical models, many of which are also motivated by a Sterile neutrinos may be produced in the early Uni- desire to account for other phenomena not explained by verse via mixing with active neutrinos [21], and this pro- the Standard Model (e.g., baryogenesis, neutrino masses, duction may be resonantly enhanced by primordial lep- the hierarchy problem, etc). ton asymmetry [22]. Considerations from big bang nucle- The techniques of indirect detection use astronomi- osynthesis (BBN) [23{25] provide an upper bound on the cal observations to search for the decay and/or anni- cosmological lepton asymmetry per unit entropy density 6 hilation of DM into Standard Model particles such as L6 10 (nν nν¯)=s 2500, which we translate into ≡ − ≤ 2 arXiv:1908.09037v3 [astro-ph.HE] 8 May 2020 electrons/positrons, (anti)protons/nuclei, neutrinos, and the constraint on the active-sterile mixing angle sin 2θ photons [8]. Because photons are not deflected by as- shown in Fig.1 using the sterile-dm code [26]. We trophysical magnetic fields, it is possible to determine note that these BBN limits are particularly sensitive to their arrival direction within the angular resolution of the treatment of neutrino opacities and the plasma equa- tion of state near the QCD phase transition, with dif- ferent calculations finding different results|for example, the limits shown in Refs. [25, 27] for the same value of L6 ∗ [email protected]; orcid.org/0000-0001-8016-2170 y are nearly an order of magnitude less constraining than [email protected]; orcid.org/0000-0001-8016-2170 those from Ref. [26], which the authors of Refs. [16, 28] z [email protected]; orcid.org/0000-0002-6404-4737 x [email protected]; orcid.org/0000-0002-0005-2631 attribute to differences in the treatment of neutrino opac- { [email protected]; orcid.org/0000-0001-6142-6556 ities in the QCD epoch. (An update to the calculation in ∗∗ [email protected]; orcid.org/0000-0003-2737-5673 Ref. [27] is presented in Ref. [41], though the latter does yy [email protected]; orcid.org/0000-0001-9110-2245 2 not present an updated constraint in the mχ sin 2θ − 2 tion mechanism, with all of the sources of uncertainty discussed in the previous paragraph [55{58]. 10 10 NuSTAR Space-based x-ray observatories such as HEAO-1 [59], Chandra [60, 61], XMM-Newton [59, 62{64], Suzaku Chandra + XMM [38, 65], Fermi-GBM [39], and INTEGRAL [40, 66] have Suzaku provided the most robust constraints on the χ ν + γ 11 ! 10 decay rate for mχ 1{100 keV. The observation of an MW Sat. ' Counts unknown x-ray line at E 3:5 keV (\the 3.5-keV line") in several analyses [33{35] has' led to much interest, as well 2 2 as many follow-up analyses using different instruments 12 n Fermi-GBM i 10 and astrophysical targets [29, 37, 38, 63, 64, 67{83]. s +INTEGRAL Some suggest that the 3.5-keV line may be a signature of sterile-neutrino DM [84] or other DM candidates [85{ 13 89]; alternatively, modeling systematics [69, 71] or novel 10 astrophysical processes [90, 91] may play a role. Future BBN Limit (Resonant Production) high-spectral-resolution x-ray instruments may also be able to investigate the DM hypothesis for the origin of the 3.5-keV signal via velocity spectroscopy [92, 93]. 10 14 6 7 8 910 20 30 40 50 Since its launch in 2012, the NuSTAR observatory, due m [ keV ] to its unique large-angle aperture for unfocused x-rays, has provided the leading constraints on sterile-neutrino DM across the mass range 10{50 keV, leveraging observa- FIG. 1. The combined impact on the νMSM parameter space tions of the Bullet Cluster [29], blank-sky fields [31], the of previous NuSTAR searches [29{32] and this work is indi- Galactic center [30], and the M31 galaxy [32]. In each cated by the green region. This work provides the leading of these cases, the NuSTAR observations were originally constraints in the 10{12 keV mass range, as shown in Fig.5. performed to study non-DM phenomena; therefore, DM The tentative E ' 3:5 keV signal [33{35] is indicated by the red point. Constraints from other x-ray instruments [36{40] searches using these data had to contend with large astro- are shown for comparison. Uncertainties associated with MW physical backgrounds and/or reduced effective areas from satellite counts [28] and BBN [26, 27] are discussed in Sec.I. masking bright point sources in the field of view (FOV). Improving upon these constraints, and extending them to the NuSTAR limit of E = 3 keV (e.g., to test the ten- plane.) This lower bound may evolve as calculations are tative 3.5-keV signal), will therefore require observations refined. with lower astrophysical backgrounds, as well as an im- proved model of the low-energy NuSTAR instrumental An additional indirect constraint on sterile-neutrino background. DM arises from comparing the observed number of Milky In this paper, we present new constraints on the decay Way (MW) satellite galaxies to the results of N-body rate of sterile-neutrino DM particles using two NuSTAR cosmological simulations. Compared to cold DM, warm observations, one 10◦ above and the other 10◦ below DM particles are expected to suppress the matter power the Galactic plane,∼ chosen to minimize astrophysical∼ spectrum at small scales, reducing the number of low- x-ray emission while still remaining near the center of mass DM subhaloes orbiting the Galaxy. In Fig.1, the Galactic DM halo. These are the first NuSTAR we adopt the result of Ref. [28] with N = 47, de- subhalo observations dedicated to DM searches. rived from SDSS data. Though a complete review of subhalo constraints on the properties of particle DM is beyond the scope of this paper, we note several impor- In Sec.II, we describe the data reduction and spectral tant points. First, the Milky Way satellite population modeling of the NuSTAR data, consistently incorporat- may not resemble that of a typical galaxy of its size ing the flux from the focused and unfocused FOVs. In and morphology, and surveys of dwarf galaxies targeting Sec. III, we combine the line flux limits from these new observations to constrain the χ ν + γ decay rate for their stellar content must be corrected for completeness ! [42]. To address the former issue, surveys such as Satel- sterile neutrinos in the mass range 10{40 keV, obtaining lites Around Galactic Analogues [43] aim to study the the strongest constraints to date in the 10{12 keV mass satellites of Milky Way analogues in the local Universe.