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.

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.

The Galactic magnetic field plays an important role in the evolution of the Galaxy, but its small-scale behaviour is still poorly known. It is also unknown whether it permeates the halo of the Galaxy or not. By using observations of pulsars in the halo globular cluster 47 Tucanae, we probed the Galactic magnetic field at arcsecond scales for the first time and discovered an unexpected large gradient in the component of the magnetic field parallel to the line of sight. This gradient is aligned with a direction perpendicular to the Galactic disk and could be explained by magnetic fields amplified to some 60 {\mu}G within the globular cluster. This scenario supports the existence of a magnetized outflow that extends from the Galactic disk to the halo and interacts with the studied globular cluster.

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.

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.

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

The 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$

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.

Here, we report on the detection and verification of Fast Radio Burst FRB 180301, which occurred on UTC 2018 March 1 during the Breakthrough Listen observations with the Parkes telescope. Full-polarization voltage data of the detection were captured--a first for non-repeating FRBs--allowing for coherent de-dispersion and additional verification tests. The coherently de-dispersed dynamic spectrum of FRB 180301 shows complex, polarized frequency structure over a small fractional bandwidth. As FRB 180301 was detected close to the geosynchronous satellite band during a time of known 1-2 GHz satellite transmissions, we consider whether the burst was due to radio interference emitted or reflected from an orbiting object. Based on the preponderance of our verification tests, we find that FRB 180301 is likely of astrophysical origin, but caution that anthropogenic sources cannot conclusively be ruled out.

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.

Studying the gravitational-wave sky with pulsar timing arrays (PTAs) is a key science goal for the Square Kilometre Array (SKA) and its pathfinder telescopes. With current PTAs reaching sub-microsecond timing precision, making accurate measurements of interstellar propagation effects and mitigating them effectively has become increasingly important to realise PTA goals. As these effects are much stronger at longer wavelengths, low-frequency observations are most appealing for characterizing the interstellar medium (ISM) along the sight lines toward PTA pulsars. The Murchison Widefield Array (MWA) and the Engineering Development Array (EDA), which utilizes MWA technologies, present promising opportunities for undertaking such studies, particularly for PTA pulsars located in the southern sky. Such pulsars are also the prime targets for PTA efforts planned with the South African MeerKAT, and eventually with the SKA. In this paper we report on observations of two bright southern millisecond pulsars PSRs J0437-4715 and J2145-0750 made with these facilities; MWA observations sampling multiple frequencies across the 80-250 MHz frequency range, while the EDA providing direct-sampled baseband data to yield a large instantaneous usable bandwidth of $\sim$200 MHz. Using these observations, we investigate various aspects relating to pulsar emission and ISM properties, such as spectral evolution of the mean pulse shape, scintillation as a function of frequency, chromaticity in interstellar dispersion, and flux density spectra at low frequencies. Systematic and regular monitoring observations will help ascertain the role of low-frequency measurements in PTA experiments, while simultaneously providing a detailed characterization of the ISM toward the pulsars, which will be useful in devising optimal observing strategies for future PTA experiments.

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.

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.

We present the measurements of Faraday rotation for 477 pulsars observed by the Parkes 64-m radio telescope and the Green Bank 100-m radio telescope. Using these results along with previous measurements for pulsars and extra-galactic sources, we analyse the structure of the large-scale magnetic field in the Galactic disk. Comparison of rotation measures of pulsars in the disk at different distances as well as with rotation measures of background radio sources beyond the disk reveals large-scale reversals of the field directions between spiral arms and interarm regions. We develop a model for the disk magnetic field, which can reproduce not only these reversals but also the distribution of observed rotation measures of background sources.

