Showing 1–50 of 61 results for author: Coles, W

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array

Authors: Daniel J. Reardon, Andrew Zic, Ryan M. Shannon, George B. Hobbs, Matthew Bailes, Valentina Di Marco, Agastya Kapur, Axl F. Rogers, Eric Thrane, Jacob Askew, N. D. Ramesh Bhat, Andrew Cameron, Małgorzata Curyło, William A. Coles, Shi Dai, Boris Goncharov, Matthew Kerr, Atharva Kulkarni, Yuri Levin, Marcus E. Lower, Richard N. Manchester, Rami Mandow, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski , et al. (4 additional authors not shown)

Abstract: Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of t… ▽ More Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^α$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $α=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $α=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $α=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2σ$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 10 figures, Accepted for publication in ApJL

Determining electron column density fluctuations in a dominant scattering region using pulsar scintillation

Authors: Daniel J. Reardon, William A. Coles

Abstract: Density fluctuations in the ionised interstellar medium have a profound effect on radio pulsar observations, through angular scattering, intensity scintillations, and small changes in time delays from dispersion. Here we show that it is possible to recover the variations in dispersive delays that originate from a dominant scattering region using measurements of the dynamic spectrum of intensity sc… ▽ More Density fluctuations in the ionised interstellar medium have a profound effect on radio pulsar observations, through angular scattering, intensity scintillations, and small changes in time delays from dispersion. Here we show that it is possible to recover the variations in dispersive delays that originate from a dominant scattering region using measurements of the dynamic spectrum of intensity scintillations, provided that the pulsar velocity and scattering region location are known. We provide a theoretical framework for the technique, which involves estimating the phase gradient from the dynamic spectra and integrating that gradient to obtain phase variations. It can be used to search for "extreme scattering events" (ESEs) in pulsars for which precision dispersion delay measurements are not otherwise possible, or to separate true dispersion variations from apparent variability caused by frequency-dependent pulse shape changes. We demonstrate that it works in practice by recovering an ESE in PSR J1603$-$7202, which is known from precision dispersion delay measurements from pulsar timing. For this pulsar, we find that the phase gradients also track the long-term variations in electron column density observed by pulsar timing, indicating that the column density variations and the scattering are dominated by the same thin scattering screen. We identify a sudden increase in the scintillation strength and magnitude of phase gradients over $\sim$days in 2010, indicating a compact structure. A decrease in the electron density in 2012 was associated with persistent phase gradients and preceded a period of decreased scintillation strength and an absence of scintillation arcs. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 9 pages, 7 figures, Accepted for publication in MNRAS

Pulsar Scintillation Studies with LOFAR: II. Dual-frequency scattering study of PSR J0826+2637 with LOFAR and NenuFAR

Authors: Ziwei Wu, William A. Coles, Joris P. W. Verbiest, Krishnakumar Moochickal Ambalappat, Caterina Tiburzi, Jean-Mathias Grießmeier, Robert A. Main, Yulan Liu, Michael Kramer, Olaf Wucknitz, Nataliya Porayko, Stefan Osłowski, Ann-Sofie Bak Nielsen, Julian Y. Donner, Matthias Hoeft, Marcus Brüggen, Christian Vocks, Ralf-Jürgen Dettmar, Gilles Theureau, Maciej Serylak, Vladislav Kondratiev, James W. McKee, Golam M. Shaifullah, Ihor P. Kravtsov, Vyacheslav V. Zakharenko , et al. (6 additional authors not shown)

Abstract: Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to st… ▽ More Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to study the properties of the IISM. Larger scales ($\sim$1-100\,AU) cause measurable changes in dispersion and these can be correlated with ISS observations to estimate the fluctuation spectrum over a very wide scale range. IISM measurements can often be modeled by a homogeneous power-law spatial spectrum of electron density with the Kolmogorov ($-11/3$) spectral exponent. Here we aim to test the validity of using the Kolmogorov exponent with PSR~J0826+2637. We do so using observations of intensity scintillation, pulse broadening and dispersion variations across a wide fractional bandwidth (20 -- 180\,MHz). We present that the frequency dependence of the intensity scintillation in the high frequency band matches the expectations of a Kolmogorov spectral exponent but the pulse broadening in the low frequency band does not change as rapidly as predicted with this assumption. We show that this behavior is due to an inhomogeneity in the scattering region, specifically that the scattering is dominated by a region of transverse size $\sim$40\,AU. The power spectrum of the electron density, however, maintains the Kolmogorov spectral exponent from spatial scales of 5$\times10^{-6}$\,AU to $\sim$100\,AU. △ Less

Submitted 25 February, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS, typo fixed

Interstellar scintillation and polarization of PSR B0656+14 in the Monogem Ring

Authors: Jumei Yao, William A. Coles, Richard N. Manchester, Daniel R. Stinebring, Michael Kramer, Na Wang, Di Li, Weiwei Zhu, Yi Feng, Jianping Yuan, Pei Wang

Abstract: High sensitivity interstellar scintillation and polarization observations of PSR~B0656+14 made at three epochs over a year using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) show that the scattering is dominated by two different compact regions. We identify the one nearer to the pulsar with the shell of the Monogem Ring, thereby confirming the association. The other is probably… ▽ More High sensitivity interstellar scintillation and polarization observations of PSR~B0656+14 made at three epochs over a year using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) show that the scattering is dominated by two different compact regions. We identify the one nearer to the pulsar with the shell of the Monogem Ring, thereby confirming the association. The other is probably associated with the Local Bubble. We find that the observed position angles of the pulsar spin axis and the space velocity are significantly different, with a separation of $19\fdg3\pm$0\fdg8, inconsistent with a previously published near-perfect alignment of $1\degr\pm 2\degr$. The two independent scattering regions are clearly defined in the secondary spectra which show two strong forward parabolic arcs. The arc curvatures imply that the scattering screens corresponding to the outer and inner arcs are located approximately 28~pc from PSR B0656+14 and 185~pc from the Earth, respectively. Comparison of the observed Doppler profiles with electromagnetic simulations shows that both scattering regions are mildly anisotropic. For the outer arc, we estimate the anisotropy $A_R$ to be approximately 1.3, with the scattering irregularities aligned parallel to the pulsar velocity. For the outer arc, we compare the observed delay profiles with delay profiles computed from a theoretical strong-scattering model. Our results suggest that the spatial spectrum of the scattering irregularities in the Monogem Ring is flatter than Kolmogorov, but further observations are required to confirm this. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 21 pages, 11 figures, accepted for publication in the Astronomical Journal

