ABSTRACT We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a ...large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of neutron stars, we place strong and quantitative constraints on the properties of the equation of state between times the nuclear saturation density. We find that around the preferred equation of state predicts radii between 10.1 and 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential.
ABSTRACT
We analysed the thermal emission from the entire surface of the millisecond pulsar PSR J0437−4715 observed in the ultraviolet and soft X-ray bands. For this, we calculated non-magnetized, ...partially ionized atmosphere models of hydrogen, helium, and iron compositions and included plasma frequency effects that may affect the emergent spectrum. This is particularly true for the coldest atmospheres composed of iron (up to a few per cent changes in the soft X-ray flux). Employing a Markov chain Monte Carlo method, we found that the spectral fits favour a hydrogen atmosphere, disfavour a helium composition, and rule out iron atmosphere and blackbody models. By using a Gaussian prior on the dust extinction, based on the latest 3D map of Galactic dust, and accounting for the presence of hot polar caps found in the previous work, we found that the hydrogen atmosphere model results in a well-constrained neutron star radius ${R_{\rm NS}}= 13.6^{+0.9}_{-0.8}{\, {\rm km}}$ and bulk surface temperature ${T_{\rm eff}^{\infty }}=\left(2.3\pm 0.1\right){\times 10^{5}}{\, {\rm K}}$. This relatively large radius favours a stiff equation of state and disfavours a strange quark composition inside neutron stars.
Abstract We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data ...by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25° with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.
We report on NICER observations of the magnetar SGR 1935+2154, covering its 2020 burst storm and long-term persistent emission evolution up to ∼90 days postoutburst. During the first 1120 s taken on ...April 28 00:40:58 UTC, we detect over 217 bursts, corresponding to a burst rate of >0.2 bursts s−1. Three hours later, the rate was 0.008 bursts s−1, remaining at a comparatively low level thereafter. The T90 burst duration distribution peaks at 840 ms; the distribution of waiting times to the next burst is fit with a lognormal with an average of 2.1 s. The 1-10 keV burst spectra are well fit by a blackbody, with an average temperature and area of kT = 1.7 keV and R2 = 53 km2. The differential burst fluence distribution over ∼3 orders of magnitude is well modeled with a power-law form dN/dF ∝ F−1.5 0.1. The source persistent emission pulse profile is double-peaked hours after the burst storm. We find that the burst peak arrival times follow a uniform distribution in pulse phase, though the fast radio burst associated with the source aligns in phase with the brighter peak. We measure the source spin-down from heavy-cadence observations covering days 21-39 postoutburst, Hz s−1, a factor of 2.7 larger than the value measured after the 2014 outburst. Finally, the persistent emission flux and blackbody temperature decrease rapidly in the early stages of the outburst, reaching quiescence 40 days later, while the size of the emitting area remains unchanged.
Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late 2017 September in a giant outburst with a peak luminosity of 2 × 1039(d/7 kpc)2 erg s−1 (0.1-10 keV), with no ...formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER) and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L ∼ 1038 erg s−1, the pulse profiles transition from single peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turnover to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ∼ 1013 G.
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 yr ...of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire time span, 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 eight new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88 × 10−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 μHz−1. For one glitch around which NICER observes 2 d 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.
Abstract
In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state ...characterizing dense matter. Two independent analyses found a mass of ∼1.3–1.4
M
⊙
and a radius of ∼13 km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyze, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and also jointly analyze XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multimodal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of ∼1.4
M
⊙
, 11.5 km and ∼1.7
M
⊙
, 14.5 km, depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them.
The Neutron Star Interior Composition Explorer (NICER) has extensively monitored the 2019 August outburst of the 401 Hz millisecond X-ray pulsar SAX J1808.4-3658. In this Letter, we report on the ...detection of a bright helium-fueled Type I X-ray burst. With a bolometric peak flux of (2.3 0.1) × 10−7 erg s−1 cm−2, this was the brightest X-ray burst among all bursting sources observed with NICER to date. The burst shows a remarkable two-stage evolution in flux, emission lines at 1.0 and 6.7 keV, and burst oscillations at the known pulsar spin frequency, with 4% fractional sinusoidal amplitude. We interpret the burst flux evolution as the detection of the local Eddington limits associated with the hydrogen and helium layers of the neutron star envelope. The emission lines are likely associated with Fe, due to reprocessing of the burst emission in the accretion disk.
The Neutron star Interior Composition Explorer observed several rotation-powered millisecond pulsars (MSPs) to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these ...pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections (≥4.7 after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614−3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241−5236, while PSR J1552+5437 and PSR J1744−1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7 and almost 3 confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected MSPs exhibit a significant phase offset between their X-ray and radio pulses.
In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the ...Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, allowing for pulsar timing based astrometry of this accreting millisecond pulsar.