The Square Kilometre Array will be an amazing instrument for pulsar astronomy. While the full SKA will be sensitive enough to detect all pulsars in the Galaxy visible from Earth, already with SKA1, pulsar searches will discover enough pulsars to increase the currently known population by a factor of four, no doubt including a range of amazing unknown sources. Real time processing is needed to deal with the 60 PB of pulsar search data collected per day, using a signal processing pipeline required to perform more than 10 POps. Here we present the suggested design of the pulsar search engine for the SKA and discuss challenges and solutions to the pulsar search venture.

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 five-component profile of the 2.7-ms pulsar J0337+1715 appears to exhibit the best example to date of a core/double-cone emission-beam structure in a millisecond pulsar (MSP). Moreover, three other MSPs, the Binary Pulsar B1913+16, B1953+29 and J1022+1001, seem to exhibit core/single-cone profiles. These configurations are remarkable and important because it has not been clear whether MSPs and slow pulsars exhibit similar emission-beam configurations, given that they have considerably smaller magnetospheric sizes and magnetic field strengths. MSPs thus provide an extreme context for studying pulsar radio emission. Particle currents along the magnetic polar flux tube connect processes just above the polar cap through the radio-emission region to the light-cylinder and the external environment. In slow pulsars radio-emission heights are typically about 500 km around where the magnetic field is nearly dipolar, and estimates of the physical conditions there point to radiation below the plasma frequency and emission from charged solitons by the curvature process. We are able to estimate emission heights for the four MSPs and carry out a similar estimation of physical conditions in their much lower emission regions. We find strong evidence that MSPs also radiate by curvature emission from charged solitons.

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.

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

A four-dimensional statistical description of electromagnetic radiation is developed and applied to the analysis of radio pulsar polarization. The new formalism provides an elementary statistical explanation of the modal broadening phenomenon in single pulse observations. It is also used to argue that the degree of polarization of giant pulses has been poorly defined in past studies. Single and giant pulse polarimetry typically involves sources with large flux densities and observations with high time resolution, factors that necessitate consideration of source-intrinsic noise and small-number statistics. Self noise is shown to fully explain the excess polarization dispersion previously noted in single pulse observations of bright pulsars, obviating the need for additional randomly polarized radiation. Rather, these observations are more simply interpreted as an incoherent sum of covariant, orthogonal, partially polarized modes. Based on this premise, the four-dimensional covariance matrix of the Stokes parameters may be used to derive mode-separated pulse profiles without any assumptions about the intrinsic degrees of mode polarization. Finally, utilizing the small-number statistics of the Stokes parameters, it is established that the degree of polarization of an unresolved pulse is fundamentally undefined; therefore, previous claims of highly polarized giant pulses are unsubstantiated. Unpublished supplementary material is appended after the bibliography.

A new method of matrix template matching is presented in the context of pulsar timing analysis. Pulse arrival times are typically measured using only the observed total intensity light curve. The new technique exploits the additional timing information available in the polarization of the pulsar signal by modeling the transformation between two polarized light curves in the Fourier domain. For a number of millisecond pulsars, arrival time estimates derived from polarimetric data are predicted to exhibit greater precision and accuracy than those derived from the total intensity alone. Furthermore, the transformation matrix produced during template matching may be used to calibrate observations of other point sources. Unpublished supplementary material is appended after the bibliography.

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.

A statistical framework is presented for the study of the orthogonally polarized modes of radio pulsar emission via the covariances between the Stokes parameters. To accommodate the typically heavy-tailed distributions of single-pulse radio flux density, the fourth-order joint cumulants of the electric field are used to describe the superposition of modes with arbitrary probability distributions. The framework is used to consider the distinction between superposed and disjoint modes with particular attention to the effects of integration over finite samples. If the interval over which the polarization state is estimated is longer than the timescale for switching between two or more disjoint modes of emission, then the modes are unresolved by the instrument. The resulting composite sample mean exhibits properties that have been attributed to mode superposition, such as depolarization. Because the distinction between disjoint modes and a composite sample of unresolved disjoint modes depends on the temporal resolution of the observing instrumentation, the arguments in favour of superposed modes of pulsar emission are revisited and observational evidence for disjoint modes is described. In principle, the four-dimensional covariance matrix that describes the distribution of sample mean Stokes parameters can be used to distinguish between disjoint modes, superposed modes, and a composite sample of unresolved disjoint modes. More comprehensive and conclusive interpretation of the covariance matrix requires more detailed consideration of various relevant phenomena, including temporally correlated subpulse modulation (e.g. jitter), statistical dependence between modes (e.g. covariant intensities and partial coherence), and multipath propagation effects (e.g. scintillation and scattering). Unpublished supplementary material is appended after the bibliography.

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.

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.

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

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

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.

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

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

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.

