Sharon M. Morsink

If dark matter (DM) exists in halos around rotating neutron stars (NSs), it will be essential to understand the effects of rotation on the distribution of DM and baryonic matter (BM) in the stars to interpret observations. In this work, we construct rapidly rotating dark matter admixed neutron stars (DANS) with DM halos using the two-fluid approximation, where the BM and DM interact only through gravity. Our goal is to describe rapidly rotating millisecond-period DANS spun up by the accretion of BM from a zero angular momentum state. We extend the Rapidly Rotating Neutron Star ( RNS ) code to compute axisymmetric configurations in which the BM rotates rigidly while the DM remains torque-free and differentially rotates through the frame-dragging of spacetime. For the first time, we examine in detail local and global definitions of mass in general relativity for two-fluid systems, showing how their differences affect the interpretation of baryonic and dark component masses. We compute energy density and frame-dragging frequency profiles for DANS with three different characteristic DM halos. We demonstrate that rapid BM rotation reduces DM halo sizes if central energy densities are kept constant between non-rotating and rotating models. We also construct sequences of DANS to create mass and radius curves and compare rotating and non-rotating cases. Finally, we quantify deviations in the spacetime metric outside the baryonic surfaces of these sequences of stars caused by the DM halos. We hypothesize that the size of this quantity could indicate whether a DM halo will significantly impact X-ray pulse profile modeling. These results provide a framework for assessing the observational consequences of DM halos around rapidly rotating NSs.

Neutron Star Interior Composition Explorer (NICER) data have been used to estimate the masses and radii of the rotation-powered millisecond pulsars PSR J0030+0451, PSR J0740+6620, PSR J0437−4715, PSR J1231−1411, and PSR J0614−3329, sometimes in joint analyses with XMM-Newton data. These measurements provide invaluable information about the properties of cold, catalyzed matter beyond nuclear saturation density. Here, we present the results of our modeling of NICER data on PSR J0437−4715 using several different models of hot thermal X-ray-emitting spots on the stellar surface. For this pulsar, previous Nuclear Spectroscopic Telescope Array observations established that there is also a modulated nonthermal component to the emission, but the previously published analysis of NICER data did not model this component. We find that the Bayesian evidence is significantly higher when the modulated nonthermal component is included and that omission of this component leads to poorer fits to the bolometric NICER data and thus risks bias in the resulting radius estimates. Our models, which we pursue to inferential convergence, therefore have modulated nonthermal emission, and our headline model has in addition three uniform-temperature thermally emitting circular spots. Using this model, the symmetric 68% credible range in the radius is 11.9–15.5 km, which at the independently measured mass of M = 1.418 ± 0.044 M _⊙ is consistent with previous reports of the radius of the ∼1.4 M _⊙ pulsar PSR J0030+0451. We discuss the implications of this measurement for the equation of state of dense matter.

The oblate Schwarzschild (OS) approximation is a method often used to compute the flux of X-rays emitted from a rapidly rotating neutron star. In this approximation, the oblate shape of the rotating star is embedded in the Schwarzschild metric, which is used to compute the redshift of photon energies as they propagate from the star to the telescope. In this paper, we demonstrate that there are small errors introduced by the standard treatment of photon redshift in the OS approximation and provide a simple method to correct these errors. These errors are constant in phase, so this results in a constant absolute reduction in the flux. For PSR J0740+6620, the most rapidly spinning of the pulsars observed by the Neutron star Interior Composition Explorer (NICER), we estimate the flux errors are less than 1%, which is an order of magnitude smaller than the uncertainty in the distance, so this does not affect the mass and radius constraints found for this pulsar. The errors for the other pulsars observed by NICER are even smaller. However, this correction should be included when analyzing data for more rapidly rotating X-ray pulsars with spin frequencies near 600 Hz.

