Search NASASearch

Engineering topics

R D Ferdman

Publications and source records attributed to R D Ferdman.

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

The International Pulsar Timing Array: Second Data Release

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

neutron

Relativistic Shapiro Delay Measurements of an Extremely Massive Millisecond Pulsar

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing—the process of accounting for every rotation of a pulsar over long time periods—can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system1. One of these, called relativistic Shapiro delay2, can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-yr data set with recent orbitalphase- specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be 2.14+0:10-0:09 M⊙ (68.3% credibility interval; the 95.4% credibility interval is 2.14+0:20-0:18 M⊙). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

H T Cromartie

High-precision x-ray timing of three millisecond pulsars with NICER: Stability estimates and comparison with radio

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars (MSPs) without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the MSPs PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821−24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the σz measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities σz≈3×10−14 for PSR B1937+21 and on the order of 10−12 for PSRs B1821−24 and J0218+4232. Within the span of our NICER data we do not yet see the characteristic break point in the slope of σz; detection of such a break would indicate that further improvement in the cumulative root-mean-square timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821−24 and PSR B1937+21 on timescales of >2 yr.

Pulsars

PSR J2234+0611: A New Laboratory for Stellar Evolution

We report the timing results for PSR J2234+0611, a 3.6 ms pulsar in a 32 day, eccentric (e = 0.13) orbit with a helium white dwarf. The precise timing and eccentric nature of the orbit allow measurements of an unusual number of parameters: (a) a precise proper motion of 27.10(3) mas yr−1 and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; enabling an estimate of the transverse velocity, 123(5) km s−1. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a 3D determination of the system’s velocity; (b) precise measurements of the rate of advance of periastron yields a total system mass of 1.6518 +0.0033-0.0035 Msun; (c) a Shapiro delay measurement, h3 = 82 ± 14 ns, despite the orbital inclination not being near 90°; combined with the measurement of the total mass yields a pulsar mass of 1.353 +0.014-0.017 Msun and a companion mass of 0.298 +0.015-0.0120 Msun; (d) we measure precisely the secular variation of the projected semimajor axis and detect a significant annual orbital parallax; together these allow a determination of the 3D orbital geometry of the system, including an unambiguous orbital inclination (i = 138.7 +2.5-2.2 deg) and a position angle for the line of nodes ( omega = 44 +5-4 deg). We discuss the component masses to investigate the hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the 3D position, motion, and orbital orientation of the system, plus the precise pulsar and WD masses and the latter’s optical detection make this system a unique test of our understanding of white dwarfs and their atmospheres.

K Stovall

Tests of Gravitational Symmetries with Pulsar Binary J1713+0747

Symmetries play a fundamental role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: (1) the gravitational constant’s constancy as part of locational invariance of gravitation; (2) the universality of free fall (UFF) for strongly self-gravitating bodies; (3) the post-Newtonian parameter ˆα3 in gravitational Lorentz invariance. Based on the pulsar timing result of the combined data set from the North American Nanohertz Gravitational Observatory and the European Pulsar Timing Array, we find G˙ /G = (−0.1 ± 0.9) × 10−12 yr−1, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain an improved test for a UFF violation by a strongly self-gravitating mass falling in the gravitational field of our Galaxy, with a limit of |Delta| < 0.002 (95 per cent C.L.). Finally, we derive an improved limit on the self-acceleration of a gravitationally bound rotating body, to a preferred reference frame in the Universe, with −3 × 10−20 < ˆα3 < 4 × 10−20 (95 per cent C.L.). These results are based on direct UFF and ˆα3 tests using pulsar binaries, and they overcome various limitations of previous tests of this kind.

W W Zhu