No Arabic abstract
Observations of thermonuclear X-ray bursts from accreting neutron stars (NSs) in low-mass X-ray binary systems can be used to constrain NS masses and radii. Most previous work of this type has set these constraints using Planck function fits as a proxy: both the models and the data are fit with diluted blackbody functions to yield normalizations and temperatures which are then compared against each other. Here, for the first time, we fit atmosphere models of X-ray bursting NSs directly to the observed spectra. We present a hierarchical Bayesian fitting framework that uses state-of-the-art X-ray bursting NS atmosphere models with realistic opacities and relativistic exact Compton scattering kernels as a model for the surface emission. We test our approach against synthetic data, and find that for data that are well-described by our model we can obtain robust radius, mass, distance, and composition measurements. We then apply our technique to Rossi X-ray Timing Explorer observations of five hard-state X-ray bursts from 4U 1702-429. Our joint fit to all five bursts shows that the theoretical atmosphere models describe the data well but there are still some unmodeled features in the spectrum corresponding to a relative error of 1-5% of the energy flux. After marginalizing over this intrinsic scatter, we find that at 68% credibility the circumferential radius of the NS in 4U 1702-429 is R = 12.4+-0.4 km, the gravitational mass is M=1.9+-0.3 Msun, the distance is 5.1 < D/kpc < 6.2, and the hydrogen mass fraction is X < 0.09.
Determining neutron star (NS) radii and masses can help to understand the properties of matter at supra-nuclear densities. Thermal emission during thermonuclear X-ray bursts from NSs in low-mass X-ray binaries provides a unique opportunity to study NS parameters, because of the high fluxes, large luminosity variations and the related changes in the spectral properties. The standard cooling tail method uses hot NS atmosphere models to convert the observed spectral evolution during cooling stages of X-ray bursts to the Eddington flux F_Edd and the stellar angular size Omega. These are then translated to the constraints on the NS mass M and radius R. Here we present the improved, direct cooling tail method that generalises the standard approach. First, we adjust the cooling tail method to account for the bolometric correction to the flux. Then, we fit the observed dependence of the blackbody normalization on flux with a theoretical model directly on the M-R plane by interpolating theoretical dependences to a given gravity, hence ensuring only weakly informative priors for M and R instead of F_Edd and Omega. The direct cooling method is demonstrated using a photospheric radius expansion burst from SAX J1810.8--2609, which has happened when the system was in the hard state. Comparing to the standard cooling tail method, the confidence regions are shifted by 1 sigma towards larger radii, giving R=11.5-13.0 km at M=1.3-1.8 M_sun for this NS.
We introduce a new, powerful method to constrain properties of neutron stars (NSs). We show that the total mass of GW170817 provides a reliable constraint on the stellar radius if the merger did not result in a prompt collapse as suggested by the interpretation of associated electromagnetic emission. The radius R_1.6 of nonrotating NSs with a mass of 1.6 M_sun can be constrained to be larger than 10.68_{-0.04}^{+0.15} km, and the radius R_max of the nonrotating maximum mass configuration must be larger than 9.60_{-0.03}^{+0.14} km. We point out that detections of future events will further improve these constraints. Moreover, we show that a future event with a signature of a prompt collapse of the merger remnant will establish even stronger constraints on the NS radius from above and the maximum mass M_max of NSs from above. These constraints are particularly robust because they only require a measurement of the chirp mass and a distinction between prompt and delayed collapse of the merger remnant, which may be inferred from the electromagnetic signal or even from the presence/absence of a ringdown gravitational-wave (GW) signal. This prospect strengthens the case of our novel method of constraining NS properties, which is directly applicable to future GW events with accompanying electromagnetic counterpart observations. We emphasize that this procedure is a new way of constraining NS radii from GW detections independent of existing efforts to infer radius information from the late inspiral phase or postmerger oscillations, and it does not require particularly loud GW events.
Eddington-limited X-ray bursts from neutron stars can be used in conjunction with other spectroscopic observations to measure neutron star masses, radii, and distances. In order to quantify some of the uncertainties in the determination of the Eddington limit, we analysed a large sample of photospheric radius-expansion thermonuclear bursts observed with the Rossi X-ray Timing Explorer. We identified the instant at which the expanded photosphere touches down back onto the surface of the neutron star and compared the corresponding touchdown flux to the peak flux of each burst. We found that for the majority of sources, the ratio of these fluxes is smaller than 1.6, which is the maximum value expected from the changing gravitational redshift during the radius expansion episodes (for a 2M_sun neutron star). The only sources for which this ratio is larger than 1.6 are high inclination sources that include dippers and Cyg X-2. We discuss two possible geometric interpretations of this effect and show that the inferred masses and radii of neutron stars are not affected by this bias. On the other hand, systematic uncertainties as large as ~50% may be introduced to the distance determination.
Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a stars equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14} M_odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.
X-ray spectral analysis of quiescent low-mass X-ray binaries (LMXBs) has been one of the most common tools to measure the radius of neutron stars (NSs) for over a decade. So far, this method has been mainly applied to NSs in globular clusters, primarily because of their well-constrained distances. Here, we study Chandra data of seven transient LMXBs in the Galactic plane in quiescence to investigate the potential of constraining the radius (and mass) of the NSs inhabiting these systems. We find that only two of these objects had X-ray spectra of sufficient quality to obtain reasonable constraints on the radius, with the most stringent being an upper limit of $Rlesssim$14.5 km for EXO 0748-676 (for assumed ranges for mass and distance). Using these seven sources, we also investigate systematic biases on the mass/radius determination; for Aql X-1 we find that omitting a power-law spectral component when it does not seem to be required by the data, results in peculiar trends in the obtained radius with changing mass and distance. For EXO 0748-676 we find that a slight variation in the lower limit of the energy range chosen for the fit leads to systematically different masses and radii. Finally, we simulated Athena spectra and found that some of the biases can be lifted when higher quality spectra are available and that, in general, the search for constraints on the equation of state of ultra-dense matter via NS radius and mass measurements may receive a considerable boost in the future.