Search

We study the effects of the time-variable properties of thermonuclear X-ray bursts on modeling their millisecond-period burst oscillations. We apply the pulse profile modeling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer (NICER) to infer masses, radii, and geometric parameters of neutron stars. By simulating and analyzing a large set of models, we show that overlooking burst time-scale variability in temperatures and sizes of the hot emitting regions can result in substantial bias in the inferred mass and radius. To adequately infer neutron star properties, it is essential to develop a model for the time variable properties or invest a substantial amount of computational time in segmenting the data into non-varying pieces. We discuss prospects for constraints from proposed future X-ray telescopes., Comment: Accepted for publication in MNRAS

We present a pulse timing analysis of NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658 during the outburst that started on 2022 August 19. Similar to previous outbursts, after decaying from a peak luminosity of $\simeq 1\times10^{36} \, \mathrm{erg \, s^{-1}}$ in about a week, the pulsar entered in a $\sim 1$ month-long reflaring stage. Comparison of the average pulsar spin frequency during the outburst with those previously measured confirmed the long-term spin derivative of $\dot{\nu}_{\textrm{SD}}=-(1.15\pm0.06)\times 10^{-15} \, \mathrm{Hz\,s^{-1}}$, compatible with the spin-down torque of a $\approx 10^{26} \, \mathrm{G \, cm^3}$ rotating magnetic dipole. For the first time in the last twenty years, the orbital phase evolution shows evidence for a decrease of the orbital period. The long-term behaviour of the orbit is dominated by a $\sim 11 \, \mathrm{s}$ modulation of the orbital phase epoch consistent with a $\sim 21 \, \mathrm{yr}$ period. We discuss the observed evolution in terms of a coupling between the orbit and variations in the mass quadrupole of the companion star., Comment: 11 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

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 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620., 29 pages, 18 figures (10 of which are figure sets), 1 animation, 10 tables. Published in ApJ

We have simulated the X-ray polarization data that can be obtained with the Imaging X-ray Polarimetry Explorer, when observing accretion-powered millisecond pulsars. We estimated the necessary exposure times for SAX J1808.4$-$3658 in order to obtain different accuracies in the measured time-dependent Stokes profiles integrated over all energy channels. We found that the measured relative errors depend strongly on the relative configuration of the observer and the emitting hotspot. The improvement in the minimum relative error in Stokes $Q$ and $U$ parameters as a function of observing time $t$ scales as $1/\sqrt{t}$, and spans the range from 30-90% with 200 ks exposure time to 20-60% with 500 ks exposure time (in case of data binned in 19 phase bins). The simulated data were also used to predict how accurate measurements of the geometrical parameters of the neutron star can be made when modelling only $Q$ and $U$ parameters, but not the flux. We found that the observer inclination and the hotspot co-latitude could be determined with better than 10 deg accuracy for most of the cases we considered. These measurements can be used to further constrain neutron star mass and radius when combined with modelling of the X-ray pulse profile., Comment: 12 pages, 11 figures, published in A&A

We have developed a complete theory for the calculation of the observed Stokes parameters for radiation emitted from the surface of a rapidly rotating neutron star (NS) using the oblate Schwarzschild approximation. We accounted for the rotation of the polarization plane due to relativistic effects along the path from the stellar surface to the observer. The results were shown to agree with those obtained by performing full numerical general relativistic ray-tracing with the \textsc{arcmancer} code. We showed that the obtained polarization angle (PA) profiles may differ substantially from those derived for a spherical star. We demonstrated that assuming incorrect shape for the star can lead to biased constraints for NS parameters when fitting the polarization data. Using a simplified model, we also made a rough estimate of how accurately the geometrical parameters of an accreting NS can be determined using the X-ray polarization measurements of upcoming polarimeters like the Imaging X-ray Polarimeter Explorer (IXPE) or the enhanced X-ray Timing and Polarimetry (eXTP) mission., 11 pages, 10 figures, accepted in A&A

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

The aim of this work is to study the possible effects and biases on the radius constraints for rotation-powered millisecond pulsars when using Thomson approximation to describe electron scattering in the atmosphere models, instead of using exact formulation for Compton scattering. We compare the differences between the two models in the energy spectrum and angular distribution of the emitted radiation. We also analyse a self-generated synthetic phase-resolved energy spectrum, based on Compton atmosphere and the most X-ray luminous rotation-powered millisecond pulsars observed by the Neutron star Interior Composition ExploreR (NICER). We derive constraints for the neutron star parameters using both the Compton and Thomson models. The results show that the method works by reproducing the correct parameters with the Compton model. However, biases are found in size and the temperature of the emitting hot spot, when using the Thomson model. The constraints on the radius are still not significantly changed, and therefore the Thomson model seems to be adequate if we are interested only in the radius measurements using NICER., Comment: 6 pages, 9 figures, published in A&A

We computed accurate atmosphere models of rotation-powered millisecond pulsars in which the polar caps of a neutron star (NS) are externally heated by magnetospheric return currents. The external ram pressure, energy losses, and stopping depth of the penetrating charged particles were computed self-consistently with the atmosphere model, instead of assuming a simplified deep-heated atmosphere in radiative equilibrium. We used exact Compton scattering formalism to model the properties of the emergent X-ray radiation. The deep-heating approximation was found to be valid only if most of the heat originates from ultra-relativistic bombarding particles with Lorentz factors of $\gamma \gtrsim 100$. In the opposite regime, the atmosphere attains a distinct two-layer structure with an overheated optically thin skin on top of an optically thick cool plasma. The overheated skin strongly modifies the emergent radiation: it produces a Compton-upscattered high-energy tail in the spectrum and alters the radiation beaming pattern from limb darkening to limb brightening for emitted hard X-rays. This kind of drastic change in the emission properties can have a significant impact on the inferred NS pulse profile parameters as performed, for example, by Neutron star Interior Composition ExploreR. Finally, the connection between the energy distribution of the return current particles and the atmosphere emission properties offers a new tool to probe the exact physics of pulsar magnetospheres., Comment: 13 pages, 10 figures, published in A&A

Thermonuclear X-ray bursts on the surface of neutron stars (NSs) can enrich photosphere with metals, which may imprint the photoionization edges on the burst spectra. We report here the discovery of absorption edges in the spectra of the type I X-ray burst from the NS low-mass X-ray binary GRS 1747$-$312 in Terzan 6 during observations by the Rossi X-ray Timing Explorer. We find that the edge energy evolves from $9.45\pm0.51$ keV to $\sim6$ keV and then back to $9.44\pm0.40$ keV during the photospheric radius expansion phase and remains at $8.06\pm0.66$ keV in the cooling tail. The photoionization absorption edges of hydrogen-like Ni, Fe or Fe/Ni mixture and the bound-bound transitions of metals may be responsible for the observed spectral features. The ratio of the measured absorption edge energy in the cooling tail to the laboratory value of hydrogen-like Ni(Fe) edge energy allows us to estimate the gravitational redshift factor $1+z=1.34\pm0.11$($1+z=1.15\pm0.09$). The evolution of the spectral parameters during the cooling tail are well-described by the metal-rich atmosphere models. The combined constraints on the NS mass and radius from the direct cooling method and the tidal deformability strongly suggest very high atmospheric abundance of the iron group elements and limit the distance to the source to $11\pm1$ kpc., matches the published version