Detection of strong scattering close to the eclipse region of PSR B1957+20

Authors: J. T. Bai, S. Dai, Q. J. Zhi, W. A. Coles, D. Li, W. W. Zhu, G. Hobbs, G. J. Qiao, N. Wang, J. P. Yuan, M. D. Filipovic, J. B. Wang, Z. C. Pan, L. H. Shang, S. J. Dang, S. Q. Wang, C. C. Miao

Abstract: We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $τ\proptoΔ{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering du… ▽ More We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $τ\proptoΔ{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 8 pages, 4 figures, MNRAS accepted

Strong-field Gravity Tests with the Double Pulsar

Authors: M. Kramer, I. H. Stairs, R. N. Manchester, N. Wex, A. T. Deller, W. A. Coles, M. Ali, M. Burgay, F. Camilo, I. Cognard, T. Damour, G. Desvignes, R. D. Ferdman, P. C. C. Freire, S. Grondin, L. Guillemot, G. B. Hobbs, G. Janssen, R. Karuppusamy, D. R. Lorimer, A. G. Lyne, J. W. McKee, M. McLaughlin, L. E. Muench, B. B. P. Perera , et al. (5 additional authors not shown)

Abstract: Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating… ▽ More Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. We observe retardation and aberrational light-bending that allow determination of the pulsar's spin direction. In total, we have detected seven post-Keplerian (PK) parameters, more than for any other binary pulsar. For some of these effects, the measurement precision is so high that for the first time we have to take higher-order contributions into account. These include contributions of A's effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and effects of the equation of state of super-dense matter on the observed PK parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars. We present the currently most precise test of general relativity's (GR's) quadrupolar description of gravitational waves, validating GR's prediction at a level of $1.3 \times 10^{-4}$ (95% conf.). We demonstrate the utility of the Double Pulsar for tests of alternative theories by focusing on two specific examples and discuss some implications for studies of the interstellar medium and models for the formation of the Double Pulsar. Finally, we provide context to other types of related experiments and prospects for the future. △ Less

Submitted 14 December, 2021; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 56 pages, 23 Figures. Published by Physical Review X. Uploaded ArXiv version is authors' (pre-proof) version with abbreviated abstract. Resubmission for updated acknowledgement information. For final published version see https://link.aps.org/doi/10.1103/PhysRevX.11.041050

Journal ref: Phys. Rev. X 11, 041050 (2021)

The first evidence for three-dimensional spin-velocity alignment in pulsars

Authors: Jumei Yao, Weiwei Zhu, Richard N. Manchester, William A. Coles, Di Li, Na Wang, Michael Kramer, Daniel R. Stinebring, Yi Feng, Wenming Yan, Chenchen Miao, Mao Yuan, Pei Wang, Jiguang Lu

Abstract: More than 50 years after the discovery of pulsars and confirmation of their association with supernova explosions, the origin of the initial spin and velocity of pulsars remains largely a mystery. The typical space velocities of several hundred km/s have been attributed to "kicks" resulting from asymmetries either in the supernova ejecta or in the neutrino emission. Observations have shown a stron… ▽ More More than 50 years after the discovery of pulsars and confirmation of their association with supernova explosions, the origin of the initial spin and velocity of pulsars remains largely a mystery. The typical space velocities of several hundred km/s have been attributed to "kicks" resulting from asymmetries either in the supernova ejecta or in the neutrino emission. Observations have shown a strong tendency for alignment of the pulsar space velocity and spin axis in young pulsars but, up to now, these comparisons have been restricted to two dimensions. We report here the first evidence for three-dimensional alignment between the spin and velocity vectors, largely based on observations made with the Five-hundred-meter Aperture Spherical radio Telescope of the pulsar PSR~J0538+2817 and its associated supernova remnant S147. Analysis of these and related observations has enabled us to determine the location of the pulsar within the supernova remnant and hence its radial velocity. Current simulations of supernova explosions have difficulty producing such three-dimensional alignment. Our results, which depend on the unprecedented sensitivity of the new observations, add another dimension to the intriguing correlation between pulsar spin-axis and birth-kick directions, thus deepening the mysteries surrounding the birth of neutron stars. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 33 pages, 8 figures, Accepted by Nature Astronomy

Searching for gravitational wave bursts from cosmic string cusps with the Parkes Pulsar Timing Array

Authors: N. Yonemaru, S. Kuroyanagi, G. Hobbs, K. Takahashi, X. -J. Zhu, W. A. Coles, S. Dai, E. Howard, R. Manchester, D. Reardon, C. Russell, R. Shannon, N. Thyagarajan, R. Spiewak, J. -B. Wang

Abstract: Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no… ▽ More Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA). △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication by MNRAS

Precision orbital dynamics from interstellar scintillation arcs for PSR J0437-4715

Authors: Daniel J. Reardon, William A. Coles, Matthew Bailes, N. D. Ramesh Bhat, Shi Dai, George B. Hobbs, Matthew Kerr, Richard N. Manchester, Stefan Oslowski, Aditya Parthasarathy, Christopher J. Russell, Ryan M. Shannon, Renee Spiewak, Lawrence Toomey, Artem V. Tuntsov, Willem van Straten, Mark A. Walker, Jingbo Wang, Lei Zhang, Xing-Jiang Zhu

