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 ${\mu}$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.

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes.

The anomalous X-ray pulsar XTE J1810$-$197 was the first magnetar found to emit pulsed radio emission. After spending almost a decade in a quiescent, radio-silent state, the magnetar was reported to have undergone a radio outburst in December, 2018. We observed radio pulsations from XTE J1810$-$197 during this early phase of its radio revival using the Ultra-Wideband Low receiver system of the Parkes radio telescope, obtaining wideband (704 MHz to 4032 MHz) polarization pulse profiles, single pulses and flux density measurements. Dramatic changes in polarization and rapid variations of the position angle of linear polarization across the main pulse and in time have been observed. The pulse profile exhibits similar structures throughout our three observations (over a week time scale), displaying a small amount of profile evolution in terms of polarization and pulse width across the wideband. We measured a flat radio spectrum across the band with a positive spectral index, in addition to small levels of flux and spectral index variability across our observing span. The observed wideband polarization properties are significantly different compared to those taken after the 2003 outburst, and therefore provide new information about the origin of radio emission.

Breakthrough Listen is a ten-year initiative to search for signatures of technologies created by extraterrestrial civilizations at radio and optical wavelengths. Here, we detail the digital data recording system deployed for Breakthrough Listen observations at the 64-m aperture CSIRO Parkes Telescope in New South Wales, Australia. The recording system currently implements two recording modes: a dual-polarization, 1.125 GHz bandwidth mode for single beam observations, and a 26-input, 308-MHz bandwidth mode for the 21-cm multibeam receiver. The system is also designed to support a 3 GHz single-beam mode for the forthcoming Parkes ultra-wideband feed. In this paper, we present details of the system architecture, provide an overview of hardware and software, and present initial performance results.

We present the spin frequency and profile evolution of the radio pulsar J1119$-$6127 following magnetar-like X-ray bursts from the system in 2016 July. Using data from the Parkes radio telescope, we observe a smooth and fast spin-down process subsequent to the X-ray bursts resulting in a net change in the pulsar rotational frequency of $\Delta\nu\approx-4\times10^{-4}$\,Hz. During the transition, a net spin-down rate increase of $\Delta\dot\nu\approx-1\times10^{-10}$\,Hz\,s$^{-1}$ is observed, followed by a return of $\dot{\nu}$ to its original value. In addition, the radio pulsations disappeared after the X-ray bursts and reappeared about two weeks later with the flux density at 1.4\,GHz increased by a factor of five. The flux density then decreased and undershot the normal flux density followed by a slow recovery back to normal. The pulsar's integrated profile underwent dramatic and short-term changes in total intensity, polarization and position angle. Despite the complex evolution, we observe correlations between the spin-down rate, pulse profile shape and radio flux density. Strong single pulses have been detected after the X-ray bursts with their energy distributions evolving with time. The peculiar but smooth spin frequency evolution of PSR~J1119$-$6127 accompanied by systematic pulse profile and flux density changes are most likely to be a result of either reconfiguration of the surface magnetic fields or particle winds triggered by the X-ray bursts. The recovery of spin-down rate and pulse profile to normal provides us the best case to study the connection between high magnetic-field pulsars and magnetars.

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.

Multi-decade observing campaigns of the globular clusters 47 Tucanae and M15 have led to an outstanding number of discoveries. Here, we report on the latest results of the long-term observations of the pulsars in these two clusters. For most of the pulsars in 47 Tucanae we have measured, among other things, their higher-order spin period derivatives, which have in turn provided stringent constraints on the physical parameters of the cluster, such as its distance and gravitational potential. For M15, we have studied the relativistic spin precession effect in PSR B2127+11C. We have used full-Stokes observations to model the precession effect, and to constrain the system geometry. We find that the visible beam of the pulsar is swiftly moving away from our line of sight and may very soon become undetectable. On the other hand, we expect to see the opposite emission beam sometime between 2041 and 2053.

Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio pulsar PSR J1119-6127 at X-ray, with XMM-Newton & NuSTAR, and at radio energies with Parkes radio telescope, during a period of magnetar-like bursts. The rotationally powered radio emission shuts off coincident with the occurrence of multiple X-ray bursts, and recovers on a time scale of ~70 seconds. These observations of related radio and X-ray phenomena further solidify the connection between radio pulsars and magnetars, and suggest that the pair plasma produced in bursts can disrupt the acceleration mechanism of radio emitting particles.