Radio telescopes are used to accurately measure the time of arrival (ToA) of radio pulses in pulsar timing experiments that target mostly millisecond pulsars (MSPs) due to their high rotational stability. This allows for detailed study of MSPs and forms the basis of experiments to detect gravitational waves. Apart from intrinsic and propagation effects, such as pulse-to-pulse jitter and dispersion variations in the interstellar medium, timing precision is limited in part by the following: polarization purity of the telescope's orthogonally polarized receptors, the signal-to-noise ratio (S/N) of the pulsar profile, and the polarization fidelity of the system. Using simulations, we present how fundamental limitations in recovering the true polarization reduce the precision of ToA measurements. Any real system will respond differently to each source observed depending on the unique pulsar polarization profile. Using the profiles of known MSPs we quantify the limits of observing system specifications that yield satisfactory ToA measurements, and we place a practical design limit beyond which improvement of the system results in diminishing returns. Our aim is to justify limits for the front-end polarization characteristics of next generation radio telescopes, leading to the Square Kilometre Array (SKA).

We report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5$\sigma$ and 6.5$\sigma$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>\alpha> -4.9$ ($S_{\rm \nu} \propto \nu^\alpha$). We present a Monte Carlo analysis supporting the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than $100\%$ correlation being due to the relative sensitivities of the two instruments and the width of the spectral index distribution. Our observations are consistent with the hypothesis that the spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 MHz.

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.

The Square Kilometre Array (SKA) will use pulsars to enable precise measurements of strong gravity effects in pulsar systems, which yield tests of gravitational theories that cannot be carried out anywhere else. The Galactic census of pulsars will discover dozens of relativistic pulsar systems, possibly including pulsar -- black hole binaries which can be used to test the "cosmic censorship conjecture" and the "no-hair theorem". Also, the SKA's remarkable sensitivity will vastly improve the timing precision of millisecond pulsars, allowing probes of potential deviations from general relativity (GR). Aspects of gravitation to be explored include tests of strong equivalence principles, gravitational dipole radiation, extra field components of gravitation, gravitomagnetism, and spacetime symmetries.

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 report on the detection of the millisecond pulsar PSR J0437-4715 with the Murchison Widefield Array (MWA) at a frequency of 192 MHz. Our observations show rapid modulations of pulse intensity in time and frequency that arise from diffractive scintillation effects in the interstellar medium (ISM), as well as prominent drifts of intensity maxima in the time-frequency plane that arise from refractive effects. Our analysis suggests that the scattering screen is located at a distance of $\sim$80-120 pc from the Sun, in disagreement with a recent claim that the screen is closer ($\sim$10 pc). Comparisons with higher frequency data from Parkes reveals a dramatic evolution of the pulse profile with frequency, with the outer conal emission becoming comparable in strength to that from the core and inner conal regions. As well as demonstrating high time resolution science capabilities currently possible with the MWA, our observations underscore the potential to conduct low-frequency investigations of timing-array millisecond pulsars, which may lead to increased sensitivity for the detection of nanoHertz gravitational waves via the accurate characterisation of ISM effects.

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.

Cyclic spectroscopy is a signal processing technique that was originally developed for engineering applications and has recently been introduced into the field of pulsar astronomy. It is a powerful technique with many attractive features, not least of which is the explicit rendering of information about the relative phases in any filtering imposed on the signal, thus making holography a more straightforward proposition. Here we present methods for determining optimum estimates of both the filter itself and the statistics of the unfiltered signal, starting from a measured cyclic spectrum. In the context of radio pulsars these quantities tell us the impulse response of the interstellar medium and the intrinsic pulse profile. We demonstrate our techniques by application to 428 MHz Arecibo data on the millisecond pulsar B1937+21, obtaining the pulse profile free from the effects of interstellar scattering. As expected, the intrinsic profile exhibits main- and inter-pulse components that are narrower than they appear in the scattered profile; it also manifests some weak, but sharp features that are revealed for the first time at low frequency. We determine the structure of the received electric-field envelope as a function of delay and Doppler-shift. Our delay-Doppler image has a high dynamic-range and displays some pronounced, low-level power concentrations at large delays. These concentrations imply strong clumpiness in the ionized interstellar medium, on AU size-scales, which must adversely affect the timing of B1937+21.

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.