Abstract: Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the l… ▽ More Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line-of-sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our $>$5000 observations and we show that they originate from scattering by thin screens located at distances $D_1 = 89.8 \pm 0.4$ pc and $D_2 = 124 \pm 3$ pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle $i = 137.1 \pm 0.3^\circ$, and longitude of ascending node, $Ω=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars, and provides parameters essential for testing theories of gravity and constraining neutron star masses. △ Less

Submitted 27 September, 2020; originally announced September 2020.

Comments: 19 pages, 10 figures. Accepted for publication in ApJ

The Parkes Pulsar Timing Array Project: Second data release

Authors: M. Kerr, D. J. Reardon, G. Hobbs, R. M. Shannon, R. N. Manchester, S. Dai, C. J. Russell, S. -B. Zhang, W. van Straten, S. Osłowski, A. Parthasarathy, R. Spiewak, M. Bailes, N. D. R. Bhat, A. D. Cameron, W. A. Coles, J. Dempsey, X. Deng, B. Goncharov, J. F Kaczmarek, M. J. Keith, P. D. Lasky, M. E. Lower, B. Preisig, J. M. Sarkissian , et al. (5 additional authors not shown)

Abstract: We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is… ▽ More We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 $μ$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files. △ Less

Submitted 21 March, 2020; originally announced March 2020.

Comments: 26 pages, 6 figures. Accepted for publication in PASA. The data release can be downloaded from https://doi.org/10.25919/5db90a8bdeb59

A pulsar-based timescale from the International Pulsar Timing Array

Authors: G. Hobbs, L. Guo, R. N. Caballero, W. Coles, K. J. Lee, R. N. Manchester, D. J. Reardon, D. Matsakis, M. L. Tong, Z. Arzoumanian, M. Bailes, C. G. Bassa, N. D. R. Bhat, A. Brazier, S. Burke-Spolaor, D. J. Champion, S. Chatterjee, I. Cognard, S. Dai, G. Desvignes, T. Dolch, R. D. Ferdman, E. Graikou, L. Guillemot, G. H. Janssen , et al. (34 additional authors not shown)

Abstract: We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties… ▽ More We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: accepted by MNRAS

PSR J1926-0652: A Pulsar with Interesting Emission Properties Discovered at FAST

Authors: Lei Zhang, Di Li, George Hobbs, Crispin H. Agar, Richard N. Manchester, Patrick Weltevrede, William A. Coles, Pei Wang, Weiwei Zhu, Zhigang Wen, Jianping Yuan, Andrew D. Cameron, Shi Dai, Kuo Liu, Qijun Zhi, Chenchen Miao, Mao Yua, Shuyun Cao, Li Feng, Hengqian Gan, Long Gao, Xuedong Gu, Minglei Guo, Qiaoli Hao, Lin Huang , et al. (37 additional authors not shown)

Abstract: We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual… ▽ More We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 13pages with 12 figures

Modelling annual and orbital variations in the scintillation of the relativistic binary PSR J1141$-$6545

Authors: D. J. Reardon, W. A. Coles, G. Hobbs, S. Ord, M. Kerr, M. Bailes, N. D. R. Bhat, V. Venkatraman Krishnan

Abstract: We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pu… ▽ More We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pulsar timing alone. These include: the orbital inclination $i$; the longitude of the ascending node $Ω$; and the pulsar system transverse velocity. We use the annual variations to resolve the previous ambiguity in the sense of the inclination angle. Using the correct sense, and a prior probability distribution given by a constraint from pulsar timing ($i=73\pm3^\circ$), we find $Ω=24.8\pm1.8^\circ$ and we estimate the pulsar distance to be $D=10^{+4}_{-3}$ kpc. This then gives us an estimate of this pulsar's proper motion of $μ_α\cosδ=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $μ_δ=1.8\pm0.6$ mas yr$^{-1}$ in declination. Finally, we obtain measurements of the spatial structure of the interstellar electron density fluctuations, including: the spatial scale and anisotropy of the diffraction pattern; the distribution of scattering material along the line of sight; and spatial variation in the strength of turbulence from epoch to epoch. We find that the scattering is dominated by a thin screen at a distance of $(0.724\pm0.008)D$, with an anisotropy axial ratio $A_{\rm r} = 2.14\pm0.11$. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: 17 pages, 8 figures, 2 tables. Accepted for publication in MNRAS

Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter

Authors: Nataliya K. Porayko, Xingjiang Zhu, Yuri Levin, Lam Hui, George Hobbs, Aleksandra Grudskaya, Konstantin Postnov, Matthew Bailes, N. D. Ramesh Bhat, William Coles, Shi Dai, James Dempsey, Michael J. Keith, Matthew Kerr, Michael Kramer, Paul D. Lasky, Richard N. Manchester, Stefan Osłowski, Aditya Parthasarathy, Vikram Ravi, Daniel J. Reardon, Pablo A. Rosado, Christopher J. Russell, Ryan M. Shannon, Renée Spiewak , et al. (5 additional authors not shown)

Abstract: It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, pr… ▽ More It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities. △ Less

Submitted 16 December, 2019; v1 submitted 7 October, 2018; originally announced October 2018.

Comments: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables

Journal ref: Phys. Rev. D 98, 102002 (2018)

Polarimetry of the eclipsing pulsar PSR J1748$-$2446A

Authors: X. P. You, R. N. Manchester, W. A. Coles, G. B. Hobbs, R. Shannon