This paper is the second in a series where we report the results of the long-term timing of the millisecond pulsars (MSPs) in 47 Tucanae with the Parkes 64-m radio telescope. We obtain improved timing parameters that provide additional information for studies of the cluster dynamics: a) the pulsar proper motions yield an estimate of the proper motion of the cluster as a whole ($\mu_{\alpha}\, = \, 5.00\, \pm \, 0.14\, \rm mas \, yr^{-1}$, $\mu_{\delta}\, = \, -2.84\, \pm \, 0.12\, \rm mas \, yr^{-1}$) and the motion of the pulsars relative to each other. b) We measure the second spin-period derivatives caused by the change of the pulsar line-of-sight accelerations; 47 Tuc H, U and possibly J are being affected by nearby objects. c) For ten binary systems we now measure changes in the orbital period caused by their acceleration in the gravitational field of the cluster. From all these measurements, we derive a cluster distance no smaller than $\sim\,$4.69 kpc and show that the characteristics of these MSPs are very similar to their counterparts in the Galactic disk. We find no evidence in favour of an intermediate mass black hole at the centre of the cluster. Finally, we describe the orbital behaviour of the four "black widow" systems. Two of them, 47 Tuc J and O, exhibit orbital variability similar to that observed in other such systems, while for 47 Tuc I and R the orbits seem to be remarkably stable. It appears, therefore, that not all "black widows" have unpredictable orbital behaviour.

We have observed the Vela pulsar for one year using a Phased Array Feed (PAF) receiver on the 12-metre antenna of the Parkes Test-Bed Facility. These observations have allowed us to investigate the stability of the PAF beam-weights over time, to demonstrate that pulsars can be timed over long periods using PAF technology and to detect and study the most recent glitch event that occurred on 12 December 2016. The beam-weights are shown to be stable to 1% on time scales on the order of three weeks. We discuss the implications of this for monitoring pulsars using PAFs on single dish telescopes.

We report on 2.4 yr of radio timing measurements of the magnetar PSR J1622$-$4950 using the Parkes telescope, between 2011 November and 2014 March. During this period the torque on the neutron star (inferred from the rotational frequency derivative) varied greatly, though much less erratically than in the 2 yr following its discovery in 2009. During the last year of our measurements the frequency derivative decreased in magnitude monotonically by 20\%, to a value of $-1.3\times10^{-13}$ s$^{-2}$, a factor of 8 smaller than when discovered. The flux density continued to vary greatly during our monitoring through 2014 March, reaching a relatively steady low level after late 2012. The pulse profile varied secularly on a similar timescale as the flux density and torque. A relatively rapid transition in all three properties is evident in early 2013. After PSR J1622$-$4950 was detected in all of our 87 observations up to 2014 March, we did not detect the magnetar in our resumed monitoring starting in 2015 January and have not detected it in any of the 30 observations done through 2016 September.

For the past couple of decades, the Parkes radio telescope has been regularly observing the millisecond pulsars in 47 Tucanae (47 Tuc). This long-term timing program was designed to address a wide range of scientific issues related to these pulsars and the globular cluster where they are located. In this paper, the first of a series, we address one of these objectives: the characterization of four previously known binary pulsars for which no precise orbital parameters were known, namely 47 Tuc P, V, W and X (pulsars 47 Tuc R and Y are discussed elsewhere). We determined the previously unknown orbital parameters of 47 Tuc V and X and greatly improved those of 47 Tuc P and W. For pulsars W and X we obtained, for the first time, full coherent timing solutions across the whole data span, which allowed a much more detailed characterization of these systems. 47 Tuc W, a well-known tight eclipsing binary pulsar, exhibits a large orbital period variability, as expected for a system of its class. 47 Tuc X turns out to be in a wide, extremely circular, 10.9-day long binary orbit and its position is ~3.8 arcmin away from the cluster center, more than three times the distance of any other pulsar in 47 Tuc. These characteristics make 47 Tuc X a very different object with respect to the other pulsars of the cluster.

We report on timing, flux density, and polarimetric observations of the transient magnetar and 5.54 s radio pulsar XTE J1810-197 using the GBT, Nancay, and Parkes radio telescopes beginning in early 2006, until its sudden disappearance as a radio source in late 2008. Repeated observations through 2016 have not detected radio pulsations again. The torque on the neutron star, as inferred from its rotation frequency derivative f-dot, decreased in an unsteady manner by a factor of 3 in the first year of radio monitoring. In contrast, during its final year as a detectable radio source, the torque decreased steadily by only 9%. The period-averaged flux density, after decreasing by a factor of 20 during the first 10 months of radio monitoring, remained steady in the next 22 months, at an average of 0.7+/-0.3 mJy at 1.4 GHz, while still showing day-to-day fluctuations by factors of a few. There is evidence that during this last phase of radio activity the magnetar had a steep radio spectrum, in contrast to earlier behavior. There was no secular decrease that presaged its radio demise. During this time the pulse profile continued to display large variations, and polarimetry indicates that the magnetic geometry remained consistent with that of earlier times. We supplement these results with X-ray timing of the pulsar from its outburst in 2003 up to 2014. For the first 4 years, XTE J1810-197 experienced non-monotonic excursions in f-dot by at least a factor of 8. But since 2007, its f-dot has remained relatively stable near its minimum observed value. The only apparent event in the X-ray record that is possibly contemporaneous with the radio shut-down is a decrease of ~20% in the hot-spot flux in 2008-2009, to a stable, minimum value. However, the permanence of the high-amplitude, thermal X-ray pulse, even after the radio demise, implies continuing magnetar activity.

