Search NASA⌕ Search

Engineering topics

P S Ray

Publications and source records attributed to P S Ray.

MAXI J1957+032: A New Accreting Millisecond X-ray Pulsar in an Ultra-Compact Binary

The detection of coherent X-ray pulsations at ∼314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterized by an orbital period of ∼1 h and a projected semimajor axis of 14 lt-ms. The neutron star binary mass function suggests a minimum donor mass of 1.7 × 10 −2 M ⊙ , assuming a neutron star mass of 1.4 M ⊙ and a binary inclination angle lower than 60 deg. This assumption is supported by the lack of eclipses or dips in the X-ray light curve of the source. We characterized the 0.5–10 keV energy spectrum of the source in outburst as the superposition of a relatively cold black-body-like thermal emission compatible with the emission from the neutron star surface and a Comptonization component with photon index consistent with a typical hard state. We did not find evidence for iron K α lines or reflection components.

X-rays:binaries↗

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

PSR J0740+6620 has a gravitational mass of 2.08±0.07Me, 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.71.52.6km(68%).We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451,the masses of two other∼2Mepulsars, 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.4Meneutron star is known to±4%(68% credibility)and the radius of a 2.08Meneutron 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.4Meneutron star and 12.35±0.75 km for a 2.08Me neutron star.

M C Miller↗

Refined Mass and Geometric Measurements of the High-mass PSR J0740+6620

We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8 ms radio pulsar in orbit with a likely ultracool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100 m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining∼1.5 yr of additional, high-cadence timing data with previous measurements confirms and improves on previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass=-+mM2.08p0.070.07 M (68.3% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of=-+d1.140.150.17kpc (68.3% credibility). This improved distance confirms the ultracool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on mp for J0740+6620 by an order of magnitude within the next few years

Neutron stars↗

Discovery and Timing of Three Millisecond Pulsars in Radio and Gamma-Rays with the Giant Metrewave Radio Telescope and Fermi Large Area Telescope

We performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope γ-ray sources using the Giant Metrewave Radio Telescope (GMRT) at 322 and 607 MHz. In this paper we report the discovery of three millisecond pulsars (MSPs), PSR J0248+4230, PSR J1207–5050, and PSR J1536–4948. We conducted follow-up timing observations for ∼5 yr with the GMRT and derived phase-coherent timing models for these MSPs. PSR J0248+4230 and J1207–5050 are isolated MSPs having periodicities of 2.60 ms and 4.84 ms. PSR J1536–4948 is a 3.07 ms pulsar in a binary system with an orbital period of ∼62 days about a companion of a minimum mass of 0.32 M⊙. We also present multifrequency pulse profiles of these MSPs from the GMRT observations. PSR J1536–4948 is an MSP with an extremely wide pulse profile having multiple components. Using the radio timing ephemeris we subsequently detected γ-ray pulsations from these three MSPs, confirming them as the sources powering the γ-ray emission. For PSR J1536–4948 we performed combined radio–γ-ray timing using ∼11.6 yr of γ-ray pulse times of arrival (TOAs) along with the radio TOAs. PSR J1536–4948 also shows evidence for pulsed γ-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a ≥10 GeV point source. The multiwavelength pulse profiles of all three MSPs offer challenges to models of radio and γ-ray emission in pulsar magnetospheres.

High energy astrophysics↗

RX J0529.8-6556: a BeXRB Pulsar with an Evolving Optical Period and Out of Phase X-ray Outbursts

We report the results of eROSITA and NICER observations of the 2020 June outburst of the Be/X-ray binary pulsar RX J0529.8−6556 in the Large Magellanic Cloud, along with the analysis of archival X-ray and optical data from this source. We find two anomalous features in the system’s behaviour. First, the pulse profile observed by NICER during maximum luminosity is similar to that observed by XMM–Newton in 2000, despite the fact that the X-ray luminosity was different by two orders of magnitude. In contrast, a modest decrease in luminosity in the 2020 observations generated a significant change in pulse profile. Secondly, we find that the historical optical outbursts are not strictly periodic, as would be expected if the outbursts were triggered by periastron passage, as is generally assumed. The optical peaks are also not coincident with the X-ray outbursts. We suggest that this behaviour may result from a misalignment of the Be star disc and the orbital plane, which might cause changes in the timing of the passage of the neutron star through the disc as it processes. We conclude that the orbital period of the source remains unclear.

H Treiber↗

Modeling the Uncertainties of Solar System Ephemerides for Robust Gravitational-wave Searches with Pular-timing Arrays

The regularity of pulsar emissions becomes apparent once we reference the pulses’ times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earthʼs position with respect to the solar system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar-timing arrays. Indeed, recent array data sets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article, we describe the motivation, construction, and application of a physical model of solar system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiterʼs orbital elements) most relevant to gravitational-wave searches with pulsar-timing arrays. This model, BAYESEPHEM, was used to derive ephemeris-robust results in NANOGravʼs 11 yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11 yr data set and that this degeneracy will vanish with improved ephemerides and with pulsar-timing data sets that extend well beyond a single Jovian orbital period.

M Vallisneri↗

The NANOGrav 11 yr Data Set: Contraints on Planetary Masses Around 45 Millisecond Pulsars

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11 yr data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11 yr data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary-mass limit for our median-sensitivity pulsar within this period range is 1 Mmoon (P 100 days)-2/3.

E A Behrens↗

The NANOGrav 11 yr data set: Evolution of Gravitational-wave Background Statistics

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a gravitational-wave strain spectral index of −2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new “time slicing” analysis of the 11 yr data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential “false-positive” signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and we show how independently modeling noise across 30 Fourier components, spanning NANOGrav’s frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally, we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red-noise process with low spectral index noise, and we use a tailored noise model to remove these effects.

J S Hazboun↗

The NANOGrav 11 yr data set: Limits on Gravitational Wave Memory

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs).While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays, the nonoscillatory GW memory effect is detectable. Further, any burst of GWs will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11 yr data set for GW memory. This data set is sensitive to very low-frequency GWs of ∼3 to 400 nHz (periods of ∼11 yr–1 month). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of 2.5 × 10−14, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at <0.4 yr−1. That strain corresponds to an SMBHB merger with reduced mass of ηM ~2 × 1010Mmoonand inclination of ι = π/3 at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9 yr data set as well. This analysis found an anomolous signal, which does not appear in the 11 yr data set. This signal is not a GW, and its origin remains unknown.

K Aggarwal↗

A NICER view of PSR J0030+0541: Millisecond Pulsar Parameter Estimation

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 Req are, respectively, -1.34+ M 0.16 0.15 and-12.71+ km 1.19 1.14 , while the compactness = -GM R c 0.156+ eq 2 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.

T E Riley↗

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 Re 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 Re = 13.02+1.24-1.06 km andM = 1.44+0.15-0.14 Mmoon (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 Re and M that are consistent with those found in this work. We show that our measurements of Re and M for PSRJ0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

M C Miller↗

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

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.

A V Bilous↗

NGC 300 ULX1: Spin Evolution, Super-Eddington Accretion, and Outflows

NGC300 ULX1 is an ultraluminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 yr, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of >~20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disc.

G Vasilopoulos↗