Abstract: Observations with the Parkes radio telescope of the eclipsing millisecond binary pulsar PSR J1748$-$2446A, which is in the globular cluster Terzan 5, are presented. These include the first observations of this pulsar in the 3\,GHz frequency band, along with simultaneous observations in the 700\,MHz band and new observations around 1400\,MHz. We show that the pulsar signal is not eclipsed in the 3\… ▽ More Observations with the Parkes radio telescope of the eclipsing millisecond binary pulsar PSR J1748$-$2446A, which is in the globular cluster Terzan 5, are presented. These include the first observations of this pulsar in the 3\,GHz frequency band, along with simultaneous observations in the 700\,MHz band and new observations around 1400\,MHz. We show that the pulsar signal is not eclipsed in the 3\,GHz band and observe the known eclipses in the lower frequency bands. We find that the observed pulse signal becomes depolarized during particular orbital phases and postulate that this depolarization occurs because of rotation-measure fluctuations resulting from turbulence in the stellar wind responsible for the eclipses. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 9 papes, 3 figures, ApJ accepted

Extreme Scattering Events Towards Two Young Pulsars

Authors: M. Kerr, W. A. Coles, C. A. Ward, S. Johnston, A. V. Tuntsov, R. M. Shannon

Abstract: We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 a… ▽ More We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 au, and strong flux density enhancement before and after the ESE may arise from refractive focusing. Long observations reveal scintillation arcs characteristic of interference between rays scattered at large angles, and the clearest arcs appear during the ESE. The arcs suggest scattering by a screen 100-200 pc from the earth, perhaps ionized filamentary structure associated with the boundary of the local bubble(s). Towards PSR J1740-3015 we observed a "double dip" in the measured flux density similar to ESEs observed towards compact extragalactic radio sources. The observed shape is consistent with that produced by a many-au scale diverging plasma lens with electron density ~500cm$^{-3}$ . The continuing ESE is at least 1500 d long, making it the longest detected event to date. These detections, with materially different observational signatures, indicate that well-calibrated pulsar monitoring is a keen tool for ESE detection and ISM diagnostics. They illustrate the strong rôle au-scale non-Kolmogorov density fluctuations and the local ISM structure play in such events and are key to understanding both their intrinsic physics and their impact on other phenomena, particularly fast radio bursts. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: MNRAS accepted

Comparison of Pulsar Positions from Timing and Very Long Baseline Astrometry

Authors: J. B. Wang, W. A. Coles, G. Hobbs, R. M. Shannon, R. N. Manchester, M. Kerr, J. P. Yuan, N. Wang, M. Bailes, N. D. R. Bhat, S. Dai J. Dempsey, M. J. Keith, P. D. Lasky, Y. Levin, S. Os lowski, V. Ravi, D. J. Reardon, P A. Rosado, C. J. Russell, R. Spiewak, W. van Straten, L. Toomey, L. Wen, X. -P. You, X. -J. Zhu

Abstract: Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radi… ▽ More Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$μ$as level between ephemerides. △ Less

Submitted 31 March, 2017; originally announced April 2017.

Comments: accepted by MNRAS,22 pages, 5 figures

Wide-band Profile Domain Pulsar Timing Analysis

Authors: L. Lentati, M. Kerr, S. Dai, M. P. Hobson, R. M. Shannon, G. Hobbs, M. Bailes, N. D. Ramesh Bhat, S. Burke-Spolaor, W. Coles, J. Dempsey, P. D. Lasky, Y. Levin, R. N. Manchester, S. Oslowski, V. Ravi, D. J. Reardon, P. A. Rosado, R. Spiewak, W. van Straten, L. Toomey, J. Wang, L. Wen, X. You, X. Zhu

Abstract: We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the… ▽ More We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $\sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band. △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: 22 pages, submitted to MNRAS

The disturbance of a millisecond pulsar magnetosphere

Authors: R. M. Shannon, L. T. Lentati, M. Kerr, M. Bailes, N. D. R. Bhat, W. A. Coles, S. Dai, J. Dempsey, G. Hobbs, M. J. Keith, P. D. Lasky, Y. Levin, R. N. Manchester, S. Oslowski, V. Ravi, D. J. Reardon, P. A. Rosado, R. Spiewak, W. van Straten, L. Toomey, J. -B. Wang, L. Wen, X. -P. You, X. -J. Zhu

Abstract: Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 month… ▽ More Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability. △ Less

Submitted 6 August, 2016; originally announced August 2016.

Comments: Accepted to ApJ Letters

Detecting Pulsars with Interstellar Scintillation in Variance Images

Authors: S. Dai, S. Johnston, M. E. Bell, W. A. Coles, G. Hobbs, R. D. Ekers, E. Lenc

Abstract: Pulsars are the only cosmic radio sources known to be sufficiently compact to show diffractive interstellar scintillations. Images of the variance of radio signals in both time and frequency can be used to detect pulsars in large-scale continuum surveys using the next generation of synthesis radio telescopes. This technique allows a search over the full field of view while avoiding the need for ex… ▽ More Pulsars are the only cosmic radio sources known to be sufficiently compact to show diffractive interstellar scintillations. Images of the variance of radio signals in both time and frequency can be used to detect pulsars in large-scale continuum surveys using the next generation of synthesis radio telescopes. This technique allows a search over the full field of view while avoiding the need for expensive pixel-by-pixel high time resolution searches. We investigate the sensitivity of detecting pulsars in variance images. We show that variance images are most sensitive to pulsars whose scintillation time-scales and bandwidths are close to the subintegration time and channel bandwidth. Therefore, in order to maximise the detection of pulsars for a given radio continuum survey, it is essential to retain a high time and frequency resolution, allowing us to make variance images sensitive to pulsars with different scintillation properties. We demonstrate the technique with Murchision Widefield Array data and show that variance images can indeed lead to the detection of pulsars by distinguishing them from other radio sources. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 8 papes, 9 figures, accepted for publication in MNRAS

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

Authors: L. Lentati, R. M. Shannon, W. A. Coles, J. P. W. Verbiest, R. van Haasteren, J. A. Ellis, R. N. Caballero, R. N. Manchester, Z. Arzoumanian, S. Babak, C. G. Bassa, N. D. R. Bhat, P. Brem, M. Burgay, S. Burke-Spolaor, D. Champion, S. Chatterjee, I. Cognard, J. M. Cordes, S. Dai, P. Demorest, G. Desvignes, T. Dolch, R. D. Ferdman, E. Fonseca , et al. (58 additional authors not shown)