In a search of the unidentified Fermi gamma-ray source 3FGL J1417.5-4402 with the Parkes radio telescope, we discovered PSR J1417-4402, a 2.66 ms pulsar having the same 5.4 day orbital period as the optical and X-ray binary identified by Strader et al. The existence of radio pulsations implies that the neutron star is currently not accreting. Substantial outflows from the companion render the radio pulsar undetectable for more than half of the orbit, and may contribute to the observed Halpha emission. Our initial pulsar observations, together with the optically inferred orbit and inclination, imply a mass ratio of 0.171+/-0.002, a companion mass of M_2=0.33+/-0.03 Msun, and a neutron star mass in the range 1.77<M_1<2.13 Msun. However, there remains a discrepancy between the distance of 4.4 kpc inferred from the optical properties of the companion and the smaller radio dispersion measure distance of 1.6 kpc. The smaller distance would reduce the inferred Roche-lobe filling factor, increase the inferred inclination angle, and decrease the masses. As a wide binary, PSR J1417-4402 differs from the radio-eclipsing black widow and redback pulsars being discovered in large numbers by Fermi. It is probably a system that began mass transfer onto the neutron star after the companion star left the main sequence. The companion should end its evolution as a He white dwarf in a 6-20 day orbit, i.e., as a typical binary millisecond pulsar companion.

In a search with the Parkes radio telescope of 56 unidentified Fermi-LAT gamma-ray sources, we have detected 11 millisecond pulsars (MSPs), 10 of them discoveries, of which five were reported in Kerr et al. (2012). We did not detect radio pulsations from another six pulsars now known in these sources. We describe the completed survey, which included multiple observations of many targets done to minimize the impact of interstellar scintillation, acceleration effects in binary systems, and eclipses. We consider that 23 of the 39 remaining sources may still be viable pulsar candidates. We present timing solutions and polarimetry for five of the MSPs, and gamma-ray pulsations for PSR J1903-7051 (pulsations for five others were reported in the second Fermi-LAT catalog of gamma-ray pulsars). Two of the new MSPs are isolated and five are in >1 d circular orbits with 0.2-0.3 Msun presumed white dwarf companions. PSR J0955-6150, in a 24 d orbit with a ~0.25 Msun companion but eccentricity of 0.11, belongs to a recently identified class of eccentric MSPs. PSR J1036-8317 is in an 8 hr binary with a >0.14 Msun companion that is probably a white dwarf. PSR J1946-5403 is in a 3 hr orbit with a >0.02 Msun companion with no evidence of radio eclipses.

"Perytons" are millisecond-duration transients of terrestrial origin, whose frequency-swept emission mimics the dispersion of an astrophysical pulse that has propagated through tenuous cold plasma. In fact, their similarity to FRB 010724 had previously cast a shadow over the interpretation of "fast radio bursts," which otherwise appear to be of extragalactic origin. Until now, the physical origin of the dispersion-mimicking perytons had remained a mystery. We have identified strong out-of-band emission at 2.3--2.5 GHz associated with several peryton events. Subsequent tests revealed that a peryton can be generated at 1.4 GHz when a microwave oven door is opened prematurely and the telescope is at an appropriate relative angle. Radio emission escaping from microwave ovens during the magnetron shut-down phase neatly explain all of the observed properties of the peryton signals. Now that the peryton source has been identified, we furthermore demonstrate that the microwaves on site could not have caused FRB 010724. This and other distinct observational differences show that FRBs are excellent candidates for genuine extragalactic transients.

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 $\lambda^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.

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.

