Andrea N. Lommen

The International Pulsar Timing Array project: using pulsars as a gravitational wave detector

The International Pulsar Timing Array project combines observations of pulsars from both northern and southern hemisphere observatories with the main aim of detecting ultra-low frequency (similar to 10(-9)-10(-8) Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

LIMITS ON THE STOCHASTIC GRAVITATIONAL WAVE BACKGROUND FROM THE NORTH AMERICAN NANOHERTZ OBSERVATORY FOR GRAVITATIONAL WAVES

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best rms timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our data set to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^(–1)) < 7 × 10^(–15) (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909–3744.

CONSTRAINING THE PROPERTIES OF SUPERMASSIVE BLACK HOLE SYSTEMS USING PULSAR TIMING: APPLICATION TO 3C 66B

General expressions for the expected timing residuals induced by gravitational wave (G-wave) emission from a slowly evolving, eccentric, binary black hole system are derived here for the first time. These expressions are used to search for the signature of G-waves emitted by the proposed supermassive binary black hole system in 3C 66B. We use data from long-term timing observations of the radio pulsar PSR B1855+09. For the case of a circular orbit, the emitted G-waves should generate clearly detectable fluctuations in the pulse-arrival times of PSR B1855+09. Since no G-waves are detected, the waveforms are used in a Monte Carlo analysis in order to place limits on the mass and eccentricity of the proposed black hole system. The analysis presented here rules out the adopted system with 95% confidence. The reported analysis also demonstrates several interesting features of a G-wave detector based on pulsar timing.General expressions for the expected timing residuals induced by gravitational wave (G-wave) emission from a slowly evolving, eccentric, binary black hole system are derived here for the first time. These expressions are used to search for the signature of G-waves emitted by the proposed supermassive binary black hole system in 3C 66B. We use data from long-term timing observations of the radio pulsar PSR B1855+09. For the case of a circular orbit, the emitted G-waves should generate clearly detectable fluctuations in the pulse-arrival times of PSR B1855+09. Since no G-waves are detected, the waveforms are used in a Monte Carlo analysis in order to place limits on the mass and eccentricity of the proposed black hole system. The analysis presented here rules out the adopted system with 95% confidence. The reported analysis also demonstrates several interesting features of a G-wave detector based on pulsar timing. Subject headings: black hole physics — gravitational waves — pulsars: general — pulsars: individual (B1855+09)

Arecibo Pulsar Survey Using ALFA. II. The Young, Highly Relativistic Binary Pulsar J1906+0746

We report the discovery of PSR J1906+0746, a young 144 ms pulsar in a highly relativistic 3.98 hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of � 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4 GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, GreenBank,andParkes,wehavemeasuredthespin-downandbinaryparametersofthepulsaranditsbasicspectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance (7N57 � 0N03 yr � 1 ), we infer a total system mass of 2:61 � 0:02 M� . While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birthrate of binaries similar to PSR J1906+0746is � 60Myr � 1 .ThisimpliesthatPSRJ1906+0746willmakeasignificantcontributiontothecomputed cosmic inspiral rate of compact binary systems. Subject headingg pulsars: general — pulsars: individual (PSR J1906+0746)

The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr^(−1) for a GWB from supermassive black hole binaries of A_(gw) < 1.5 x 10^(-15). We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, Ω_(gw)(f)h^2 < 4.2 x 10^(-10). Our limit on the cosmic string GWB, Ω_(gw)(f)h^2 < 2.2 x 10^(-10), translates to a conservative limit on the cosmic string tension with Gµ < 3.3 x 10^(-8), a factor of four better than the joint Planck and high-l cosmic microwave background data from other experiments.

MASSES, PARALLAX, AND RELATIVISTIC TIMING OF THE PSR J1713+0747 BINARY SYSTEM

We report on 12 years of observations of PSR J1713+0747, a pulsar in a 68 day orbit with a white dwarf. Pulse times of arrival were measured with uncertainties as small as 200 ns. The timing data yielded measurements of the relativistic Shapiro delay, perturbations of pulsar orbital elements due to secular and annual motion of the Earth, and the pulsars parallax, as well as pulse spin-down, astrometric, and Keplerian measurements. The observations constrain the masses of the pulsar and secondary star to be m1 = 1.3 ± 0.2 M☉ and m2 = 0.28 ± 0.03 M☉, respectively (68% confidence). Combining the theoretical orbital period-core mass relation with the observational constraints yields a somewhat higher pulsar mass, m1 = 1.53 M☉. The parallax is π = 0.89 ± 0.08 mas, corresponding to a distance of 1.1 ± 0.1 kpc; the precision of the parallax measurement is limited by uncertainties in the electron content of the solar wind. The transverse velocity is unusually small, 33 ± 3 km s-1. We find significant timing noise on timescales of several years, but no more than expected by extrapolating timing noise statistics from the slow pulsar population. With the orientation of the binary orbit fully measured, we are able to improve on previous tests of equivalence principle violations.

Constraining the properties of the proposed supermassive black hole system in 3c66b: Limits from pulsar timing

General expressions for the expected timing residuals induced by gravitational wave (G-wave) emission from a slowly evolving, eccentric, binary black hole system are derived here for the first time. These expressions are used to search for the signature of G-waves emitted by the proposed supermassive binary black hole system in 3C 66B. We use data from long-term timing observations of the radio pulsar PSR B1855+09. For the case of a circular orbit, the emitted G-waves should generate clearly detectable fluctuations in the pulse-arrival times of PSR B1855+09. Since no G-waves are detected, the waveforms are used in a Monte Carlo analysis in order to place limits on the mass and eccentricity of the proposed black hole system. The analysis presented here rules out the adopted system with 95% confidence. The reported analysis also demonstrates several interesting features of a G-wave detector based on pulsar timing.General expressions for the expected timing residuals induced by gravitational wave (G-wave) emission from a slowly evolving, eccentric, binary black hole system are derived here for the first time. These expressions are used to search for the signature of G-waves emitted by the proposed supermassive binary black hole system in 3C 66B. We use data from long-term timing observations of the radio pulsar PSR B1855+09. For the case of a circular orbit, the emitted G-waves should generate clearly detectable fluctuations in the pulse-arrival times of PSR B1855+09. Since no G-waves are detected, the waveforms are used in a Monte Carlo analysis in order to place limits on the mass and eccentricity of the proposed black hole system. The analysis presented here rules out the adopted system with 95% confidence. The reported analysis also demonstrates several interesting features of a G-wave detector based on pulsar timing. Subject headings: black hole physics — gravitational waves — pulsars: general — pulsars: individual (B1855+09)

Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits: Limits from the North American Nanohertz Observatory for Gravitational Waves

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of

MEASURING THE MASS OF SOLAR SYSTEM PLANETS USING PULSAR TIMING

h_0\lesssim 3.8\times 10^{-14}

New pulsars from an arecibo drift scan search

at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the