Abstract: We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single… ▽ More We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: 29 pages. 16 figures. Accepted for publication in MNRAS

Broadband Meter-Wavelength Observations of Ionospheric Scintillation

Authors: R. A. Fallows, W. A. Coles, D. McKay, J. Vierinen, I. I. Virtanen, M. Postila, Th. Ulich, C-F. Enell, A. Kero, T. Iinatti, M. Lehtinen, M. Orispää, T. Raita, L. Roininen, E. Turunen, M. Brentjens, N. Ebbendorf, M. Gerbers, T. Grit, P. Gruppen, H. Meulman, M. Norden, J-P. de Reijer, A. Schoenmakers, K. Stuurwold

Abstract: Intensity scintillations of cosmic radio sources are used to study astrophysical plasmas like the ionosphere, the solar wind, and the interstellar medium. Normally these observations are relatively narrow band. With Low Frequency Array (LOFAR) technology at the Kilpisjärvi Atmospheric Imaging Receiver Array (KAIRA) station in northern Finland we have observed scintillations over a 3 octave bandwid… ▽ More Intensity scintillations of cosmic radio sources are used to study astrophysical plasmas like the ionosphere, the solar wind, and the interstellar medium. Normally these observations are relatively narrow band. With Low Frequency Array (LOFAR) technology at the Kilpisjärvi Atmospheric Imaging Receiver Array (KAIRA) station in northern Finland we have observed scintillations over a 3 octave bandwidth. ``Parabolic arcs'', which were discovered in interstellar scintillations of pulsars, can provide precise estimates of the distance and velocity of the scattering plasma. Here we report the first observations of such arcs in the ionosphere and the first broad-band observations of arcs anywhere, raising hopes that study of the phenomenon may similarly improve the analysis of ionospheric scintillations. These observations were made of the strong natural radio source Cygnus-A and covered the entire 30-250\,MHz band of KAIRA. Well-defined parabolic arcs were seen early in the observations, before transit, and disappeared after transit although scintillations continued to be obvious during the entire observation. We show that this can be attributed to the structure of Cygnus-A. Initial results from modeling these scintillation arcs are consistent with simultaneous ionospheric soundings taken with other instruments, and indicate that scattering is most likely to be associated more with the topside ionosphere than the F-region peak altitude. Further modeling and possible extension to interferometric observations, using international LOFAR stations, are discussed. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: 11 pages, 17 figures

Journal ref: Fallows, R A et al, Broadband meter-wavelength observations of ionospheric scintillation, Journal of Geophysical Research Space Physics, 119(12), 10-544, 2014

Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays

Authors: D. R. Madison, X. -J. Zhu, G. Hobbs, W. Coles, R. M. Shannon, J. Wang, C. Tiburzi, R. N. Manchester, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, J. Dempsey, M. Keith, M. Kerr, P. Lasky, Y. Levin, S. Oslowski, V. Ravi, D. Reardon, P. Rosado, R. Spiewak, W. van Straten, L. Toomey, L. Wen , et al. (1 additional authors not shown)

Abstract: By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "… ▽ More By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise. △ Less

Submitted 27 October, 2015; originally announced October 2015.

Comments: Accepted for publication by MNRAS

Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array

Authors: D. J. Reardon, G. Hobbs, W. Coles, Y. Levin, M. J. Keith, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, M. Kerr, P. D. Lasky, R. N. Manchester, S. Osłowski, V. Ravi, R. M. Shannon, W. van Straten, L. Toomey, J. Wang, L. Wen, X. P. You, X. -J. Zhu

Abstract: We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual… ▽ More We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date. △ Less

Submitted 15 October, 2015; originally announced October 2015.

Comments: 21 pages, 5 figures, 7 tables. Accepted for publication in MNRAS

A study of spatial correlations in pulsar timing array data

Authors: Caterina Tiburzi, George Hobbs, Matthew Kerr, William Coles, Shi Dai, Richard Manchester, Andrea Possenti, Ryan Shannon, Xiaopeng You

Abstract: Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, err… ▽ More Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, errors in the planetary ephemeris, the solar wind and instrumentation issues. All these effects can induce significant false detections of gravitational waves. We test mitigation routines to account for clock errors, ephemeris errors and the solar wind. We demonstrate that it is non-trivial to find an effective mitigation routine for the planetary ephemeris and emphasise that other spatially correlated signals may be present in the data. △ Less

Submitted 8 October, 2015; originally announced October 2015.

Comments: 13 pages, 7 figures. Accepted for publication in MNRAS

Gravitational waves from binary supermassive black holes missing in pulsar observations

Authors: R. M. Shannon, V. Ravi, L. T. Lentati, P. D. Lasky, G. Hobbs, M. Kerr, R. N. Manchester, W. A. Coles, Y. Levin, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, M. J. Keith, S. Osłowski, D. J. Reardon, W. van Straten, L. Toomey, J. -B. Wang, L. Wen, J. S. B. Wyithe, X. -J. Zhu

Abstract: Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background,… ▽ More Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves. △ Less

Submitted 24 September, 2015; originally announced September 2015.

Comments: Published in Science 25 September 2015

Pulsar Observations of Extreme Scattering Events

Authors: W. A. Coles, M. Kerr, R. M. Shannon, G. Hobbs, R. N. Manchester, X. P. You, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, M. J. Keith, Y. Levin, S. Oslowski, V. Ravi, D. Reardon, L. Toomey, W. van Straten, J. B. Wang, L. Wen, X. J. Zhu

Abstract: Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the elect… ▽ More Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin. △ Less

Submitted 25 June, 2015; originally announced June 2015.

Comments: Accepted by ApJ

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

