Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fast Folding Algorithm searches. We present 19 new discoveries, bringing the total to 27 discoveries in the programme, and we report the results of follow-up timing observations at Parkes for 26 of these pulsars, including the millisecond pulsar PSR J1421-4409; the faint, highly-modulated, slow pulsar PSR J1646-1910; and the nulling pulsar PSR J1337-4441. We present new timing solutions for 23 pulsars, and we report flux densities, modulation indices, and polarization properties.

The soft-gamma repeater Swift J1818.0$-$1607 is only the fifth magnetar found to exhibit pulsed radio emission. Using the Ultra-Wideband Low receiver system of the Parkes radio telescope, we conducted a 3 h observation of Swift J1818.0$-$1607. Folding the data at a rotation period of $P = 1.363$ s, we obtained wideband polarization profiles and flux density measurements covering radio frequencies between 704-4032 MHz. After measuring, and then correcting for the pulsar's rotation measure of $1442.0 \pm 0.2$ rad m$^{-2}$, we find the radio profile is between 80-100 per cent linearly polarized across the wide observing band, with a small amount of depolarization at low frequencies that we ascribe to scatter broadening. We also measure a steep spectral index of $\alpha = -2.26^{+0.02}_{-0.03}$ across our large frequency range, a significant deviation from the flat or inverted spectra often associated with radio-loud magnetars. The steep spectrum and temporal rise in flux density bears some resemblance to the behaviour of the magnetar-like, rotation-powered pulsar PSR J1119$-$6127. This leads us to speculate that Swift J1818.0$-$1607 may represent an additional link between rotation-powered pulsars and magnetars.

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

The current gravitational wave detectors have identified a surprising population of heavy stellar mass black holes, and an even larger population of coalescing neutron stars. The first observations have led to many dramatic discoveries and the confirmation of general relativity in very strong gravitational fields. The future of gravitational wave astronomy looks bright, especially if additional detectors with greater sensitivity, broader bandwidth, and better global coverage can be implemented. The first discoveries add impetus to gravitational wave detectors designed to detect in the nHz, mHz and kHz frequency bands. This paper reviews the century-long struggle that led to the recent discoveries, and reports on designs and possibilities for future detectors. The benefits of future detectors in the Asian region are discussed, including analysis of the benefits of a detector located in Australia.

In Paper I of this series, we detected a significant value of the braking index ($n$) for 19 young, high-$\dot{E}$ radio pulsars using $\sim$ 10 years of timing observations from the 64-m Parkes radio telescope. Here we investigate this result in more detail using a Bayesian pulsar timing framework to model timing noise and to perform selection to distinguish between models containing exponential glitch recovery and braking index signatures. We show that consistent values of $n$ are maintained with the addition of substantial archival data, even in the presence of glitches. We provide strong arguments that our measurements are unlikely due to exponential recovery signals from unseen glitches even though glitches play a key role in the evolution of a pulsar's spin frequency. We conclude that, at least over decadal time scales, the value of $n$ can be significantly larger than the canonical 3 and discuss the implications for the evolution of pulsars.

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.

