Engineering topics

Morsink, S. M.

Publications and source records attributed to Morsink, S. M..

5 records

DOE OSTI2021

The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

PSR J0740+6620 has a gravitational mass of 2.08 ± 0.07 M ⊙ , which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 13.7$^{+2.6}_{–1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451, the masses of two other ~2 M ⊙ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at ~1.5-5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4 M ⊙ neutron star is known to ±4% (68% credibility) and the radius of a 2.08 M ⊙ neutron star is known to ±5%. The full radius range that spans the ±1σ credible intervals of all the radius estimates in the three frameworks is 12.45 ± 0.65 km for a 1.4 M ⊙ neutron star and 12.35 ± 0.75 km for a 2.08 M ⊙ neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation

Here, we report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius R eq are, respectively, ${1.34}_{-0.16}^{+0.15}\,{M}_{\odot }$ and ${12.71}_{-1.19}^{+1.14}\,\mathrm{km}$, while the compactness ${GM}/{R}_{\mathrm{eq}}{c}^{2}={0.156}_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

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 (EoS) of this dense matter is to measure both a star's 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. In this work, 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.06}^{+1.24}$ km and $M={1.44}_{-0.14}^{+0.15}\,{M}_{\odot }$ (68%). The independent analysis reported in the companion paper by Riley et al. 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 EoS of cold, catalyzed matter above nuclear saturation density.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

NEUTRON STAR MASS–RADIUS CONSTRAINTS USING EVOLUTIONARY OPTIMIZATION

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 modeling pulse profiles of thermonuclear X-ray burst oscillations from hot spots 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 a mission similar to 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 observer’s line of sight, with an accuracy of ≤3% compared to the true value and with ≤5% 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 ≤5%, with respective uncertainties of ≤7% and ≤10%. The accuracy and precision depend on the observer inclination and spot colatitude, with values of ∼1% achievable in mass and radius if both the inclination and colatitude are ≳60°.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE IMPACT OF SURFACE TEMPERATURE INHOMOGENEITIES ON QUIESCENT NEUTRON STAR RADIUS MEASUREMENTS

Fitting the thermal X-ray spectra of neutron stars (NSs) in quiescent X-ray binaries can constrain the masses and radii of NSs. The effect of undetected hot spots on the spectrum, and thus on the inferred NS mass and radius, has not yet been explored for appropriate atmospheres and spectra. A hot spot would harden the observed spectrum, so that spectral modeling tends to infer radii that are too small. However, a hot spot may also produce detectable pulsations. We simulated the effects of a hot spot on the pulsed fraction and spectrum of the quiescent NSs X5 and X7 in the globular cluster 47 Tucanae, using appropriate spectra and beaming for hydrogen atmosphere models, incorporating special and general relativistic effects, and sampling a range of system angles. We searched for pulsations in archival Chandra HRC-S observations of X5 and X7, placing 90% confidence upper limits on their pulsed fractions below 16%. We use these pulsation limits to constrain the temperature differential of any hot spots, and to then constrain the effects of possible hot spots on the X-ray spectrum and the inferred radius from spectral fitting. We find that hot spots below our pulsation limit could bias the spectroscopically inferred radius downward by up to 28%. For Cen X-4 (which has deeper published pulsation searches), an undetected hot spot could bias its inferred radius downward by up to 10%. Improving constraints on pulsations from quiescent LMXBs may be essential for progress in constraining their radii.

79 ASTRONOMY AND ASTROPHYSICS↗