Authors: S. Dai, G. Hobbs, R. N. Manchester, M. Kerr, R. M. Shannon, W. van Straten, A. Mata, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, W. A. Coles, S. Johnston, M. J. Keith, Y. Levin, S. Oslowski, D. Reardon, V. Ravi, J. M. Sarkissian, C. Tiburzi, L. Toomey, H. G. Wang, J. -B. Wang, L. Wen, R. X. Xu, W. M. Yan , et al. (1 additional authors not shown)

Abstract: We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximate… ▽ More We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile. △ Less

Submitted 5 March, 2015; originally announced March 2015.

Comments: 41 pages, 28 figures, accepted for publication in MNRAS

Searching for gravitational wave memory bursts with the Parkes Pulsar Timing Array

Authors: J. B. Wang, G. Hobbs, W. Coles, R. M. Shannon, X. J. Zhu, D. R. Madison, M. Kerr, V. Ravi, M. J. Keith, R. N. Manchester, Y. Levin, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, S. Oslowski, W. van Straten, L. Toomey, N. Wang, L. Wen

Abstract: Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simu… ▽ More Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected. △ Less

Submitted 13 October, 2014; originally announced October 2014.

Comments: Accepted by MNRAS

An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set

Authors: X. -J. Zhu, G. Hobbs, L. Wen, W. A. Coles, J. -B. Wang, R. M. Shannon, R. N. Manchester, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, M. J. Keith, M. Kerr, Y. Levin, D. R. Madison, S. Osłowski, V. Ravi, L. Toomey, W. van Straten

Abstract: We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each p… ▽ More We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$. △ Less

Submitted 18 September, 2014; v1 submitted 21 August, 2014; originally announced August 2014.

Comments: A few typos corrected, to match the published version, 13 pages, 11 figures

Journal ref: MNRAS 2014 444 (4): 3709-3720

Limitations in timing precision due to single-pulse shape variability in millisecond pulsars

Authors: R. M. Shannon, S. Osłowski, S. Dai, M. Bailes, G. Hobbs, R. N. Manchester, W. van Straten, C. A. Raithel, V. Ravi, L. Toomey, N. D. R. Bhat, S. Burke-Spolaor, W. A. Coles, M. J. Keith, M. Kerr, Y. Levin, J. M. Sarkissian, J. -B. Wang, L. Wen, X. -J. Zhu

Abstract: High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short… ▽ More High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities. △ Less

Submitted 18 June, 2014; originally announced June 2014.

Comments: 21 pages, 23 figures, accepted by MNRAS

Interstellar Scintillation of the Double Pulsar J0737$-$3039

Authors: B. J. Rickett, W. A. Coles, C. F. Nava, M. A. McLaughlin, S. M. Ransom, F. Camilo, R. D. Ferdman, P. C. C. Freire, M. Kramer, A. G. Lyne, I. H. Stairs

Abstract: We report here a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737$-$3039 over the course of 18 months. As in earlier work (Coles et al., 2005) the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be… ▽ More We report here a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737$-$3039 over the course of 18 months. As in earlier work (Coles et al., 2005) the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be determined both by VLBI and timing observations. The orbital velocity and inclination is almost completely determined from timing observations, but the direction of the orbital angular momentum is not known. Since the Earth's velocity is known, and can be compared with the orbital velocity by its effect on the timescale of the ISS, we can determine the orientation $Ω$ of the pulsar orbit with respect to equatorial coordinates ($Ω= 65\pm2$ deg). We also resolve the ambiguity ($i= 88.7$ or $91.3$ deg) in the inclination of the orbit deduced from the measured Shapiro delay by our estimate $i=88.1\pm0.5$ deg. This relies on analysis of the ISS over both frequency and time and provides a model for the location, anisotropy, turbulence level and transverse phase gradient of the IISM. We find that the IISM can be well-modeled during each observation, typically of a few orbital periods, but its turbulence level and mean velocity vary significantly over the 18 months. △ Less

Submitted 3 April, 2014; originally announced April 2014.

Comments: Accepted for publication in Astrophysical Journal

Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

Authors: R. M. Shannon, V. Ravi, W. A. Coles, G. Hobbs, M. J. Keith, R. N. Manchester, J. S. B. Wyithe, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, J. Khoo, Y. Levin, S. Osłowski, J. M. Sarkissian, W. van Straten, J. P. W. Verbiest, J-B. Wang

Abstract: The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% conf… ▽ More The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be $Ω_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability. △ Less

Submitted 16 October, 2013; originally announced October 2013.

Comments: Published in Science 18 October 2013

Interplanetary spacecraft navigation using pulsars

Authors: X. P. Deng, G. Hobbs, X. P. You, M. T. Li, M. J. Keith, R. M. Shannon, W. Coles, R. N. Manchester, J. H. Zheng, X. Z. Yu, D. Gao, X. Wu, D. Chen

Abstract: We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond… ▽ More We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s. △ Less

Submitted 20 July, 2013; originally announced July 2013.

Comments: accepted by Advances in Space Research

Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing

Authors: M. J. Keith, W. Coles, R. M. Shannon, G. B. Hobbs, R. N. Manchester, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. J. Champion, A. Chaudhary, A. W. Hotan, J. Khoo, J. Kocz, S. Oslowski, V. Ravi, J. E. Reynolds, J. Sarkissian, W. van Straten, D. R. B. Yardley

Abstract: Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the det… ▽ More Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 16 pages, 8 figures. Accepted for publication in MNRAS

The Parkes Pulsar Timing Array Project

Authors: R. N. Manchester, G. Hobbs, M. Bailes, W. A. Coles, W. van Straten, M. J. Keith, R. M. Shannon, N. D. R. Bhat, A. Brown, S. G. Burke-Spolaor, D. J. Champion, A. Chaudhary, R. T. Edwards, G. Hampson, A. W. Hotan, A. Jameson, F. A. Jenet, M. J. Kesteven, J. Khoo, J. Kocz, K. Maciesiak, S. Oslowski, V. Ravi, J. R. Reynolds, J. M. Sarkissian , et al. (6 additional authors not shown)