We analyse dispersion measure (DM) variations in six years of radio observations of more than 160 young pulsars, all gamma ray candidates for the Fermi gamma ray telescope mostly located close to the Galactic plane. DMs were fit across 256 MHz of bandwidth for observations centred at 1.4 GHz and across three frequencies -- 0.7 GHz, 1.4 GHz, and 3.1 GHz -- where multifrequency observations were available. Changes in dispersion measure, dDM/dt, were calculated using a weighted linear fit across all epochs of available DMs. DM variations were detected at a 3\sigma level in 11 pulsars, four of which were above 5\sigma: PSRs J0835-4510, J0908-4913, J1824-1945, and J1833-0827. We find that after 28 years of gradual decline, the DM of PSR J0835-4510 is now increasing. The magnitude of variations in three of the four (PSRs J0835-4510, J0908-4913, and J1833-0827) are above what models would predict for an ISM dominated by Kolmogorov turbulence. We attribute this excess as likely due to the pulsar's local environment - the supernova remnants near PSRs J0835-4510 and J1833-0827 and the pulsar wind nebula around PSR J0908-4913. Upper limits were determined for all pulsars without detectable values of dDM/dt, most limits were found to lie above the levels of variations predicted by ISM theory. We find our results to be consistent with scattering estimates from the NE2001 model along these lines of sight.

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.

A new method of polarimetric calibration is presented in which the instrumental response is derived from regular observations of PSR J0437-4715 based on the assumption that the mean polarized emission from this millisecond pulsar remains constant over time. The technique is applicable to any experiment in which high-fidelity polarimetry is required over long time scales; it is demonstrated by calibrating 7.2 years of high-precision timing observations of PSR J1022+1001 made at the Parkes Observatory. Application of the new technique followed by arrival time estimation using matrix template matching yields post-fit residuals with an uncertainty-weighted standard deviation of 880 ns, two times smaller than that of arrival time residuals obtained via conventional methods of calibration and arrival time estimation. The precision achieved by this experiment yields the first significant measurements of the secular variation of the projected semi-major axis, the precession of periastron, and the Shapiro delay; it also places PSR J1022+1001 among the ten best pulsars regularly observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is shown that the timing accuracy of a large fraction of the pulsars in the PPTA is currently limited by the systematic timing error due to instrumental polarization artifacts. More importantly, long-term variations of systematic error are correlated between different pulsars, which adversely affects the primary objectives of any pulsar timing array experiment. These limitations may be overcome by adopting the techniques presented in this work, which relax the demand for instrumental polarization purity and thereby have the potential to reduce the development cost of next-generation telescopes such as the Square Kilometre Array.

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.

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.

PSRCHIVE is an open-source, object-oriented, scientific data analysis software library and application suite for pulsar astronomy. It implements an extensive range of general-purpose algorithms for use in data calibration and integration, statistical analysis and modeling, and visualisation. These are utilised by a variety of applications specialised for tasks such as pulsar timing, polarimetry, radio frequency interference mitigation, and pulse variability studies. This paper presents a general overview of PSRCHIVE functionality with some focus on the integrated interfaces developed for the core applications.

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