It is thought that neutron stars in low-mass binary systems can accrete matter and angular momentum from the companion star and be spun-up to millisecond rotational periods. During the accretion stage, the system is called a low-mass X-ray binary, and bright X-ray emission is observed. When the rate of mass transfer decreases in the later evolutionary stages, these binaries host a radio millisecond pulsar whose emission is powered by the neutron star's rotating magnetic field. This evolutionary model is supported by the detection of millisecond X-ray pulsations from several accreting neutron stars and also by the evidence for a past accretion disc in a rotation-powered millisecond pulsar. It has been proposed that a rotation-powered pulsar may temporarily switch on during periods of low mass inflow in some such systems. Only indirect evidence for this transition has hitherto been observed. Here we report observations of accretion-powered, millisecond X-ray pulsations from a neutron star previously seen as a rotation-powered radio pulsar. Within a few days after a month-long X-ray outburst, radio pulses were again detected. This not only shows the evolutionary link between accretion and rotation-powered millisecond pulsars, but also that some systems can swing between the two states on very short timescales.

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.

We have conducted a new search for radio pulsars in compact binary systems in the Parkes multi-beam pulsar survey (PMPS) data, employing novel methods to remove the Doppler modulation from binary motion. This has yielded unparalleled sensitivity to pulsars in compact binaries. The required computation time of approximately 17000 CPU core years was provided by the distributed volunteer computing project Einstein@Home, which has a sustained computing power of about 1 PFlop/s. We discovered 24 new pulsars in our search, of which 18 were isolated pulsars, and six were members of binary systems. Despite the wide filterbank channels and relatively slow sampling time of the PMPS data, we found pulsars with very large ratios of dispersion measure (DM) to spin period. Among those is PSR J1748-3009, the millisecond pulsar with the highest known DM (approximately 420 pc/cc). We also discovered PSR J1840-0643, which is in a binary system with an orbital period of 937 days, the fourth largest known. The new pulsar J1750-2536 likely belongs to the rare class of intermediate-mass binary pulsars. Three of the isolated pulsars show long-term nulling or intermittency in their emission, further increasing this growing family. Our discoveries demonstrate the value of distributed volunteer computing for data-driven astronomy and the importance of applying new analysis methods to extensively searched data.

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.

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.

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.