Abstract: A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA conc… ▽ More A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars. △ Less

Submitted 23 October, 2012; originally announced October 2012.

Comments: 36 pages, 13 figures, accepted by Publications of the Astronomical Society of Australia (PASA)

Development of a pulsar-based timescale

Authors: G. Hobbs, W. Coles, R. N. Manchester, M. J. Keith, R. M. Shannon, D. Chen, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. Champion, A. Chaudhary, A. Hotan, J. Khoo, J. Kocz, Y. Levin, S. Oslowski, B. Preisig, V. Ravi, J. E. Reynolds, J. Sarkissian, W. van Straten, J. P. W. Verbiest, D. Yardley, X. P. You

Abstract: Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic… ▽ More Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale. △ Less

Submitted 17 August, 2012; originally announced August 2012.

Comments: Accepted for publication in MNRAS

Optimal Interpolation and Prediction in Pulsar Timing

Authors: X. P. Deng, W. Coles, G. Hobbs, M. J. Keith, R. N. Manchester, R. M. Shannon, J. H. Zheng

Abstract: For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of t… ▽ More For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty. △ Less

Submitted 27 April, 2012; originally announced April 2012.

Comments: 9 pages, accepted by MNRAS

Measurement of the electron density and magnetic field of the solar wind using millisecond pulsars

Authors: X. P. You, W. A. Coles, G. B. Hobbs, R. N. Manchester

Abstract: The magnetic field of the solar wind near the Sun is very difficult to measure directly. Measurements of Faraday rotation of linearly polarized radio sources occulted by the solar wind provide a unique opportunity to estimate this magnetic field, and the technique has been widely used in the past. However Faraday rotation is a path integral of the product of electron density and the projection of… ▽ More The magnetic field of the solar wind near the Sun is very difficult to measure directly. Measurements of Faraday rotation of linearly polarized radio sources occulted by the solar wind provide a unique opportunity to estimate this magnetic field, and the technique has been widely used in the past. However Faraday rotation is a path integral of the product of electron density and the projection of the magnetic field on the path. The electron density near the Sun can be measured by several methods, but it is quite variable. Here we show that it is possible to measure the path integrated electron density and the Faraday rotation simultaneously at 6-10 $R_\odot$ using millisecond pulsars as the linearly polarized radio source. By analyzing the Faraday rotation measurements with and without the simultaneous electron density observations we show that these observations significantly improve the accuracy of the magnetic field estimates. △ Less

Submitted 10 February, 2012; originally announced February 2012.

Comments: 6 pages, 5 figures, accepted by MNRAS

Pulsar timing analysis in the presence of correlated noise

Authors: W. Coles, G. Hobbs, D. J. Champion, R. N. Manchester, J. P. W. Verbiest

Abstract: Pulsar timing observations are usually analysed with least-square-fitting procedures under the assumption that the timing residuals are uncorrelated (statistically "white"). Pulsar observers are well aware that this assumption often breaks down and causes severe errors in estimating the parameters of the timing model and their uncertainties. Ad hoc methods for minimizing these errors have been dev… ▽ More Pulsar timing observations are usually analysed with least-square-fitting procedures under the assumption that the timing residuals are uncorrelated (statistically "white"). Pulsar observers are well aware that this assumption often breaks down and causes severe errors in estimating the parameters of the timing model and their uncertainties. Ad hoc methods for minimizing these errors have been developed, but we show that they are far from optimal. Compensation for temporal correlation can be done optimally if the covariance matrix of the residuals is known using a linear transformation that whitens both the residuals and the timing model. We adopt a transformation based on the Cholesky decomposition of the covariance matrix, but the transformation is not unique. We show how to estimate the covariance matrix with sufficient accuracy to optimize the pulsar timing analysis. We also show how to apply this procedure to estimate the spectrum of any time series with a steep red power-law spectrum, including those with irregular sampling and variable error bars, which are otherwise very difficult to analyse. △ Less

Submitted 26 July, 2011; originally announced July 2011.

Comments: Accepted by MNRAS

Rotation measure variations for 20 millisecond pulsars

Authors: Wenming Yan, R. N. Manchester, Willem van Straten, John Reynolds, George Hobbs, Na Wang, Matthew Bailes, Ramesh Bhat, Sarah Burke-Spolaor, David Champion, Ankur Chaudhary, William Coles, Aidan Hotan, Jonathan Khoo, Stefan Oslowski, John Sarkissian, Daniel Yardley

Abstract: We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the Inter… ▽ More We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight. △ Less

Submitted 21 May, 2011; originally announced May 2011.

Comments: Accepted for publication in Astrophysics & Space Science

Polarization observations of 20 millisecond pulsars

Authors: Wenming Yan, R. N. Manchester, Willem van Straten, John Reynolds, George Hobbs, Na Wang, Matthew Bailes, Ramesh Bhat, Sarah Burke-Spolaor, David Champion, William Coles, Aidan Hotan, Jonathan Khoo, Stefan Oslowski, John Sarkissian, Joris Verbiest, Daniel Yardley

Abstract: Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/no… ▽ More Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: 15pages, 21 figures, 3 tables, accepted for publication in MNRAS

On detection of the stochastic gravitational-wave background using the Parkes pulsar timing array

Authors: D. R. B. Yardley, W. A. Coles, G. B. Hobbs, J. P. W. Verbiest, R. N. Manchester, W. van Straten, F. A. Jenet, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. J. Champion, A. W. Hotan, S. Oslowski, J. E. Reynolds, J. M. Sarkissian

