Showing papers by "Wynn C. G. Ho published in 2019"

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

Abstract: 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. 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 km and (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.

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

Abstract: 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, and , while the compactness is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

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

Abstract: 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 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. 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.

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

Abstract: 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 (NICER) 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_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst 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.

A Nicer View of PSR J0030+0451: Implications for the Dense Matter Equation of State

Abstract: Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's Neutron Star Interior Composition Explorer (NICER) mission. In this Letter we study the implications of the mass–radius inference reported for this source by Riley et al. for the dense matter equation of state (EoS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EoS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EoS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EoS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass–radius likelihood function derived by Riley et al. using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

A NICER view of PSR J0030+0451: Implications for the dense matter equation of state

Abstract: Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-Ray Data Set

Abstract: We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437−4715, J0030+0451, J1231−1411, and J2124−3358, selected as targets for constraining the mass–radius relation of neutron stars and the dense matter equation of state (EoS) via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass–radius relation and the dense matter EoS. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437−4715, for which multiple temperatures are required.

A $NICER$ View of PSR J0030+0451: Evidence for a Global-scale Multipolar Magnetic Field

Abstract: Recent modeling of Neutron Star Interior Composition Explorer observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar’s surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451’s magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

Abstract: GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

Abstract: We present a search for gravitational waves from 222 pulsars with rotation frequencies ≳10 Hz. We use advanced LIGO data from its first and second observing runs spanning 2015–2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the l = m = 2 mass quadrupole mode, with a frequency at twice that of the pulsar’s rotation, and the l = 2, m = 1 mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the l = m = 2 mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars’ spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711−6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of 1038 kg m2 for the star’s moment of inertia, and imply a gravitational-wave-derived upper limit on the star’s ellipticity of 1.2 × 10−8. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

NICER X-Ray Observations of Seven Nearby Rotation-powered Millisecond Pulsars

Abstract: The Neutron star Interior Composition Explorer observed several rotation-powered millisecond pulsars (MSPs) to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections (≥4.7σ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614−3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241−5236, while PSR J1552+5437 and PSR J1744−1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7σ and almost 3σ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass–radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected MSPs exhibit a significant phase offset between their X-ray and radio pulses.

Diffusive nuclear burning in cooling simulations and application to new temperature data of the Cassiopeia A neutron star

Abstract: A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns yH ≳ 107 g cm-2, while DNB can make helium envelopes more transparent to heat flux for surface temperatures Ts ≳ 2x106 K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the long-term cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of ~2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.

Anti-glitches in the Ultraluminous Accreting Pulsar NGC 300 ULX-1 Observed with NICER

Abstract: We report the discovery of three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). We find evidence for three sudden spin-down events of magnitudes $\\Delta\ u = -23, -30,$ and $-43 \\,\\mu$Hz (fractional amplitudes $\\Delta\ u/\ u = -4.4, -5.5,$ and $-7.7 \\times 10^{-4}$). These are larger in magnitude (and opposite in sign) than any radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfliud components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse.

A NICER view of PSR J0030+0451: evidence for a global-scale multipolar magnetic field

Abstract: Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-ray Data Set.

Abstract: We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

A systematic study of soft X-ray pulse profiles of magnetars in quiescence

Abstract: Magnetars are neutron stars with extremely high surface magnetic fields. They show diverse X-ray pulse profiles in the quiescent state. We perform a systematic Fourier analysis of their soft X-ray pulse profiles. We find that most magnetars have a single-peaked profile and hence have low amplitudes of the second Fourier harmonics (A2). On the other hand, the pulsed fraction (PF) spreads over a wide range. We compared the results with theoretical profiles assuming various surface hotspot asymmetries, viewing geometries, and beaming functions. We found that a single value of the intensity ratio r between two antipodal hotspots is unable to reproduce the observed distribution of A2 and PF for all magnetars. The inferred r is probably anticorrelated with the thermal luminosity, implying that high-luminosity magnetars tend to have two symmetric hotspots. Our results are consistent with theoretical predictions, for which the existence of an evolving toroidal magnetic field breaks the symmetry of the surface temperature.

A joint NICER and XMM-Newton view of the “Magnificent” thermally emitting X-ray Isolated Neutron Star RX J1605.3+3249

Abstract: Thermally emitting X-ray isolated neutron stars represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer physical parameters of the neutron star. Among the seven known members of this class, RX J1605.3+3249 is the only one that still lacks confirmation of its spin period. Here we analyze NICER and XMM-Newton observations of RX J1605.3+3249, in order to address its timing and spectral behavior. Contrary to a previous tentative detection, but in agreement with the recent work by Pires et al. (2019), we find no significant pulsation with pulsed fraction higher than 1.3% (3{\sigma}) for periods above 150 ms. We also find a limit of 2.6% for periods above 2 ms, despite searches in different energy bands. The X-ray spectrum can be fit by either a double-blackbody model or by a single-temperature magnetized atmosphere model, both modified by a Gaussian absorption line at ~0.44 keV. The origin of the absorption feature as a proton cyclotron line or as an atomic transition in the neutron star atmosphere is discussed. The predictions of the best-fit X-ray models extended to IR, optical and UV bands are compared with archival data. Our results are interpreted in the framework of a fallback disk scenario.

Targeted search for young radio pulsars in the SMC : discovery of two new pulsars

Abstract: We report the first rotation powered pulsars discovered in the Small Magellanic Cloud (SMC) in more than a decade. PSR J0043-73 and PSR J0052-72 were discovered during a Parkes Multi-Beam (PMB) survey of MCSNR J0127-7332, and five new, optically selected, supernova remnant (SNR) candidates identified by the XMM Newton survey. In addition to the candidates, we adjusted the PMB rotation to include an additional nine SNRs and pulsar wind nebulae. We searched for young pulsars (1 - 200 ms) employing a Fourier analysis with PRESTO, as well as a search for longer period pulsars (200 ms - 360 s) with a fast folding algorithm. Our targeted survey had a limiting flux density of 0.039 mJy for periods greater than 50 ms. Although not the main target of this search it was also sensitive to millisecond pulsars. PSR J0043-73 has a period and dispersion measure of 937.42937 (26) ms and 115.1 (3.4) pc cm$^{-3}$ respectively, and PSR J0052-72 has a period of 191.444328 (46) ms and a DM of 158.6 (1.6) pc cm$^{-3}$.

NICER X-ray Observations of Seven Nearby Rotation-Powered Millisecond Pulsars

Abstract: NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

Modelling decretion discs in Be/X-ray binaries

Abstract: As the largest population of high mass X-ray binaries, Be/X-ray binaries provide an excellent laboratory to investigate the extreme physics of neutron stars. It is generally accepted that Be stars possess a circumstellar disc, providing an additional source of accretion to the stellar winds present around young hot stars. Interaction between the neutron star and the disc is often the dominant accretion mechanism. A large amount of work has gone into modelling the properties of these circumstellar discs, allowing for the explanation of a number of observable phenomena. In this paper, smoothed particle hydrodynamics simulations are performed whilst varying the model parameters (orbital period, eccentricity, the mass ejection rate of the Be star and the viscosity and orientation of the disc). The relationships between the model parameters and the disc's characteristics (base gas density, the accretion rate of the neutron star and the disc's size) are presented. The observational evidence for a dependency of the size of the Be star's circumstellar disc on the orbital period (and semi-major axis) is support ed by the simulations.

