We detect a weak unidentified emission line at E = (3.55−3.57)±0.03 keV in a stacked XMM-Newton spectrum of 73 galaxy clusters spanning a redshift range 0.01 − 0.35. MOS and PN observations independently show the presence of the line at consistent energies. When the full sample is divided into three subsamples (Perseus, Centaurus+Ophiuchus+Coma, and all others), the line is seen at > 3σ statistical significance in all three independent MOS spectra and the PN "all others" spectrum. The line is also detected at the same energy in the Chandra ACIS-S and ACIS-I spectra of the Perseus cluster, with a flux consistent with XMM-Newton (however, it is not seen in the ACIS-I spectrum of Virgo). The line is present even if we allow maximum freedom for all the known thermal emission lines. However, it is very weak (with an equivalent width in the full sample of only ∼ 1 eV) and located within 50-110 eV of several known faint lines; the detection is at the limit of the current instrument capabilities and subject to significant modeling uncertainties. On the origin of this line, we argue that there should be no atomic transitions in thermal plasma at this energy. An intriguing possibility is the decay of sterile neutrino, a long-sought dark matter particle candidate. Assuming that all dark matter is in sterile neutrinos with m s = 2E = 7.1 keV, our detection in the full sample corresponds to a neutrino decay mixing angle sin 2 (2θ) ≈ 7 × 10 −11 , below the previous upper limits. However, based on the cluster masses and distances, the line in Perseus is much brighter than expected in this model, significantly deviating from other subsamples. This appears to be because of an anomalously bright line at E = 3.62 keV in Perseus, which could be an Ar xvii dielectronic recombination line, although its emissivity would have to be 30 times the expected value and physically difficult to understand. In principle, such an anomaly might explain our line detection in other subsamples as well, though it would stretch the line energy uncertainties. Another alternative is the above anomaly in the Ar line combined with the nearby 3.51 keV K line also exceeding expectation by a factor 10-20. Confirmation with Chandra and Suzaku, and eventually Astro-H, are required to determine the nature of this new line. (APJ HAS THE ABRIDGED ABSTRACT) Γ γ (m s , θ) = 1.38 × 10 −29 s −1 sin 2 2θ 10 −7 m s 1 keV 5 ,(1) where the particle mass m s and the "mixing angle" θ

PSR J0740+6620 has a gravitational mass of 2.08 ± 0.07 M ⊙, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 13.7 − 1.5 + 2.6 km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451, the masses of two other ∼2 M ⊙ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at ∼1.5–5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4 M ⊙ neutron star is known to ±4% (68% credibility) and the radius of a 2.08 M ⊙ neutron star is known to ±5%. The full radius range that spans the ±1σ credible intervals of all the radius estimates in the three frameworks is 12.45 ± 0.65 km for a 1.4 M ⊙ neutron star and 12.35 ± 0.75 km for a 2.08 M ⊙ neutron star.

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be -+ 10 0.06 0.05 ( [ ])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples. Unified Astronomy Thesaurus concepts: Millisecond pulsars (1062); Rotation powered pulsars (1408); Pulsars (1306); Radio pulsars (1353); X-ray astronomy (1810); Neutron stars (1108)