Abstract: We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and n… ▽ More We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and not previously presented in the literature are discussed. These include the effects of quadratic removal, irregular sampling, and variable errors which exacerbate the spectral leakage inherent in estimating the steep red spectrum of the gravitational-wave background. These observations are found to be consistent with the null hypothesis, that no gravitational-wave background is present, with 76 percent confidence. We show that the detection statistic is dominated by the contributions of only a few pulsars because of the inhomogeneity of this data set. The issues of detecting the signature of a gravitational-wave background with future observations are discussed. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: 12 pages, 8 figures, 7 tables, accepted for publication in MNRAS

Developing a pulsar-based timescale

Authors: G. Hobbs, W. Coles, R. Manchester, D. Chen

Abstract: We show how pulsar observations may be used to construct a time standard that is independent of terrestrial time standards. The pulsar time scale provides a method to determine the stability of terrestrial time standards over years to decades. Here, we summarise the method, provide initial results and discuss the possibilities and limitations of our pulsar time scale. We show how pulsar observations may be used to construct a time standard that is independent of terrestrial time standards. The pulsar time scale provides a method to determine the stability of terrestrial time standards over years to decades. Here, we summarise the method, provide initial results and discuss the possibilities and limitations of our pulsar time scale. △ Less

Submitted 23 November, 2010; originally announced November 2010.

Comments: 6 page conference proceedings from Journees conference in Paris (20-22 Sept 2010)

Interferometric imaging with the 32 element Murchison Wide-field Array

Authors: S. M. Ord, D. A. Mitchell, R. B. Wayth, L. J. Greenhill, G. Bernardi, S. Gleadow, R. G. Edgar, M. A. Clark, G. Allen, W. Arcus, L. Benkevitch, J. D. Bowman, F. H. Briggs, J. D. Bunton, S. Burns, R. J. Cappallo, W. A. Coles, B. E. Corey, L. deSouza, S. S. Doeleman, M. Derome, A. Deshpande, D. Emrich, R. Goeke, M. R. Gopalakrishna , et al. (34 additional authors not shown)

Abstract: The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 M… ▽ More The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 MHz to 300 MHz with an instantaneous field of view that is tens of degrees wide and a resolution of a few arcminutes. A 32-station prototype of the MWA has been recently commissioned and a set of observations taken that exercise the whole acquisition and processing pipeline. We present Stokes I, Q, and U images from two ~4 hour integrations of a field 20 degrees wide centered on Pictoris A. These images demonstrate the capacity and stability of a real-time calibration and imaging technique employing the weighted addition of warped snapshots to counter extreme wide field imaging distortions. △ Less

Submitted 8 October, 2010; originally announced October 2010.

Comments: Accepted for publication in PASP. This is the draft before journal typesetting corrections and proofs so does contain formatting and journal style errors, also has with lower quality figures for space requirements

Measuring the mass of solar system planets using pulsar timing

Authors: D. J. Champion, G. B. Hobbs, R. N. Manchester, R. T. Edwards, D. C. Backer, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, W. Coles, P. B. Demorest, R. D. Ferdman, W. M. Folkner, A. W. Hotan, M. Kramer, A. N. Lommen, D. J. Nice, M. B. Purver, J. M. Sarkissian, I. H. Stairs, W. van Straten, J. P. W. Verbiest, D. R. B. Yardley

Abstract: High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period a… ▽ More High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets. △ Less

Submitted 21 August, 2010; originally announced August 2010.

Comments: Accepted for publication in ApJ

Scattering of pulsar radio emission by the interstellar plasma

Authors: W. A. Coles, B. J. Rickett, J. J. Gao, G. Hobbs, J. P. W. Verbiest

Abstract: We present simulations of scattering phenomena which are important in pulsar observations, but which are analytically intractable. The simulation code, which has also been used for solar wind and atmospheric scattering problems, is available from the authors. These simulations reveal an unexpectedly important role of dispersion in combination with refraction. We demonstrate the effect of analyzing… ▽ More We present simulations of scattering phenomena which are important in pulsar observations, but which are analytically intractable. The simulation code, which has also been used for solar wind and atmospheric scattering problems, is available from the authors. These simulations reveal an unexpectedly important role of dispersion in combination with refraction. We demonstrate the effect of analyzing observations which are shorter than the refractive scale. We examine time-of-arrival fluctuations in detail: showing their correlation with intensity and dispersion measure; providing a heuristic model from which one can estimate their contribution to pulsar timing observations; and showing that much of the effect can be corrected making use of measured intensity and dispersion. Finally, we analyze observations of the millisecond pulsar J0437$-$4715, made with the Parkes radio telescope, that show timing fluctuations which are correlated with intensity. We demonstrate that these timing fluctuations can be corrected, but we find that they are much larger than would be expected from scattering in a homogeneous turbulent plasma with isotropic density fluctuations. We do not have an explanation for these timing fluctuations. △ Less

Submitted 26 May, 2010; originally announced May 2010.

The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves

Authors: D. R. B. Yardley, G. B. Hobbs, F. A. Jenet, J. P. W. Verbiest, Z. L. Wen, R. N. Manchester, W. A. Coles, W. van Straten, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. J. Champion, A. W. Hotan, J. M. Sarkissian

Abstract: We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complem… ▽ More We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed. △ Less

Submitted 25 August, 2010; v1 submitted 10 May, 2010; originally announced May 2010.

Comments: fixed error in equation (4). [13 pages, 6 figures, 1 table, published in MNRAS]

Status Update of the Parkes Pulsar Timing Array

Authors: J. P. W. Verbiest, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. J. Champion, W. Coles, G. B. Hobbs, A. W. Hotan, F. Jenet, J. Khoo, K. J. Lee, A. Lommen, R. N. Manchester, S. Oslowski, J. Reynolds, J. Sarkissian, W. van Straten, D. R. B. Yardley, X. P. You

Abstract: The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a succ… ▽ More The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers. △ Less

Submitted 14 December, 2009; originally announced December 2009.

Comments: 10 pages, 3 figures, Amaldi 8 conference proceedings, accepted by Classical & Quantum Gravity

Journal ref: Class.Quant.Grav.27:084015,2010