Neutron Star Interior Composition Explorer X-ray Timing of the Radio and Gamma-ray Quiet Pulsars PSR J1412+7922 AND PSR J1849-0001

Abstract: We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a power law. Our refined timing measurements enable deeper searches for pulsations at other wavelengths and sensitive targeted searches by LIGO/Virgo for continuous gravitational waves from these neutron stars. Using the sensitivity of LIGO's first observing run, we estimate constraints that a gravitational wave search of these pulsars would obtain on the size of their mass deformation and r-mode fluid oscillation.

Soft excess in the quiescent Be/X-ray pulsar RX J0812.4–3114

Abstract: We report a 72 ks XMM-Newton observation of the Be/X-ray pulsar (BeXRP) RX J0812.4-3114 in quiescence ($L_X \approx 1.6 \times 10^{33}~\mathrm{erg~s^{-1}}$). Intriguingly, we find a two component spectrum, with a hard power-law ($\Gamma \approx 1.5$) and a soft blackbody-like excess below $\approx 1~\mathrm{keV}$. The blackbody component is consistent in $kT$ with a prior quiescent Chandra observation reported by Tsygankov et al. and has an inferred blackbody radius of $\approx 10~\mathrm{km}$, consistent with emission from the entire neutron star (NS) surface. There is also mild evidence for an absorption line at $\approx 1~\mathrm{keV}$ and/or $\approx 1.4~\mathrm{keV}$. The hard component shows pulsations at $P \approx 31.908~\mathrm{s}$ (pulsed fraction $0.84 \pm 0.10$), agreeing with the pulse period seen previously in outbursts, but no pulsations were found in the soft excess (pulsed fraction $\lesssim 31\%$). We conclude that the pulsed hard component suggests low-level accretion onto the neutron star poles, while the soft excess seems to originate from the entire NS surface. We speculate that, in quiescence, the source switches between a soft thermal-dominated state (when the propeller effect is at work) and a relatively hard state with low-level accretion, and use the propeller cutoff to estimate the magnetic field of the system to be $\lesssim 8.4 \times 10^{11}~\mathrm{G}$. We compare the quiescent thermal $L_X$ predicted by the standard deep crustal heating model to our observations and find that RX J0812.4-3114 has a high thermal $L_X$, at or above the prediction for minimum cooling mechanisms. This suggests that RX J0812.4-3114 either contains a relatively low-mass NS with minimum cooling, or that the system may be young enough that the NS has not fully cooled from the supernova explosion.

