An Investigation Into Possible Systematic Effects on Neutron Star Radius Estimates using NICER-like Synthetic Data
Neutron star cores contain the densest matter in the observable universe. The state of this matter is of interest in numerous fields, but laboratory experiments cannot explore this matter. Although the composition of the matter would be of great interest, macroscopic observables such as the neutron star mass-radius relation depend primarily on the equation of state (EOS). As a result, precise and reliable radius measurements would be valuable in constraining the EOS. However, most attempts at radius measurements are susceptible to systematic errors, meaning the inferred radius can be significantly biased even though the fit to the data appears to be statistically good. Previous studies suggested that radii inferred using X-ray data provided by NASA’s NICER mission may be more immune from such systematic errors. This is in part because, compared with previous measurements that obtained averaged spectra and fluxes, NICER observes millisecond pulsars by timing the photons so precisely, it is possible to obtain the spectrum as a function of rotational phase and see variations such as heated regions on the star rotate into and out of view hundreds of times per second. This extra information seems promising to break degeneracies and mitigate systematic errors, but a more in-depth study is necessary. We report the first steps of that study, in which we generate NICER-like synthetic data and determine the quality of fit and bias in radius obtained when we fit the data using a model different from the model used to generate the data.