Timing pulses of pulsars has proved to be a most powerful technique useful to a host of research areas in astronomy and physics. Importantly, the precision of this timing is not only affected by radiometer noise, but also by intrinsic pulse shape changes, interstellar medium (ISM) evolution, instrumental distortions, etc. In this paper we review the known causes of pulse shape variations and assess their effect on the precision and accuracy of a single measurement of pulse arrival time with current instrumentation. Throughout this analysis we use the brightest and most precisely timed millisecond pulsar (MSP), PSR J0437-4715, as a case study, and develop a set of diagnostic tools to evaluate profile stability in timing observations. We conclude that most causes of distortion can be either corrected by state-of-the-art techniques or taken into account in the estimation of time-of-arrival (TOA) uncertainties. The advent of a new generation of radio telescopes (e.g. the Square Kilometre Array, SKA), and their increase in collecting area has sparked speculation about the timing precision achievable through increases in gain. Based on our analysis of current data, we predict that for normal-brightness MSPs a TOA precision of between 80 and 230 ns can be achieved at 1.4 GHz with 10-minute integrations by the SKA. The actual rms timing residuals for each pulsar will be approximately at the same level only if all the other influences on timing precision (e.g. ISM, spin noise) are either corrected, modelled, or negligible.

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

The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.

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.

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 Crab pulsar is well-known for its anomalous giant radio pulse emission. Past studies have concentrated only on the very bright pulses or were insensitive to the faint end of the giant pulse luminosity distribution. With our new instrumentation offering a large bandwidth and high time resolution combined with the narrow radio beam of the Westerbork Synthesis Radio Telescope (WSRT), we seek to probe the weak giant pulse emission regime. The WSRT was used in a phased array mode, resolving a large fraction of the Crab nebula. The resulting pulsar signal was recorded using the PuMa II pulsar backend and then coherently dedispersed and searched for giant pulse emission. After careful flux calibration, the data were analysed to study the giant pulse properties. The analysis includes the distributions of the measured pulse widths, intensities, energies, and scattering times. The weak giant pulses are shown to form a separate part of the intensity distribution. The large number of giant pulses detected were used to analyse scattering and scintillation in giant pulses. We report for the first time the detection of giant pulse emission at both the main- and interpulse phases within a single rotation period. The rate of detection is consistent with the appearance of pulses at either pulse phase as being independent. These pulse pairs were used to examine the scintillation timescales within a single pulse period.

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.

Observations of the binary pulsar PSR J1141-6545 using the Parkes radio telescope over 9.3 years show clear time-variations in pulse width, shape and polarization. We interpret these variations in terms of relativistic precession of the pulsar spin axis about the total angular momentum vector of the system. Over the nine years, the pulse width at the 50% level has changed by more than a factor of three. Large variations have also been observed in the 1400-MHz mean flux density. The pulse polarization has been monitored since 2004 April using digital filterbank systems and also shows large and systematic variations in both linear and circular polarization. Position angle variations, both across the pulse profile and over the data span, are complex, with major differences between the central and outer parts of the pulse profile. Modelling of the observed position angle variations by relativistic precession of the pulsar spin axis shows that the spin-orbit misalignment angle is about 110 deg and that the precessional phase has passed through 180 deg during the course of our observations. At the start of our observations, the line-of-sight impact parameter was about 4 deg in magnitude and it reached a minimum very close to 0 deg around early 2007, consistent with the observed pulse width variations. We have therefore mapped approximately one half of the emission beam, showing that it is very asymmetric with respect to the magnetic axis. The derived precessional parameters imply that the pre-supernova star had a mass of about 2 Msun and that the supernova recoil kick velocity was relatively small. With the reversal in the rate of change of the impact parameter, we predict that over the next decade we will see a reversed "replay" of the variations observed in the past decade.

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.

This paper defines the mathematical convention adopted to describe an electromagnetic wave and its polarisation state, as implemented in the PSRCHIVE software and represented in the PSRFITS definition. Contrast is made between the convention that has been widely accepted by pulsar astronomers and the IAU/IEEE definitions of the Stokes parameters. The former is adopted as the PSR/IEEE convention, and a set of useful parameters are presented for describing the differences between the PSR/IEEE standard and the conventions (either implicit or explicit) that form part of the design of observatory instrumentation. To aid in the empirical determination of instrumental convention parameters, well-calibrated average polarisation profiles of PSR J0304+1932 and PSR J0742-2822 are presented at radio wavelengths of approximately 10, 20, and 40 cm.

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.

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from two of them. PSR J1741-2054, with period P=413ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM=4.7pc/cc, implying a distance of ~0.4kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1GeV corresponds to 25% of the spin-down luminosity of dE/dt=9.4e33erg/s. The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks with the first peak lagging the radio pulse by delta=0.29P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P=143ms, and its DM=115pc/cc suggests a distance of ~3.6kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by delta=0.15P the first of two narrow gamma-ray peaks that are separated by Delta=0.50P. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA source TeV J2032+4130.

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