XMM-Newton detection and spectrum of the second fastest spinning pulsar PSR J0952-0607

Abstract: With a spin frequency of 707 Hz, PSR J0952-0607 is the second fastest spinning pulsar known. It was discovered in radio by LOFAR in 2017 at an estimated distance of either 0.97 or 1.74 kpc and has a low-mass companion with a 6.42 hr orbital period. We report discovery of the X-ray counterpart of PSR J0952-0607 using XMM-Newton. The X-ray spectra can be well-fit by a single power law model (Gamma = 2.5) or by a thermal plus power law model (kTeff = 40 eV and Gamma = 1.4). We do not detect evidence of variability, such as that due to orbital modulation from pulsar wind and companion star interaction. Because of its fast spin rate, PSR J0952-0607 is a crucial source for understanding the r-mode instability, which can be an effective mechanism for producing gravitational waves. Using the high end of our measured surface temperature, we infer a neutron star core temperature of ~10^7 K, which places PSR J0952-0607 within the window for the r-mode to be unstable unless an effect such as superfluid mutual friction damps the fluid oscillation. The measured luminosity limits the dimensionless r-mode amplitude to be less than ~1x10^-9.

Neutron Star Interior Composition Explorer X-ray Timing of the Radio and $\gamma$-ray Quiet Pulsars PSR J1412+7922 AND PSR J1849-0001

Abstract: We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a power law. Our refined timing measurements enable deeper searches for pulsations at other wavelengths and sensitive targeted searches by LIGO/Virgo for continuous gravitational waves from these neutron stars. Using the sensitivity of LIGO's first observing run, we estimate constraints that a gravitational wave search of these pulsars would obtain on the size of their mass deformation and r-mode fluid oscillation.

The semicentennial binary system PSR J2032+4127 at periastron: X-ray photometry, optical spectroscopy, and SPH modelling

Abstract: X-ray photometry and optical spectra are presented covering the periastron passage of the highly-eccentric, ~50 year binary system PSR J2032+4127 in November 2017. This system consists of a 143 ms pulsar in orbit around a massive OB star, MT 91-213. The data show dramatic changes during the encounter as the pulsar wind collided with the stellar wind. The X-ray flux rose on the approach to periastron, then underwent a major dip in the few days around periastron, and then gradually declined over the next few weeks. The optical spectroscopy revealed a steady decline in the H? line strength on the approach to periastron (from an Equivalent Width of -15A to -7A) implying a truncation of the OB star's circumstellar disk by the approaching neutron star. Smooth Particle Hydrodynamic (SPH) modelling is used here to model the system within the context of the observed behaviour and predict the geometrical configuration of the circumstellar disk with respect to the pulsar's orbit.

Modelling the observable behaviour of SXP 5.05

