{"id":"d23c3dfe-91c4-4c6f-a1da-436916c32dd1","arxiv_id":"1908.09037","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A 190 kilosecond NuSTAR search at high Galactic latitude finds no sterile-neutrino dark matter decay line and improves the leading upper limit by a factor of about two for 10-12 keV masses.","lead":"Astronomers used the NuSTAR space telescope to search for the faint X-ray line that sterile-neutrino dark matter would produce when it decays. They found no line, and in the 10-12 keV mass range they set the world's best limits on this decay, about twice as strong as before.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Frozen internal background power-law calibration from Earth-occulted data is the pivotal assumption; freeing it in science fits would test whether the 10–12 keV factor-of-2 improvement is robust.","rationale":"The paper's headline claim is the factor-of-2 improvement in sterile-neutrino limits at 10–12 keV. The sensitivity in this band rests on the accuracy of the internal background model, especially the low-energy power-law that is calibrated on Earth-occulted data and frozen in science fits (Sec. II D). If that power-law normalization or index were biased, it would shift the continuum and the derived line limits. The paper reports residuals (8–9 keV, 15–20 keV) that demonstrate the model is not perfect, and the MC simulations in Sec. III B are generated from the same model, so they validate the procedure rather than the absolute background. However, the paper partially mitigates this: residuals are localized and appear in only two of four spectra; the 10–12 keV negative residuals are narrower than the detector resolution and at different energies in three modules, pointing to statistics; and a flat 7.5% systematic weakens the 10–12 keV limit by only ~1.5, leaving it still ~1.5 stronger than Ref. [31]. A correlated error in the power-law shape is not fully covered by the flat systematic test, so the proposed freeing test is the right check. Given the good fit quality and the existing robustness tests, I do not expect the concern to overturn the central claim; the verdict remains ACCEPT.","tokens_in":24702,"tokens_out":9933,"duration_ms":103692,"concrete_test":"Free the internal power-law spectral index and normalization in the science-mode fits (per FPM/obsID), with the 95–110 keV band anchoring the internal continuum, and recompute the 95% upper limits on Γ for 10–40 keV. Compare the 10–12 keV limits to the frozen-calibration result. If the limits weaken by more than 50%, the factor-of-2 improvement is not robust to the calibration assumption; if they shift by less than 30%, the frozen calibration is adequate. As a secondary cross-check, repeat using power-law parameters drawn from the covariance of the Earth-occultation fit to propagate its uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that these observations provide the leading constraints on sterile-neutrino dark matter in the 10–12 keV mass range depends on the absolute accuracy of the low-energy internal background model. The internal power-law that dominates below ~10 keV is not fit to the science data; its spectral index and relative normalization are frozen to values obtained from Earth-occulted data (Sec. II D). This assumes the internal background is identical between occulted and science modes and that a single power-law accurately describes the residual instrumental component across 5–20 keV. The paper itself reports residuals at 8–9 keV and 15–20 keV (Sec. II D, Figs. 3–4), demonstrating the model is imperfect. A normalization or index bias in this frozen power-law would create a broad continuum error and directly shift the derived DM line limits; the specific 10–12 keV improvement is partly driven by negative residuals in three of the four spectra (Sec. III B). The Monte Carlo validation in Sec. III B generates mock spectra from the same best-fit model, so it tests the statistical coverage of the fitting procedure under the assumed model, not the absolute accuracy of that model. The flat 7.5% systematic test bounds the effect to a factor ~1.5, which is helpful, but a flat systematic in quadrature does not mimic a correlated shape error of the power-law normalization or index. Thus, the headline improvement is robust only if the occultation-calibrated power-law is valid at the ~10% level in the 5–6 keV band.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a search for monoenergetic x-ray lines from the radiative decay of sterile-neutrino dark matter using two dedicated NuSTAR observations of the Galactic bulge region, with a combined cleaned exposure of approximately 190 ks. The fields are located approximately 10 degrees above and below the Galactic plane, chosen to reduce astrophysical backgrounds while remaining close to the dense dark-matter halo. The analysis models the instrumental and astrophysical backgrounds with a six-component spectral model, including a low-energy internal power-law calibrated from Earth-occulted data and frozen during science-mode fits. No significant line is found in the 5-20 keV band, and upper limits are set on the decay rate for sterile-neutrino masses 10-40 keV. The authors report that their constraints are the strongest to date in the 10-12 keV mass range, improving on previous NuSTAR limits by a factor of about two. Monte Carlo simulations are used to show that the observed limits are consistent with statistical fluctuations, and a flat 7.5% systematic is found to weaken the limits by a factor of about 1.5.","tokens_in":24944,"tokens_out":6894,"duration_ms":70134,"significance":"If the systematic treatment is adequate, this result provides a meaningful improvement in the sterile-neutrino dark-matter parameter space, specifically in the 10-12 keV mass range where the leading constraints are currently set by NuSTAR. The paper is careful in several respects: independent fits are performed for each focal-plane module, the line-search procedure conservatively allows a dark-matter line to absorb background lines, the internal power-law is calibrated on Earth-occulted data, and the Monte Carlo validation demonstrates that the observed limits are consistent with statistical expectations. The explicit discussion of model residuals and systematic tests, including a flat 7.5% systematic, is commendable and provides some confidence in the robustness of the no-detection claim. The result has implications for the νMSM and for future x-ray searches for decaying dark matter, and the paper clearly identifies the need for an improved NuSTAR background model, especially at lower energies.","major_comments":[{"comment":"The internal power-law spectral index and normalization, which dominate the instrumental background below about 10 keV, are frozen to values derived from Earth-occulted data and are not allowed to vary in the science-mode fits. The paper itself reports residuals at 8-9 keV and 15-20 keV (Figs. 3 and 4), demonstrating that the model is imperfect. The Monte Carlo simulations in Sec. III B generate mock spectra from the same best-fit model, so they test the statistical coverage of the fitting procedure under the assumed background, not the absolute accuracy of that background. The flat 7.5% systematic test is a global normalization uncertainty and does not mimic a correlated shape error in the power-law normalization or index. Because the headline improvement in the 10-12 keV range is partly driven by negative residuals in three of the four spectra, a correlated background-shape bias could in principle alter the derived limits. I therefore request an additional robustness test in which the internal power-law normalization and index are allowed to vary, for example with Gaussian priors centered on the occultation values, and the resulting change in the 10-12 keV limits is reported. This would directly address whether the factor-of-two improvement is robust against plausible background-model uncertainties and would substantially strengthen the central claim.","section":"II D, III B"}],"minor_comments":[{"comment":"The Monte Carlo validation uses 100 mock spectra generated from a single best-fit spectrum (FPMA of obsID 40410001002) to set the expected limit bands, with only 10 full realizations used for a cross-check. It would be useful to state the statistical uncertainty on the quoted 68% and 95% bands arising from the finite number of simulations, and to clarify whether the full-realization cross-check was sufficient to ensure that the simplified procedure does not bias the expected bands.","section":"III B"},{"comment":"The GRXE component is described as absorbing contribution from un-modeled point sources, reflected Earth x-rays, and any low-energy instrumental backgrounds not captured by the default model. This is a sensible catch-all, but it would be clearer to the reader if the caption of Table II explicitly stated that the quoted GRXE flux is therefore an effective (background-inclusive) value and not a pure measurement of Galactic ridge emission, even though this is mentioned in the text.","section":"II D and Table II"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid, well-written analysis, but the frozen internal power-law calibration is a central assumption whose effect on the headline improvement is not fully quantified. The requested test—freeing the power-law parameters with priors—is straightforward and would convert a potential concern into a demonstrated robustness. I believe the authors can address this with a modest addition to the analysis. If they do so, the paper would be acceptable. No other issues are of concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper before the field moves on: it is the first NuSTAR observation dedicated to a dark-matter line search, and it delivers the leading constraint on sterile-neutrino decay in the 10–12 keV mass range, a factor ~2 better than the previous record, with only ~190 ks of exposure. The improvement is real, though part of it comes from downward statistical fluctuations, which the authors say plainly. That honesty is typical of the paper: it reports residuals, tests its own assumptions, and compares its limits to Monte Carlo expectations. It is a well-executed extension of an established 0-bounce technique, not a methodological breakthrough, but the new observations and the factor-of-two gain are genuinely new. What the paper does well: the line search is conservative, letting the DM line absorb background lines, so background excesses weaken limits rather than create false signals; the four spectra are fit independently and combined only at the chi-squared level; J-factor uncertainties are checked against several density profiles; and the MC simulations show the observed limits sit inside the expected statistical band. The authors also disclose that the 10–12 keV improvement is driven partly by negative residuals that appear at different energies in three modules, with bin widths much narrower than the detector resolution—evidence that those residuals are statistical rather than a coherent background feature. The soft spots are real but proportionate. The pivotal assumption is that the low-energy internal background power-law, calibrated on Earth-occulted data, is valid in science mode at the ~10% level across 5–20 keV. The paper itself shows residuals at 8–9 and 15–20 keV, and the Monte Carlo validation is conditional on the same background model, so it tests statistical coverage, not absolute background accuracy. The flat 7.5% systematic test is useful but does not mimic a correlated shape error in the power-law index or normalization. I would have liked a test where the power-law index is freed in science fits to see how much the 10–12 keV limit moves. Still, these are caveats, not a fatal flaw: the residuals are mostly excesses, the conservative scanning procedure pushes limits in the safe direction, and the factors involved are consistent with the claimed ~2x improvement being robust at the factor ~1.5 level. Bottom line: this is a solid, citable result for anyone working on sterile-neutrino or keV-mass dark matter. It deserves a serious referee and, in my view, publication. The main message—that low-background off-plane fields are far more efficient than deep pointings near the plane—is well supported and likely to shape future observing proposals.","headline":"A careful, transparent null result that improves the world-best sterile-neutrino DM limit at 10–12 keV by ~2x; the main caveat is the frozen internal background power-law, but the paper's own checks keep the central claim credible.","tokens_in":714,"tokens_out":1082,"would_cite":true,"duration_ms":23046,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"No sterile-neutrino dark-matter decay lines appear in 190 ks of dedicated NuSTAR observations, tightening the best limits at 10–12 keV masses by roughly a factor of two.","keywords":["sterile neutrino dark matter","x-ray line search","NuSTAR","0-bounce aperture","radiative decay","Galactic bulge","νMSM","indirect dark matter detection"],"falsifier":"An independent re-analysis of the same observations with an updated NuSTAR background model should either reproduce the 8–9 and 15–20 keV excesses as statistical fluctuations or identify them as instrumental lines; if they are instrumental, the limits in the affected mass ranges are weaker than quoted, and a true dark-matter line in those regions could be hidden.","tokens_in":24461,"feed_emoji":"🔭","tokens_out":6198,"duration_ms":57959,"temperature":0.7,"pith_summary":"This paper tries to establish whether sterile neutrinos can be the dark matter by looking for the x-ray line their radiative decay would produce. It analyzes two dedicated NuSTAR observations, totaling about 190 ks, aimed roughly ten degrees above and below the Galactic plane to keep astrophysical backgrounds low while staying near the dark-matter-dense bulge. No anomalous line is found between 5 and 20 keV, which corresponds to sterile-neutrino masses of 10–40 keV. Interpreted through neutrino-mixing production, the nondetection gives the strongest constraints yet in the 10–12 keV mass range, improving the previous best limits by a factor of about two.","feed_headline":"NuSTAR sees no sterile-neutrino decay lines in 5-20 keV","feed_subtitle":"Two dedicated 190 ks exposures tighten leading sterile-neutrino limits at 10-12 keV by a factor of two.","key_machinery":"The mechanism that does the work is NuSTAR's 0-bounce aperture: unfocused x-rays reach the detector array over a large effective field of view of roughly 4.5 square degrees, and the predicted dark-matter signal is the line flux times the response $E_{\\rm Be}(E) A_{0b}\\Delta\\Omega_{0b} J$, where $J$ is the line-of-sight integral of the dark-matter density in each field. The analysis chain scans a narrow Gaussian line at each mass bin, re-fits the full six-component spectral model, and combines the four spectra's $χ^2$ curves to set a 95% upper limit on the decay rate. The internal background's low-energy power-law is calibrated on Earth-occulted data and frozen in science-mode fits.","core_discovery":"The central claim is that the radiative decay channel $χ \\to \\nu + \\gamma$ does not produce a detectable monochromatic x-ray line in the 5–20 keV band of these observations, and the resulting upper limits on the decay rate are the leading constraints for sterile-neutrino dark matter in the 10–12 keV mass range. The two fields were modeled individually, combining 0-bounce and 2-bounce apertures, and the line search allowed a putative dark-matter line to absorb any background excess, a conservative choice. The limits are comparable to earlier NuSTAR searches across 10–40 keV despite a factor of roughly fifty less exposure, and Monte Carlo simulations indicate the observed limit fluctuations sit inside the expected statistical band, apart from known excesses at 8–9 and 15–20 keV that the authors treat as consistent with statistics.","pith_inferences":["A repeat observation of these two fields, or a search that splits the exposure into time bins, could test whether the strong 10–12 keV limit is partly a downward statistical fluctuation, as the paper itself suggests.","If the 8–9 and 15–20 keV excesses prove to be unmodeled instrumental lines, an updated background model could either strengthen the limits in those mass ranges or reveal that the current conservative procedure hides a real line.","The same off-plane strategy could be extended to other high-$J$-factor, low-background targets or to lower energies once the 3–5 keV background is modeled, potentially testing the 3.5-keV line interpretation."],"forward_implications":["Sterile neutrinos in the 10–12 keV mass range with mixing angles above the new limit are excluded, narrowing the parameter space left for the νMSM.","Dedicated low-background fields near the bulge reach sensitivity comparable to searches with several megaseconds of exposure in only ~190 ks, validating the field-selection strategy.","In the full 10–40 keV mass range the null result keeps the decay-rate limit comparable to previous NuSTAR blank-sky, Galactic-center, and M31 searches.","Closing the remaining νMSM window above 10 keV will require roughly a fourfold sensitivity gain, or about 4 Ms of similar exposures, and an improved instrumental background model for energies below 5 keV."],"supporting_citations":[{"why":"Previous NuSTAR blank-sky search; supplies the leading 10–12 keV limits that this paper improves by a factor of about two.","marker":"[31]"},{"why":"Previous NuSTAR Galactic-center search whose line-search and 0-bounce methods this analysis follows.","marker":"[30]"},{"why":"Previous NuSTAR M31 analysis; source of the 0-bounce technique, the 3–5 keV exclusion, and the χ2-combination procedure.","marker":"[32]"},{"why":"Default NuSTAR instrumental background model, including the internal continuum, lines, and effective-area calculations.","marker":"[98]"},{"why":"Radiative decay relation that converts decay-rate limits to mixing-angle limits for sterile-neutrino dark matter.","marker":"[17, 18]"},{"why":"HEAO-1 and INTEGRAL measurements fixing the cosmic X-ray background flux and spectral index.","marker":"[105, 106]"},{"why":"NFW dark-matter density profile and Galactic bulge density model used to compute J-factors.","marker":"[121, 122]"},{"why":"Dynamical constraints on the Milky Way dark-matter distribution used for the conservative sNFW profile and local density.","marker":"[123]"}],"fun_headline_variants":["NuSTAR finds no sterile-neutrino decay lines in new bulge data","NuSTAR tightens sterile-neutrino limits by factor of two","No sterile-neutrino line, but NuSTAR improves low-mass bounds","NuSTAR's 190-ks stares rule out sterile-neutrino lines","Sterile-neutrino dark matter eludes NuSTAR's deep search"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes the NuSTAR instrumental background in the 5–20 keV search band is described by the default model with the low-energy power-law fixed from Earth-occulted data, and the residual excesses at 8–9 and 15–20 keV show that this assumption is only approximate.","fun_headline_variants_meta":{"raw":{"variants":["NuSTAR finds no sterile-neutrino decay lines in new bulge data","NuSTAR tightens sterile-neutrino limits by factor of two","No sterile-neutrino line, but NuSTAR improves low-mass bounds","NuSTAR's 190-ks stares rule out sterile-neutrino lines","Sterile-neutrino dark matter eludes NuSTAR's deep search"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000309,"raw_usage":{"total_tokens":1753,"prompt_tokens":923,"completion_tokens":830,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":729}},"tokens_in":539,"tokens_out":830,"duration_ms":8442,"temperature":1.0,"reasoning_tokens":729,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:23:55.187564+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent re-analysis of the same observations with an updated NuSTAR background model should either reproduce the 8–9 and 15–20 keV excesses as statistical fluctuations or identify them as instrumental lines; if they are instrumental, the limits in the affected mass ranges are weaker than quoted, and a true dark-matter line in those regions could be hidden.","supporting_citations":[],"review_version":1}