An RXTE Search for the Sterile Neutrino Decay in Galaxy Clusters

We have used long observations of galaxy clusters obtained with the Rossi X-ray Timing Explorer to search for the 3.55 keV line from sterile neutrino decay. If a lepton-number asymmetry exists in one or more types of active neutrinos in the early Universe, sterile neutrinos can be produced via the Shi–Fuller mechanism. The data consist of 11 clusters observed for a total of 3.1 megaseconds using the Proportional Counter Array. A 2.5σ excess of emission over a thermal model is found over the energy span of the 3.55 keV line in the combined spectra of the eight clusters that individually have an excess. These residuals are added to increase the signal to noise ratio of the excess, which is then modeled with a Gaussian to simulate the instrumental spectral response. We find a significant correlation (r = 0.76) for a line centered at 3.6 keV with a model flux of 3.07 × 10⁻⁵ ph cm⁻² s⁻¹. Mixing angle for detected clusters ranges from 2.0 to 21.6 × 10⁻¹⁰. The decay rate inferred from the line flux is strongly correlated (r = 0.87) with cluster temperature, which is due to hotter, more massive clusters having a larger amount of dark matter. Approximately half of the total flux comes from the Coma cluster. The mixing angle for Coma is calculated to be 6.2 × 10⁻¹⁰. We fit the Coma cluster spectrum with two different three-component models. The first includes a Gaussian fixed at 3.55 keV to model soft emission. The flux of the Gaussian is 5.6 × 10⁻¹² ph cm⁻² s⁻¹ or 1.3% of the total flux. The second three-component model uses a second thermal component to model soft emission. This model gives a temperature of 0–17 keV for the second thermal component and a lower temperature for the hot component. This indicates that the second thermal component is modeling high-energy residuals rather than low ones, where the Gaussian is. Though our line fluxes exceed most reported detections and upper limits, they do not overproduce the dark matter. We conclude that some fraction of the marginally detected excess could be attributed to the decay line since low-temperature thermal emission and systematics fail to model it completely.