Abstract: SXP 5.05 is a Be/X-ray binary with a neutron star companion located in the Small Magellanic Cloud. It was first detected in 2013, and later that year, SXP 5.05 underwent a massive optical and X-ray outburst. This outburst dwarfs any other optical event that has been observed for this system during the last 5 years. The large increase in optical brightness of the system implies an increase in the size and density of the Be star's circumstellar disc. The X-ray data show two occultations of the neutron star per orbit and is not consistent with a neutron star passing behind the Be star, and hence the disc is responsible for these occultations. In this paper, we model the outburst of Be/neutron star binary SXP 5.05 as being due to a large increase in mass ejection by the Be star. The neutron star passes directly through the growing disc, and it is shown that the resulting obscuration can qualitatively explain the observed X-ray behaviour of the system. We find the only way to reproduce the timescales of the observed optical behaviour of the system is to increase the mass ejection substantially for a short time (<2 orbits) and to decrease the viscosity during the event. The general behaviour of the observed X-ray and Halpha line emission are also reproduced by the model. However, the inferred mass ejection and viscosity needed to produce a sufficiently rapid increase of disc size are both higher than suggested by previous works.

XMM-Newton detection and spectrum of the second fastest spinning pulsar PSR J0952-0607

Abstract: With a spin frequency of 707 Hz, PSR J0952-0607 is the second fastest spinning pulsar known. It was discovered in radio by LOFAR in 2017 at an estimated distance of either 0.97 or 1.74 kpc and has a low-mass companion with a 6.42 hr orbital period. We report discovery of the X-ray counterpart of PSR J0952-0607 using XMM-Newton. The X-ray spectra can be well-fit by a single power law model (Gamma = 2.5) or by a thermal plus power law model (kTeff = 40 eV and Gamma = 1.4). We do not detect evidence of variability, such as that due to orbital modulation from pulsar wind and companion star interaction. Because of its fast spin rate, PSR J0952-0607 is a crucial source for understanding the r-mode instability, which can be an effective mechanism for producing gravitational waves. Using the high end of our measured surface temperature, we infer a neutron star core temperature of ~10^7 K, which places PSR J0952-0607 within the window for the r-mode to be unstable unless an effect such as superfluid mutual friction damps the fluid oscillation. The measured luminosity limits the dimensionless r-mode amplitude to be less than ~1x10^-9.

Erratum: Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data (Astrophysical Journal (2019) 879 (10) DOI: 10.3847/1538-4357/Ab20cb)

Abstract: Two analysis errors have been identified that affect the results for a handful of the high-value pulsars given in Table 1 of Abbott et al. (2019). One affects the Bayesian analysis for the five pulsars that glitched during the analysis period, and the other affects the 5n-vector analysis for J0711-6830. Updated results after correcting the errors are shown in Table 1, which now supersedes the results given for those pulsars in Table 1 of Abbott et al. (2019). Updated versions of figures can be seen in Figures 1-4. Bayesian analysis.-For the glitching pulsars, the signal phase evolution caused by the glitch was wrongly applied twice and was therefore not consistent with our expected model of the pulsar phase. This error did not affect the F/G-statistic or 5n-vector analysis. Analyses of the five pulsars PSR J0205+6449, PSR J0534+2200, PSR J0835-4510, PSR J1028-5819, and PSR J1718-3825 have been repeated after correcting for the error. There are small quantitative differences in the results, but the changes do not affect the main conclusions of the paper. The largest differences are for PSR J0835-4510 (the Vela pulsar), for which the updated upper limits from the Bayesian method are found to be between 1.1 and 2 times larger than those obtained when the error was present. This appears primarily to be due to the error leading to the decohering of a strong spectral line in the LIGO Livingston detector and thus lowering the amplitude limit. 5n-vector analysis.-An error was also identified in the settings of the 5n-vector analysis, which affected the upper limit computation at the rotation frequency for C21 95% of J0711-6830. Specifically, we found an incorrect choice for the range of amplitudes used to inject simulated signals in the O2 data. The updated upper limit is about 2.5 times worse than that obtained when the error was present. This error did not affect the Bayesian or F/G-statistic results. (Table Presented) (Figure Presented).

Cooling of the Cassiopeia A neutron star and the effect of diffusive nuclear burning

Abstract: The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star.

Cooling of the Cassiopeia A neutron star and the effect of diffusive nuclear burning

Abstract: The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star.