The PULSE@Parkes project has been designed to monitor the rotation of radio pulsars over time spans of days to years. The observations are obtained using the Parkes 64-m and 12-m radio telescopes by Australian and international high school students. These students learn the basis of radio astronomy and undertake small projects with their observations. The data are fully calibrated and obtained with the state-of-the-art pulsar hardware available at Parkes. The final data sets are archived and are currently being used to carry out studies of 1) pulsar glitches, 2) timing noise, 3) pulse profile stability over long time scales and 4) the extreme nulling phenomenon. The data are also included in other projects such as gamma-ray observatory support and for the Parkes Pulsar Timing Array project. In this paper we describe the current status of the project and present the first scientific results from the Parkes 12-m radio telescope. We emphasise that this project offers a straightforward means to enthuse high school students and the general public about radio astronomy while obtaining scientifically valuable data sets.

We present the first polarization measurements of the radio emission from RRAT J1819$-$1458. Our observations, conducted in parallel to regular timing sessions, have yielded a small number of bright and polarized pulses. The polarization characteristics and integrated profile resemble those of normal pulsars with average spin-down energy (Edot): moderate to low linear polarization in the integrated profile despite relatively high polarization in the individual pulses. On average, a small degree of circular polarization is also observed. The polarization position angle executes a remarkably smooth, steep S-shaped curve, interrupted by two orthogonal jumps. Based on the shape of the PA swing, we place some constraints on the emission geometry. We compare these polarization properties to those of other radio emitting neutron star populations, including young pulsars, pulsars with a high surface magnetic field and radio emitting magnetars. From the polarization measurements, the Faraday rotation measure of this RRAT is derived.

We report on polarimetric observations of 100 pulsars centered on 774 MHz, made using the Green Bank Telescope, presenting their polarization profiles and polarized flux densities and comparing them with previous observations when possible. For 67 pulsars, these are the first such measurements made. Polarization profiles of 8 millisecond pulsars in our sample show wide profiles and flat position-angle curves. Strong linear polarization, sometimes approaching 100% of the total intensity, has been detected in all or a part of the average pulse profiles of some pulsars. In general, circular polarization is very weak, although it is observed to be extremely strong in the leading component of PSR J1920+2650. Sense reversal of circular polarization as a function of pulse phase has been detected from both core and other components of more than 20 pulsars. Any relationship between the spin-down luminosity and the percentage of linear polarization is not evident in our data at this frequency.

Since 2004 we have been carrying out a pulsar survey of the Cygnus region with the Westerbork Synthesis Radio Telescope (WSRT) at a frequency of 328 MHz. The survey pioneered a novel interferometric observing mode, termed 8gr8 (eight-grate), whereby multiple simultaneous digital beams provide high sensitivity over a large field of view. Since the Cygnus region is known to contain OB associations, it is likely that pulsars are formed here. Simulations have shown that this survey could detect 70 pulsars, which would increase our understanding of the radio pulsar population in this region. We also aim to expand the known population of intermittent and rotating radio transient (RRAT)-like pulsars. In this paper we describe our methods of observation, processing and data analysis, and we present the first results. Our observing method exploits the way a regularly spaced, linear array of telescopes yields a corresponding regularly spaced series of so-called ``grating'' beams on the sky. By simultaneously forming a modest number (eight) of offset digital beams, we can utilize the entire field of view of each WSRT dish, but retain the coherently summed sensitivity of the entire array. For the processing we performed a large number of trial combinations of period and dispersion measure (DM) using a computer cluster. In the first processing cycle of the WSRT 8gr8 Cygnus Survey, we have discovered three radio pulsars, with spin periods of 1.657, 1.099 and 0.445 seconds. These pulsars have been observed on a regular basis since their discovery, both in a special follow-up programme as well as in the regular timing programme. The timing solutions are presented in this paper. We also discuss this survey method in the context of the SKA and its pathfinders.

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.