The past decade has seen Parkes once again involved in a wide range of space tracking activities that have added to its illustrious legacy. This contribution is a personal recollection of those tracking efforts - both real and celluloid. We begin in a light-hearted vein with some behind-the-scenes views of the popular film, "The DISH", and then turn to more serious contributions; discussing the vital role of the telescope in alleviating the great "traffic jam" at Mars in 2003/04 and salvaging the Doppler Wind Experiment as the Huygens probe descended though the atmosphere of Saturn's largest moon, Titan, in mid-decade. We cap off the decade with a discussion of the search for the missing Apollo 11 slow-scan TV tapes.

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.

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.

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified Fermi-LAT sources in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived from timing observations of PSR J0101-6422 (P =2.57 ms, DM=12 pc cm-3), we have detected {\gamma}-ray pulsations and measured its proper motion. Its {\gamma}-ray spectrum (a power law of {\Gamma} = 0.9 with a cutoff at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and {\gamma}-ray light curves challenge simple geometric models of emission. The high success rate of this survey-enabled by selecting {\gamma}-ray sources based on their detailed spectral characteristics-and other similarly successful searches indicate that a substantial fraction of the local population of MSPs may soon be known.

The Parkes pulsar data archive currently provides access to 144044 data files obtained from observations carried out at the Parkes observatory since the year 1991. Around 10^5 files are from surveys of the sky, the remainder are observations of 775 individual pulsars and their corresponding calibration signals. Survey observations are included from the Parkes 70cm and the Swinburne Intermediate Latitude surveys. Individual pulsar observations are included from young pulsar timing projects, the Parkes Pulsar Timing Array and from the PULSE@Parkes outreach program. The data files and access methods are compatible with Virtual Observatory protocols. This paper describes the data currently stored in the archive and presents ways in which these data can be searched and downloaded.

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.

SGR 1833-0832 was discovered on 2010 March 19 thanks to the Swift detection of a short hard X-ray burst and follow-up X-ray observations. Since then, it was repeatedly observed with Swift, Rossi X-ray Timing Explorer, and XMM-Newton. Using these data, which span about 225 days, we studied the long-term spectral and timing characteristics of SGR 1833-0832. We found evidence for diffuse emission surrounding SGR 1833-0832, which is most likely a halo produced by the scattering of the point source X-ray radiation by dust along the line of sight, and we show that the source X-ray spectrum is well described by an absorbed blackbody, with temperature kT=1.2 keV and absorbing column nH=(10.4+/-0.2)E22 cm^-2, while different or more complex models are disfavoured. The source persistent X-ray emission remained fairly constant at about 3.7E-12 erg/cm^2/s for the first 20 days after the onset of the bursting episode, then it faded by a factor 40 in the subsequent 140 days, following a power-law trend with index alpha=-0.5. We obtained a phase-coherent timing solution with the longest baseline (225 days) to date for this source which, besides period P=7.5654084(4) s and period derivative dP/dt=3.5(3)E-12 s/s, includes higher order period derivatives. We also report on our search of the counterpart to the SGR at radio frequencies using the Australia Telescope Compact Array and the Parkes radio telescope. No evidence for radio emission was found, down to flux densities of 0.9 mJy (at 1.5 GHz) and 0.09 mJy (at 1.4 GHz) for the continuum and pulsed emissions, respectively, consistently with other observations at different epochs.

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.

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.

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.

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.

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.

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

The faint radio supernova remnant SNR G315.9-0.0 is notable for a long and thin trail that extends outward perpendicular from the edge of its approximately circular shell. In a search with the Parkes telescope we have found a young and energetic pulsar that is located at the tip of this collimated linear structure. PSR J1437-5959 has period P = 61 ms, characteristic age tau_c = 114 kyr, and spin-down luminosity dE/dt = 1.4e36 erg/s. It is very faint, with a flux density at 1.4 GHz of about 75 uJy. From its dispersion measure of 549 pc/cc, we infer d ~ 8 kpc. At this distance and for an age comparable to tau_c, the implied pulsar velocity in the plane of the sky is V_t = 300 km/s for a birth at the center of the SNR, although it is possible that the SNR/pulsar system is younger than tau_c and that V_t > 300 km/s. The highly collimated linear feature is evidently the pulsar wind trail left from the supersonic passage of PSR J1437-5959 through the interstellar medium surrounding SNR G315.9-0.0.

Analysis of high-precision timing observations of an array of approx. 20 millisecond pulsars (a so-called "timing array") may ultimately result in the detection of a stochastic gravitational-wave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a five-year timescale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within five years to a decade are good.

The first direct detection of gravitational waves may be made through observations of pulsars. The principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low frequency gravitational waves (f ~ 10^-9 to 10^-8 Hz). Such waves are expected to be caused by coalescing supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable signal could have been produced in the inflationary era or by cosmic strings. In this paper we review the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection such as ground- and space-based interferometer systems.

Analysis of ten years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, Pbdot, determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 Solar System ephemerides. The precise measurement of Pbdot allows a limit on the variation of Newton's gravitational constant, |Gdot/G| < 23 x 10^{-12} 1/yr. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a(Sun)/c| < 1.5 x 10^{-18} 1/s at 95% confidence, and derive a pulsar mass, m(psr) = 1.76 +/- 0.20 M, one of the highest estimates so far obtained.

We present a novel experiment with the capacity to independently measure both the electron density and the magnetic field of the solar corona. We achieve this through measurement of the excess Faraday rotation due to propagation of the polarised emission from a number of pulsars through the magnetic field of the solar corona. This method yields independent measures of the integrated electron density, via dispersion of the pulsed signal and the magnetic field, via the amount of Faraday rotation. In principle this allows the determination of the integrated magnetic field through the solar corona along many lines of sight without any assumptions regarding the electron density distribution. We present a detection of an increase in the rotation measure of the pulsar J1801$-$2304 of approximately 160 \rad at an elongation of 0.95$^\circ$ from the centre of the solar disk. This corresponds to a lower limit of the magnetic field strength along this line of sight of $> 393\mu\mathrm{G}$. The lack of precision in the integrated electron density measurement restricts this result to a limit, but application of coronal plasma models can further constrain this to approximately 20mG, along a path passing 2.5 solar radii from the solar limb. Which is consistent with predictions obtained using extensions to the Source Surface models published by Wilcox Solar Observatory

We present an analysis of the variations seen in the dispersion measures (DMs) of 20 millisecond pulsars observed as part of the Parkes Pulsar Timing Array project. We carry out a statistically rigorous structure function analysis for each pulsar and show that the variations seen for most pulsars are consistent with those expected for an interstellar medium characterised by a Kolmogorov turbulence spectrum. The structure functions for PSRs J1045-4509 and J1909-3744 provide the first clear evidence for a large inner scale, possibly due to ion-neutral damping. We also show the effect of the solar wind on the DMs and show that the simple models presently implemented into pulsar timing packages cannot reliably correct for this effect. For the first time we clearly show how DM variations affect pulsar timing residuals and how they can be corrected in order to obtain the highest possible timing precision. Even with our presently limited data span, the residuals (and all parameters derived from the timing) for six of our pulsars have been significantly improved by correcting for the DM variations.

The double pulsar system, PSR J0737-3039A/B, is unique in that both neutron stars are detectable as radio pulsars. This, combined with significantly higher mean orbital velocities and accelerations when compared to other binary pulsars, suggested that the system would become the best available testbed for general relativity and alternative theories of gravity in the strong-field regime. Here we report on precision timing observations taken over the 2.5 years since its discovery and present four independent strong-field tests of general relativity. Use of the theory-independent mass ratio of the two stars makes these tests uniquely different from earlier studies. By measuring relativistic corrections to the Keplerian description of the orbital motion, we find that the ``post-Keplerian'' parameter s agrees with the value predicted by Einstein's theory of general relativity within an uncertainty of 0.05%, the most precise test yet obtained. We also show that the transverse velocity of the system's center of mass is extremely small. Combined with the system's location near the Sun, this result suggests that future tests of gravitational theories with the double pulsar will supersede the best current Solar-system tests. It also implies that the second-born pulsar may have formed differently to the usually assumed core-collapse of a helium star.

Exploiting a five-year span of data, we present improved timing solutions for the five millisecond pulsars known in the globular cluster NGC 6752. They include proper motion determinations for the two outermost pulsars in the cluster, PSR J1910-5959A and PSR J1910-5959C. The values of the proper motions are in agreement with each other within current uncertainties, but do not match (at 4 sigma and 2 sigma level respectively) with the value of the proper motion of the entire globular cluster derived in the optical band. Implications of these results for the cluster membership of the two pulsars are investigated. Prospects for the detection of the Shapiro delay in the binary system J1910-5959A are also discussed.

Using a statistically rigorous analysis method, we place limits on the existence of an isotropic stochastic gravitational wave background using pulsar timing observations. We consider backgrounds whose characteristic strain spectra may be described as a power-law dependence with frequency. Such backgrounds include an astrophysical background produced by coalescing supermassive black-hole binary systems and cosmological backgrounds due to relic gravitational waves and cosmic strings. Using the best available data, we obtain an upper limit on the energy density per unit logarithmic frequency interval of \Omega^{\rm SMBH}_g(1/8yr) h^2 <= 1.9 x 10^{-8} for an astrophysical background which is five times more stringent than the earlier Kaspi et al. (1994) limit of 1.1 x 10^{-7}. We also provide limits on a background due to relic gravitational waves and cosmic strings of \Omega^{\rm relic}_g(1/8yr) h^2 <= 2.0 x 10^{-8} and \Omega^{\rm cs}_g(1/8yr) h^2 <= 1.9 x 10^{-8} respectively. All of the quoted upper limits correspond to a 0.1% false alarm rate together with a 95% detection rate. We discuss the physical implications of these results and highlight the future possibilities of the Parkes Pulsar Timing Array project. We find that our current results can 1) constrain the merger rate of supermassive binary black hole systems at high red shift, 2) rule out some relationships between the black hole mass and the galactic halo mass, 3) constrain the rate of expansion in the inflationary era and 4) provide an upper bound on the dimensionless tension of a cosmic string background.

