No Arabic abstract
The Neutron star Interior Composition ExploreR (NICER) is expected to launch in early 2017 and will gather X-ray data on neutron stars and other high-energy sources from a berth on the International Space Station. Its prime scientific goal is to measure the masses and radii of non-accreting neutron stars via fits to the energy-dependent waveforms produced by the rotation of hot spots on their surfaces. These measurements will provide valuable input to theoretical models of cold matter beyond nuclear density. Here we propose that PSR J1614$-$2230, despite its low count rate, is a promising source to observe with NICER. The reason is that XMM-Newton observations suggest that the fractional oscillation amplitude from PSR J1614$-$2230 could be high enough that this star cannot be very compact. We show that if we analyze 0.5 Ms of NICER data and 0.1 Ms of nearby off-source data and combine that analysis with the known mass of this star, we would find a robust lower limit to the radius with a statistical uncertainty of only $sim 0.5-0.7$ km. We also show that even if there is an unmodeled nonthermal component modulated at the pulsation frequency, good statistical fits could rule out significant biases. The low count rate will make reliable upper limits on the radius difficult, but the lower limit could rule out some equations of state that are currently being discussed. This analysis would require a good estimate of the non-source background, so Chandra observations of the vicinity of PSR J1614$-$2230 would be helpful.
We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X-ray fields. PSR J0740+6620 exhibits a broad double-peaked profile with a separation of ~0.4 in phase. PSR J1614-2230, on the other hand, has a broad single-peak profile. The broad modulations with soft X-ray spectra of both pulsars are indicative of thermal radiation from one or more small regions of the stellar surface. We show the NICER detections of X-ray pulsations for both pulsars and also discuss the phase relationship to their radio pulsations. In the case of PSR J0740+6620, this paper documents the data reduction performed to obtain the pulsation detection and prepare for pulse profile modeling analysis.
Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASAs NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.
We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsars surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16%$ and $84%$ quantiles in marginal posterior mass.
PSR J0740$+$6620 has a gravitational mass of $2.08pm 0.07~M_odot$, 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^{+2.6}_{-1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030$+$0451, the masses of two other $sim 2~M_odot$ 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 $sim 1.5-3$ times nuclear saturation density. For a given framework, when all measurements are included the radius of a $1.4~M_odot$ neutron star is known to $pm 4$% (68% credibility) and the radius of a $2.08~M_odot$ neutron star is known to $pm 5$%. The full radius range that spans the $pm 1sigma$ credible intervals of all the radius estimates in the three frameworks is $12.45pm 0.65$ km for a $1.4~M_odot$ neutron star and $12.35pm 0.75$ km for a $2.08~M_odot$ neutron star.
PSR J0537-6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 years of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire timespan, which overlaps with the third observing run of LIGO/Virgo, thus enabling the most sensitive gravitational wave searches of this potentially strong gravitational wave-emitting pulsar. We find that the short-term braking index between glitches decreases towards a value of 7 or lower at longer times since the preceding glitch. By combining NICER and RXTE data, we measure a long-term braking index n=-1.25+/-0.01. Our analysis reveals 8 new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88x10^-7 yr^-1. The new glitches follow the seemingly unique time-to-next-glitch---glitch-size correlation established previously using RXTE data, with a slope of 5 d microHz^-1. For one glitch around which NICER observes two days on either side, we search for but do not see clear evidence of spectral nor pulse profile changes that may be associated with the glitch.