The Pulsar Machine II (PuMa II) is the new flexible pulsar processing backend system at the Westerbork Synthesis Radio Telescope (WSRT), specifically designed to take advantage of the upgraded WSRT. The instrument is based on a computer cluster running the Linux operating system, with minimal custom hardware. A maximum of 160 MHz analogue bandwidth sampled as 8X20 MHz subbands with 8-bit resolution can be recorded on disks attached to separate computer nodes. Processing of the data is done in the additional 32-nodes allowing near real time coherent dedispersion for most pulsars observed at the WSRT. This has doubled the bandwidth for pulsar observations in general, and has enabled the use of coherent dedispersion over a bandwidth eight times larger than was previously possible at the WSRT. PuMa II is one of the widest bandwidth coherent dedispersion machines currently in use and has a maximum time resolution of 50ns. The system is now routinely used for high precision pulsar timing studies, polarization studies, single pulse work and a variety of other observational work.

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 dynamic spectrum of a radio pulsar is an in-line digital hologram of the ionised interstellar medium. It has previously been demonstrated that such holograms permit image reconstruction, in the sense that one can determine an approximation to the complex electric field values as a function of Doppler-shift and delay, but to date the quality of the reconstructions has been poor. Here we report a substantial improvement in the method which we have achieved by simultaneous optimisation of the thousands of coefficients that describe the electric field. For our test spectrum of PSR B0834+06 we find that the model provides an accurate representation of the data over the full 63 dB dynamic range of the observations: residual differences between model and data are noise-like. The advent of interstellar holography enables detailed quantitative investigation of the interstellar radio-wave propagation paths for a given pulsar at each epoch of observation; we illustrate this using our test data which show the scattering material to be structured and highly anisotropic. The temporal response of the medium exhibits a scattering tail out to beyond 100 microsec and a pulse arrival time measurement at this frequency and this epoch of observation would be affected by a mean delay of 15 microsec due to multipath propagation in the interstellar medium.

We have used the Parkes radio telescope to study the polarized emission from the anomalous X-ray pulsar XTE J1810-197 at frequencies of 1.4, 3.2, and 8.4 GHz. We find that the pulsed emission is nearly 100% linearly polarized. The position angle of linear polarization varies gently across the observed pulse profiles, varying little with observing frequency or time, even as the pulse profiles have changed dramatically over a period of 7 months. In the context of the standard pulsar "rotating vector model," there are two possible interpretations of the observed position angle swing coupled with the wide profile. In the first, the magnetic and rotation axes are substantially misaligned and the emission originates high in the magnetosphere, as seen for other young radio pulsars, and the beaming fraction is large. In the second interpretation, the magnetic and rotation axes are nearly aligned and the line of sight remains in the emission zone over almost the entire pulse phase. We deprecate this possibility because of the observed large modulation of thermal X-ray flux. We have also measured the Faraday rotation caused by the Galactic magnetic field, RM = +77 rad/m^2, implying an average magnetic field component along the line of sight of 0.5 microG.

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.

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.

The large-scale magnetic field of our Galaxy can be probed in three dimensions using Faraday rotation of pulsar signals. We report on the determination of 223 rotation measures from polarization observations of relatively distant southern pulsars made using the Parkes radio telescope. Combined with previously published observations these data give clear evidence for large-scale counterclockwise fields (viewed from the north Galactic pole) in the spiral arms interior to the Sun and weaker evidence for a counterclockwise field in the Perseus arm. However, in interarm regions, including the Solar neighbourhood, we present evidence that suggests that large-scale fields are clockwise. We propose that the large-scale Galactic magnetic field has a bisymmetric structure with reversals on the boundaries of the spiral arms. Streaming motions associated with spiral density waves can directly generate such a structure from an initial inwardly directed radial field. Large-scale fields increase toward the Galactic Center, with a mean value of about 2~$\mu$G in the Solar neighbourhood and 4~$\mu$G at a Galactocentric radius of 3 kpc.

