Search NASASearch

Engineering topics

Paul S Ray

Publications and source records attributed to Paul S Ray.

Discovery, Timing, and Multiwavelength Observations of the Black Widow Millisecond Pulsar PSR J1555–2908

We report the discovery of PSR J1555−2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 M⊙. This fast and energetic ((E =3 ×10(exp 35)) ̇ erg s(exp -1) millisecond pulsar was first detected as a gamma-ray point source in Fermi Large Area Telescope (LAT) sky survey observations. Guided by a steep-spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the discovery and timing of radio and gamma-ray pulsations, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ≈10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system.

Millisecond pulsars

A Comprehensive X-ray Report on AT2019wey

The Galactic low-mass X-ray binary AT2019wey (ATLAS19bcxp, SRGAJ043520.9+552226, SRGEJ043523.3+552234, ZTF19acwrvzk) was discovered as a new optical transient in Dec 2019, and independently as an X-ray transient in Mar 2020. In this paper, we present comprehensive NICER, NuSTAR, Chandra, Swift, and MAXI observations of AT2019wey from ~1 year prior to the discovery to the end of September 2020. AT2019wey appeared as a ~ 1 mCrab source and stayed at this flux density for several months, displaying a hard X-ray spectrum that can be modeled as a power-law with photon index Γ ~ 1.8. In June 2020 it started to brighten and reached ~20 mCrab in ~2 months. The inclination of this system can be constrained to I ≲ 30° by modelling the reflection spectrum. Starting from late-August (~ 59082 MJD), AT2019wey entered into the hard-intermediate state (HIMS), and underwent a few week-long timescale outbursts, where the brightening in soft X-rays is correlated with the enhancement of a thermal component. Low-frequency quasi-periodic oscillation (QPO) was observed in the HIMS. We detect no pulsation and in timing analysis of the NICER and NuSTAR data. The X-ray states and power spectra of AT2019wey are discussed against the landscape of low-mass X-ray binaries.

Individual X-rays

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars

Enhanced X-Ray Emission Coinciding With Giant Radio Pulses From the Crab Pulsar

Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars that last a few microseconds and are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside of radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percentage points during GRPs. We observed the Crab Pulsar simultaneously at x-ray and radio wavelengths, finding enhancement of the x-ray emission by 3.8 ± 0.7% (a 5.4sdetection)coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts.

Teruaki Enoto

Multimessenger Gravitational-wave Searches with Pulsar Timing Arrays: Application to 3C 66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65±0.02)×109Me using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over “blind” pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

Zaven Arzoumanian

Return of the Big Glitcher: NICER Timing and Glitches of PSR J0537-6910

PSR J0537−6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 yr of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire time span, which overlaps with the third observing run of LIGO/Virgo, thus enabling the most sensitive gravitational wave searches of this potentially strong gravitational wave-emitting pulsar. We find that the short-term braking index between glitches decreases towards a value of 7 or lower at longer times since the preceding glitch. By combining NICER and RXTE data, we measure a long-term braking index n = −1.25 ± 0.01. Our analysis reveals eight new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88 × 10−7 yr−1. The new glitches follow the seemingly unique time-to-next-glitch–glitch-size correlation established previously using RXTE data, with a slope of 5 d μHz−1. For one glitch around which NICER observes 2 d on either side, we search for but do not see clear evidence of spectral nor pulse profile changes that may be associated with the glitch.

Wynn C G Ho

Timing the Pulsations of The Accreting Millisecond Pulsar SAX J1808.4-3658 During Its 2019 Outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

Peter Bult

Timing of the Accreting Millisecond Pulsar IGR J17591-2342: Evidence of Spin-down during Accretion

We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591–2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ( ˙ν ∼ −7 × 10−14 Hz s−1), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field Beq = 2.8(3) × 108 G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy is compatible with the down-scattering of hard X-ray photons in the disc or the neutron star surface.

A Sanna

A NICER View of Spectral and Profile Evolution for Three X-ray-Emitting Millisecond Pulsars

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821−24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at >3σ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821−24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821−24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

Dominick M Rowan

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

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.

Sebastien Guillot

A Joint NICER and XMM-Newton View of the "Magnificent" Thermally Emitting X-ray Isolated Neutron Star RX J1605.3+3249

Thermally emitting X-ray isolated neutron stars (NSs) represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer the physical parameters of the NS. 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., we find no significant pulsation with a pulsed fraction higher than 1.3% (3σ) 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 NS 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.

Christian Malacaria

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

We present evidence for 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. Our timing analysis reveals three sudden spin-down events of magnitudes Δν=−23, −30, and −43 μHz (fractional amplitudes Δν/ν=−4.4, −5.5, and −7.7×10−4). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known 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-superfluid 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.

Paul S Ray

Discovery of Soft X-ray Pulsations from PSR J1231-1411 using NICER

We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231−1411 using the Neutron Star Interior Composition Explorer (NICER). The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s−1 in the 0.35–1.5 keV band. We analyzed Fermi Large Area Telescope (LAT) data to update the pulse timing model covering 10 yr of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and γ rays provides insight into the geometry of the system.

gamma rays

Using Pulsars For Clock Steering and Time Transfer

As we prepare for extended missions on the Moon and in cislunar space, precise time synchronization between lunar and terrestrial clocks becomes increasingly critical for navigation, communication, and science measurements. One promising technique is the use of millisecond pulsars, which are highly magnetized, rapidly rotating neutron stars emitting beams of electromagnetic radiation. These celestial clocks offer extraordinary stability and precision, with rotational periods ranging from about 1 to 10 milliseconds and timing stability comparable to laboratory atomic clocks on Earth. Pulsars can provide a common time reference for clock comparisons, allowing time transfer independent of a communications link.

Spacecraft Communications, Command and Tracking