While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($\nu$) and spin-down frequency ($\dot{\nu}$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $\nu^{a} |\dot{\nu}|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position.

We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 \mu$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution.

Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring precession of the orbit resulting from rapid rotation of the WD. Lense-Thirring precession, an effect of relativistic frame-dragging, is a prediction of general relativity. This detection is consistent with the evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period $< 200$ seconds.

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project.

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.

The past four years have seen a scientific revolution through the birth of a new field: gravitational-wave astronomy. The first detection of gravitational waves---recognised by the 2017 Nobel Prize in Physics---provided unprecedented tests of general relativity while unveiling a previously unknown class of massive black holes, thirty times more massive than the Sun. The subsequent detection of gravitational waves from a merging binary neutron star confirmed the hypothesised connection between binary neutron stars and short gamma-ray bursts while providing an independent measurement of the expansion of the Universe. The discovery enabled precision measurement of the speed of gravity while shedding light on the origin of heavy elements. At the time of writing, the Laser Interferometer Gravitational-wave Observatory (LIGO) and its European partner, Virgo, have published the detection of eleven gravitational-wave events. New, not-yet-published detections are announced on a nearly weekly basis. This fast-growing catalogue of gravitational-wave transients is expected to yield insights into a number of topics, from the equation of state of matter at supra-nuclear densities to the fate of massive stars. The science potential of 3G observatories is enormous, enabling measurements of gravitational waves from the edge of the Universe and precise determination of the neutron star equation of state. Australia is well-positioned to help develop the required technology. The Mid-term Review for the Decadal plan for Australian astronomy 2016-2025 should consider investment in a scoping study for an Australian Gravitational-Wave Pathfinder that develops and validates core technologies required for the global 3G detector network.

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.

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.

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

The smooth spin-down of young pulsars is perturbed by two non-deterministic phenomenon, glitches and timing noise. Although the timing noise provides insights into nuclear and plasma physics at extreme densities, it acts as a barrier to high-precision pulsar timing experiments. An improved methodology based on Bayesian inference is developed to simultaneously model the stochastic and deterministic parameters for a sample of 85 high-$\dot{E}$ radio pulsars observed for $\sim$ 10 years with the 64-m Parkes radio telescope. Timing noise is known to be a red process and we develop a parametrization based on the red-noise amplitude ($A_{\rm red}$) and spectral index ($\beta$). We measure the median $A_{\rm red}$ to be $-10.4^{+1.8}_{-1.7}$ yr$^{3/2}$ and $\beta$ to be $-5.2^{+3.0}_{-3.8}$ and show that the strength of timing noise scales proportionally to $\nu^{1}|\dot{\nu}|^{-0.6\pm0.1}$, where $\nu$ is the spin frequency of the pulsar and $\dot{\nu}$ its spin-down rate. Finally, we measure significant braking indices for 19 pulsars, proper motions for two pulsars and discuss the presence of periodic modulation in the arrival times of five pulsars.

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.

We study the polarization properties of 22 known rotating radio transients (RRATs) with the 64-m Parkes radio telescope and present the Faraday rotation measures (RMs) for the 17 with linearly polarized flux exceeding the off-pulse noise by 3$\sigma$. Each RM was estimated using a brute-force search over trial RMs that spanned the maximum measurable range $\pm1.18 \times 10^5 \, \mathrm{rad \, m^2}$ (in steps of 1 $\mathrm{rad \, m^2}$), followed by an iterative refinement algorithm. The measured RRAT RMs are in the range |RM| $\sim 1$ to $\sim 950$ rad m$^{-2}$ with an average linear polarization fraction of $\sim 40$ per cent. Individual single pulses are observed to be up to 100 per cent linearly polarized. The RMs of the RRATs and the corresponding inferred average magnetic fields (parallel to the line-of-sight and weighted by the free electron density) are observed to be consistent with the Galactic plane pulsar population. Faraday rotation analyses are typically performed on accumulated pulsar data, for which hundreds to thousands of pulses have been integrated, rather than on individual pulses. Therefore, we verified the iterative refinement algorithm by performing Monte Carlo simulations of artificial single pulses over a wide range of S/N and RM. At and above a S/N of 17 in linearly polarized flux, the iterative refinement recovers the simulated RM value 100 per cent of the time with a typical mean uncertainty of $\sim5$ rad m$^{-2}$. The method described and validated here has also been successfully used to determine reliable RMs of several fast radio bursts (FRBs) discovered at Parkes.

We describe the ongoing `Survey for Magnetars, Intermittent pulsars, Rotating radio transients and Fast radio bursts' (SMIRF), performed using the newly refurbished UTMOST telescope. SMIRF repeatedly sweeps the southern Galactic plane performing real-time periodicity and single-pulse searches, and is the first survey of its kind carried out with an interferometer. SMIRF is facilitated by a robotic scheduler which is capable of fully autonomous commensal operations. We report on the SMIRF observational parameters, the data analysis methods, the survey's sensitivities to pulsars, techniques to mitigate radio frequency interference and present some early survey results. UTMOST's wide field of view permits a full sweep of the Galactic plane to be performed every fortnight, two orders of magnitude faster than previous surveys. In the six months of operations from January to June 2018, we have performed $\sim 10$ sweeps of the Galactic plane with SMIRF. Notable blind re-detections include the magnetar PSR J1622$-$4950, the RRAT PSR J0941$-$3942 and the eclipsing pulsar PSR J1748$-$2446A. We also report the discovery of a new pulsar, PSR J1705$-$54. Our follow-up of this pulsar with the UTMOST and Parkes telescopes at an average flux limit of $\leq 20$ mJy and $\leq 0.16$ mJy respectively, categorizes this as an intermittent pulsar with a high nulling fraction of $< 0.002$

We detail a new fast radio burst (FRB) survey with the Molonglo Radio Telescope, in which six FRBs were detected between June 2017 and December 2018. By using a real-time FRB detection system, we captured raw voltages for five of the six events, which allowed for coherent dedispersion and very high time resolution (10.24 $\mu$s) studies of the bursts. Five of the FRBs show temporal broadening consistent with interstellar and/or intergalactic scattering, with scattering timescales ranging from 0.16 to 29.1 ms. One burst, FRB181017, shows remarkable temporal structure, with 3 peaks each separated by 1 ms. We searched for phase-coherence between the leading and trailing peaks and found none, ruling out lensing scenarios. Based on this survey, we calculate an all-sky rate at 843 MHz of $98^{+59}_{-39}$ events sky$^{-1}$ day$^{-1}$ to a fluence limit of 8 Jy-ms: a factor of 7 below the rates estimated from the Parkes and ASKAP telescopes at 1.4 GHz assuming the ASKAP-derived spectral index $\alpha=-1.6$ ($F_{\nu}\propto\nu^{\alpha}$). Our results suggest that FRB spectra may turn over below 1 GHz. Optical, radio and X-ray followup has been made for most of the reported bursts, with no associated transients found. No repeat bursts were found in the survey.

Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these observed $\dot P$ measurements to acceleration of these pulsars along their lines-of-sight (LOSs) due to the Galactic gravitational field. Using models of the Galactic mass distribution and pulsar velocities, we constrain the distances to these pulsars, placing them on the far-side of the Galaxy, providing the first accurate distance measurements to pulsars located in this region and allowing us to consider the electron density along these LOSs. We find the new electron density model YMW16 to be more consistent with these observations than the previous model NE2001. The LOS dynamics further constrain the model-dependent intrinsic $\dot P$ values for these pulsars and they are consistent with measurements for other known pulsars. In the future, the independent distance measurements to these and other pulsars near the GC would allow us to constrain the Galactic gravitational potential more accurately.

Future GW detector networks and EM observatories will provide a unique opportunity to observe the most luminous events in the Universe involving matter in extreme environs. They will address some of the key questions in physics and astronomy: formation and evolution of compact binaries, sites of formation of heavy elements and the physics of jets.

We present the earliest X-ray observations of the 2018 outburst of XTE J1810-197, the first outburst since its 2003 discovery as the prototypical transient and radio-emitting anomalous X-ray pulsar (AXP). The Monitor of All-sky X-ray Image (MAXI) detected XTE J1810-197 immediately after a November 20-26 visibility gap, contemporaneous with its reactivation as a radio pulsar, first observed on December 8. On December 13 the Nuclear Spectroscopic Telescope Array (NUSTAR) detected X-ray emission up to at least 30 keV, with a spectrum well-characterized by a blackbody plus power-law model with temperature kT = 0.74+/-0.02 keV and photon index Gamma = 4.4+/-0.2 or by a two-blackbody model with kT = 0.59+/-0.04 keV and kT = 1.0+/-0.1 keV, both including an additional power-law component to account for emission above 10 keV, with Gamma_h = -0.2+/-1.5 and Gamma_h = 1.5+/-0.5, respectively. The latter index is consistent with hard X-ray flux reported for the non-transient magnetars. In the 2-10 keV bandpass, the absorbed flux is 2E-10 erg/s/cm^2, a factor of 2 greater than the maximum flux extrapolated for the 2003 outburst. The peak of the sinusoidal X-ray pulse lags the radio pulse by approx. 0.13 cycles, consistent with their phase relationship during the 2003 outburst. This suggests a stable geometry in which radio emission originates on magnetic field lines containing currents that heat a spot on the neutron star surface. However, a measured energy-dependent phase shift of the pulsed X-rays suggests that all X-ray emitting regions are not precisely co-aligned.

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.

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 $\Omega$; 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 $\Omega=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 $\mu_{\alpha}\cos{\delta}=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $\mu_{\delta}=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$.

PSR J1141$-$6545 is a precessing binary pulsar that has the rare potential to reveal the two-dimensional structure of a non-recycled pulsar emission cone. It has undergone $\sim 25 \deg$ of relativistic spin precession in the $\sim18$ years since its discovery. In this paper, we present a detailed Bayesian analysis of the precessional evolution of the width of the total intensity profile, to understand the changes to the line-of-sight impact angle ($\beta$) of the pulsar using four different physically motivated prior distribution models. Although we cannot statistically differentiate between the models with confidence, the temporal evolution of the linear and circular polarisations strongly argue that our line-of-sight crossed the magnetic pole around MJD 54000 and that only two models remain viable. For both these models, it appears likely that the pulsar will precess out of our line-of-sight in the next $3-5$ years, assuming a simple beam geometry. Marginalising over $\beta$ suggests that the pulsar is a near-orthogonal rotator and provides the first polarization-independent estimate of the scale factor ($\mathbb{A}$) that relates the pulsar beam opening angle ($\rho$) to its rotational period ($P$) as $\rho = \mathbb{A}P^{-0.5}$ : we find it to be $> 6 \rm~deg~s^{0.5}$ at 1.4 GHz with 99\% confidence. If all pulsars emit from opposite poles of a dipolar magnetic field with comparable brightness, we might expect to see evidence of an interpulse arising in PSR J1141$-$6545, unless the emission is patchy.

We report on the latest six pulsars discovered through our standard pipeline in the intermediate-latitude region (|b| < 15 deg) of the Parkes High Time Resolution Universe Survey (HTRU). We also present timing solutions for the new discoveries and for 19 further pulsars for which only discovery parameters were previously published. Highlights of the presented sample include the isolated millisecond pulsar J1826-2415, the long-period binary pulsar J1837-0822 in a mildly eccentric 98-day orbit with a > 0.27 M_sun companion, and the nulling pulsar J1638-4233, detected only 10% of the time. Other interesting objects are PSR J1757-1500, exhibiting sporadic mode changes, and PSR J1635-2616 showing one glitch over 6 years. The new discoveries bring the total count of HTRU intermediate-latitude pulsars to 113, 25% of which are recycled pulsars. This is the higest ratio of recycled over ordinary pulsars discoveries of all recent pulsar surveys in this region of the sky. Among HTRU recycled pulsars, four are isolated objects. Comparing the characteristics of Galactic fully-recycled isolated MSPs with those of eclipsing binaries ('spiders'), from which the former are believed to have formed, we highlight a discrepancy in their spatial distribution. This may reflect a difference in the natal kick, hence, possibly, a different formation path. On the other hand, however, isolated fully-recycled MSPs spin periods are, on average, longer than those of spiders, in line with what one would expect, from simple magnetic-dipole spin-down, if the former were indeed evolved from the latter.

The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type.

We present an overview and the first results from a large-scale pulsar timing programme that is part of the UTMOST project at the refurbished Molonglo Observatory Synthesis Radio Telescope (MOST) near Canberra, Australia. We currently observe more than 400 mainly bright southern radio pulsars with up to daily cadences. For 205 (8 in binaries, 4 millisecond pulsars) we publish updated timing models, together with their flux densities, flux density variability, and pulse widths at 843 MHz, derived from observations spanning between 1.4 and 3 yr. In comparison with the ATNF pulsar catalogue, we improve the precision of the rotational and astrometric parameters for 123 pulsars, for 47 by at least an order of magnitude. The time spans between our measurements and those in the literature are up to 48 yr, which allows us to investigate their long-term spin-down history and to estimate proper motions for 60 pulsars, of which 24 are newly determined and most are major improvements. The results are consistent with interferometric measurements from the literature. A model with two Gaussian components centred at 139 and $463~\text{km} \: \text{s}^{-1}$ fits the transverse velocity distribution best. The pulse duty cycle distributions at 50 and 10 per cent maximum are best described by log-normal distributions with medians of 2.3 and 4.4 per cent, respectively. We discuss two pulsars that exhibit spin-down rate changes and drifting subpulses. Finally, we describe the autonomous observing system and the dynamic scheduler that has increased the observing efficiency by a factor of 2-3 in comparison with static scheduling.

We have performed a new search for radio pulsars in archival data of the intermediate and high Galactic latitude parts of the Southern High Time Resolution Universe pulsar survey. This is the first time the entire dataset has been searched for binary pulsars, an achievement enabled by GPU-accelerated dedispersion and periodicity search codes nearly 50 times faster than the previously used pipeline. Candidate selection was handled entirely by a Machine Learning algorithm, allowing for the assessment of 17.6 million candidates in a few person-days. We have also introduced an outlier detection algorithm for efficient radio-frequency interference (RFI) mitigation on folded data, a new approach that enabled the discovery of pulsars previously masked by RFI. We discuss implications for future searches, particularly the importance of expanding work on RFI mitigation to improve survey completeness. In total we discovered 23 previously unknown sources, including 6 millisecond pulsars and at least 4 pulsars in binary systems. We also found an elusive but credible redback candidate that we have yet to confirm.

Fast radio bursts (FRBs) are millisecond pulses of radio emission of seemingly extragalactic origin. More than 50 FRBs have now been detected, with only one seen to repeat. Here we present a new FRB discovery, FRB 110214, which was detected in the high latitude portion of the High Time Resolution Universe South survey at the Parkes telescope. FRB 110214 has one of the lowest dispersion measures of any known FRB (DM = 168.9$\pm$0.5 pc cm$^{-3}$), and was detected in two beams of the Parkes multi-beam receiver. A triangulation of the burst origin on the sky identified three possible regions in the beam pattern where it may have originated, all in sidelobes of the primary detection beam. Depending on the true location of the burst the intrinsic fluence is estimated to fall in the range of 50 -- 2000 Jy ms, making FRB 110214 one of the highest-fluence FRBs detected with the Parkes telescope. No repeating pulses were seen in almost 100 hours of follow-up observations with the Parkes telescope down to a limiting fluence of 0.3 Jy ms for a 2-ms pulse. Similar low-DM, ultra-bright FRBs may be detected in telescope sidelobes in the future, making careful modeling of multi-beam instrument beam patterns of utmost importance for upcoming FRB surveys.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

We report the detection of a glitch event in the pulsar J1709$-$4429 (also known as B1706$-$44) during regular monitoring observations with the Molonglo Observatory Synthesis Telescope (UTMOST). The glitch was found during timing operations, in which we regularly observe over 400 pulsars with up to daily cadence, while commensally searching for Rotating Radio Transients, pulsars, and FRBs. With a fractional size of $\Delta\nu/\nu \approx 52.4 \times10^{-9}$, the glitch reported here is by far the smallest known for this pulsar, attesting to the efficacy of glitch searches with high cadence using UTMOST.

We report a new Fast Radio Burst (FRB) discovered in real-time as part of the UTMOST project at the Molonglo Observatory Synthesis Radio Telescope (MOST). FRB170827 is the first detected with our low-latency ($< 24$ s), machine-learning-based FRB detection system. The FRB discovery was accompanied by the capture of voltage data at the native time and frequency resolution of the observing system, enabling coherent dedispersion and detailed off-line analysis, which have unveiled fine temporal and frequency structure. The dispersion measure (DM) of 176.80 $\pm$ 0.04 pc cm$^{-3}$, is the lowest of the FRB population. The Milky Way contribution along the line of sight is $\sim$ 40 pc cm$^{-3}$, leaving an excess DM of $\sim$ 145 pc cm$^{-3}$. The FRB has a fluence $>$ 20 $\pm$ 7 Jy ms, and is narrow, with a width of $\sim$ 400 $\mu$s at 10$\%$ of its maximum amplitude. However, the burst shows three temporal components, the narrowest of which is $\sim$ 30 $\mu$s, and a scattering timescale of $4.1 \pm 2.7$ $\mu$s. The FRB shows spectral modulations on frequency scales of 1.5 MHz and 0.1 MHz. Both are prominent in the dynamic spectrum, which shows a very bright region of emission between 841 and 843 MHz, and weaker, patchy emission across the entire band. We show the fine spectral structure could arise in the FRB host galaxy, or its immediate vicinity.

We report on the polarization properties of two fast radio bursts (FRBs): 151230 and 160102 discovered in the SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes radio telescope. FRB 151230 is observed to be 6 $\pm$ 11% circularly polarized and 35 $\pm$ 13 % linearly polarized with a rotation measure (RM) consistent with zero. Conversely, FRB 160102 is observed to have a circular polarization fraction of 30 $\pm$ 11 %, linear polarization fraction of 84 $\pm$ 15 % for RM = -221(6) rad m$^{-2}$ and the highest measured DM (2596.1 $\pm$ 0.3 pc cm$^{-3}$) for an FRB to date. We examine possible progenitor models for FRB 160102 in extragalactic, non-cosmological and cosmological scenarios. After accounting for the Galactic foreground contribution, we estimate the intrinsic RM to be -256(9) rad m$^{-2}$ in the low-redshift case and $\sim 2.4 \times 10^{2}$ rad m$^{-2}$ in the high-redshift case. We assess the relative likeliness of these scenarios and how each can be tested. We also place constraints on the scattering measure and study the impact of scattering on the signal's polarization position angle.

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.

Radio pulsars have been responsible for many astonishing astrophysical and fundamental physics breakthroughs since their discovery 50 years ago. In this review I will discuss many of the highlights, most of which were only possible because of the provision of large-scale observing facilities. The next 50 years of pulsar astronomy can be very bright, but only if our governments properly plan and fund the infrastructure necessary to enable future discoveries. Being a small sub-field of astronomy places an onus on the pulsar community to have an open-source/open access approach to data, software, and major observing facilities to enable new groups to emerge to keep the field vibrant.

The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results.

Pulsars in relativistic binary systems have emerged as fantastic natural laboratories for testing theories of gravity, the most prominent example being the double pulsar, PSR J0737$-$3039. The HTRU-South Low Latitude pulsar survey represents one of the most sensitive blind pulsar surveys taken of the southern Galactic plane to date, and its primary aim has been the discovery of new relativistic binary pulsars. Here we present our binary pulsar searching strategy and report on the survey's flagship discovery, PSR J1757$-$1854. A 21.5-ms pulsar in a relativistic binary with an orbital period of 4.4 hours and an eccentricity of 0.61, this double neutron star (DNS) system is the most accelerated pulsar binary known, and probes a relativistic parameter space not yet explored by previous pulsar binaries.

We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757$-$1854. With simulations suggesting 3-$\sigma$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories.

We present the discovery of a binary millisecond pulsar (MSP), PSR J2322$-$2650, found in the Southern section of the High Time Resolution Universe survey. This system contains a 3.5-ms pulsar with a $\sim10^{-3}$ M$_{\odot}$ companion in a 7.75-hour circular orbit. Follow-up observations at the Parkes and Lovell telescopes have led to precise measurements of the astrometric and spin parameters, including the period derivative, timing parallax, and proper motion. PSR J2322$-$2650 has a parallax of $4.4\pm1.2$ mas, and is thus at an inferred distance of $230^{+90}_{-50}$ pc, making this system a candidate for optical studies. We have detected a source of $R\approx26.4$ mag at the radio position in a single $R$-band observation with the Keck Telescope, and this is consistent with the blackbody temperature we would expect from the companion if it fills its Roche lobe. The intrinsic period derivative of PSR J2322$-$2650 is among the lowest known, $4.4(4)\times10^{-22}$ s s$^{-1}$, implying a low surface magnetic field strength, $4.0(4)\times10^7$ G. Its mean radio flux density of 160 $\mu$Jy combined with the distance implies that its radio luminosity is the lowest ever measured, $0.008(5)$ mJy kpc$^2$. The inferred population of these systems in the Galaxy may be very significant, suggesting that this is a common MSP evolutionary path.

The initial results from timing observations of PSR J1141-6545, a relativistic pulsar white-dwarf binary system, are presented. Predictions from the timing baseline hint at the most stringent test of gravity by an asymmetric binary yet. The timing precision has been hindered by the dramatic variations of the pulse profile due to geodetic precession, a pulsar glitch and red timing noise. Methods to overcome such timing irregularities are briefly presented along with preliminary results from the test of the General Theory of Relativity (GR) from this pulsar

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We present a study of the spectral properties of 441 pulsars observed with the Parkes radio telescope near the centre frequencies of 728, 1382 and 3100 MHz. The observations at 728 and 3100 MHz were conducted simultaneously using the dual-band 10-50cm receiver. These high-sensitivity, multi-frequency observations provide a systematic and uniform sample of pulsar flux densities. We combine our measurements with spectral data from the literature in order to derive the spectral properties of these pulsars. Using techniques from robust regression and information theory we classify the observed spectra in an objective, robust and unbiased way into five morphological classes: simple or broken power law, power law with either low or high-frequency cut-off and log-parabolic spectrum. While about $79 \%$ of the pulsars that could be classified have simple power law spectra, we find significant deviations in 73 pulsars, 35 of which have curved spectra, 25 with a spectral break and 10 with a low-frequency turn-over. We identify 11 gigahertz-peaked spectrum (GPS) pulsars, with 3 newly identified in this work and 8 confirmations of known GPS pulsars; 3 others show tentative evidence of GPS, but require further low-frequency measurements to support this classification. The weighted mean spectral index of all pulsars with simple power law spectra is $-1.60 \pm 0.03$. The observed spectral indices are well described by a shifted log-normal distribution. The strongest correlations of spectral index are with spin-down luminosity, magnetic field at the light-cylinder and spin-down rate. We also investigate the physical origin of the observed spectral features and determine emission altitudes for three pulsars.

The Molonglo Observatory Synthesis Telescope (MOST) is an 18,000 square meter radio telescope situated some 40 km from the city of Canberra, Australia. Its operating band (820-850 MHz) is now partly allocated to mobile phone communications, making radio astronomy challenging. We describe how the deployment of new digital receivers (RX boxes), Field Programmable Gate Array (FPGA) based filterbanks and server-class computers equipped with 43 GPUs (Graphics Processing Units) has transformed MOST into a versatile new instrument (the UTMOST) for studying the dynamic radio sky on millisecond timescales, ideal for work on pulsars and Fast Radio Bursts (FRBs). The filterbanks, servers and their high-speed, low-latency network form part of a hybrid solution to the observatory's signal processing requirements. The emphasis on software and commodity off-the-shelf hardware has enabled rapid deployment through the re-use of proven 'software backends' for its signal processing. The new receivers have ten times the bandwidth of the original MOST and double the sampling of the line feed, which doubles the field of view. The UTMOST can simultaneously excise interference, make maps, coherently dedisperse pulsars, and perform real-time searches of coherent fan beams for dispersed single pulses. Although system performance is still sub-optimal, a pulsar timing and FRB search programme has commenced and the first UTMOST maps have been made. The telescope operates as a robotic facility, deciding how to efficiently target pulsars and how long to stay on source, via feedback from real-time pulsar folding. The regular timing of over 300 pulsars has resulted in the discovery of 7 pulsar glitches and 3 FRBs. The UTMOST demonstrates that if sufficient signal processing can be applied to the voltage streams it is possible to perform innovative radio science in hostile radio frequency environments.

We describe the Survey for Pulsars and Extragalactic Radio Bursts (SUPERB), an ongoing pulsar and fast transient survey using the Parkes radio telescope. SUPERB involves real-time acceleration searches for pulsars and single-pulse searches for pulsars and fast radio bursts. We report on the observational setup, data analysis, multi-wavelength/messenger connections, survey sensitivities to pulsars and fast radio bursts and the impact of radio frequency interference. We further report on the first 10 pulsars discovered in the project. Among these is PSR~J1306$-$40, a millisecond pulsar in a binary system where it appears to be eclipsed for a large fraction of the orbit. PSR~J1421$-$4407 is another binary millisecond pulsar; its orbital period is $30.7$ days. This orbital period is in a range where only highly eccentric binaries are known, and expected by theory; despite this its orbit has an eccentricity of $10^{-5}$.

We report the detection of an ultra-bright fast radio burst (FRB) from a modest, 3.4-day pilot survey with the Australian Square Kilometre Array Pathfinder. The survey was conducted in a wide-field fly's-eye configuration using the phased-array-feed technology deployed on the array to instantaneously observe an effective area of $160$ deg$^2$, and achieve an exposure totaling $13200$ deg$^2$ hr. We constrain the position of FRB 170107 to a region $8'\times8'$ in size (90% containment) and its fluence to be $58\pm6$ Jy ms. The spectrum of the burst shows a sharp cutoff above $1400$ MHz, which could be either due to scintillation or an intrinsic feature of the burst. This confirms the existence of an ultra-bright ($>20$ Jy ms) population of FRBs.

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.

HIPSR (HI-Pulsar) is a digital signal processing system for the Parkes 21-cm Multibeam Receiver that provides larger instantaneous bandwidth, increased dynamic range, and more signal processing power than the previous systems in use at Parkes. The additional computational capacity enables finer spectral resolution in wideband HI observations and real-time detection of Fast Radio Bursts during pulsar surveys. HIPSR uses a heterogeneous architecture, consisting of FPGA-based signal processing boards connected via high-speed Ethernet to high performance compute nodes. Low-level signal processing is conducted on the FPGA-based boards, and more complex signal processing routines are conducted on the GPU-based compute nodes. The development of HIPSR was driven by two main science goals: to provide large bandwidth, high-resolution spectra suitable for 21-cm stacking and intensity mapping experiments; and to upgrade the Berkeley-Parkes-Swinburne Recorder (BPSR), the signal processing system used for the High Time Resolution Universe (HTRU) Survey and the Survey for Pulsars and Extragalactic Radio Bursts (SUPERB).

We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and latitude of 24.66^{\circ}, 5.28^{\circ}, 25 degrees away from the Galactic Center. The burst was found to be 43$\pm$5% linearly polarized with a rotation measure (RM) in the range -9 < RM < 12 rad m^{-2} (95% confidence level), consistent with zero. The burst was followed-up with 11 telescopes to search for radio, optical, X-ray, gamma-ray and neutrino emission. Neither transient nor variable emission was found to be associated with the burst and no repeat pulses have been observed in 17.25 hours of observing. The sightline to the burst is close to the Galactic plane and the observed physical properties of FRB 150215 demonstrate the existence of sight lines of anomalously low RM for a given electron column density. The Galactic RM foreground may approach a null value due to magnetic field reversals along the line of sight, a decreased total electron column density from the Milky Way, or some combination of these effects. A lower Galactic DM contribution might explain why this burst was detectable whereas previous searches at low latitude have had lower detection rates than those out of the plane.

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$\mu$as level between ephemerides.

We present the first interferometric detections of Fast Radio Bursts (FRBs), an enigmatic new class of astrophysical transient. In a 180-day survey of the Southern sky we discovered 3 FRBs at 843 MHz with the UTMOST array, as part of commissioning science during a major ongoing upgrade. The wide field of view of UTMOST ($\approx 9$ deg$^{2}$) is well suited to FRB searches. The primary beam is covered by 352 partially overlapping fan-beams, each of which is searched for FRBs in real time with pulse widths in the range 0.655 to 42 ms, and dispersion measures $\leq$2000 pc cm$^{-3}$. Detections of FRBs with the UTMOST array places a lower limit on their distances of $\approx 10^4$ km (limit of the telescope near-field) supporting the case for an astronomical origin. Repeating FRBs at UTMOST or an FRB detected simultaneously with the Parkes radio telescope and UTMOST, would allow a few arcsec localisation, thereby providing an excellent means of identifying FRB host galaxies, if present. Up to 100 hours of follow-up for each FRB has been carried out with the UTMOST, with no repeating bursts seen. From the detected position, we present 3$\sigma$ error ellipses of 15 arcsec x 8.4 deg on the sky for the point of origin for the FRBs. We estimate an all-sky FRB rate at 843 MHz above a fluence $\cal F_\mathrm{lim}$ of 11 Jy ms of $\sim 78$ events sky$^{-1}$ d$^{-1}$ at the 95 percent confidence level. The measured rate of FRBs at 843 MHz is of order two times higher than we had expected, scaling from the FRB rate at the Parkes radio telescope, assuming that FRBs have a flat spectral index and a uniform distribution in Euclidean space. We examine how this can be explained by FRBs having a steeper spectral index and/or a flatter log$N$-log$\mathcal{F}$ distribution than expected for a Euclidean Universe.

Advances in astronomy are intimately linked to advances in digital signal processing (DSP). This special issue is focused upon advances in DSP within radio astronomy. The trend within that community is to use off-the-shelf digital hardware where possible and leverage advances in high performance computing. In particular, graphics processing units (GPUs) and field programmable gate arrays (FPGAs) are being used in place of application-specific circuits (ASICs); high-speed Ethernet and Infiniband are being used for interconnect in place of custom backplanes. Further, to lower hurdles in digital engineering, communities have designed and released general-purpose FPGA-based DSP systems, such as the CASPER ROACH board, ASTRON Uniboard and CSIRO Redback board. In this introductory article, we give a brief historical overview, a summary of recent trends, and provide an outlook on future directions.

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.

Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.

We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $\Gamma= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase-resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsar's phase-folded lightcurve with the pulsar's well-defined mass and distance from radio timing observations.

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.

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

High time resolution radio surveys over the last few years have discovered a population of millisecond-duration transient bursts called Fast Radio Bursts (FRBs), which remain of unknown origin. FRBs exhibit dispersion consistent with propagation through a cold plasma and dispersion measures indicative of an origin at cosmological distances. In this paper we perform Monte Carlo simulations of a cosmological population of FRBs, based on assumptions consistent with observations of their energy distribution, their spatial density as a function of redshift and the properties of the interstellar and intergalactic media. We examine whether the dispersion measures, fluences, inferred redshifts, signal-to-noises and effective widths of known FRBs are consistent with a cosmological population. Statistical analyses indicate that at least 50 events at Parkes are required to distinguish between a constant co-moving FRB density, and a FRB density that evolves with redshift like the cosmological star formation rate density.

We report the first radio interferometric search at 843 MHz for fast transients, particularly Fast Radio Bursts (FRBs). The recently recommissioned Swinburne University of Technology's digital backend for the Molonglo Observatory Synthesis Telescope array (the UTMOST) with its large collecting area (18,000 $\mathrm{m^2}$) and wide instantaneous field of view (7.80 $\mathrm{deg^2}$) is expected to be an efficient tool to detect FRBs. As an interferometer it will be capable of discerning whether the FRBs are truly a celestial population. We show that UTMOST at full design sensitivity can detect an event approximately every few days. We report on 2 preliminary FRB surveys at about 7% and 14% respectively of the array's final sensitivity. Several pulsars have been detected via single pulses and no FRBs were discovered with pulse widths ($W$), in the range 655.36 $\mu$s $< W < 41.9$ ms and dispersion measures (DMs) in the range $100 < $DM$< 2000$ $\mathrm{pc\,cm^{-3}}$. This non-detection sets a 2$\sigma$ upper limit of the sky rate of not more than 1000 events $\mathrm{sky^{-1}}$ $\mathrm{day^{-1}}$ at 843 MHz down to a flux limit of 11 Jy for 1 ms FRBs. We show that this limit is consistent with previous survey limits at 1.4 GHz and 145 MHz and set a lower limit on the mean spectral index of FRBs of $\alpha > -3.2$.

Here we present a catalogue of known Fast Radio Burst (FRB) sources in the form of an online catalogue, FRBCAT. The catalogue includes information about the instrumentation used for the observations for each detected burst, the measured quantities from each observation, and model-dependent quantities derived from observed quantities. To aid in consistent comparisons of burst properties such as width and signal-to-noise ratios we have reprocessed all the bursts for which we have access to the raw data, with software which we make available. The originally derived properties are also listed for comparison. The catalogue is hosted online as a MySQL database which can also be downloaded in tabular or plain text format for off-line use. This database will be maintained for use by the community for studies of the FRB population as it grows.

The detection of five new fast radio bursts (FRBs) found in the High Time Resolution Universe high latitude survey is presented. The rate implied is 6$^{+4}_{-3}\times~10^3$ (95%) FRBs sky$^{-1}$ day$^{-1}$ above a fluence of between 0.13 and 5.9 Jy ms for FRBs between 0.128 and 262 ms in duration. One of these FRBs has a clear two-component profile, each component is similar to the known population of single component FRBs and are separated by 2.4(4) ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intra-channel smearing. The two-component FRB also has the highest dispersion measure (1629 pc cm$^{-3}$) of any FRB to-date. Many of the proposed models to explain FRBs use a single high energy event involving compact objects (such as neutron star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing or orbital events however could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

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.

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.

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.

Several theories exist to explain the source of the bright, millisecond duration pulses known as fast radio bursts (FRBs). If the progenitors of FRBs are non-cataclysmic, such as giant pulses from pulsars, pulsar-planet binaries, or magnetar flares, FRB emission may be seen to repeat. We have undertaken a survey of the fields of eight known FRBs from the High Time Resolution Universe survey to search for repeating pulses. Although no repeat pulses were detected the survey yielded the detection of a new FRB, described in Petroff et al. (2015a). From our observations we rule out periodic repeating sources with periods P $\leq$ 8.6 hours and rule out sources with periods 8.6 < P < 21 hours at the 90% confidence level. At P $\geq$ 21 hours our limits fall off as ~1/P. Dedicated and persistent observations of FRB source fields are needed to rule out repetition on longer timescales, a task well-suited to next generation wide-field transient detectors.

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.

We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.

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.

Globular clusters are highly efficient radio pulsar factories. These pulsars can be used as precision probes of the clusters' structure, gas content, magnetic field, and formation history; some of them are also highly interesting in their own right because they probe exotic stellar evolution scenarios as well as the physics of dense matter, accretion, and gravity. Deep searches with SKA1-MID and SKA1-LOW will plausibly double to triple the known population. Such searches will only require one to a few tied-array beams, and can be done during early commissioning of the telescope - before an all-sky pulsar survey using hundreds to thousands of tied-array beams is feasible. With SKA2 it will be possible to observe most of the active radio pulsars within a large fraction of the Galactic globular clusters, an estimated population of 600 - 3700 observable pulsars (those beamed towards us). This rivals the total population of millisecond pulsars that can be found in the Galactic field; fully characterizing it will provide the best-possible physical laboratories as well as a rich dynamical history of the Galactic globular cluster system.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.

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.

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}$.

We describe SPINN (Straightforward Pulsar Identification using Neural Networks), a high-performance machine learning solution developed to process increasingly large data outputs from pulsar surveys. SPINN has been cross-validated on candidates from the southern High Time Resolution Universe (HTRU) survey and shown to identify every known pulsar found in the survey data while maintaining a false positive rate of 0.64%. Furthermore, it ranks 99% of pulsars among the top 0.11% of candidates, and 95% among the top 0.01%. In conjunction with the PEASOUP pipeline (Barr et al., in prep.), it has already discovered four new pulsars in a re-processing of the intermediate Galactic latitude area of HTRU, three of which have spin periods shorter than 5 milliseconds. SPINN's ability to reduce the amount of candidates to visually inspect by up to four orders of magnitude makes it a very promising tool for future large-scale pulsar surveys. In an effort to provide a common testing ground for pulsar candidate selection tools and stimulate interest in their development, we also make publicly available the set of candidates on which SPINN was cross-validated.

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.

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$\sigma$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

Single pulses from radio pulsars contain a wealth of information about emission and propagation in the magnetosphere and insight into their timing properties. It was recently demonstrated that single-pulse emission is responsible for limiting the timing stability of the brightest of millisecond pulsars. We report on an analysis of more than a million single-pulses from PSR J0437-4715 and present various statistical properties such as the signal-to-noise ratio (S/N) distribution, timing and polarimetry of average profiles integrated from subpulses with chosen S/N cut-offs, modulation properties of the emission, phase-resolved statistics of the S/N, and two dimensional spherical histograms of the polarization vector orientation. The last of these indicates the presence of orthogonally polarised modes (OPMs). Combined with the dependence of the polarisation fraction on the S/N and polarimetry of the brightest pulses, the existence of OPMs constrains pulsar emission mechanisms and models for the plasma physics in the magnetosphere.

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

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 present a polarimetric analysis of 49 long-period pulsars discovered as part of the High Time Resolution Universe (HTRU) southern survey. The sources exhibit the typical characteristics of "old" pulsars, with low fractional linear and circular polarisation and narrow, multicomponent profiles. Although the position angle swings are generally complex, for two of the analysed pulsars (J1622-3751 and J1710-2616) we obtained an indication of the geometry via the rotating vector model. We were able to determine a value of the rotation measure (RM) for 34 of the sources which, when combined with their dispersion measures (DM), yields an integrated magnetic field strength along the line of sight. With the data presented here, the total number of values of RM associated to pulsars discovered during the HTRU southern survey sums to 51. The RMs are not consistent with the hypothesis of a counter-clockwise direction of the Galactic magnetic field within an annulus included between 4 and 6 kpc from the Galactic centre. A partial agreement with a counter-clockwise sense of the Galactic magnetic field within the spiral arms is however found in the area of the Carina-Sagittarius arm.

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

Searches for transient astrophysical sources often reveal unexpected classes of objects that are useful physical laboratories. In a recent survey for pulsars and fast transients we have uncovered four millisecond-duration radio transients all more than 40{\deg} from the Galactic plane. The bursts' properties indicate that they are of celestial rather than terrestrial origin. Host galaxy and intergalactic medium models suggest that they have cosmological redshifts of 0.5 to 1, and distances of up to 3 gigaparsecs. No temporally coincident x- or gamma-ray signature was identified in association with the bursts. Characterization of the source population and identification of host galaxies offers an opportunity to determine the baryonic content of the Universe.

We have used millisecond pulsars (MSPs) from the southern High Time Resolution Universe (HTRU) intermediate latitude survey area to simulate the distribution and total population of MSPs in the Galaxy. Our model makes use of the scale factor method, which estimates the ratio of the total number of MSPs in the Galaxy to the known sample. Using our best fit value for the z-height, z=500 pc, we find an underlying population of MSPs of 8.3(\pm 4.2)*10^4 sources down to a limiting luminosity of L_min=0.1 mJy kpc^2 and a luminosity distribution with a steep slope of d\log N/d\log L = -1.45(\pm 0.14). However, at the low end of the luminosity distribution, the uncertainties introduced by small number statistics are large. By omitting very low luminosity pulsars, we find a Galactic population above L_min=0.2 mJy kpc^2 of only 3.0(\pm 0.7)*10^4 MSPs. We have also simulated pulsars with periods shorter than any known MSP, and estimate the maximum number of sub-MSPs in the Galaxy to be 7.8(\pm 5.0)*10^4 pulsars at L=0.1 mJy kpc^2. In addition, we estimate that the high and low latitude parts of the southern HTRU survey will detect 68 and 42 MSPs respectively, including 78 new discoveries. Pulsar luminosity, and hence flux density, is an important input parameter in the model. Some of the published flux densities for the pulsars in our sample do not agree with the observed flux densities from our data set, and we have instead calculated average luminosities from archival data from the Parkes Telescope. We found many luminosities to be very different than their catalogue values, leading to very different population estimates. Large variations in flux density highlight the importance of including scintillation effects in MSP population studies.

Next-generation radio arrays, including the Square Kilometre Array (SKA) and its pathfinders, will open up new avenues for exciting transient science at radio wavelengths. Their innovative designs, comprising a large number of small elements, pose several challenges in digital processing and optimal observing strategies. The Giant Metrewave Radio Telescope (GMRT) presents an excellent test-bed for developing and validating suitable observing modes and strategies for transient experiments with future arrays. Here we describe the first phase of the ongoing development of a transient detection system for GMRT that is planned to eventually function in a commensal mode with other observing programs. It capitalizes on the GMRT's interferometric and sub-array capabilities, and the versatility of a new software backend. We outline considerations in the plan and design of transient exploration programs with interferometric arrays, and describe a pilot survey that was undertaken to aid in the development of algorithms and associated analysis software. This survey was conducted at 325 and 610 MHz, and covered 360 deg$^2$ of the sky with short dwell times. It provides large volumes of real data that can be used to test the efficacies of various algorithms and observing strategies applicable for transient detection. We present examples that illustrate the methodologies of detecting short-duration transients, including the use of sub-arrays for higher resilience to spurious events of terrestrial origin, localization of candidate events via imaging and the use of a phased array for improved signal detection and confirmation. In addition to demonstrating applications of interferometric arrays for fast transient exploration, our efforts mark important steps in the roadmap toward SKA-era science.

At the highest levels of pulsar timing precision achieved to date, experiments are limited by noise intrinsic to the pulsar. This stochastic wideband impulse modulated self-noise (SWIMS) limits pulsar timing precision by randomly biasing the measured times of arrival and thus increasing the root mean square (rms) timing residual. We discuss an improved methodology of removing this bias in the measured times of arrival by including information about polarized radiation. Observations of J0437-4715 made over a one-week interval at the Parkes Observatory are used to demonstrate a nearly 40 per cent improvement in the rms timing residual with this extended analysis. In this way, based on the observations over a 64 MHz bandwidth centred at 1341 MHz with integrations over 16.78 s we achieve a 476 ns rms timing residual. In the absence of systematic error, these results lead to a predicted rms timing residual of 30 ns in one hour integrations; however the data are currently limited by variable Faraday rotation in the Earth's ionosphere. The improvement demonstrated in this work provides an opportunity to increase the sensitivity in various pulsar timing experiments, for example pulsar timing arrays that pursue the detection of the stochastic background of gravitational waves. The fractional improvement is highly dependent on the properties of the pulse profile and the stochastic wideband impulse modulated self-noise of the pulsar in question.

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.

We present the results from nearly three years of monitoring of the variations in dispersion measure (DM) along the line-of-sight to 11 millisecond pulsars using the Giant Metrewave Radio Telescope (GMRT). These results demonstrate accuracies of single epoch DM estimates of the order of 5x10^(-4) cm^(-3) pc. A preliminary comparison with the Parkes Pulsar Timing Array (PPTA) data shows that the measured DM fluctuations are comparable. We show effects of DM variations due to the solar wind and solar corona and compare with the existing 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.

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.

Pulsars in close binary systems have provided some of the most stringent tests of strong-field gravity to date. The pulsar--white-dwarf binary system J1141-6545 is specifically interesting due to its gravitational asymmetry which makes it one of the most powerful probes of tensor-scalar theories of gravity. We give an overview of current gravitational tests provided by the J1141-6545 binary system and comment on how anomalous accelerations, geodetic precession and timing instabilities may be prevented from limiting future tests of gravity to come from this system.

We present 75 pulsars discovered in the mid-latitude portion of the High Time Resolution Universe survey, 54 of which have full timing solutions. All the pulsars have spin periods greater than 100 ms, and none of those with timing solutions are in binaries. Two display particularly interesting behaviour; PSR J1054-5944 is found to be an intermittent pulsar, and PSR J1809-0119 has glitched twice since its discovery. In the second half of the paper we discuss the development and application of an artificial neural network in the data-processing pipeline for the survey. We discuss the tests that were used to generate scores and find that our neural network was able to reject over 99% of the candidates produced in the data processing, and able to blindly detect 85% of pulsars. We suggest that improvements to the accuracy should be possible if further care is taken when training an artificial neural network; for example ensuring that a representative sample of the pulsar population is used during the training process, or the use of different artificial neural networks for the detection of different types of pulsars.

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.

(abridged) We report the results of a 10-year timing campaign on PSR J1738+0333, a 5.85-ms pulsar in a low-eccentricity 8.5-hour orbit with a low-mass white dwarf companion (...) The measurements of proper motion and parallax allow for a precise subtraction of the kinematic contribution to the observed orbital decay; this results in a significant measurement of the intrinsic orbital decay: (-25.9 +/- 3.2) \times 10^{-15} s/s. This is consistent with the orbital decay from the emission of gravitational waves predicted by general relativity, (-27.7 +1.5/-1.9) \times 10^{-15} s/s (...). This agreement introduces a tight upper limit on dipolar gravitational wave emission, a prediction of most alternative theories of gravity for asymmetric binary systems such as this. We use this limit to derive the most stringent constraints ever on a wide class of gravity theories, where gravity involves a scalar field contribution. When considering general scalar-tensor theories of gravity, our new bounds are more stringent than the best current solar-system limits over most of the parameter space, and constrain the matter-scalar coupling constant {\alpha}_0^2 to be below the 10^{-5} level. For the special case of the Jordan-Fierz-Brans-Dicke, we obtain the one-sigma bound {\alpha}_0^2 < 2 \times 10^{-5}, which is within a factor two of the Cassini limit. We also use our limit on dipolar gravitational wave emission to constrain a wide class of theories of gravity which are based on a generalization of Bekenstein's Tensor-Vector-Scalar gravity (TeVeS), a relativistic formulation of Modified Newtonian Dynamics (MOND).

Here we report on observations of the radio magnetar PSR J1622-4950 at frequencies from 1.4 to 17 GHz. We show that although its flux density is varying up to a factor of ~10 within a few days, it has on average decreased by a factor of 2 over the last 700 days. At the same time, timing analysis indicates a trend of decreasing spin-down rate over our entire data set, again of about a factor of 2 over 700 days, but also an erratic variability in the spin-down rate within this time span. Integrated pulse profiles are often close to 100 per cent linearly polarized, but large variations in both the profile shape and fractional polarization are regularly observed. Furthermore, the behaviour of the position angle of the linear polarization is very complex - offsets in both the absolute position angle and the phase of the position angle sweep are often seen and the occasional presence of orthogonal mode jumps further complicates the picture. However, model fitting indicates that the magnetic and rotation axes are close to aligned. Finally, a single pulse analysis has been carried out at four observing frequencies, demonstrating that the wide pulse profile is built up of narrow spikes of emission, with widths that scale inversely with observing frequency. All three of the known radio magnetars seem to have similar characteristics, with highly polarized emission, time-variable flux density and pulse profiles, and with spectral indices close to zero.

We report on the pulse-to-pulse energy distributions and phase-resolved modulation properties for catalogued pulsars in the southern High Time Resolution Universe intermediate-latitude survey. We selected the 315 pulsars detected in a single-pulse search of this survey, allowing a large sample unbiased regarding any rotational parameters of neutron stars. We found that the energy distribution of many pulsars is well-described by a log-normal distribution, with few deviating from a small range in log-normal scale and location parameters. Some pulsars exhibited multiple energy states corresponding to mode changes, and implying that some observed "nulling" may actually be a mode-change effect. PSRJ1900-2600 was found to emit weakly in its previously-identified "null" state. We found evidence for another state-change effect in two pulsars, which show bimodality in their nulling time scales; that is, they switch between a continuous-emission state and a single-pulse-emitting state. Large modulation occurs in many pulsars across the full integrated profile, with increased sporadic bursts at leading and trailing sub-beam edges. Some of these high-energy outbursts may indicate the presence of "giant pulse" phenomena. We found no correlation with modulation and pulsar period, age, or other parameters. Finally, the deviation of integrated pulse energy from its average value was generally quite small, despite the significant phase-resolved modulation in some pulsars; we interpret this as tenuous evidence of energy regulation between distinct pulsar sub-beams.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

Incoherent dedispersion is a computationally intensive problem that appears frequently in pulsar and transient astronomy. For current and future transient pipelines, dedispersion can dominate the total execution time, meaning its computational speed acts as a constraint on the quality and quantity of science results. It is thus critical that the algorithm be able to take advantage of trends in commodity computing hardware. With this goal in mind, we present analysis of the 'direct', 'tree' and 'sub-band' dedispersion algorithms with respect to their potential for efficient execution on modern graphics processing units (GPUs). We find all three to be excellent candidates, and proceed to describe implementations in C for CUDA using insight gained from the analysis. Using recent CPU and GPU hardware, the transition to the GPU provides a speed-up of 9x for the direct algorithm when compared to an optimised quad-core CPU code. For realistic recent survey parameters, these speeds are high enough that further optimisation is unnecessary to achieve real-time processing. Where further speed-ups are desirable, we find that the tree and sub-band algorithms are able to provide 3-7x better performance at the cost of certain smearing, memory consumption and development time trade-offs. We finish with a discussion of the implications of these results for future transient surveys. Our GPU dedispersion code is publicly available as a C library at: http://dedisp.googlecode.com/

Exploration of the time-domain radio sky has huge potential for advancing our knowledge of the dynamic universe. Past surveys have discovered large numbers of pulsars, rotating radio transients and other transient radio phenomena; however, they have typically relied upon off-line processing to cope with the high data and processing rate. This paradigm rules out the possibility of obtaining high-resolution base-band dumps of significant events or of performing immediate follow-up observations, limiting analysis power to what can be gleaned from detection data alone. To overcome this limitation, real-time processing and detection of transient radio events is required. By exploiting the significant computing power of modern graphics processing units (GPUs), we are developing a transient-detection pipeline that runs in real-time on data from the Parkes radio telescope. In this paper we discuss the algorithms used in our pipeline, the details of their implementation on the GPU and the challenges posed by the presence of radio frequency interference.

We present the discovery of six millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) survey for pulsars and fast transients carried out with the Parkes radio telescope. All six are in binary systems with approximately circular orbits and are likely to have white dwarf companions. PSR J1017-7156 has a high flux density and a narrow pulse width, making it ideal for precision timing experiments. PSRs J1446-4701 and J1125-5825 are coincident with gamma-ray sources, and folding the high-energy photons with the radio timing ephemeris shows evidence of pulsed gamma-ray emission. PSR J1502-6752 has a spin period of 26.7 ms, and its low period derivative implies that it is a recycled pulsar. The orbital parameters indicate it has a very low mass function, and therefore a companion mass much lower than usually expected for such a mildly recycled pulsar. In addition we present polarisation profiles for all 12 MSPs discovered in the HTRU survey to date. Similar to previous observations of MSPs, we find that many have large widths and a wide range of linear and circular polarisation fractions. Their polarisation profiles can be highly complex, and although the observed position angles often do not obey the rotating vector model, we present several examples of those that do. We speculate that the emission heights of MSPs are a substantial fraction of the light cylinder radius in order to explain broad emission profiles, which then naturally leads to a large number of cases where emission from both poles is observed.

We demonstrate that the sensitivity of high-precision pulsar timing experiments will be ultimately limited by the broadband intensity modulation that is intrinsic to the pulsar's stochastic radio signal. That is, as the peak flux of the pulsar approaches that of the system equivalent flux density, neither greater antenna gain nor increased instrumental bandwidth will improve timing precision. These conclusions proceed from an analysis of the covariance matrix used to characterise residual pulse profile fluctuations following the template matching procedure for arrival time estimation. We perform such an analysis on 25 hours of high-precision timing observations of the closest and brightest millisecond pulsar, PSR J0437-4715. In these data, the standard deviation of the post-fit arrival time residuals is approximately four times greater than that predicted by considering the system equivalent flux density, mean pulsar flux and the effective width of the pulsed emission. We develop a technique based on principal component analysis to mitigate the effects of shape variations on arrival time estimation and demonstrate its validity using a number of illustrative simulations. When applied to our observations, the method reduces arrival time residual noise by approximately 20%. We conclude that, owing primarily to the intrinsic variability of the radio emission from PSR J0437-4715 at 20 cm, timing precision in this observing band better than 30 - 40 ns in one hour is highly unlikely, regardless of future improvements in antenna gain or instrumental bandwidth. We describe the intrinsic variability of the pulsar signal as stochastic wideband impulse modulated self-noise (SWIMS) and argue that SWIMS will likely limit the timing precision of every millisecond pulsar currently observed by Pulsar Timing Array projects as larger and more sensitive antennae are built in the coming decades.

Millisecond pulsars are thought to be neutron stars that have been spun-up by accretion of matter from a binary companion. Although most are in binary systems, some 30% are solitary, and their origin is therefore mysterious. PSR J1719-1438, a 5.7 ms pulsar, was detected in a recent survey with the Parkes 64m radio telescope. We show that it is in a binary system with an orbital period of 2.2 h. Its companion's mass is near that of Jupiter, but its minimum density of 23 g cm$^{-3}$ suggests that it may be an ultra-low mass carbon white dwarf. This system may thus have once been an Ultra Compact Low-Mass X-ray Binary, where the companion narrowly avoided complete destruction.

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.

We present observations of PSRs J0437-4715, J0738-4042, J0835-4510, J0908-4913, J1048-5832, J1622-4950, J1644-4559, J1721-3532 and J1740-3015 at 17 GHz using the Parkes radio telescope. All 9 were detected at 17 GHz, additionally, we detected PSR J0835-4510 and J1622-4950 at 24 GHz. Polarisation profiles of each pulsar and the variation with frequency are discussed. In general, we find that the highly polarised edge components of young pulsars continue to dominate their profiles at 17 GHz. Older pulsars (>10^5 years) appear to be almost completely depolarised. Our detection of PSR J0437-4715 is the highest frequency observation of a millisecond pulsar to date, and implies a luminosity at 17 GHz of 14 {\mu}Jy kpc^2, and a mean spectral index of 2.2. We find that the spectral index of the magnetar PSR J1622-4950 is flat between 1.4 and 24 GHz, similar to the other known radio magnetars XTE J1810-197 and 1E 1547.0-5408. The profile is similar to that at 3.1 GHz, and is highly linearly polarised. Analysis of the frequency evolution of the profile of PSR J0835-4510 show that the profile is made of four components that vary with frequency only in their amplitude. The width and separation of the components remains fixed and the spectral index of each component can be determined independently.

We present the search methods and initial results for transient radio signals in the High Time Resolution Universe (HTRU) Survey. The HTRU survey's single-pulse search, the software designed to perform the search, and a determination of the HTRU survey's sensitivity to single pulses are described. Initial processing of a small fraction of the survey has produced 11 discoveries, all of which are sparsely-emitting neutron stars, as well as provided confirmation of two previously unconfirmed neutron stars. Most of the newly discovered objects lie in regions surveyed previously, indicating both the improved sensitivity of the HTRU survey observing system and the dynamic nature of the radio sky. The cycles of active and null states in nulling pulsars, rotating radio transients (RRATs), and long-term intermittent pulsars are explored in the context of determining the relationship between these populations, and of the sensitivity of a search to the various radio-intermittent neutron star populations. This analysis supports the case that many RRATs are in fact high-null-fraction pulsars (i. e. with null fraction >~0.95), and indicates that intermittent pulsars appear distinct from nulling pulsars in their activity cycle timescales. We find that in the measured population, there is a deficit of pulsars with typical emission timescales greater than ~300 s that is not readily explained by selection effects. The HTRU low-latitude survey will be capable of addressing whether this deficit is physical. We predict that the HTRU survey will explore pulsars with a broad range of nulling fractions (up to and beyond 0.999), and at its completion is likely to increase the currently known RRATs by a factor of more than two.

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 discovery of 5 millisecond pulsars found in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods from ~2.3 to ~7.5 ms, and all are in binary systems with orbital periods ranging from ~0.3 to ~150 d. In four of these systems, the most likely companion is a white dwarf, with minimum masses of ~0.2 Solar Masses. The other pulsar, J1731-1847, has a very low mass companion and exhibits eclipses, and is thus a member of the "black widow" class of pulsar binaries. These eclipses have been observed in bands centred near frequencies of 700, 1400 and 3000 MHz, from which measurements have been made of the electron density in the eclipse region. These measurements have been used to examine some possible eclipse mechanisms. The eclipse and other properties of this source are used to perform a comparison with the other known eclipsing and "black widow" pulsars. These new discoveries occupy a short-period and high-dispersion measure (DM) region of parameter space, which we demonstrate is a direct consequence of the high time and frequency resolution of the HTRU survey. The large implied distances to our new discoveries makes observation of their companions unlikely with both current optical telescopes and the Fermi Gamma-ray Space Telescope. The extremely circular orbits make any advance of periastron measurements highly unlikely. No relativistic Shapiro delays are obvious in any of the systems, although the low flux densities would make their detection difficult unless the orbits were fortuitously edge-on.

We have embarked on a survey for pulsars and fast transients using the 13-beam Multibeam receiver on the Parkes radio telescope. Installation of a digital backend allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels and sampled every 64 us. Limits of the receiver package restrict us to a 340 MHz observing band centred at 1352 MHz. The factor of eight improvement in frequency resolution over previous multibeam surveys allows us to probe deeper into the Galactic plane for short duration signals such as the pulses from millisecond pulsars. We plan to survey the entire southern sky in 42641 pointings, split into low, mid and high Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond pulsars. With ~30% of the mid-latitude survey complete, we have re-detected 223 previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. The newly discovered millisecond pulsars tend to have larger dispersion measures than those discovered in previous surveys, as expected from the improved time and frequency resolution of our instrument.

A survey of the Galactic plane in the region $-60\degree \leq l \leq 30\degree$, $|b| \leq 0.25\degree$ was carried out using the seven-beam Parkes Methanol Multibeam (MMB) receiver, which operates at a frequency of 6.5 GHz. Three pulsars were discovered, and 16 previously known pulsars detected. In this paper we present two previously-unpublished discoveries, both with extremely high dispersion measures, one of which is very close, in angular distance, to the Galactic centre. The survey data also contain the first known detection, at radio frequencies, of the radio magnetar PSR J1550-5418. Our survey observation was made 46 days prior to that previously published and places constraints on the beginning of pulsed radio emission from the source. The detection of only three previously undiscovered pulsars argues that there are few pulsars in the direction of the inner Galaxy whose flux density spectrum is governed by a flat power law. However, these pulsars would be likely to remain undetected at lower frequencies due to the large amount of scatter broadening which affects pulsars with high values of dispersion measure. Surveys with future telescopes at high observing frequencies will, therefore, play an important role in the discovery of pulsars at the Galactic centre. By simulating pulsar surveys of the Galaxy with Phase 1 SKA at frequencies of 1.4 GHz and 10 GHz, we find that high-frequency observations are the only way to discover and observe the Galactic-centre pulsar population.

Three years ago, the report of a solitary radio burst was thought to be the first discovery of a rare, impulsive event of unknown extragalactic origin (Lorimer et al. 2007). The extragalactic interpretation was based on the swept-frequency nature of the event, which followed the dispersive delay expected from an extragalactic pulse. We report here on the detection of 16 pulses, the bulk of which exhibit a frequency sweep with a shape and magnitude resembling the Lorimer Burst. These new events were detected in a sidelobe of the Parkes Telescope and are of clearly terrestrial origin, with properties unlike any known sources of terrestrial broad-band radio emission. The new detections cast doubt on the extragalactic interpretation of the original burst, and call for further sophistication in radio-pulse survey techniques to identify the origin of the anomalous terrestrial signals and definitively distinguish future extragalactic pulse detections from local signals. The ambiguous origin of these seemingly dispersed, swept-frequency signals suggest that radio-pulse searches using multiple detectors will be the only experiments able to provide definitive information about the origin of new swept-frequency radio burst detections.

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.

DSPSR is a high-performance, open-source, object-oriented, digital signal processing software library and application suite for use in radio pulsar astronomy. Written primarily in C++, the library implements an extensive range of modular algorithms that can optionally exploit both multiple-core processors and general-purpose graphics processing units. After over a decade of research and development, DSPSR is now stable and in widespread use in the community. This paper presents a detailed description of its functionality, justification of major design decisions, analysis of phase-coherent dispersion removal algorithms, and demonstration of performance on some contemporary microprocessor architectures.

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.

As part of a survey for radio pulsars with the Parkes 64-m telescope we have discovered PSR J1622-4950, a pulsar with a 4.3-s rotation period. Follow-up observations show that the pulsar has the highest inferred surface magnetic field of the known radio pulsars (B ~ 3e14 G), exhibits significant timing noise and appears to have an inverted spectrum. Unlike the vast majority of the known pulsar population, PSR J1622-4950 appears to switch off for many hundreds of days and even in its on-state exhibits extreme variability in its flux density. Furthermore, the integrated pulse profile changes shape with epoch. All of these properties are remarkably similar to the only two magnetars previously known to emit radio pulsations. The position of PSR J1622-4950 is coincident with an X-ray source that, unlike the other radio pulsating magnetars, was found to be in quiescence. We conclude that our newly discovered pulsar is a magnetar - the first to be discovered via its radio emission.

The work presented here examines populations of double compact binary systems and tidally enhanced collapsars. We make use of BINPOP and BINKIN, two components of a recently developed population synthesis package. Results focus on correlations of both binary and spatial evolutionary population characteristics. Pulsar and long duration gamma-ray burst observations are used in concert with our models to draw the conclusions that: double neutron star binaries can merge rapidly on timescales of a few million years (much less than that found for the observed double neutron star population), common envelope evolution within these models is a very important phase in double neutron star formation, and observations of long gamma-ray burst projected distances are more centrally concentrated than our simulated coalescing double neutron star and collapsar Galactic populations. Better agreement is found with dwarf galaxy models although the outcome is strongly linked to the assumed birth radial distribution. The birth rate of the double neutron star population in our models range from 4-160 Myr^-1 and the merger rate ranges from 3-150 Myr^-1. The upper and lower limits of the rates results from including electron capture supernova kicks to neutron stars and decreasing the common envelope efficiency respectively. Our double black hole merger rates suggest that black holes should receive an asymmetric kick at birth.

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

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.

We present the results of a search for transient radio bursts of between 0.125 and 32 millisecond duration in two archival pulsar surveys of intermediate galactic latitudes with the Parkes multibeam receiver. Fourteen new neutron stars have been discovered, seven of which belong to the recently identified "rotating radio transients" (RRATs) class. Here we describe our search methodology, and discuss the new detections in terms of how the RRAT population relates to the general population of pulsars. The new detections indicate (1) that the galactic z-distribution of RRATs in the surveys closely resembles the distribution of pulsars, with objects up to 0.86 kpc from the galactic plane; (2) where measurable, the RRAT pulse widths are similar to that of individual pulses from pulsars of similar period, implying a similar beaming fraction; and (3) our new detections span a variety of nulling fractions, and thus we postulate that the RRATs may simply be nulling pulsars that are only "on" for less than a pulse period. Finally, the newly discovered object PSR J0941-39 may represent a link between pulsars and RRATs. This bizarre object was discovered as an RRAT, but in follow-up observations often appeared as a bright (~10 mJy) pulsar with a low nulling fraction. It is obvious therefore that a neutron star can oscillate between being an RRAT and a pulsar. Crucially, the sites of the RRAT pulses are coincident with the pulsar's emission, implying that the two emission mechanisms are linked, and that RRATs are not just pulsars observed from different orientations.

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.

Using Chandra data taken on 2008 June, we detected pulsations at 2.59439(4) s in the soft gamma-ray repeater SGR 1627-41. This is the second measurement of the source spin period and allows us to derive for the first time a long-term spin-down rate of (1.9 +/- 0.4)E-11 s/s. From this value we infer for SGR 1627-41 a characteristic age of 2.2 kyr, a spin-down luminosity of 4E+34 erg/s (one of the highest among sources of the same class), and a surface dipole magnetic field strength of 2E+14 G. These properties confirm the magnetar nature of SGR 1627-41; however, they should be considered with caution since they were derived on the basis of a period derivative measurement made using two epochs only and magnetar spin-down rates are generally highly variable. The pulse profile, double-peaked and with a pulsed fraction of (13 +/- 2)% in the 2-10 keV range, closely resembles that observed by XMM-Newton in 2008 September. Having for the first time a timing model for this SGR, we also searched for a pulsed signal in archival radio data collected with the Parkes radio telescope nine months after the previous X-ray outburst. No evidence for radio pulsations was found, down to a luminosity level 10-20 times fainter (for a 10% duty cycle and a distance of 11 kpc) than the peak luminosity shown by the known radio magnetars.

Accurate measurement of pulsar distances via astrometry using very long baseline interferometry enables the improvement of Galactic electron density distribution models, improving distance estimates for the vast majority of pulsars for which parallax measurements are unavailable. However, pulsars at southern declinations have been under-represented in previous interferometric astrometry campaigns. In order to redress this imbalance, we have conducted a two-year astrometric campaign targeting eight southern pulsars with the Australian Long Baseline Array. The program summarized in this paper has resulted in the measurement of seven new pulsar parallaxes, with success on objects down to a mean flux density of 0.8 mJy at 1600 MHz. Our results highlight the substantial uncertainties that remain when utilizing free electron density models for individual pulsar distances. Until this study, PSR J0630-2834 was believed to convert 16% of its spin-down energy into x-rays, but our measured parallax distance of 332 (+52 -40) pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to be almost a factor of two more distant than previously thought, making its conversion of spin-down energy to x-rays the most efficient known (>1%). The 8.5 second radio pulsar J2144-3933 was found to be closer than previously predicted, making its apparent 1400 MHz radio luminosity the lowest of any known pulsar (20 microJy kpc^2). We have examined the growing population of neutron stars with accurate parallaxes to determine the effect of distance errors on the underlying neutron star velocity distribution, and find that typical distance errors may be biasing the estimated mean pulsar velocity upwards by 5%, and are likely to exaggerate the distribution's high-velocity tail.

The double pulsar J0737-3039A/B is a unique system with which to test gravitational theories in the strong-field regime. However, the accuracy of such tests will be limited by knowledge of the distance and relative motion of the system. Here we present very long baseline interferometry observations which reveal that the distance to PSR J0737-3039A/B is 1150+220-160 pc, more than double previous estimates, and confirm its low transverse velocity (~9 km/s). Combined with a decade of pulsar timing, these results will allow tests of gravitational radiation emission theories at the 0.01% level, putting stringent constraints on theories which predict dipolar gravitational radiation. They also allow insight into the system's formation and the source of its high-energy emission.

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.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

Using the recently upgraded Long Baseline Array, we have measured the trigonometric parallax of PSR J0437-4715 to better than 1% precision, the most precise pulsar distance determination made to date. Comparing this VLBI distance measurement to the kinematic distance obtained from pulsar timing, which is calculated from the pulsar's proper motion and apparent rate of change of orbital period, gives a precise limit on the unmodeled relative acceleration between the Solar System and PSR J0437-4715, which can be used in a variety of applications. Firstly, it shows that Newton's gravitational constant G is stable with time (\dot{G}/G = (-5 +- 26) x 10^{-13} yr^{-1}, 95% confidence). Secondly, if a stochastic gravitational wave background existed at the currently quoted limit, this null result would fail ~50% of the time. Thirdly, it excludes Jupiter-mass planets within 226 AU of the Sun in 50% of the sky (95% confidence). Finally, the ~1% agreement of the parallax and orbital period derivative distances provides a fundamental confirmation of the parallax distance method upon which all astronomical distances are based.

We present the discovery of PSR J1410-6132, a 50-ms pulsar found during a high-frequency survey of the Galactic plane, using a 7-beam 6.3-GHz receiver on the 64-m Parkes radio telescope. The pulsar lies within the error box of the unidentified EGRET source 3EG J1410-6147, has a characteristic age of 26 kyr and a spin-down energy of 10^37 erg s^-1. It has a very high dispersion measure of 960+/-10 cm^-3 pc and the largest rotation measure of any pulsar, RM=+2400 +/- 30 rad m^-2. The pulsar is very scatter-broadened at frequencies of 1.4 GHz and below, making pulsed emission almost impossible to detect. Assuming a distance of 15 kpc, the pulsar's spin-down energy and a gamma-ray efficiency factor of ~10 per cent is sufficient to power the gamma-ray source. We therefore believe we have identified the nature of 3EG J1410-6147. This new discovery suggests that deep targeted high-frequency surveys of inner-galaxy EGRET sources could uncover further young, energetic pulsars.

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad-hoc manner. Thus no complete and direct comparison with observations of the pulsar population within the Galactic disk has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed \textit{typical} MSP formation process. Results also focus on reproducing the observed $P\dot{P}$ diagram for Galactic pulsars and how this precludes short timescales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.

Pulsars in close binary orbit around another neutron star or a massive white dwarf make ideal laboratories for testing the predictions of gravitational radiation and self-gravitational effects. We report new timing measurements of the pulsar-white-dwarf binary PSR J1141-6545, providing strong evidence that such asymmetric systems have gravitational wave losses that are consistent with general relativity. The orbit is found to be decaying at a rate of $1.04\pm0.06$ times the general relativistic prediction and the Shapiro delay is consistent with the orbital inclination angle derived from scintillation measurements. The system provides a unique test-bed for tensor-scalar theories of gravity; our current measurements place stringent constraints in the theory space, with a limit of $\alpha_0^2 < 2.1 \times 10^{-5}$ for weakly non-linear coupling and an asymptotic limit of $\alpha_0^2 < 3.4 \times 10^{-6}$ for strongly non-linear coupling, where $\alpha_0$ is the linear coupling strength of matter to an underlying scalar field. This asymptotic limit is nearly three times smaller than the Cassini bound ($\alpha_0^2 \approx 10^{-5}$).

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.

The majority of pulsar population synthesis studies performed to date have focused on isolated pulsar evolution. Those that have incorporated pulsar evolution within binary systems have tended to either treat binary evolution poorly of evolve the pulsar population in an ad-hoc manner. Here we present the first model of the Galactic field pulsar population that includes a comprehensive treatment of both binary and pulsar evolution. Synthetic observational surveys mimicking a variety of radio telescopes are then performed on this population. As such, a complete and direct comparison of model data with observations of the pulsar population within the Galactic disk is now possible. The tool used for completing this work is a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Here we give a brief overview of the method and assumptions involved with each component. Some preliminary results are also presented as well as plans for future applications of the code.

We describe in this paper a new binary millisecond pulsar, PSR J1738+0333. Using Arecibo, we have achieved good timing accuracy for this object, about 220 ns for 1-hour integrations over 100 MHz. This allowed us to measure a precise proper motion, parallax and orbital parameters for this system. We highlight the system's potential for constraining alternative theories of gravitation.

Pulsar surveys offer one of the few opportunities to monitor even a small fraction (~0.00001) of the radio sky for impulsive burst-like events with millisecond durations. In analysis of archival survey data, we have discovered a 30-Jy dispersed burst of duration <5 ms located three degrees from the Small Magellanic Cloud. The burst properties argue against a physical association with our Galaxy or the Small Magellanic Cloud. Current models for the free electron content in the Universe imply a distance to the burst of <1 Gpc No further bursts are seen in 90-hr of additional observations, implying that it was a singular event such as a supernova or coalescence of relativistic objects. Hundreds of similar events could occur every day and act as insightful cosmological probes.

Binary and Millisecond pulsars have a great deal to teach us about stellar evolution and are invaluable tools for tests of relativistic theories of gravity. Our understanding of these objects has been transformed by large-scale surveys that have uncovered a great deal of new objects, exquisitely timed by ever-improving instrumentation. Here we argue that there exists a fundamental relation between the spin period of a pulsar and its companion mass, and that this determines many of the observable properties of a binary pulsar. No recycled pulsars exist in which the minimum companion mass exceeds (P/10 ms) Solar Masses. Furthermore, the three fastest disk millisecond pulsars are either single, or possess extremely low-mass companions Mc~0.02 Mo, consistent with this relation. Finally, the four relativistic binaries for which we have actual measurements of neutron star masses, suggest that not only are their spin periods related to the companion neutron star mass, but that the kick imparted to the system depends upon it too, leading to a correlation between orbital eccentricity and spin period. The isolation of the relativistic binary pulsars in the magnetic field-period diagram is used to argue that this must be because the kicks imparted to proto-relativistic systems are usually small, leading to very few if any isolated runaway mildly-recycled pulsars. This calls into question the magnitude of supernova kicks in close binaries, which have been usually assumed to be similar to those imparted to the bulk of the pulsar population. Finally, we review some of the highlights of the Parkes precision timing efforts, which suggest 10 nanosecond timing is obtainable on PSR J1909-3744 that will aid us in searching for cosmological sources of gravitational waves.

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.

We describe the development of an FX style correlator for Very Long Baseline Interferometry (VLBI), implemented in software and intended to run in multi-processor computing environments, such as large clusters of commodity machines (Beowulf clusters) or computers specifically designed for high performance computing, such as multi-processor shared-memory machines. We outline the scientific and practical benefits for VLBI correlation, these chiefly being due to the inherent flexibility of software and the fact that the highly parallel and scalable nature of the correlation task is well suited to a multi-processor computing environment. We suggest scientific applications where such an approach to VLBI correlation is most suited and will give the best returns. We report detailed results from the Distributed FX (DiFX) software correlator, running on the Swinburne supercomputer (a Beowulf cluster of approximately 300 commodity processors), including measures of the performance of the system. For example, to correlate all Stokes products for a 10 antenna array, with an aggregate bandwidth of 64 MHz per station and using typical time and frequency resolution presently requires of order 100 desktop-class compute nodes. Due to the effect of Moore's Law on commodity computing performance, the total number and cost of compute nodes required to meet a given correlation task continues to decrease rapidly with time. We show detailed comparisons between DiFX and two existing hardware-based correlators: the Australian Long Baseline Array (LBA) S2 correlator, and the NRAO Very Long Baseline Array (VLBA) correlator. In both cases, excellent agreement was found between the correlators. Finally, we describe plans for the future operation of DiFX on the Swinburne supercomputer, for both astrophysical and geodetic science.

We present five recycled pulsars discovered during a 21-cm survey of approximately 4,150 deg^2 between 15 deg and 30 deg from the galactic plane using the Parkes radio telescope. One new pulsar, PSR J1528-3146, has a 61 ms spin period and a massive white dwarf companion. Like many recycled pulsars with heavy companions, the orbital eccentricity is relatively high (~0.0002), consistent with evolutionary models that predict less time for circularization. The four remaining pulsars have short spin periods (3 ms < P < 6 ms); three of these have probable white dwarf binary companions and one (PSR J2010-1323) is isolated. PSR J1600-3053 is relatively bright for its dispersion measure of 52.3 pc cm^-3 and promises good timing precision thanks to an intrinsically narrow feature in its pulse profile, resolvable through coherent dedispersion. In this survey, the recycled pulsar discovery rate was one per four days of telescope time or one per 600 deg^2 of sky. The variability of these sources implies that there are more millisecond pulsars that might be found by repeating this survey.

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.

We have searched for microsecond bursts of emission from millisecond pulsars in the globular cluster M28 using the Parkes radio telescope. We detected a total of 27 giant pulses from the known emitter PSR J1824-2452A. At wavelengths around 20 cm the giant pulses are scatter-broadened to widths of around 2 microseconds and follow power-law statistics. The pulses occur in two narrow phase-windows which correlate in phase with X-ray emission and trail the peaks of the integrated radio pulse-components. Notably, the integrated radio emission at these phase windows has a steeper spectral index than other emission. The giant pulses exhibit a high degree of polarization, with many being 100% elliptically polarized. Their position angles appear random. Although the integrated emission of PSR J1824-2452A is relatively stable for the frequencies and bandwidths observed, the intensities of individual giant pulses vary considerably across our bands. Two pulses were detected at both 2700 and 3500 MHz. The narrower of the two pulses is 20 ns wide at 3500 MHz. At 2700 MHz this pulse has an inferred brightness temperature at maximum of 5 x 10^37 K. Our observations suggest the giant pulses of PSR J1824-2452A are generated in the same part of the magnetosphere as X-ray emission through a different emission process to that of ordinary pulses.

We have conducted a search for giant pulses from four millisecond pulsars using the 100m Green Bank Telescope. Coherently dedispersed time-series from PSR J0218+4232 were found to contain giant pulses of very short intrinsic duration whose energies follow power-law statistics. The giant pulses are in phase with the two minima of the radio integrated pulse profile but are phase aligned with the peaks of the X-ray profile. Historically, individual pulses more than 10-20 times the mean pulse energy have been deemed to be ``giant pulses''. As only 4 of the 155 pulses had energies greater than 10 times the mean pulse-energy, we argue the emission mechanism responsible for giant pulses should instead be defined through: (a) intrinsic timescales of microsecond or nanosecond duration; (b) power-law energy statistics; and (c) emission occurring in narrow phase-windows coincident with the phase windows of non-thermal X-ray emission. Four short-duration pulses with giant-pulse characteristics were also observed from PSR B1957+20. As the inferred magnetic fields at the light cylinders of the millisecond pulsars that emit giant pulses are all very high, this parameter has previously been considered to be an indicator of giant pulse emissivity. However, the frequency of giant pulse emission from PSR~B1957+20 is significantly lower than for other millisecond pulsars that have similar magnetic fields at their light cylinders. This suggests that the inferred magnetic field at the light cylinder is a poor indicator of the rate of emission of giant pulses.

We report on nearly two years of timing observations of the low-mass binary millisecond pulsar, PSR J1909-3744 with the Caltech-Parkes-Swinburne Recorder II (CPSR2), a new instrument that gives unprecedented timing precision. Daily observations give a weighted rms residual of 74 ns, indicating an extremely low level of systematic error. We have greatly improved upon the previous parallax and proper motion measurements of PSR J1909-3744, yielding a distance of $1.14^{+0.04}_{-0.03}$ kpc and transverse velocity of $(200^{+7}_{-6})$ km s$^{-1}$. The system's orbital eccentricity is just 1.35(12)$\times10^{-7}$, the smallest yet recorded. Since their discovery, the masses of the rapidly rotating millisecond pulsars have remained a mystery, with the recycling hypothesis arguing for heavy objects, and the accretion-induced collapse of a white dwarf more consistent with neutron stars less than the Chandrashkar limit. Fortuitously, PSR J1909-3744 is an edge-on system, and our data have allowed the measurement of the range and shape of the Shapiro delay to high accuracy, giving the first precise determination of a millisecond pulsar mass to date, $m_p = (1.438 \pm 0.024) M_\odot$. The mass of PSR J1909-3744 is at the upper edge of the range observed in mildly recycled pulsars in double neutron star systems, consistent with the the recycling hypothesis. It appears that the production of millisecond pulsars is possible with the accretion of less than $0.2 M_\odot$.

We demonstrate that a natural consequence of an asymmetric kick imparted to neutron stars at birth is that the majority of double neutron star binaries should possess highly eccentric orbits. This leads to greatly accelerated orbital decay, due to the enormous increase in the emission of gravitational radiation at periastron as originally demonstrated by Peters (1964). A uniform distribution of kick velocities constrained to the orbital plane would result in ~24% of surviving binaries coalescing at least 10,000 times faster than an unperturbed circular system. Even if the planar kick constraint is lifted, ~6% of bound systems still coalesce this rapidly. In a non-negligible fraction of cases it may even be possible that the system could coalesce within 10 years of the final supernova, resulting in what we might term a "double supernova''. For systems resembling the progenitor of PSR J0737-3039A, this number is as high as \~9% (in the planar kick model). Whether the kick velocity distribution extends to the range required to achieve this is still unclear. We do know that the observed population of binary pulsars has a deficit of highly eccentric systems at small orbital periods. In contrast, the long-period systems favour large eccentricities, as expected. We argue that this is because the short-period highly eccentric systems have already coalesced and are thus selected against by pulsar surveys. This effect needs to be taken into account when using the scale-factor method to estimate the coalescence rate of double neutron star binaries. We therefore assert that the coalesence rate of such binaries is underestimated by a factor of several.

We have searched for microsecond-timescale broadband emission from a sample of eighteen millisecond pulsars. Our study places strong limits on such emission from several millisecond pulsars and shows that it is only present in a small subset of millisecond pulsars. Giant pulses of up to 64 times the mean pulse energy were detected from PSR J1823-3021A in the globular cluster NGC 6624. In contrast to the giant pulses of PSR B1937+21, nearly all of the giant pulses from PSR J1823-3021A were distributed within the trailing half of the main-pulse component of the integrated pulse profile. The fact that no giant pulses were observed on the leading side of the main-pulse component suggests that giant pulses are preferentially emitted closer to the last open field line than ordinary emission. The correlation between giant pulse emissivity and spin-down luminosity in millisecond pulsars suggests that the high period derivative of PSR J1823-3021A is intrinsic and is not just an artifact of its acceleration in the gravitational potential of NGC 6624.

We present observations that show dramatic evolution of the mean pulse profile of the relativistic binary pulsar J1141--6545 over a period of five years. This is consistent with precession of the pulsar spin axis due to relativistic spin-orbit coupling. Observations made between 1999 and 2004 with a number of instruments at the Parkes radio telescope demonstrate a steady, secular evolution of the mean total intensity profile, which increases in width by more than 50 percent during the five year period. Analysis of the changing position angle of the linearly polarised component of the mean profile suggests that our line of sight is shifting closer to the core of the emission cone. We find that the slope of the position angle swing across the centre of the pulse steepens with time and use a simplified version of the rotating vector model to constrain the magnitude and direction of the change in our line of sight angle relative to the pulsar magnetic axis. The fact that we appear to be moving deeper into the emission cone is consistent with the non-detection of this pulsar in previous surveys.

We present an investigation of the morphology and arrival times of integrated radio pulses from the binary millisecond pulsar PSR J1022+1001. This pulsar is renowned for its poor timing properties, which have been postulated to originate from variability in its average pulse profile. Although a sub-class of long-period pulsars are known to exhibit mode changes that give rise to very large deviations in their integrated profiles, this was the first millisecond pulsar thought to have an unstable mean profile. As part of a precision timing program at the Parkes radio telescope we observed this pulsar between January 2003 and March 2004 using a coherent de-dispersion system (CPSR2). A study of morphological variability during our brightest observations suggests that the pulse profile varies by at most a few percent, similar to the uncertainty in our calibration. Unlike previous authors, we find that this pulsar times extremely well. In five minute integrations of 64 MHz bands we obtain a weighted RMS residual of just 2.27 microseconds. The reduced chi-squared of our best fit is 1.43, which suggests that this pulsar can be timed to high accuracy with standard cross-correlation techniques. Combining relativistic constraints with the pulsar mass function and consideration of the Chandrasekhar mass limit on the white dwarf companion, we can constrain the inclination angle of the system to lie within the range 37 < i < 56 degrees. For reasonable pulsar masses, this suggests that the white dwarf is at least 0.9 solar masses. We also find evidence for secular evolution of the projected semi-major axis.

We report the discovery of PSR J1909-3744, a 2.95 millisecond pulsar in a nearly circular 1.53 day orbit. Its narrow pulse width of 43 microseconds allows pulse arrival times to be determined with great accuracy. We have spectroscopically identified the companion as a moderately hot (T = 8500 K) white dwarf with strong absorption lines. Radial velocity measurements of the companion will yield the mass ratio of the system. Our timing data suggest the presence of Shapiro delay; we expect that further timing observations, combined with the mass ratio, will allow the first accurate determination of a millisecond pulsar mass. We have measured the timing parallax and proper motion for this pulsar which indicate a transverse velocity of 140 (+80/-40) km/s. This pulsar's stunningly narrow pulse profile makes it an excellent candidate for precision timing experiments that attempt to detect low frequency gravitational waves from coalescing supermassive black hole binaries.

Here we report timing measurements of the relativistic binary pulsar PSR J1141-6545 that constrain the component masses and demonstrate that the orbital period derivative \dot Pb = (-4+/-1)x10^-13 is consistent with gravitational wave emission as described by the general theory of relativity. The mass of the neutron star and its companion are 1.30+/-0.02 Mo and 0.986+/-0.020 Mo respectively, suggesting a white dwarf companion, and extending the range of systems for which general relativity provides a correct description. On evolutionary grounds, the progenitor mass of PSR J1141-6545 should be near the minimum for neutron star production. Its mass is two standard deviations below the mean of the other neutron stars, suggesting a relationship between progenitor and remnant masses.

We have obtained an HI absorption spectrum of the relativistic binary PSR J1141-6545 and used it to constrain the distance to the system. The spectrum suggests that the pulsar is at, or beyond, the tangent point, estimated to be at 3.7 kpc. PSR J1141-6545 offers the promise of stringent tests of General Relativity (GR) by comparing its observed orbital period derivative with that derived from other relativistic observables. At the distance of PSR J1141-6545 it should be possible to verify GR to an accuracy of just a few percent, as contributions to the observed orbital period derivative from kinematic terms will be a small fraction of that induced by the emission of gravitational radiation. PSR J1141-6545 will thus make an exceptional gravitational laboratory.

We report a dramatic orbital modulation in the scintillation timescale of the relativistic binary pulsar J1141--6545 that both confirms the validity of the scintillation speed methodology and enables us to derive important physical parameters. We have determined the space velocity, the orbital inclination and even the longitude of periastron of the binary system, which we find to be in good agreement with that obtained from pulse timing measurements. Our data permit two equally-significant physical interpretations of the system. The system is either an edge-on binary with a high space velocity ($\sim 115$ km s$^{-1}$) or is more face-on with a much slower velocity ($\sim 45$ km s$^{-1}$). We favor the former, as it is more consistent with pulse timing and the distribution of known neutron star masses. Under this assumption, the runaway velocity of 115 km s$^{-1}$ is much greater than is expected if pulsars do not receive a natal kick at birth. The derived inclination of the binary system is (76\pm 2.5^{\circ}) degrees, implying a companion mass of 1.01 (\pm )~0.02 M(_{\odot}) and a pulsar mass of 1.29 (\pm)~0.02 M(_{\odot}). Our derived physical parameters indicate that this pulsar should prove to be an excellent laboratory for tests of gravitational wave emission.

We have detected the rare phenomenon of stable, drifting sub-pulse behaviour in two pulsars discovered in the recent Swinburne intermediate latitude pulsar survey. The pulsars; PSR J1231-47 and PSR J1919+0134, have approximate periods (P) of 1.873 and 1.6039 seconds respectively. Both pulsars have multi-component profiles, and distinct drifting is observed across them. We have identified a single drift mode in both pulsars: the drift rate for PSR J1231-47 being 5.4(1) ms P^{-1} and 5.8(2) ms P^{-1} for PSR 1919+0134. The drifting is linear across the profile with no departure from linearity at the edges within the sensitivity of our observations.

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 have conducted a search of 19 southern Galactic globular clusters for sub-millisecond pulsars at 660 MHz with the Parkes 64-m radio telescope. To minimize dispersion smearing we used the CPSR baseband recorder, which samples the 20 MHz observing band at the Nyquist rate. By possessing a complete description of the signal we could synthesize an optimal filterbank in software, and in the case of globular clusters of known dispersion measure, much of the dispersion could be removed using coherent techniques. This allowed for very high time resolution (25.6 us in most cases), making our searches in general sensitive to sub-millisecond pulsars with flux densities greater than about 3 mJy at 50 cm. No new pulsars were discovered, placing important constraints on the proportion of pulsars with very short spin periods in these clusters.

We have conducted a survey of intermediate Galactic latitudes using the 13-beam 21-cm multibeam receiver of the Parkes 64-m radio telescope. The survey covered the region enclosed by 5 deg < |b| < 15 deg and -100 deg < l < 50 deg with 4,702 processed pointings of 265 s each, for a total of 14.5 days of integration time. Thirteen 2x96-channel filterbanks provided 288 MHz of bandwidth at a centre frequency of 1374 MHz, one-bit sampled every 125 microsec and incurring ~DM/13.4 cm^-3 pc samples of dispersion smearing. The system was sensitive to slow and most millisecond pulsars in the region with flux densities greater than approximately 0.3--1.1 mJy. Offline analysis on the 64-node Swinburne workstation cluster resulted in the detection of 170 pulsars of which 69 were new discoveries. Eight of the new pulsars, by virtue of their small spin periods and period derivatives, may be recycled and have been reported elsewhere. The slow pulsars discovered are typical of those already known in the volume searched, being of intermediate to old age. Several pulsars experience pulse nulling and two display very regular drifting sub-pulses. We discuss the new discoveries and provide timing parameters for the 48 slow pulsars for which we have a phase-connnected solution.

We present the timing parameters of nine pulsars discovered in a survey of intermediate Galactic latitudes at 1400 MHz with the Parkes radio telescope. Eight of these pulsars possess small pulse periods and period derivatives thought to be indicative of ``recycling''. Six of the pulsars are in circular binary systems, including two with relatively massive white dwarf companions. We discuss the implications of these new systems for theories of binary formation and evolution. One long-period pulsar (J1410-7404) has a moderately weak magnetic field and an exceedingly narrow average pulse profile, similar to other recycled pulsars.

We present high time resolution observations of single pulses from the Vela pulsar (PSR B0833-45) made with a baseband recording system at observing frequencies of 660 and 1413 MHz. We have discovered two startling features in the 1413 MHz single pulse data. The first is the presence of giant micro-pulses which are confined to the leading edge of the pulse profile. One of these pulses has a peak flux density in excess of 2500 Jy, more than 40 times the integrated pulse peak. The second new result is the presence of a large amplitude gaussian component on the trailing edge of the pulse profile. This component can exceed the main pulse in intensity but is switched on only relatively rarely. Fluctutation spectra reveal a possible periodicity in this feature of 140 pulse periods. Unlike the rest of the profile, this component has low net polarization and emits predominantly in the orthogonal mode. This feature appears to be unique to the Vela pulsar. We have also detected microstructure in the Vela pulsar for the first time. These same features are present in the 660 MHz data. We suggest that the full width of the Vela pulse profile might be as large as 10 ms but that the conal edges emit only rarely.

We have discovered two recycled pulsars in relativistic orbits as part of the first high-frequency survey of intermediate Galactic latitudes. PSR J1157-5112 is a 44 ms pulsar and the first recycled pulsar with an ultra-massive (M > 1.14 Mo) white dwarf companion. Millisecond pulsar J1757-5322 is a relativistic circular-orbit system which will coalesce due to the emission of gravitational radiation in less than 9.5 Gyr. Of the ~40 known circular orbit pulsars, J1757-5322 and J1157-5112 have the highest projected orbital velocities. There are now three local neutron-star/white-dwarf binaries that will coalesce in less than a Hubble time, implying a large coalescence rate for these objects in the local Universe. Systems such as J1141-6545 (Kaspi et al. 2000) are potential gamma-ray burst progenitors and dominate the coalescence rate, whilst lighter systems make excellent progenitors of millisecond pulsars with planetary or ultra-low mass companions.

Timing observations of 40 mostly young pulsars using the ATNF Parkes radio telescope between 1990 January and 1998 December are reported. In total, 20 previously unreported glitches and ten other glitches were detected in 11 pulsars. These included 12 glitches in PSR J1341$- $6220, corresponding to a glitch rate of 1.5 glitches per year. We also detected the largest known glitch, in PSR J1614$-$5047, with $\Delta\nu_g/\nu \approx 6.5 \times 10^{-6}$ where $\nu = 1/P$ is the pulse frequency. Glitch parameters were determined both by extrapolating timing solutions to inter-glitch intervals and by phase-coherent timing fits across the glitch(es). Analysis of glitch parameters, both from this work and from previously published results, shows that most glitches have a fractional amplitude $\Delta\nu_g/\nu$ of between $10^{-8}$ and $10^{-6}$. There is no consistent relationship between glitch amplitude and the time since the previous glitch or the time to the following glitch, either for the ensemble or for individual pulsars. As previously recognised, the largest glitch activity is seen in pulsars with ages of order 10$^4$ years, but for about 30 per cent of such pulsars, no glitches were detected in the 8-year data span. There is some evidence for a new type of timing irregularity in which there is a significant increase in pulse frequency over a few days, accompanied by a decrease in the magnitude of the slowdown rate. Fits of an exponential recovery to post-glitch data show that for most older pulsars, only a small fraction of the glitch decays. In some younger pulsars, a large fraction of the glitch decays, but in others, there is very little decay.

We present the annual trigonometric parallax of PSR J1744-1134 derived from an analysis of pulse times of arrival. The measured parallax, pi = 2.8+/-0.3 mas ranks among the most precisely determined distances to any pulsar. The parallax distance of 357+/-39 pc is over twice that derived from the dispersion measure using the Taylor & Cordes model for the Galactic electron distribution. The mean electron density in the path to the pulsar, n_e = (0.0088 +/- 0.0009) cm^{-3}, is the lowest for any disk pulsar. We have compared the n_e for PSR J1744-1134 with those for another 11 nearby pulsars with independent distance estimates. We conclude that there is a striking asymmetry in the distribution of electrons in the local interstellar medium. The electron column densities for pulsars in the third Galactic quadrant are found to be systematically higher than for those in the first. The former correlate with the position of the well known local HI cavity in quadrant three. The excess electrons within the cavity may be in the form of HII clouds marking a region of interaction between the local hot bubble and a nearby superbubble.

We present results of Compton Gamma-Ray Observatory EGRET observations of the unidentified high-energy gamma-ray sources 2EG J1049-5847 (GEV J1047-5840, 3EG J1048-5840) and 2EG J1103-6106 (3EG J1102-6103). These sources are spatially coincident with the young, energetic radio pulsars PSRs B1046-58 and J1105-6107, respectively. We find evidence for an association between PSR B1046-58 and 2EG J1049-5847. The gamma-ray pulse profile, obtained by folding time-tagged photons having energies above 400 MeV using contemporaneous radio ephemerides, has probability of arising by chance of 1.2E-4 according to the binning-independent H-test. A spatial analysis of the on-pulse photons reveals a point source of equivalent significance 10.2 sigma. Off-pulse, the significance drops to 5.8 sigma. Archival ASCA data show that the only hard X-ray point source in the 95% confidence error box of the gamma-ray source is spatially coincident with the pulsar within the 1' uncertainty (Pivovaroff, Kaspi & Gotthelf 1999). The double peaked gamma-ray pulse morphology and leading radio pulse are similar to those seen for other gamma-ray pulsars and are well-explained in models in which the gamma-ray emission is produced in charge-depleted gaps in the outer magnetosphere. The inferred pulsed gamma-ray flux above 400 MeV, (2.5 +/- 0.6) x 10E-10 erg/cm^2/s, represents 0.011 +/- 0.003 of the pulsar's spin-down luminosity, for a distance of 3 kpc and 1 sr beaming. For PSR J1105-6107, light curves obtained by folding EGRET photons using contemporaneous radio ephemerides show no significant features. We conclude that this pulsar converts less than 0.014 of its spin-down luminosity into E > 100 MeV gamma-rays beaming in our direction (99% confidence), assuming a distance of 7 kpc, 1 sr beaming and a duty cycle of 0.5.

We report on the first firm detection of pulsed gamma-ray emission from PSR B1509-58 in the 0.75-30 MeV energy range in CGRO COMPTEL data collected over more than 6 years. The modulation significance in the 0.75-30 MeV pulse-phase distribution is 5.4 sigma and the lightcurve is similar to the lightcurves found earlier between 0.7 and 700 keV: a single broad asymmetric pulse reaching its maximum 0.38 +/- 0.03 in phase after the radio peak, compared to the offset of 0.30 found in the CGRO BATSE soft gamma-ray data, and 0.27 +/- 0.01 for RXTE (2-16 keV), compatible with ASCA (0.7-2.2 keV). Spectral analysis based on the excess counts in the broad pulse of the lightcurve shows that extrapolation of the OSSE power-law spectral fit with index -1.68 describes our data well up to 10 MeV. Above 10 MeV the spectrum breaks abruptly. The precise location of the break/bend between 10 and 30 MeV depends on the interpretation of the structure in the lightcurve measured by COMPTEL and EGRET above 10 MeV.

We report the COMPTEL detection of pulsed $\gamma$-emission from PSR B1509-58 up to at least 10 MeV using data collected over more than 6 years. The 0.75-10 MeV lightcurve is broad and reaches its maximum near radio-phase 0.38, slightly beyond the maximum found at hard X-rays/ soft $\gamma$-rays. In the 10-30 MeV energy range a strong source is present in the skymap positionally consistent with the pulsar, but we do not detect significant pulsed emission. However, the lightcurve is consistent with the pulse shape changing from a single broad pulse into a double-peak morphology. Our results significantly constrain pulsar modelling.

We present improved timing parameters for 13 millisecond pulsars (MSPs) including 9 new proper motion measurements. These new proper motions bring to 23 the number of MSPs with measured transverse velocities. In light of these new results we present and compare the kinematic properties of MSPs with those of ordinary pulsars. The mean transverse velocity of MSPs was found to be 85+/-13 km/s; a value consistent with most models for the origin and evolution of MSPs and approximately a factor of four lower than that of ordinary pulsars. We also find that, in contrast to young ordinary pulsars, the vast majority of which are moving away from the Galactic plane, almost half of the MSPs are moving towards the plane. This near isotropy would be expected of a population that has reached dynamic equilibrium. Accurate measurements of MSP velocities have allowed us to correct their measured spin-down rates for Doppler acceleration effects, and thereby derive their intrinsic magnetic field strengths and characteristic ages. We find that close to half of our sample of MSPs have a characteristic age comparable to or greater than the age of the Galaxy.

We compare the spectral properties of the millisecond and slow pulsars detected in the Parkes 70 cm survey. The mean spectral index for the millisecond pulsars (MSPs) is -1.9 +/- 0.1 whereas the mean spectral index for the slow pulsars is a surprisingly steep -1.72 +/- 0.04. A Kolmogorov-Smirnov test indicates that there is only a 72% probability that the two distributions differ. As a class, MSPs are therefore only fractionally steeper-spectrum objects than slow pulsars, as recent literature would suggest. We then model the expected distribution of millisecond pulsars in the Galaxy and find that high-frequency surveys, with sensitivities similar to the current Parkes multibeam survey, are likely to detect MSPs in large numbers. The observed distribution of MSPs will be much less isotropic than that resulting from low-frequency surveys, with 50% of detectable MSPs residing within 11 degrees of the Galactic plane in an all-sky survey.

We have carried out high-precision timing observations of the eclipsing binary PSR J2051-0827 in the 3.3 years since its discovery. These data indicate that the orbital period is decreasing at a rate of dPb/dt = (-11+-1)X10^-12. If secular, this orbital period derivative implies a decay time for the orbit of only 25 Myr which is much shorter than the expected timescale for ablation of the companion. We have also measured the proper motion of the pulsar to be 5+-3 mas/yr. Assuming the pulsar is at the dispersion-measure distance this implies a very slow transverse velocity vt=(30+-20) km/s. This combination of low velocity and short orbital period argue against formation of the system in the standard manner and we discuss the implications for its evolutionary history.

The pulsar PSR B1259-63 is in a highly eccentric 3.4-yr orbit with the Be star SS 2883. Timing observations of this pulsar, made over a 7-yr period using the Parkes 64-m radio telescope, cover two periastron passages, in 1990 August and 1994 January. The timing data cannot be fitted by the normal pulsar and Keplerian binary parameters. A timing solution including a (non-precessing) Keplerian orbit and timing noise (represented as a polynomial of fifth order in time) provide a satisfactory fit to the data. However, because the Be star probably has a significant quadrupole moment, we prefer to interpret the data by a combination of timing noise, dominated by a cubic phase term, and $\dot\omega$ and $\dot x$ terms. We show that the $\dot\omega$ and $\dot x$ are likely to be a result of a precessing orbit caused by the quadrupole moment of the tilted companion star. We further rule out a number of possible physical effects which could contribute to the timing data of PSR B1259-63 on a measurable level.

We report the discovery and follow-up timing observations of the 63 ms radio pulsar, PSR J1105-6107. The pulsar is young, having a characteristic age of only 63 kyr and, from its dispersion measure, is estimated to be at a distance of ~7 kpc from the Sun. We consider its possible association with the nearby supernova remnant G290.1-0.8 (MSH 11-61A); an association requires that the pulsar's proper motion be ~22 mas/yr (corresponding to ~650 km/s for a distance of 7 kpc) directed away from the remnant center, assuming that the characteristic age is the true age. The pulsar's spin-down luminosity, 2.5 x 10^{36} erg/s, is in the top 1% of all known pulsar spin-down luminosities. Given its estimated distance, PSR J1105-6107 is therefore likely to be observable at high energies. Indeed, it is coincident with the known CGRO/EGRET source 2EG J1103-6106; we consider the possible association and conclude that it is likely.

We present timing observations of four millisecond pulsars, using data obtained over three years at the ATNF Parkes and NRAL Jodrell Bank radio telescopes. Astrometric, spin, and binary parameters are updated, and substantially improved for three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured the time variation of the projected semi-major axis of the PSR J0437-4715 orbit due to its proper motion, and use it to constrain the inclination of the orbit and the mass of the companion. Some evidence is found for changes in the dispersion measures of PSRs J1045-4509 and J1643-1224. Limits are placed on the existence of planetary mass companions, ruling out companions with masses and orbits similar to the terrestrial planets of the solar system for eight pulsars.

We demonstrate how measuring orbital period derivatives can lead to more accurate distance estimates and transverse velocities for some nearby binary pulsars. In many cases this method will estimate distances more accurately than is possible by annual parallax, as the relative error decreases as t^-5/2. Unfortunately, distance uncertainties limit the degree to which nearby relativistic binary pulsars can be used for testing the general relativistic prediction of orbital period decay to a few percent. Nevertheless, the measured orbital period derivative of PSR B1534+12 agrees within the observational uncertainties with that predicted by general relativity if the proper-motion contribution is accounted for.

We report on optical observations of the low-mass binary millisecond pulsar systems J0034-0534 and J2145-0750. A faint (I=23.5) object was found to be coincident with the timing position of PSR J2145-0750. While a galaxy or distant main-sequence star cannot be ruled out, its magnitude is consistent with an ancient white dwarf, as expected from evolutionary models. For PSR J0034-0534 no objects were detected to a limiting magnitude of R=25.0, suggesting that the white dwarf in this system is cold. Using white dwarf cooling models, the limit on the magnitude of the PSR J0034-0534 companion suggests that at birth the pulsar in this system may have rotated with a period as short as 0.6 ms. These observations provide further evidence that the magnetic fields of millisecond pulsars do not decay on time scales shorter than 1 Gyr.

Binary radio pulsars are superb tools for mapping binary orbits, because of the precision of the pulse timing method (Taylor and Weisberg 1989). To date, all orbital parameters for binary pulsars have been derived from observations of the pulsar alone. We present the first observations of the radial velocity variations due to the binary motion of a companion to a radio pulsar. Our results demonstrate that the companion to the Small Magellanic Cloud pulsar PSR J0045-7319 is the B1V star identified by Kaspi et al. (1994). The mass ratio of the system is 6.3 +/- 1.2, which, for a neutron star mass of 1.4 Mo, implies a mass of 8.8 +/- 1.8 Mo for the companion, consistent with the mass expected for a B1V star. The inclination angle for the binary system is therefore 44 +/- 5 degrees, and the projected rotational velocity of the companion is 113 +/- 10 km/s. The heliocentric radial velocity of the binary system is consistent with that of other stars and gas in the same region of the Small Magellanic Cloud.