Search

We report the real-time discovery of a fast radio burst (FRB 131104) with the Parkes radio telescope in a targeted observation of the Carina dwarf spheroidal galaxy. The dispersion measure of the burst is 779 cm$^{-3}$ pc, exceeding predictions for the maximum line-of-sight Galactic contribution by a factor of 11. The temporal structure of the burst is characterized by an exponential scattering tail with a timescale of 2.0$^{+0.8}_{-0.5}$ ms at 1582 MHz that scales as frequency to the power $-$4.4$^{+1.6}_{-1.8}$ (all uncertainties represent 95% confidence intervals). We bound the intrinsic pulse width to be $<0.64$ ms due to dispersion smearing across a single spectrometer channel. Searches in 78 hours of follow-up observations with the Parkes telescope reveal no additional sporadic emission and no evidence for associated periodic radio emission. We hypothesize that the burst is associated with the Carina dwarf galaxy. Follow-up observations at other wavelengths are necessary to test this hypothesis., Comment: 5 pages, 3 figures, 1 table, accepted to ApJ Letters. Updated to match published version

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., Comment: Accepted by MNRAS

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}$., Comment: A few typos corrected, to match the published version, 13 pages, 11 figures

Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range $5.1\times10^{-16}

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., Comment: 21 pages, 23 figures, accepted by MNRAS

We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought., Comment: 15 pages, 7 figures, 1 table, accepted by MNRAS

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 ${\Omega}_{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., Comment: Published in Science 18 October 2013

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., Comment: 16 pages, 8 figures. Accepted for publication in MNRAS

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes \Delta\nu_g/\nu range between 10^-10 and 10^-5, where \nu = 1/P is the pulse frequency. We confirm that the distribution of \Delta\nu_g/\nu is bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dot\nu in the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dot\nu. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddot\nu values characterising the linear changes in \dot\nu are almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dot\nu| and the inter-glitch interval, as predicted by vortex-creep models., Comment: Accepted for publication in MNRAS

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., Comment: 36 pages, 13 figures, accepted by Publications of the Astronomical Society of Australia (PASA)