A new set of software applications and libraries for use in the archival and analysis of pulsar astronomical data is introduced. Known collectively as the PSRCHIVE scheme, the code was developed in parallel with a new data storage format called PSRFITS, which is based on the Flexible Image Transport System (FITS). Both of these projects utilise a modular, object-oriented design philosophy. PSRCHIVE is an open source development environment that incorporates an extensive range of C++ object classes and pre-built command line and graphical utilities. These deal transparently and simultaneously with multiple data storage formats, thereby enhancing data portability and facilitating the adoption of the PSRFITS file format. Here, data are stored in a series of modular header-data units that provide flexibility and scope for future expansion. As it is based on FITS, various standard libraries and applications may be used for data input/output and visualisation. Both PSRCHIVE and PSRFITS are made publicly available to the academic community in the hope that this will promote their widespread use and acceptance.

A new method is presented for modeling the transformation between two polarimetric pulse profiles in the Fourier domain. In practice, one is a well-determined standard with high signal-to-noise ratio and the other is an observation that is to be fitted to the standard. From this fit, both the longitudinal shift and the polarimetric transformation between the two profiles are determined. Arrival time estimates derived from the best-fit longitudinal shift are shown to exhibit greater precision than those derived from the total intensity profile alone. In addition, the polarimetric transformation obtained through this method may be used to completely calibrate the instrumental response in observations of other sources.

A mathematical framework is presented for use in the experimental determination of the polarimetric response of observatory instrumentation. Elementary principles of linear algebra are applied to model the full matrix description of the polarization measurement equation by least-squares estimation of non-linear, scalar parameters. The formalism is applied to calibrate the center element of the Parkes Multibeam receiver using observations of the millisecond pulsar, PSR J0437-4715, and the radio galaxy, 3C 218 (Hydra A).

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 present high-resolution multi-frequency single-pulse observations of the Vela pulsar, PSR B0833-45, aimed at studying micro-structure, phase-resolved intensity fluctuations and energy distributions at 1.41 and 2.30 GHz. We show that the micro-pulse width in pulsars has a period dependence. Like individual pulses, Vela's micro-pulses are highly elliptically polarized. There is a strong correlation between Stokes parameters V and I in the micro-structure. We show that the V/I distribution is Gaussian with a narrow width and that this width appears to be constant as a function of pulse phase. The phase-resolved intensity distributions of I are best fitted with log-normal statistics. Extra emission components, i.e.``bump'' and ``giant micro-pulses'', discovered by Johnston et al.(2001) are also present at 2.3 GHz. The bump component seems to be an extra component superposed on the main pulse profile but does not appear periodically. The giant micro-pulses are time-resolved and have significant jitter in their arrival times. Their flux density distribution is best fitted by a power-law, indicating a link between these features and ``classical'' giant pulses as observed for the Crab pulsar, (PSR B0531+21), PSR B1937+21 and PSR B1821-24. We find that Vela contains a mixture of emission properties representing both ``classical'' properties of radio pulsars (e.g. micro-structure, high degree of polarization, S-like position angle swing, orthogonal modes) and features which are most likely related to high-energy emission (e.g. extra profile components, giant micro-pulses). It hence represents an ideal test case to study the relationship between radio and high-energy emission in significant detail.

A new phase-coherent technique for the calibration of polarimetric data is presented. Similar to the one-dimensional form of convolution, data are multiplied by the response function in the frequency domain. Therefore, the system response may be corrected with arbitrarily high spectral resolution, effectively treating the problem of bandwidth depolarization. As well, the original temporal resolution of the data is retained. The method is therefore particularly useful in the study of radio pulsars, where high time resolution and polarization purity are essential requirements of high-precision timing. As a demonstration of the technique, it is applied to full-polarization baseband recordings of the nearby millisecond pulsar, PSR J0437-4715.

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