Anomalous X-ray pulsars (AXPs) are slowly rotating neutron stars with very bright and highly variable X-ray emission that are believed to be powered by ultra-strong magnetic fields of >1e14 G, according to the 'magnetar' model. The radio pulsations that have been observed from more than 1,700 neutron stars with weaker magnetic fields have never been detected from any of the dozen known magnetars. The X-ray pulsar XTE J1810-197 was revealed (in 2003) as the first AXP with transient emission when its luminosity increased 100-fold from the quiescent level; a coincident radio source of unknown origin was detected one year later. Here we show that XTE J1810-197 emits bright, narrow, highly linearly polarized radio pulses, observed at every rotation, thereby establishing that magnetars can be radio pulsars. There is no evidence of radio emission before the 2003 X-ray outburst (unlike ordinary pulsars, which emit radio pulses all the time), and the flux varies from day to day. The flux at all radio frequencies is approximately equal -- and at >20 GHz XTE J1810-197 is currently the brightest neutron star known. These observations link magnetars to ordinary radio pulsars, rule out alternative accretion models for AXPs, and provide a new window into the coronae of magnetars.

We present radio observations of comet 9P/Tempel 1 associated with the Deep Impact spacecraft collision of 2005 July 4. Weak 18-cm OH emission was detected with the Parkes 64-m telescope, in data averaged over July 4 to 6, at a level of 12 +/- 3 mJy km/s, corresponding to OH production rate 2.8 x 10^{28} molecules/second (Despois et al. inversion model, or 1.0 x 10^{28} /s for the Schleicher & A'Hearn model). We did not detect the HCN 1-0 line with the Mopra 22-m telescope over the period July 2 to 6. The 3 sigma limit of 0.06 K km/s for HCN on July 4 after the impact gives the limit to the HCN production rate of < 1.8 x 10^{25} /s. We did not detect the HCN 1-0 line, 6.7 GHz CH_3OH line or 3.4-mm continuum with the Australia Telescope Compact Array (ATCA) on July 4, giving further limits on any small-scale structure due to an outburst. The 3 sigma limit on HCN emission of 2.5 K km/s from the ATCA around impact corresponds to limit < 4 x 10^{29} HCN molecules released by the impact.

Analysis of 20 months of observations at the Parkes radio telescope shows secular changes in the pulsed emission from J0737-3039B, the 2.77 s pulsar of the double-pulsar system. Pulse profiles are becoming single-peaked in both bright phases of the orbital modulation although there is no clear variation in overall pulse width. The shape of the orbital modulation is also varying systematically, with both bright phases shrinking in longitude by ~ 7 deg/yr. However, the combined span of the two bright phases is relatively constant and together they are shifting to higher longitudes at a rate of ~ 3 deg/yr. We discuss the possible contributions of geodetic precession and periastron advance to the observed variations.

This first ever double pulsar system consists of two pulsars orbiting the common center of mass in a slightly eccentric orbit of only 2.4-hr duration. The pair of pulsars with pulse periods of 22 ms and 2.8 sec, respectively, confirms the long-proposed recycling theory for millisecond pulsars and provides an exciting opportunity to study the works of pulsar magnetospheres by a very fortunate geometrical alignment of the orbit relative to our line-of-sight. In particular, this binary system represents a truly unique laboratory for relativistic gravitational physics. This contribution serves as an update on the currently obtained results and their consequences for the test of general relativity in the strong-field regime. A complete and more up-to-date report of the timing results will be presented elsewhere shortly.

General relativity predicts that the spin axes of the pulsars in the double-pulsar system (PSR J0737-3039A/B) will precess rapidly, in general leading to a change in the observed pulse profiles. We have observed this system over a one-year interval using the Parkes 64-m radio telescope at three frequencies: 680, 1390 and 3030 MHz. These data, combined with the short survey observation made two years earlier, show no evidence for significant changes in the pulse profile of PSR J0737-3039A, the 22-ms pulsar. The limit on variations of the profile 10% width is about 0.5 deg per year. These results imply an angle delta between the pulsar spin axis and the orbit normal of <~ 60 deg, consistent with recent evolutionary studies of the system. Although a wide range of system parameters remain consistent with the data, the model proposed by Jenet & Ransom (2004) can be ruled out. A non-zero ellipticity for the radiation beam gives slightly but not significantly improved fits to the data, so that a circular beam describes the data equally well within the uncertainties.

This manuscript is an updated version of Kalogera et al. (2004) published in ApJ Letters to correct our calculation of the Galactic DNS in-spiral rate. The details of the original erratum submitted to ApJ Letters are given in page 6 of this manuscript. We report on the newly increased event rates due to the recent discovery of the highly relativistic binary pulsar J0737--3039 (Burgay et al. 2003). Using a rigorous statistical method, we present the calculations reported by Burgay et al., which produce a in-spiral rate for Galactic double neutron star (DNS) systems that is higher by a factor of 5-7 compared to estimates made prior to the new discovery. Our method takes into account known pulsar-survey selection effects and biases due to small-number statistics. This rate increase has dramatic implications for gravitational wave detectors. For the initial Laser Interferometer Gravitational-Wave Observatory (LIGO) detectors, the most probable detection rates for DNS in-spirals are one event per 10-630 yr; at 95% confidence, we obtain rates up to one per 3 yr. For the advanced LIGO detectors, the most probable rates are 10-500 events per year. These predictions, for the first time, bring the expectations for DNS detections by the initial LIGO detectors to the astrophysically relevant regime. We also use our models to predict that the large-scale Parkes Multibeam pulsar survey with acceleration searches could detect an average of four binary pulsars similar to those known at present.

PSR J0737-3039A is a millisecond pulsar with a spin period of 22.7 ms included in a double-neutron star system with an orbital period of 2.4 hrs. Its companion has also been detected as a radio pulsar, making this binary the first known double-pulsar system. Its discovery has important implications for relativistic gravity tests, gravitational wave detection and plasma physics. Here we will shortly describe the discovery of the first pulsar in this unique system and present the first results obtained by follow-up studies.

The merger of close binary systems containing two neutron stars should produce a burst of gravitational waves, as predicted by the theory of general relativity. A reliable estimate of the double-neutron-star merger rate in the Galaxy is crucial in order to predict whether current gravity wave detectors will be successful in detecting such bursts. Present estimates of this rate are rather low, because we know of only a few double-neutron-star binaries with merger times less than the age of the Universe. Here we report the discovery of a 22-ms pulsar, PSR J0737-3039, which is a member of a highly relativistic double-neutron-star binary with an orbital period of 2.4 hours. This system will merge in about 85 Myr, a time much shorter than for any other known neutron-star binary. Together with the relatively low radio luminosity of PSR J0737-3039, this timescale implies an order-of-magnitude increase in the predicted merger rate for double-neutron-star systems in our Galaxy (and in the rest of the Universe).

We present rotational and astrometric parameters of three millisecond pulsars located near the center of the globular cluster NGC 6266 (M62) resulting from timing observations with the Parkes radio telescope. Their accelerations toward the cluster center yield values of the cluster central density and mass-to-light ratio consistent with those derived from optical data. The three pulsars are in binary systems. One (spin period P=5.24 ms) is in a 3.5-day orbit around a companion of minimum mass 0.2 Msun. The other two millisecond pulsars (P=3.59 ms and 3.81 ms) have shorter orbital periods (3.4 hr and 5.0 hr) and lighter companions (minimum mass 0.12 Msun and 0.07 Msun respectively). The pulsar in the closest system is the fifth member of an emerging class of millisecond pulsars displaying irregular radio eclipses and having a relatively massive companion. This system is a good candidate for optical identification of the companion star. The lack of known isolated pulsars in NGC 6266 is also discussed.

Using pulse timing observations we have obtained precise parameters, including positions with about 20 mas accuracy, of five millisecond pulsars in NGC 6752. Three of them, located relatively close to the cluster center, have line-of-sight accelerations larger than the maximum value predicted by the central mass density derived from optical observation, providing dynamical evidence for a central mass-to-light ratio >~ 10, much higher than for any other globular cluster. It is likely that the other two millisecond pulsars have been ejected out of the core to their present locations at 1.4 and 3.3 half-mass radii, respectively, suggesting unusual non-thermal dynamics in the cluster core.

We report the discovery with the Parkes radio telescope of a pulsar associated with the ~ 1700 yr-old oxygen-rich composite supernova remnant G292.0+1.8. PSR J1124-5916 has period 135 ms and period derivative 7.4e-13, implying characteristic age 2900 yr, spin-down luminosity 1.2e37 erg/s, and surface magnetic field strength 1.0e13 G. Association between the pulsar and the synchrotron nebula previously identified with Chandra within this supernova remnant is confirmed by the subsequent detection of X-ray pulsations by Hughes et al. The pulsar's flux density at 1400 MHz is very small, S ~ 0.08 mJy, but the radio luminosity of Sd^2 ~ 2 mJy kpc^2 is not especially so, although it is one order of magnitude smaller than that of the least luminous young pulsar previously known. This discovery suggests that very deep radio searches should be done for pulsations from pulsar wind nebulae in which the central pulsed source is yet to be detected and possibly other more exotic neutron stars.

We have discovered 12 new millisecond pulsars in 6 globular clusters in which no pulsars were previously known, in the first two years of a search at 1.4 GHz in progress at the Parkes radio telescope. Here we briefly describe the motivation, the new hardware and software systems adopted for this survey, and we present the results obtained thus far.

Binary pulsars provide an excellent system for testing general relativity because of their intrinsic rotational stability and the precision with which radio observations can be used to determine their orbital dynamics. Measurements of the rate of orbital decay of two pulsars have been shown to be consistent with the emission of gravitational waves as predicted by general relativity, providing the most convincing evidence for the self-consistency of the theory to date. However, independent verification of the orbital geometry in these systems was not possible. Such verification may be obtained by determining the orientation of a binary pulsar system using only classical geometric constraints, permitting an independent prediction of general relativistic effects. Here we report high-precision timing of the nearby binary millisecond pulsar PSR J0437-4715, which establish the three-dimensional structure of its orbit. We see the expected retardation of the pulse signal arising from the curvature of space-time in the vicinity of the companion object (the `Shapiro delay'), and we determine the mass of the pulsar and its white dwarf companion. Such mass determinations contribute to our understanding of the origin and evolution of neutron stars.

We present the results of one year of pulse timing observations of PSR J1740-5340, an eclipsing millisecond pulsar located in the globular cluster NGC 6397. We have obtained detailed orbital parameters and a precise position for the pulsar. The radio pulsar signal shows frequent interactions with the atmosphere of the companion, and suffers significant and strongly variable delays and intensity variations over a wide range of orbital phases. These characteristics and the binary parameters indicate that the companion may be a bloated main-sequence star or the remnant (still filling its Roche lobe) of the star that spun up the pulsar. In both cases, this would be the first binary millisecond pulsar system with such a companion.

A new search of globular clusters for millisecond pulsars is in progress at Parkes. In this paper we describe the motivation, the new hardware and software systems adopted, the survey plan and the preliminary results. So far, we have discovered ten new millisecond pulsars in four clusters for which no associated pulsars were previously known.