No Arabic abstract
In this paper we show that X-ray spectral observations of the ATHENA mission, which is planned to launch in 2031, can constrain the equation of state of superdense matter. We use our well-constrained continuum fitting method for mass and radius determination of the neutron star. Model spectra of the emission from a neutron star were calculated using the atmosphere code ATM24. In the next step, those models were fitted to a simulated spectra of the neutron star calculated for ATHENAs WFI detector, using the satellite calibration files. To simulate the spectra we assumed three different values of effective temperatures, surface gravities and gravitational redshifts. There cases are related to the three different neutron star masses and radii. This analysis allows us to demonstrate the precision of our method and demonstrate the need for a fast detector onboard of ATHENA. A large grid of theoretical spectra was calculated with various parameters and a hydrogen-helium-iron composition of solar proportion. These spectra were fitted to the simulated spectrum to estimate the precision of mass and radius determination. In each case, we obtained very precise mass and radius values with errors in the range 3--10% for mass and in the range 2--8% for radius within the 1-sigma confidence error. We show here that with the ATHENA WFI detector, such a determination could be used to constrain the equation of state of superdense neutron star matter.
The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modelling pulse profiles of thermonuclear X-ray burst oscillations from hotspots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT, NICER, or LOFT. In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty, and to explore how the goodness-of-fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observers line-of-sight, with an accuracy of $le3$% compared to the true value and with $le5$% statistical uncertainty. The next best-determined are the mass and radius; for a neutron star with a spin frequency of 600,Hz, the best-fit mass and radius are accurate to $le5$%, with respective uncertainties of $le7$% and $le10$%. The accuracy and precision depend on the observer inclination and spot co-latitude, with values of $sim1$ % achievable in mass and radius if both the inclination and co-latitude are > 60 degrees.
The Rapid Burster (MXB 1730-335) is a unique object, showing both type I and type II X-ray bursts. A type I burst of the Rapid Burster was observed with Swift/XRT on 2009 March 5, showing photospheric radius expansion for the first time in this source. We report here on the mass and radius determination from this photospheric radius expansion burst using a Bayesian approach. After marginalization over the likely distance of the system (5.8-10 kpc) we obtain M=1.1+/-0.3 M_sun and R=9.6+/-1.5 km (1-sigma uncertainties) for the compact object, ruling out the stiffest equations of state for the neutron star. We study the sensitivity of the results to the distance, the color correction factor, and the hydrogen mass fraction in the envelope. We find that only the distance plays a crucial role.
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.
The dependence of the spin frequency derivative $dot{ u}$ of accreting neutron stars with a strong magnetic field (X-ray pulsars) on the mass accretion rate (bolometric luminosity, $L_{bol}$) has been investigated for eight transient pulsars in binary systems with Be stars. Using data from the Fermi/GBM and Swift/BAT telescopes, we have shown that for seven of the eight systems the dependence $dot{ u}$ can be fitted by the model of angular momentum transfer through an accretion disk, which predicts the relation $dot{ u}sim L^{6/7}_{bol}$. Hysteresis in the dependence $dot{ u}(L_{bol})$ has been confirmed in the system V 0332+53 and has been detected for the first time in the systems KS 1947+300, GRO J1008-57, and 1A 0535+26. The radius of the neutron star magnetosphere in all of the investigated systems have been estimated. We show that this quantity varies from pulsar to pulsar and depends strongly on the analytical model and the estimates for the neutron star and binary system parameters.
We analyze observations of eight quiescent low-mass X-ray binaries in globular clusters and combine them to determine the neutron star mass-radius curve and the equation of state of dense matter. We determine the effect that several uncertainties may have on our results, including uncertainties in the distance, the atmosphere composition, the neutron star maximum mass, the neutron star mass distribution, the possible presence of a hotspot on the neutron star surface, and the prior choice for the equation of state of dense matter. We find that the radius of a 1.4 solar mass neutron star is most likely from 10 to 14 km and that tighter constraints are only possible with stronger assumptions about the nature of the neutron stars, the systematics of the observations, or the nature of dense matter. Strong phase transitions are preferred over other models and interpretations of the data with a Bayes factor of 8 or more, and in this case, the radius is likely smaller than 12 km. However, radii larger than 12 km are preferred if the neutron stars have uneven temperature distributions.