Recent constraints on neutron star mass and radius have advanced our understanding of the equation of state (EOS) of cold dense matter. Some of them have been obtained by modeling the pulses of three millisecond X-ray pulsars observed by the Neutron Star Interior Composition Explorer (NICER). Here, we present a Bayesian parameter inference for a fourth pulsar, PSR J1231−1411, using the same technique with NICER and XMM-Newton data. When applying a broad mass-inclination prior from radio timing measurements and the emission region geometry model that can best explain the data, we find likely converged results only when using a limited radius prior. If limiting the radius to be consistent with the previous observational constraints and EOS analyses, we infer the radius to be 12.6 ± 0.3 km and the mass to be 1.04 − 0.03 + 0.05 M _⊙, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. If using an uninformative prior but limited between 10 and 14 km, we find otherwise similar results, but R eq = 13.5 − 0.5 + 0.3 km for the radius. In both cases, we find a nonantipodal hot region geometry where one emitting spot is at the equator or slightly above, surrounded by a large colder region, and where a noncircular hot region lies close to southern rotational pole. If using a wider radius prior, we only find solutions that fit the data significantly worse. We discuss the challenges in finding the better fitting solutions, possibly related to the weak interpulse feature in the pulse profile.

Pulse profile modeling (PPM), the technique used to infer mass, radius, and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross tests for simple hot spot geometries, focusing primarily on the implementation of the spacetime model. In this paper, we present verification for test problems that explore the more complex hot spot geometries that are now being employed in the NICER PPM analyses. We conclude that the accuracy of our computed waveforms is in general sufficiently high for analyses of current NICER data sets. We have however identified some extreme configurations where extra care may be needed.

The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using Neutron Star Interior Composition Explorer (NICER) data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio-timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gravitational mass are 12.49 − 0.88 + 1.28 km and 2.073 − 0.069 + 0.069 M _⊙, respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, with an estimated systematic error ≲ 0.1 km. This result was obtained using the best computationally feasible sampler settings providing a strong radius lower limit but a slightly more uncertain radius upper limit. The inferred radius interval is also close to the R = 12.76 − 1.02 + 1.49 km obtained by Dittmann et al., when they require the radius to be less than 16 km as we do. The results continue to disfavor very soft equations of state for dense matter, with R < 11.15 km for this high-mass pulsar excluded at the 95% probability. The results do not depend significantly on the assumed cross-calibration uncertainty between NICER and XMM-Newton. Using simulated data that resemble the actual observations, we also show that our pipeline is capable of recovering parameters for the inferred models reported in this paper.

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be 2.08 ± 0.07 M _⊙. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. and Riley et al. reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 2020 April 17, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 2021 December 28 was presented in Salmi et al. Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 April 21, totaling ∼1.1 Ms additional exposure compared to the data set analyzed in Miller et al. and Riley et al., and to data from XMM-Newton observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 12.92 − 1.13 + 2.09 km (68% credibility), a fractional uncertainty ∼83% the width of that reported in Miller et al., in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al., then our 68% credible region would become R = 12.76 − 1.02 + 1.49 km, which is close to the headline result of Salmi et al. Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

We report Bayesian inference of the mass, radius, and hot X-ray emitting region properties—using data from the Neutron Star Interior Composition ExploreR (NICER)—for the brightest rotation-powered millisecond X-ray pulsar, PSR J0437−4715. Our modeling is conditional on informative tight priors on mass, distance, and binary inclination obtained from radio pulsar timing using the Parkes Pulsar Timing Array (PPTA; Reardon et al.), and we use NICER background models to constrain the nonsource background, cross-checking with data from XMM-Newton. We assume two distinct hot emitting regions and various parameterized hot region geometries that are defined in terms of overlapping circles; while simplified, these capture many of the possibilities suggested by detailed modeling of return current heating. For the preferred model identified by our analysis, we infer a mass of M = 1.418 ± 0.037 M _⊙ (largely informed by the PPTA mass prior) and an equatorial radius of R = 11.36 − 0.63 + 0.95 km, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. This radius favors softer dense matter equations of state and is highly consistent with constraints derived from gravitational wave measurements of neutron star binary mergers. The hot regions are inferred to be nonantipodal and hence inconsistent with a pure centered dipole magnetic field.