We have used the Nanshan 25-m Radio Telescope at Xinjiang Astronomical Observatory to obtain timing observations of 87 pulsars from 2002 July to 2014 March. Using the "Cholesky" timing analysis method we have determined positions and proper motions for 48 pulsars, 24 of which are improved positions compared to previously published values. We also present the first published proper motions for nine pulsars and improved proper motions for 21 pulsars using pulsar timing and position comparison method. The pulsar rotation parameters are derived and are more accurate than previously published values for 36 pulsars. Glitches are detected in three pulsars: PSRs J1722$-$3632, J1852$-$0635 and J1957+2831. For the first two, the glitches are large, with $\Delta\nu_g/\nu > 10^{-6}$, and they are the first detected glitches in these pulsars. PSR J1722$-$3632 is the second oldest pulsar with large glitch. For the middle-age pulsars ($\tau_c > 10^5$~yr), the calculated braking indices, $|n|$, are strongly correlated with $\tau_c$ and the numbers of positive and negative values of $n$ are almost equal. For young pulsars ($\tau_c < 10^5$~yr), there is no correlation between $|n|$ and $\tau_c$ and most have $n>0$.

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

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

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

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.

By analysing the data acquired from the Parkes 64-m radio telescope at 1369 MHz, we report on the phase-stationary non-drift amplitude modulation observed in PSR J1048-5832. The high-sensitivity observations revealed that the central and trailing components of the pulse profile of this pulsar switch between a strong mode and a weak mode periodically. However, the leading component remains unchanged. Polarization properties of the strong and weak modes are investigated. Considering the similarity to mode changing, we argue that the periodic amplitude modulation in PSR J1048$-$5832 is periodic mode changing. The fluctuation spectral analysis showed that the modulation period is very short (~2.1 s or 17 P1), where P1 is the rotation period of the pulsar. We find that this periodic amplitude modulation is hard to explain by existing models that account for the periodic phenomena in pulsars like subpulse drifting.

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

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

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

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

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

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

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

PSR B1828-11 is a young pulsar once thought to be undergoing free precession and recently found instead to be switching magnetospheric states in tandem with spin-down changes. Here we show the two extreme states of the mode-changing found for this pulsar and comment briefly on its interpretation.

PSR J1825$-$0935 (PSR B1822$-$09) switches between radio-quiet (Q-mode) and radio-bright (B-mode) modes. The Q-mode is known to have a periodic fluctuation that modulates both the interpulse and the main pulse with the same period. Earlier investigators argued that the periodic Q-mode modulation is associated with drifting subpulses. We report on single-pulse observations of PSR J1825$-$0935 that were made using the Parkes 64-m radio telescope with a central frequency of 1369 MHz. The high-sensitivity observations revealed that the periodic Q-mode modulation is in fact a periodic longitude-stationary intensity modulation occurring in the interpulse and the main pulse. The fluctuation spectral analysis showed that the modulation period is about $43 P_1$, where $P_1$ is the rotation period of the pulsar. Furthermore, we confirm that the modulation patterns in the interpulse and the main pulse are phase-locked. Specifically, the intensities of the interpulse and the immediately following main pulse are more highly correlated than for the main pulse and interpulse at any other lag. Polarization properties of the strong and weak Q-mode states are different, even for the trailing part of the main pulse which does not show the periodic intensity modulation.

The young pulsar PSR B1828-11 has long been known to show correlated shape and spin-down changes with timescales of roughly 500 and 250 days, perhaps associated with large-scale magnetospheric switching. Here we present multi-hour observations with the Parkes and Green Bank Telescopes at multiple phases across the roughly 500-day cycle and show that the pulsar undergoes mode-changing between two stable, extreme profile states. The fraction of time spent in each profile state naturally accounts for the observed overall "shape parameter" (defined to be 0 for wide profiles and 1 for narrow ones); this and the variable rate of the mode transitions are directly related to the spin-down changes. We observe that the mode transition rate could plausibly function as an additional parameter governing the chaotic behaviour in this object which was proposed earlier by Seymour and Lorimer. Free precession is not needed to account for the variations.

We present integrated pulse profiles at 5~GHz for 71 pulsars, including eight millisecond pulsars (MSPs), obtained using the Shanghai Tian Ma Radio Telescope (TMRT). Mean flux densities and pulse widths are measured. For 19 normal pulsars and one MSP, these are the first detections at 5~GHz and for a further 19, including five MPSs, the profiles have a better signal-to-noise ratio than previous observations. Mean flux density spectra between 400~MHz and 9~GHz are presented for 27 pulsars and correlations of power-law spectral index are found with characteristic age, radio pseudo-luminosity and spin-down luminosity. Mode changing was detected in five pulsars. The separation between the main pulse and interpulse is shown to be frequency independent for six pulsars but a frequency dependence of the relative intensity of the main pulse and interpulse is found. The frequency dependence of component separations is investigated for 20 pulsars and three groups are found: in seven cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping; in eleven cases the separation is nearly constant; in the remain two cases the separation between the outmost components increases with frequency. We obtain the correlations of pulse widths with pulsar period and estimate the core widths of 23 multi-component profiles and conal widths of 17 multi-component profiles at 5.0~GHz using Gaussian fitting and discuss the width-period relationship at 5~GHz compared with the results at at 1.0~GHz and 8.6~GHz.

Identifying as many gamma-ray pulsars as possible in the Fermi Large Area Telescope (LAT) data helps test pulsar emission models by comparing predicted and observed properties for a large, varied sample with as little selection bias as possible. It also improves extrapolations from the observed population to estimate the contribution of unresolved pulsars to the diffuse gamma-ray emission. We use a recently developed method to determine the probability that a given gamma-ray photon comes from a known position in the sky, convolving the photon's energy with the LAT's energy-dependent point-spread-function (PSF), without the need for an accurate spatial and spectral model of the gamma-ray sky around the pulsar. The method is simple and fast and, importantly, provides probabilities, or _weights_, for gamma rays from pulsars too faint for phase-integrated detection. We applied the method to over a thousand pulsars for which we obtained rotation ephemerides from radio observations, and discovered gamma-ray pulsations from 16 pulsars, 12 young and 4 recycled. PSR J2208+4056 has spindown power Edot = 8e32 erg/s, about three times lower than the previous observed gamma-ray emission "deathline". PSRs J2208+4056 and J1816-0755 have radio interpulses, constraining their geometry and perhaps enhancing their gamma-ray luminosity. We discuss whether the deathline is an artifact of selection bias due to the pulsar distance.

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.

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

Despite considerations of mass loss from stellar evolution suggesting otherwise, the content of gas in globular clusters seems poor and hence its measurement very elusive. One way of constraining the presence of ionized gas in a globular cluster is through its dispersive effects on the radiation of the millisecond pulsars included in the cluster. This effect led Freire et al. in 2001 to the first detection of any kind of gas in a globular cluster in the case of 47 Tucanae. By exploiting the results of 12 additional years of timing, as well as the observation of new millisecond pulsars in 47 Tucanae, we revisited this measurement: we first used the entire set of available timing parameters in order to measure the dynamical properties of the cluster and the three-dimensional position of the pulsars. Then we applied and tested various gas distribution models: assuming a constant gas density, we confirmed the detection of ionized gas with a number density of $n= 0.23\pm 0.05$ cm$^{-3}$, larger than the previous determination (at 2$\sigma$ uncertainty). Models predicting a decreasing density or following the stellar distribution density are highly disfavoured. We are also able to investigate the presence of an intermediate mass black hole in the centre of the cluster, showing that is not required by the available data, with an upper limit for the mass at $\sim 4000$ M$_{\odot}$.

We have observed the remnant of supernova SN~1987A (SNR~1987A), located in the Large Magellanic Cloud (LMC), to search for periodic and/or transient radio emission with the Parkes 64\,m-diameter radio telescope. We found no evidence of a radio pulsar in our periodicity search and derived 8$\sigma$ upper bounds on the flux density of any such source of $31\,\mu$Jy at 1.4~GHz and $21\,\mu$Jy at 3~GHz. Four candidate transient events were detected with greater than $7\sigma$ significance, with dispersion measures (DMs) in the range 150 to 840\,cm$^{-3}\,$pc. For two of them, we found a second pulse at slightly lower significance. However, we cannot at present conclude that any of these are associated with a pulsar in SNR~1987A. As a check on the system, we also observed PSR~B0540$-$69, a young pulsar which also lies in the LMC. We found eight giant pulses at the DM of this pulsar. We discuss the implications of these results for models of the supernova remnant, neutron star formation and pulsar evolution.

We report on single-pulse observations of the Galactic Center magnetar PSR J1745$-$2900 that were made using the Parkes 64-m radio telescope with a central frequency of 3.1 GHz at five observing epochs between 2013 July and August. The shape of the integrated pulse profiles was relatively stable across the five observations, indicating that the pulsar was in a stable state between MJDs 56475 and 56514. This extends the known stable state of this pulsar to 6.8 months. Short term pulse shape variations were also detected. It is shown that this pulsar switches between two emission modes frequently and that the typical duration of each mode is about ten minutes. No giant pulses or subpulse drifting were observed. Apparent nulls in the pulse emission were detected on MJD 56500. Although there are many differences between the radio emission of magnetars and normal radio pulsars, they also share some properties. The detection of mode changing and pulse nulling in PSR J1745$-$2900 suggests that the basic radio emission process for magnetars and normal pulsars is the same.

We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.

The first known pulsar glitch was discovered in the Vela pulsar at both Parkes and Goldstone in March 1969. Since then the number of known glitches has grown enormously, with more than 520 glitches now known in more than 180 pulsars. Details of glitch parameters and post-glitch recoveries are described and some implications for the physics of neutron stars are discussed.

Pulsars are remarkably precise "celestial clocks" that can be used to explore many different aspects of physics and astrophysics. In this article I give a brief summary of pulsar properties and describe some of the applications of pulsar timing, including tests of theories of gravitation, efforts to detect low-frequency gravitational waves using pulsar timing arrays and establishment a "pulsar timescale".

A brief, personal, and very incomplete account of 50 years of pulsar astronomy presented at the Conference Dinner for "Physics of Neutron Stars - 2017 - 50 Years After", held in Saint Petersburg, July 2017.

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.

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

Millisecond pulsars (MSPs) are short-period pulsars that are distinguished from "normal" pulsars, not only by their short period, but also by their very small spin-down rates and high probability of being in a binary system. These properties are consistent with MSPs having a different evolutionary history to normal pulsars, viz., neutron-star formation in an evolving binary system and spin-up due to accretion from the binary companion. Their very stable periods make MSPs nearly ideal probes of a wide variety of astrophysical phenomena. For example, they have been used to detect planets around pulsars, to test the accuracy of gravitational theories, to set limits on the low-frequency gravitational-wave background in the Universe, and to establish pulsar-based timescales that rival the best atomic-clock timescales in long-term stability. MSPs also provide a window into stellar and binary evolution, often suggesting exotic pathways to the observed systems. The X-ray accretion-powered MSPs, and especially those that transition between an accreting X-ray MSP and a non-accreting radio MSP, give important insight into the physics of accretion on to highly magnetised neutron stars.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

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

We derive the Sun's offset from the local mean Galactic plane($z_\odot$) using the observed $z$ distribution of young pulsars. Pulsar distances are obtained from measurements of annual parallax, HI absorption spectra or associations where available and otherwise from the observed pulsar dispersion and a model for the distribution of free electrons in the Galaxy. We fit the cumulative distribution function for a ${\rm sech}^2(z)$ distribution function, representing an isothermal self-gravitating disk, with uncertainties being estimated using the bootstrap method. We take pulsars having characteristic age $\tau_c<10^{6.5}$~yr and located within 4.5~kpc of the Sun, omitting those within the local spiral arm and those significantly affected by the Galactic warp, and solve for $z_\odot$ and the scale height, $H$, for different cutoffs in $\tau_c$. We compute these quantities using just the independently determined distances, and these together with DM-based distances separately using the YMW16 and NE2001 Galactic electron density models. We find that an age cutoff at $10^{5.75}$~yr with YMW16 DM-distances gives the best results with a minimum uncertainty in $z_\odot$ and an asymptotically stable value for $H$ showing that, at this age and below, the observed pulsar $z$-distribution is dominated by the dispersion in their birth locations. From this sample of 115 pulsars, we obtain $z_\odot=13.4\pm$4.4~pc and $H=56.9\pm$6.5~pc, similar to estimated scale heights for OB stars and open clusters. Consistent results are obtained using the independent-only distances and using the NE2001 model for the DM-based distances.

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

Pulse arrival time measurements over fourteen years with the Nanshan 25-m and Parkes 64-m radio telescopes have been used to determine the average profile and timing properties for PSR J1757$-$2421 (B1754$-$24). Analysis of the radio profile data shows a large variation of spectral index across the profile and an unusual increase in component separation with increasing frequency. The timing observations show that PSR B1754$-$24 underwent a large glitch with a fractional increase in spin frequency of $\Delta\nu_{\rm g}/\nu\sim 7.8\times10^{-6}$ in May 2011. At the time of the glitch there was a large permanent jump in the rotational frequency accompanied by two exponential recovery terms with timescales of 15(6) and 97(15) days, respectively.

We present a new model for the distribution of free electrons in the Galaxy, the Magellanic Clouds and the intergalactic medium (IGM) that can be used to estimate distances to real or simulated pulsars and fast radio bursts (FRBs) based on their dispersion measure (DM). The Galactic model has an extended thick disk representing the so-called warm interstellar medium, a thin disk representing the Galactic molecular ring, spiral arms based on a recent fit to Galactic HII regions, a Galactic Center disk and seven local features including the Gum Nebula, Galactic Loop I and the Local Bubble. An offset of the Sun from the Galactic plane and a warp of the outer Galactic disk are included in the model. Parameters of the Galactic model are determined by fitting to 189 pulsars with independently determined distances and DMs. Simple models are used for the Magellanic Clouds and the IGM. Galactic model distances are within the uncertainty range for 86 of the 189 independently determined distances and within 20\% of the nearest limit for a further 38 pulsars. We estimate that 95\% of predicted Galactic pulsar distances will have a relative error of less than a factor of 0.9. The predictions of YMW16 are compared to those of the TC93 and NE2001 models showing that YMW16 performs significantly better on all measures. Timescales for pulse broadening due to interstellar scattering are estimated for (real or simulated) Galactic and Magellanic Cloud pulsars and FRBs.

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.

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

A pulsar timing array is a Galactic-scale detector of nanohertz gravitational waves (GWs). Its target signals contain two components: the `Earth term' and the `pulsar term' corresponding to GWs incident on the Earth and pulsar respectively. In this work we present a Frequentist method for the detection and localization of continuous waves that takes into account the pulsar term and is significantly faster than existing methods. We investigate the role of pulsar terms by comparing a full-signal search with an Earth-term-only search for non-evolving black hole binaries. By applying the method to synthetic data sets, we find that (i) a full-signal search can slightly improve the detection probability (by about five percent); (ii) sky localization is biased if only Earth terms are searched for and the inclusion of pulsar terms is critical to remove such a bias; (iii) in the case of strong detections (with signal-to-noise ratio $\gtrsim$ 30), it may be possible to improve pulsar distance estimation through GW measurements.

We present proper motions for 15 pulsars which are observed regularly by the Nanshan 25-m radio telescope. Two methods, the frequentist method (Coles et al.2011) and the Bayesian (Lentati et al. 2014) method, are used and the results are compared. We demonstrate that the two methods can be applied to young pulsar data sets that exhibit large amounts of timing noise with steep spectral exponents and give consistent results. The measured positions also agree with very-long-baseline interferometric positions. Proper motions for four pulsars are obtained for the first time, and improved values are obtained for five pulsars.

Aims. To make a detailed study of the nulling and subpulse drifting in PSR J1727$-$2739 for investigation of its radiation properties. Methods. The observations were carried out on 20 March, 2004 using the Parkes 64-m radio telescope, with a central frequency of 1518 MHz. A total of 5568 single pulses were analysed. Results. This pulsar shows well defined nulls with lengths lasting from 6 to 281 pulses and separated by burst phases ranging from 2 to 133 pulses. We estimate a nulling fraction of around 68\%. No emission in the average pulse profile integrated over all null pulses is detected with significance above $3\sigma$. Most transitions from nulls to bursts are within a few pulses, whereas the transitions from bursts to nulls exhibit two patterns of decay: decrease gradually or rapidly. In the burst phase, we find that there are two distinct subpulse drift modes with vertical spacing between the drift bands of $9.7 \pm 1.6$ and $5.2 \pm 0.9$ pulse periods, while sometimes there is a third mode with no subpulse drifting. Some mode transitions occur within a single burst, while others are separated by nulls. Different modes have different average pulse profiles. Possible physical mechanisms are discussed.

Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $\approx$ 1300 d time scale, and that a larger fraction ( $\gtrsim 25\%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $\approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars.

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.

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.

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.

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.

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

In this paper, we report radio observations of the Galactic Center magnetar PSR J1745-2900 at six epochs between June and October, 2014. These observations were carried out using the new Shanghai Tian Ma Radio Telescope at a frequency of 8.6 GHz. Both the flux density and integrated profile of PSR J1745-2900 show dramatic changes from epoch to epoch showing that the pulsar was in its "erratic" phase. On MJD 56836, the flux density of this magnetar was about 8.7 mJy, which was ten times large than that reported at the time of discovery, enabling a single-pulse analysis. The emission is dominated by narrow "spiky" pulses which follow a log-normal distribution in peak flux density. From 1913 pulses, we detected 53 pulses whose peak flux density is ten times greater than that of the integrated profile. They are concentrated in pulse phase at the peaks of the integrated profile. The pulse widths at the 50% level of these bright pulses was between 0.2 to 0.9 deg, much narrower than that of integrated profile (~12 deg). The observed pulse widths may be limited by interstellar scattering. No clear correlation was found between the widths and peak flux density of these pulses and no evidence was found for subpulse drifting. Relatively strong spiky pulses are also detected in the other five epochs of observation, showing the same properties as that detected in MJD 56836. These strong spiky pulses cannot be classified as "giant" pulses but are more closely related to normal pulse emission.

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.

Complementary to ground-based laser interferometers, pulsar timing array experiments are being carried out to search for nanohertz gravitational waves. Using the world's most powerful radio telescopes, three major international collaborations have collected $\sim$10-year high precision timing data for tens of millisecond pulsars. In this paper we give an overview on pulsar timing experiments, gravitational wave detection in the nanohertz regime, and recent results obtained by various timing array projects.

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.

Pulsars are wonderful gravitational probes. Their tiny size and stellar mass give their rotation periods a stablility comparable to that of atomic frequency standards. This is especially true of the rapidly rotating "millisecond pulsars" (MSPs). Many of these rapidly rotating pulsars are in orbit with another star, allowing pulsar timing to probe relativistic perturbations to the orbital motion. Pulsars have provided the most stringent tests of theories of relativistic gravitation, especially in the strong-field regime, and have shown that Einstein's general theory of relativity is an accurate description of the observed motions. Many other gravitational theories are effectively ruled out or at least severely constrained by these results. MSPs can also be used to form a "Pulsar Timing Array" (PTA). PTAs are Galactic-scale interferometers that have the potential to directly detect nanohertz gravitational waves from astrophysical sources. Orbiting super-massive black holes in the cores of distant galaxies are the sources most likely to be detectable. Although no evidence for gravitational waves has yet been found in PTA data sets, the latest limits are seriously constraining current ideas on galaxy and black-hole evolution in the early Universe.

We present timing observations of four millisecond pulsars discovered in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. PSRs J1552-4937 and J1843-1448 are isolated objects with spin periods of 6.28 and 5.47 ms respectively. PSR J1727-2946 is in a 40-day binary orbit and has a spin period of 27 ms. The 4.43-ms pulsar J1813-2621 is in a circular 8.16-day binary orbit around a low-mass companion star with a minimum companion mass of 0.2 solar masses. Combining these results with detections from five other Parkes multibeam surveys, gives a well-defined sample of 56 pulsars with spin periods below 20 ms. We develop a likelihood analysis to constrain the functional form which best describes the underlying distribution of spin periods for millisecond pulsars. The best results were obtained with a log-normal distribution. A gamma distribution is less favoured, but still compatible with the observations. Uniform, power-law and Gaussian distributions are found to be inconsistent with the data. Galactic millisecond pulsars being found by current surveys appear to be in agreement with a log-normal distribution which allows for the existence of pulsars with periods below 1.5 ms.

Pulsar timing arrays (PTAs) can be used to search for very low frequency ($10^{-9}$--$10^{-7}$ Hz) gravitational waves (GWs). In this paper we present a general method for the detection and localization of single-source GWs using PTAs. We demonstrate the effectiveness of this new method for three types of signals: monochromatic waves as expected from individual supermassive binary black holes in circular orbits, GWs from eccentric binaries and GW bursts. We also test its implementation in realistic data sets that include effects such as uneven sampling and heterogeneous data spans and measurement precision. It is shown that our method, which works in the frequency domain, performs as well as published time-domain methods. In particular, we find it equivalent to the $\mathcal{F}_{e}$-statistic for monochromatic waves. We also discuss the construction of null streams -- data streams that have null response to GWs, and the prospect of using null streams as a consistency check in the case of detected GW signals. Finally, we present sensitivities to individual supermassive binary black holes in eccentric orbits. We find that a monochromatic search that is designed for circular binaries can efficiently detect eccentric binaries with both high and low eccentricities, while a harmonic summing technique provides greater sensitivities only for binaries with moderate eccentricities.

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.

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

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.

In a series of three previous papers, formulation and specifics of the representation of World Coordinate Transformations in FITS data have been presented. This fourth paper deals with encoding time. Time on all scales and precisions known in astronomical datasets is to be described in an unambiguous, complete, and self-consistent manner. Employing the well--established World Coordinate System (WCS) framework, and maintaining compatibility with the FITS conventions that are currently in use to specify time, the standard is extended to describe rigorously the time coordinate. World coordinate functions are defined for temporal axes sampled linearly and as specified by a lookup table. The resulting standard is consistent with the existing FITS WCS standards and specifies a metadata set that achieves the aims enunciated above.

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

The Five-hundred-meter Aperture Spherical Telescope (FAST) will become one of the world-leading telescopes for pulsar timing array (PTA) research. The primary goals for PTAs are to detect (and subsequently study) ultra-low-frequency gravitational waves, to develop a pulsar-based time standard and to improve solar system planetary ephemerides. FAST will have the sensitivity to observe known pulsars with significantly improved signal-to-noise ratios and will discover a large number of currently unknown pulsars. We describe how FAST will contribute to PTA research and show that jitter- and timing-noise will be the limiting noise processes for FAST data sets. Jitter noise will limit the timing precision achievable over data spans of a few years while timing noise will limit the precision achievable over many years.

The pulsar PSR J1756$-$2251 resides in a relativistic double neutron star (DNS) binary system with a 7.67-hr orbit. We have conducted long-term precision timing on more than 9 years of data acquired from five telescopes, measuring five post-Keplerian parameters. This has led to several independent tests of general relativity (GR), the most constraining of which shows agreement with the prediction of GR at the 4% level. Our measurement of the orbital decay rate disagrees with that predicted by GR, likely due to systematic observational biases. We have derived the pulsar distance from parallax and orbital decay measurements to be 0.73$_{-0.24}^{+0.60}$ kpc (68%) and < 1.2 kpc (95% upper limit), respectively; these are significantly discrepant from the distance estimated using Galactic electron density models. We have found the pulsar mass to be 1.341$\pm$0.007 M$_\odot$, and a low neutron star (NS) companion mass of 1.230$\pm$0.007 M$_\odot$. We also determined an upper limit to the spin-orbit misalignment angle of 34{\deg} (95%) based on a system geometry fit to long-term profile width measurements. These and other observed properties have led us to hypothesize an evolution involving a low mass loss, symmetric supernova progenitor to the second-formed NS companion, as is thought to be the case for the double pulsar system PSR J0737$-$3039A/B. This would make PSR J1756$-$2251 the second compact binary system providing concrete evidence for this type of NS formation channel.

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.

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

We present an analysis of $23~$yr of pulse arrival times for PSR B1259$-$63. The pulsar is in a binary orbit about its approximately $20~M_{\odot}$ companion LS$~$2883. Our best-fitting timing solution has none of the pulse-number ambiguities that have plagued previous attempts to model the binary orbit. We measure significant first and second time derivatives of the projected semimajor axis and longitude of the periastron of the orbit. These variations are found to be consistent with the precession of the orbital plane due to classical spin-orbit coupling. The derived moment of inertia and spin of the companion are consistent with the companion rotating at near-breakup velocity. The system configuration is also consistent with the geometry derived from both the polarisation of the radio emission and the eclipse of the pulsar by the equatorial disc of the companion. We find strong evidence for orbital period decay that can be attributed to mass loss from the companion star. We also measure a significant proper motion that locates the birth of the system in the Centaurus OB1 association. By combining proper motion of the pulsar with radial velocity measurements of the companion, we measure the three-dimensional velocity of the system. This velocity is used to constrain the masses of the stars prior to the supernova explosion and the kick the pulsar received at or immediately after the explosion.

As the remnant of Supernova (SN) 1987A has been getting brighter over time, new observations at high frequencies have allowed imaging of the radio emission at unprecedented detail. We present a new radio image at 44 GHz of the supernova remnant (SNR), derived from observations performed with the Australia Telescope Compact Array (ATCA) in 2011. The diffraction-limited image has a resolution of $349\times225$ mas, which is the highest achieved to date in high-dynamic range images of the SNR. We also present a new image at 18 GHz, also derived from ATCA observations performed in 2011, which is super-resolved to $0./!/!^{\prime\prime}25$. The new 44 and 18 GHz images yield the first high-resolution spectral index map of the remnant. The comparison of the 44 GHz image with contemporaneous X-ray and H$\alpha$ observations allows further investigations of the nature of the remnant asymmetry and sheds more light into the progenitor hypotheses and SN explosion. In light of simple free-free absorption models, we discuss the likelihood of detecting at 44 GHz the possible emission originating from a pulsar wind nebula (PWN) or a compact source in the centre of the remnant.

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

We present radio imaging observations of supernova remnant 1987A at 9 GHz, taken with the Australia Telescope Compact Array over 21 years from 1992 to 2013. By employing a Fourier modeling technique to fit the visibility data, we show that the remnant structure has evolved significantly since day 7000 (mid-2006): the emission latitude has gradually decreased, such that the overall geometry has become more similar to a ring structure. Around the same time, we find a decreasing trend in the east-west asymmetry of the surface emissivity. These results could reflect the increasing interaction of the forward shock with material around the circumstellar ring, and the relative weakening of the interaction with the lower-density material at higher latitudes. The morphological evolution caused an apparent break in the remnant expansion measured with a torus model, from a velocity of 4600+150-200 km/s between day 4000 and 7000 to 2400+100-200 km/s after day 7000. However, we emphasize that there is no conclusive evidence for a physical slowing of the shock at any given latitude in the expanding remnant, and that a change of radio morphology alone appears to dominate the evolution. This is supported by our ring-only fits which show a constant expansion of 3890+/-50 km/s without deceleration between days 4000 and 9000. We suggest that once the emission latitude no longer decreases, the expansion velocity obtained from the torus model should return to the same value as that measured with the ring model.

The International Pulsar Timing Array (IPTA) is an organisation whose raison d'etre is to facilitate collaboration between the three main existing PTAs (the EPTA in Europe, NANOGrav in North America and the PPTA in Australia) in order to realise the benefits of combined PTA data sets in reaching the goals of PTA projects. Currently, shared data sets for 39 pulsars are available for IPTA-based projects. Operation of the IPTA is administered by a Steering Committee consisting of six members, two from each PTA, plus the immediate past Chair in a non-voting capacity. A Constitution and several Agreements define the framework for the collaboration. Web pages provide information both to members of participating PTAs and to the general public. With support from an NSF PIRE grant, the IPTA facilitates the organisation of annual Student Workshops and Science Meetings. These are very valuable both in training new students and in communicating current results from IPTA-based research.

We present a study of PSR J1723-2837, an eclipsing, 1.86 ms millisecond binary radio pulsar discovered in the Parkes Multibeam survey. Radio timing indicates that the pulsar has a circular orbit with a 15 hr orbital period, a low-mass companion, and a measurable orbital period derivative. The eclipse fraction of ~15% during the pulsar's orbit is twice the Roche lobe size inferred for the companion. The timing behavior is significantly affected by unmodeled systematics of astrophysical origin, and higher-order orbital period derivatives are needed in the timing solution to account for these variations. We have identified the pulsar's (non-degenerate) companion using archival ultraviolet, optical, and infrared survey data and new optical photometry. Doppler shifts from optical spectroscopy confirm the star's association with the pulsar and indicate a pulsar-to-companion mass ratio of 3.3 +/- 0.5, corresponding to a companion mass range of 0.4 to 0.7 Msun and an orbital inclination angle range of between 30 and 41 degrees, assuming a pulsar mass range of 1.4-2.0 Msun. Spectroscopy indicates a spectral type of G for the companion and an inferred Roche-lobe-filling distance that is consistent with the distance estimated from radio dispersion. The features of PSR J1723-2837 indicate that it is likely a "redback" system. Unlike the five other Galactic redbacks discovered to date, PSR J1723-2837 has not been detected as a gamma-ray source with Fermi. This may be due to an intrinsic spin-down luminosity that is much smaller than the measured value if the unmeasured contribution from proper motion is large.

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

We present new imaging observations of the remnant of Supernova (SN) 1987A at 44 GHz, performed in 2011 with the Australia Telescope Compact Array (ATCA). The $0\farcs35\times0\farcs23$ resolution of the diffraction-limited image is the highest achieved to date in high-dynamic range. We also present a new ATCA image at 18 GHz derived from 2011 observations, which is super-resolved to $0\farcs25$. The flux density is 40$\pm$2 mJy at 44 GHz and 81$\pm$6 mJy at 18 GHz. At both frequencies, the remnant exhibits a ring-like emission with two prominent lobes, and an east-west brightness asymmetry that peaks on the eastern lobe. A central feature of fainter emission appears at 44 GHz. A comparison with previous ATCA observations at 18 and 36 GHz highlights higher expansion velocities of the remnant eastern side. The 18-44 GHz spectral index is $\alpha=-0.80$ ($S_{\nu}\propto\nu^{\alpha}$). The spectral index map suggests slightly steeper values at the brightest sites on the eastern lobe, whereas flatter values are associated with the inner regions. The remnant morphology at 44 GHz generally matches the structure seen with contemporaneous X-ray and H$\alpha$ observations. Unlike the H$\alpha$ emission, both the radio and X-ray emission peaks on the eastern lobe. The regions of flatter spectral index align and partially overlap with the optically-visible ejecta. Simple free-free absorption models suggest that emission from a pulsar wind nebula or a compact source inside the remnant may now be detectable at high frequencies, or at low frequencies if there are holes in the ionised component of the ejecta.

The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

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.

The study of high-magnetic-field pulsars is important for examining the relationships between radio pulsars, magnetars, and X-ray-isolated neutron stars (XINSs). Here we report on X-ray observations of three such high-magnetic-field radio pulsars. We first present the results of a deep XMM-Newton observation of PSR J1734-3333, taken to follow up on its initial detection in 2009. The pulsar's spectrum is well fit by a blackbody with a temperature of 300 +/- 60 eV, with bolometric luminosity L_bb = 2.0(+2.2 -0.7)e+32 erg/s = 0.0036E_dot for a distance of 6.1 kpc. We detect no X-ray pulsations from the source, setting a 1 sigma upper limit on the pulsed fraction of 60% in the 0.5-3 keV band. We compare PSR J1734-3333 to other rotation-powered pulsars of similar age and find that it is significantly hotter, supporting the hypothesis that the magnetic field affects the observed thermal properties of pulsars. We also report on XMM-Newton and Chandra observations of PSRs B1845-19 and J1001-5939. We do not detect either pulsar, setting 3 sigma upper limits on their blackbody temperatures of 48 and 56 eV, respectively. Despite the similarities in rotational properties, these sources are significantly cooler than all but one of the XINSs, which we attribute to the two groups having been born with different magnetic fields and hence evolving differently.

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

We report on an X-ray observation of the 166 Myr old radio pulsar J0108-1431 with XMM-Newton. The X-ray spectrum can be described by a power-law model with a relatively steep photon index Gamma~3 or by a combination of thermal and non-thermal components, e.g., a power-law component with fixed photon index Gamma~2 plus a blackbody component with a temperature of kT=0.11 keV. The two-component model appears more reasonable considering different estimates for the hydrogen column density. The non-thermal X-ray efficiency in the single power-law model is eta^PL (1-10 keV) = L^PL (1-10 keV) / Edot ~ 0.003, higher than in most other X-ray detected pulsars. In the case of the combined model, the non-thermal and thermal X-ray efficiencies are even higher, eta^PL (1-10 keV) ~ eta^bb ~ 0.006. We detected X-ray pulsations at the radio period of P=0.808s with significance of 7sigma. The pulse shape in the folded X-ray lightcurve (0.15-2 keV) is asymmetric, with statistically significant contributions from up to 5 leading harmonics. Pulse profiles at two different energy ranges differ slightly: the profile is asymmetric at low energies, 0.15-1 keV, while at higher energies, 1-2 keV, it has a nearly sinusodial shape. The radio pulse peak leads the 0.15-2 keV X-ray pulse peak by delta phi = 0.06 +/- 0.03.

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.

We report on the detection of sporadic, strong single pulses co-existing with a periodic weak emission in the duration of weak mode of PSR B0826-34. The intensities and durations of these pulses are comparable with that of the sub-pulses in the strong mode, and these pulses are distributed within the phase ranges of the main-pulse and interpulse of the strong-mode average profile. These results suggest that there are most possibly sporadic, very short timescale turn-on of strong-mode emission during the weak-mode state of the pulsar. The emission features of the bursts of strong pulses of PSR B0826-34 during its weak-mode state are similar to those of the rotating radio transients (RRATs). PSR B0826-34 is the second pulsar known which oscillates between pulsar-like and RRAT-like modes.

We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

The first pulsar observations were made at Parkes on March 8, 1968, just 13 days after the publication of the discovery paper by Hewish and Bell. Since then, Parkes has become the world's most successful pulsar search machine, discovering nearly two thirds of the known pulsars, among them many highly significant objects. It has also led the world in pulsar polarisation and timing studies. In this talk I will review the highlights of pulsar work at Parkes from those 1968 observations to about 2006 when the Parkes Multibeam Pulsar Survey was essentially completed and the Parkes Pulsar Timing Array project was established.

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

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

PSR J1740-3052 is a young pulsar in orbit around a companion that is most likely a B-type main-sequence star. Since its discovery more than a decade ago, data have been taken at several frequencies with instruments at the Green Bank, Parkes, Lovell, and Westerbork telescopes. We measure scattering timescales in the pulse profiles and dispersion measure changes as a function of binary orbital phase and present evidence that both of these vary as would be expected due to a wind from the companion star. Using pulse arrival times that have been corrected for the observed periodic dispersion measure changes, we find a timing solution spanning 1997 November to 2011 March. This includes measurements of the advance of periastron and the change in the projected semimajor axis of the orbit and sets constraints on the orbital geometry. From these constraints, we estimate that the pulsar received a kick of at least ~50 km/s at birth. A quasi-periodic signal is present in the timing residuals with a period of 2.2 times the binary orbital period. The origin of this signal is unclear.

We report on radiation properties of extreme nulling pulsar J1502-5653, by analyzing the data acquired from the Parkes 64-m telescope at 1374 MHz. The radio emission from this pulsar exhibits sequences of several tens to several hundreds consecutive burst pulses, separated by null pulses, and the appearance of the emission seems quasi-periodic. The null fraction from the data is estimated to be 93.6%. No emission is detected in the integrated profile of all null pulses. Systematic modulations of pulse intensity and phase are found at the beginning of burst-pulse sequences just after null. The intensity usually rises to a maximum for the first few pulses, then declines exponentially afterwards, and becomes stable after few tens of pulse periods. The peak phase appears at later longitudes for the first pulse, then drifts to earlier longitudes rapidly, and then systematic drifting gradually vanishes while the intensity becomes stable. In this pulsar, the intensity variation and phase modulation of pulses are correlated in a short duration after the emission starts following a null. Observed properties of the pulsar are compared with other nulling pulsars published previously, and the possible explanation for phase modulation is discussed.

The Parkes Pulsar Data Archive currently provides access to 165,755 data files obtained from observations carried out at the Parkes Observatory since the year 1991. Data files and access methods are compliant with the Virtual Observatory protocol. This paper provides a tutorial on how to make use of the Parkes Pulsar Data Archive and provides example queries using on-line interfaces.

We investigate the changes in polarization position angle in radiation from pulsar A around the eclipse in the Double Pulsar system PSR J0737-3039A/B at the 20 cm and 50 cm wavelengths using the Parkes 64-m telescope. The changes are ~2\sigma\ during and shortly after the eclipse at 20 cm but less significant at 50 cm. We show that the changes in position angle during the eclipse can be modelled by differential synchrotron absorption in the eclipse regions. Position angle changes after the eclipse are interpreted as Faraday rotation in the magnetotail of pulsar B. Implied charge densities are consistent with the Goldreich-Julian density, suggesting that the particle energies in the magnetotail are mildly relativistic.

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

We present 35 ks Chandra ACIS observations of the 42 Myr old radio pulsar PSR B1451-68. A point source is detected 0.32" +/- 0.73" from the expected radio pulsar position. It has ~200 counts in the 0.3-8 keV energy range. We identify this point source as the X-ray counterpart of the radio pulsar. PSR B1451-68 is located close to a 2MASS point source, for which we derive 7% as the upper limit on the flux contribution to the measured pulsar X-ray flux. The pulsar spectrum can be described by either a power-law model with photon index Gamma=2.4 (+0.4/-0.3) and a unrealistically high absorbing column density N(H)= (2.5 (+1.2/-1.3)) * 10^(21) cm^-2, or by a combination of a kT=0.35 (+0.12/-0.07) keV blackbody and a Gamma = 1.4 +/- 0.5 power-law component for N(H)[DM]= 2.6 * 10^(20) cm^-2, estimated from the pulsar dispersion measure. At the parallactic, Lutz-Kelker bias corrected distance of 480 pc, the non-thermal X-ray luminosities in the 0.3-8 keV energy band are either Lx(nonth)= (11.3 +/- 1.7) * 10^(29) erg/s or Lx(nonth)= (5.9 (+4.9/-5.0)) * 10^(29) erg/s, respectively. This corresponds to non-thermal X-ray efficiencies of either eta(nonth)= Lx(nonth) / (dE/dt) ~ 0.005 or 0.003, respectively.

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

We present the first optical observations of the unique system J0737-3039 (composed of two pulsars, hereafter PSR-A and PSR-B). Ultra-deep optical observations, performed with the High Resolution Camera of the Advanced Camera for Surveys on board the Hubble Space Telescope could not detect any optical emission from the system down to m_F435W=27.0 and m_F606W=28.3. The estimated optical flux limits are used to constrain the three-component (two thermal and one non-thermal) model recently proposed to reproduce the XMM-Newton X-ray spectrum. They suggest the presence of a break at low energies in the non-thermal power law component of PSR-A and are compatible with the expected black-body emission from the PSR-B surface. The corresponding efficiency of the optical emission from PSR-A's magnetosphere would be comparable to that of other Myr-old pulsars, thus suggesting that this parameter may not dramatically evolve over a time-scale of a few Myr.

We present 13 cm and 20 cm radio observations of the magnetic chemically peculiar star CU Virginis taken with the Australia Telescope Compact Array. We detect two circularly polarised radio pulses every rotation period which confirm previous detections. In the first pulse, the lower frequency emission arrives before the higher frequency emission and the ordering reverses in the second pulse. In order to explain the frequency dependence of the time between the two pulses, we construct a geometric model of the magnetosphere of CU Virginis, and consider various emission angles relative to the magnetic field lines. A simple electron cyclotron maser emission model, in which the emission is perpendicular to the magnetic field lines, is not consistent with our data. A model in which the emission is refracted through cold plasma in the magnetosphere is shown to have the correct pulse arrival time frequency dependence.

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than seven years; four pulsars are in binary systems and one is isolated. We are able to measure the pulsars' proper motions and derive an estimate for their space velocities. The measured two-dimensional velocities are in the range 70-210 km/s, consistent with those measured for other MSPs. We also use all the available proper motion information for isolated and binary MSPs to update the known velocity distribution for these populations. As found by earlier works, we find that the velocity distribution of binary and isolated MSPs are indistinguishable with the current data. Four of the pulsars in our observing program are highly recycled with low-mass white dwarf companions and we are able to derive accurate binary parameters for these systems. For three of these binary systems we are able to place initial constraints on the pulsar masses with best-fit values in the range 1.0-1.6 M_sun. The implications of the results presented here to our understanding of binary pulsar evolution are discussed. The updated parameters for the binary systems studied here, together with recently discovered similar systems, allowed us to update previous limits on the the violation of the strong equivalence principle through the parameter |Delta| to 4.6x10^-3 (95% confidence) and the violation of Lorentz-invariance/momentum-conservation through the parameter |hat{alpha}_3| to 5.5x10^-20 (95% confidence).

PSR J1734-3333 is a radio pulsar rotating with a period P=1.17 s and slowing down with a period derivative Pdot=2.28 x 10^{-12}, the third largest among rotation-powered pulsars. These properties are midway between those of normal rotation-powered pulsars and magnetars, two populations of neutron stars that are notably different in their emission properties. Here we report on the measurement of the second period derivative of the rotation of PSR J1734-3333 and calculate a braking index n=0.9 +- 0.2. This value is well below 3, the value expected for an electromagnetic braking due to a constant magnetic dipole, and indicates that this pulsar may soon have the rotational properties of a magnetar. While there are several mechanisms which could lead to such a low braking index, we discuss this observation, together with the properties exhibited by some other high-Pdot rotation-powered pulsars, and interpret it as evidence of a possible evolutionary route for magnetars through a radio-pulsar phase, supporting a unified description of the two classes of object.

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

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

We have analysed the long- and short-term time dependence of the pulse arrival times and the pulse detection rates for eight Rotating Radio Transient (RRAT) sources from the Parkes Multi-beam Pulsar Survey (PMPS). We find significant periodicities in the individual pulse arrival times from six RRATs. These periodicities range from 30 minutes to 2100 days and from one to 16 independent (i.e. non-harmonically related) periodicities are detected for each RRAT. In addition, we find that pulse emission is a random (i.e. Poisson) process on short (hour-long) time scales but that most of the objects exhibit longer term (months-years) non-random behaviour. We find that PSRs J1819-1458 and J1317-5759 emit more doublets (two consecutive pulses) and triplets (three consecutive pulses) than is expected in random pulse distributions. No evidence for such an excess is found for the other RRATs. There are several different models for RRAT emission depending on both extrinsic and intrinsic factors which are consistent with these properties.

High-magnetic-field pulsars represent an important class of objects for studying the relationship between magnetars and radio pulsars. Here we report on four Chandra observations of the high-magnetic-field pulsar J1718-3718 ($B=7.4\times10^{13}$\,G) taken in 2009 as well as on a re-analysis of 2002 Chandra observations of the region. We also report an improved radio position for this pulsar based on ATCA observations. We detect X-ray pulsations at the pulsar's period in the 2009 data, with pulsed fraction 52%$\pm$13% in the 0.8--2.0 keV band. We find that the X-ray pulse is aligned with the radio pulse. The data from 2002 and 2009 show consistent spectra and fluxes; a merged overall spectrum is well fit by a blackbody of temperature $186^{+19}_{-18}$\,eV, slightly higher than predicted by standard cooling models, however, the best-fit neutron star atmosphere model is consistent with standard cooling. We find the bolometric luminosity $L^{\infty}_{\rm bb}=4^{+5}_{-2} \times10^{32}$\,erg\,s$^{-1}\sim0.3\dot{E}$, for a distance of 4.5 kpc. We compile measurements of the temperatures of all X-ray detected high-$B$ pulsars as well as those of low-$B$ radio pulsars and find evidence for the former being on average hotter than the latter.

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.

Pulsar timing observations are used to place constraints on the rate of coalescence of supermassive black-hole (SMBH) binaries as a function of mass and redshift. In contrast to the indirect constraints obtained from other techniques, pulsar timing observations provide a direct constraint on the black-hole merger rate. This is possible since pulsar timing is sensitive to the gravitational waves (GWs) emitted by these sources in the final stages of their evolution. We find that upper bounds calculated from the recently published Parkes Pulsar Timing Array data are just above theoretical predictions for redshifts below 10. In the future, with improved timing precision and longer data spans, we show that a non-detection of GWs will rule out some of the available parameter space in a particular class of SMBH binary merger models. We also show that if we can time a set of pulsars to 10ns timing accuracy, for example, using the proposed Square Kilometre Array, it should be possible to detect one or more individual SMBH binary systems.

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 describe our studies of the radio and high-energy properties of Rotating Radio Transients (RRATs). We find that the radio pulse intensity distributions are log-normal, with power-law tails evident in two cases. For the three RRATs with coverage over a wide range of frequency, the mean spectral index is -1.7\pm0.1, roughly in the range of normal pulsars. We do not observe anomalous magnetar-like spectra for any RRATs. Our 94-ks XMM-Newton observation of the high magnetic field RRAT J1819-1458 reveals a blackbody spectrum (kT ~130 eV) with an unusual absorption feature at ~1 keV. We find no evidence for X-ray bursts or other X-ray variability. We performed a correlation analysis of the X-ray photons with radio pulses detected in concurrent observations with the Green Bank, Effelsberg, and Parkes telescopes. We find no evidence for any correlations between radio pulse emission and X-ray photons, perhaps suggesting that sporadicity is not due to variations in magnetospheric particle density but to changes in beaming or coherence.

Millisecond pulsars are intrinsically very stable clocks and precise measurement of their observed pulse periods can be used to study a wide variety of astrophysical phenomena. In particular, observations of a large sample of millisecond pulsars at regular intervals, constituting a Pulsar Timing Array (PTA), can be used as a detector of low-frequency gravitational waves and to establish a standard of time independent of terrestrial atomic timescales. Three major timing array projects have been established: The European Pulsar Timing Array (EPTA), the North American pulsar timing array (NANOGrav) and the Parkes Pulsar Timing Array (PPTA). Results from the PPTA project are described in some detail and future prospects for PTA projects are discussed.

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 show how pulsar observations may be used to construct a time standard that is independent of terrestrial time standards. The pulsar time scale provides a method to determine the stability of terrestrial time standards over years to decades. Here, we summarise the method, provide initial results and discuss the possibilities and limitations of our pulsar time scale.

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

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

The magnetic chemically peculiar star CU Virginis is a unique astrophysical laboratory for stellar magnetospheres and coherent emission processes. It is the only known main sequence star to emit a radio pulse every rotation period. Here we report on new observations of the CU Virginis pulse profile in the 13 and 20\,cm radio bands. The profile is known to be characterised by two peaks of 100$\%$ circularly polarised emission that are thought to arise in an electron-cyclotron maser mechanism. We find that the trailing peak is stable at both 13 and 20\,cm, whereas the leading peak is intermittent at 13\,cm. Our measured pulse arrival times confirm the discrepancy previously reported between the putative stellar rotation rates measured with optical data and with radio observations. We suggest that this period discrepancy might be caused by an unknown companion or by instabilities in the emission region. Regular long-term pulse timing and simultaneous multi-wavelength observations are essential to clarify the behaviour of this emerging class of transient radio source.

We present the evolution of the radio emission from the 2.8-s pulsar of the double pulsar system PSR J0737-3039A/B. We provide an update on the Burgay et al. (2005) analysis by describing the changes in the pulse profile and flux density over five years of observations, culminating in the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases by 0.177 mJy/yr at the brightest orbital phases and the pulse profile evolves from a single to a double peak, with a separation rate of 2.6 deg/yr. The pulse profile changes are most likely caused by relativistic spin precession, but can not be easily explained with a circular hollow-cone beam as in the model of Clifton & Weisberg (2008). Relativistic spin precession, coupled with an elliptical beam, can model the pulse profile evolution well. This particular beam shape predicts geometrical parameters for the two bright orbital phases which are consistent and similar to those derived by Breton et al. (2008). However, the observed decrease in flux over time and B's eventual disappearance cannot be easily explained by the model and may be due to the changing influence of A on B.

Seven years of pulse time-of-arrival measurements have been collected from observations of the young pulsar PSR B2334+61 using the Nanshan radio telescope of Urumqi Observatory. A phase-connected solution has been obtained over the whole data span, 2002 August to 2009 August. This includes a very large glitch that occurred between 2005 August 26 and September 8 (MJDs 53608 and 53621). The relative increase in rotational frequency for this glitch, $\Delta\nu_{g}/\nu~\sim~20.5\times10^{-6}$, is the largest ever seen. Although accounting for less than 1\% of the glitch, there were two well-defined exponential decay terms with time constants of 21 and 147 days respectively. There was also a large long-term increase in the spindown rate with $\Delta\dot\nu_p/\dot\nu \sim 0.011$ at the time of the glitch. A highly significant oscillation with a period of close to one year is seen in the post-glitch residuals. It is very unlikely that this can be accounted for by a pulsar position error or proper motion -- it appears to result from effects interior to the neutron star. Implications of these results for pulsar glitch models are discussed.

We investigate the radio and gamma-ray beaming properties of normal and millisecond pulsars by selecting two samples from the known populations. The first, Sample G, contains pulsars which are detectable in blind searches of gamma-ray data from the Fermi Large Area Telescope. The second, Sample R, contains pulsars detectable in blind radio searches which have spin-down luminosities Edot > 10^{34} erg/s. We analyse the fraction of the gamma-ray-selected Sample G which have detectable radio pulses and the fraction of the radio-selected Sample R which have detectable gamma-ray pulses. Twenty of our 35 Sample G pulsars have already observed radio pulses. This rules out low-altitude polar-cap beaming models if, as is currently believed, gamma-ray beams are generated in the outer magnetosphere and are very wide. We further find that, for the highest-Edot pulsars, the radio and gamma-ray beams have comparable beaming factors, i.e., the beams cover similar regions of the sky as the star rotates. For lower-Edot gamma-ray emitting pulsars, the radio beams have about half of the gamma-ray sky coverage. These results suggest that, for high-Edot young and millisecond pulsars, the radio emission originates in wide beams from regions high in the pulsar magnetosphere, probably close to the null-charge surface and to the gamma-ray emitting regions. Furthermore, it suggests that for these high-Edot pulsars, as in the gamma-ray case, features in the radio profile represent caustics in the emission beam pattern.

Pulsars are very stable clocks in space which have many applications to problems in physics and astrophysics. Observations of double-neutron-star binary systems have given the first observational evidence for the existence of gravitational waves (GWs) and shown that Einstein's general theory of relativity is an accurate description of gravitational interactions in the regime of strong gravity. Observations of a large sample of pulsars spread across the celestial sphere forming a "Pulsar Timing Array" (PTA), can in principle enable a positive detection of the GW background in the Galaxy. The Parkes Pulsar Timing Array (PPTA) is making precise timing measurements of 20 millisecond pulsars at three radio frequencies and is approaching the level of timing precision and data spans which are needed for GW detection. These observations will also allow us to establish a "Pulsar Timescale" and to detect or limit errors in the Solar System ephemerides used in pulsar timing analyses. Combination of PPTA data with that of other groups to form an International Pulsar Timing Array (IPTA) will enhance the sensitivity to GWs and facilitate reaching other PTA goals. The principal source of GWs at the nanoHertz frequencies to which PTAs are sensitive is believed to be super-massive binary black holes in the cores of distant galaxies. Current results do not signficantly limit models for formation of such black-hole binary systems, but in a few years we expect that PTAs will either detect GWs or seriously constrain current ideas about black-hole formation and galaxy mergers. Future instruments such as the Square Kilometre Array (SKA) should not only detect GWs from astrophysical sources but also enable detailed studies of the sources and the gravitational theories used to account for the GW emission.

We report the discovery of 31.18 ms pulsations from the INTEGRAL source IGR J14003-6326 using the Rossi X-ray Timing Explorer (RXTE). This pulsar is most likely associated with the bright Chandra X-ray point source lying at the center of G310.6-1.6, a previously unrecognised Galactic composite supernova remnant with a bright central non-thermal radio and X-ray nebula, taken to be the pulsar wind nebula (PWN). PSR J1400-6325 is amongst the most energetic rotation-powered pulsars in the Galaxy, with a spin-down luminosity of Edot = 5.1E+37 erg.s-1. In the rotating dipole model, the surface dipole magnetic field strength is B_s = 1.1E+12 G and the characteristic age tau_c = P/2Pdot = 12.7 kyr. The high spin-down power is consistent with the hard spectral indices of the pulsar and the nebula of 1.22 +/- 0.15 and 1.83 +/- 0.08, respectively, and a 2-10 keV flux ratio F_PWN/F_PSR ~ 8. Follow-up Parkes observations resulted in the detection of radio emission at 10 and 20 cm from PSR J1400-6325 at a dispersion measure of ~ 560 cm-3 pc, which implies a relatively large distance of 10 +/- 3 kpc. However, the resulting location off the Galactic Plane of ~ 280 pc would be much larger than the typical thickness of the molecular disk, and we argue that G310.6-1.6 lies at a distance of ~ 7 kpc. There is no gamma-ray counterpart to the nebula or pulsar in the Fermi data published so far. A multi-wavelength study of this new composite supernova remnant, from radio to very-high energy gamma-rays, suggests a young (< 1000 yr) system, formed by a sub-energetic (~ 1E+50 ergs), low ejecta mass (M_ej ~ 3 Msun) SN explosion that occurred in a low-density environment (n_0 ~ 0.01 cm-3).

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

PSR J1802-2124 is a 12.6-ms pulsar in a 16.8-hour binary orbit with a relatively massive white dwarf (WD) companion. These properties make it a member of the intermediate-mass class of binary pulsar (IMBP) systems. We have been timing this pulsar since its discovery in 2002. Concentrated observations at the Green Bank Telescope, augmented with data from the Parkes and Nancay observatories, have allowed us to determine the general relativistic Shapiro delay. This has yielded pulsar and white dwarf mass measurements of 1.24(11) and 0.78(4) solar masses (68% confidence), respectively. The low mass of the pulsar, the high mass of the WD companion, the short orbital period, and the pulsar spin period may be explained by the system having gone through a common-envelope phase in its evolution. We argue that selection effects may contribute to the relatively small number of known IMBPs.

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.

Glitches detected in pulsar timing observations at the Nanshan radio telescope of Urumqi Observatory between 2002 July and 2008 December are presented. In total, 29 glitches were detected in 19 young pulsars, with this being the first detection of a glitch in 12 of these pulsars. Fractional glitch amplitudes range from a few parts in 10^{-11} to 3.9 x 10^{-6}. Three "slow" glitches are identified in PSRs J0631+1036, B1822-09 and B1907+10. Post-glitch recoveries differ greatly from pulsar to pulsar and for different glitches in the same pulsar. Most large glitches show some evidence for exponential post-glitch recovery on timescales of 100 -- 1000 days, but in some cases, e.g., PSR B1758-23, there is little or no recovery. Beside exponential recoveries, permanent increases in slowdown rate are seen for the two large glitches in PSRs B1800-21 and B1823-13. These and several other pulsars also show a linear increase in nudot following the partial exponential recovery, which is similar to the Vela pulsar post-glitch behaviour. Analysis of the whole sample of known glitches show that fractional glitch amplitudes are correlated with characteristic age with a peak at about 10^5 years, but there is a spread of two or three orders of magnitude at all ages. Glitch activity is positively correlated with spin-down rate, again with a wide spread of values.

We present extensive observations of the radio emission from the remnant of SN 1987A made with the Australia Telescope Compact Array (ATCA), since the first detection of the remnant in 1990. The radio emission has evolved in time providing unique information on the interaction of the supernova shock with the circumstellar medium. We particularly focus on the monitoring observations at 1.4, 2.4, 4.8 and 8.6 GHz, which have been made at intervals of 4-6 weeks. The flux density data show that the remnant brightness is now increasing exponentially, while the radio spectrum is flattening. The current spectral index value of -0.68 represents an 18+/-3% increase over the last 8 years. The exponential trend in the flux is also found in the ATCA imaging observations at 9 GHz, which have been made since 1992, approximately twice a year, as well as in the 843 MHz data set from the Molonglo Observatory Synthesis Telescope from 1987 to March 2007. Comparisons with data at different wavelengths (X-ray, H\alpha) are made. The rich data set that has been assembled in the last 22 years forms a basis for a better understanding of the evolution of the supernova remnant.

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.

We report observations of the remnant of Supernova 1987A with the High Resolution Camera (HRC) onboard the Chandra X-ray Observatory. A direct image from the HRC resolves the annular structure of the X-ray remnant, confirming the morphology previously inferred by deconvolution of lower resolution data from the Advanced CCD Imaging Spectrometer. Detailed spatial modeling shows that the a thin ring plus a thin shell gives statistically the best description of the overall remnant structure, and suggests an outer radius 0.96" +/- 0.05" +/- 0.03" for the X-ray-emitting region, with the two uncertainties corresponding to the statistical and systematic errors, respectively. This is very similar to the radius determined by a similar modeling technique for the radio shell at a comparable epoch, in contrast to previous claims that the remnant is 10-15% smaller at X-rays than in the radio band. The HRC observations put a flux limit of 0.010 cts/s (99% confidence level, 0.08-10 keV range) on any compact source at the remnant center. Assuming the same foreground neutral hydrogen column density as towards the remnant, this allows us to rule out an unobscured neutron star with surface temperature T^\infty > 2.5MK observed at infinity, a bright pulsar wind nebula or a magnetar.

We report here on two years of timing of 168 pulsars using the Parkes radio telescope. The vast majority of these pulsars have spin-down luminosities in excess of 10^34 erg/s and are prime target candidates to be detected in gamma-rays by the Fermi Gamma-Ray Space Telescope. We provide the ephemerides for the ten pulsars being timed at Parkes which have been detected by Fermi in its first year of operation. These ephemerides, in conjunction with the publicly available photon list, can be used to generate gamma-ray profiles from the Fermi archive. We will make the ephemerides of any pulsars of interest available to the community upon request. In addition to the timing ephemerides, we present the parameters for 14 glitches which have occurred in 13 pulsars, seven of which have no previously known glitch history. The Parkes timing programme, in conjunction with Fermi observations, is expected to continue for at least the next four years.

The aftermath of supernova (SN) 1987A continues to provide spectacular insights into the interaction between a SN blastwave and its circumstellar en- vironment. We here present 36 GHz observations from the Australia Telescope Compact Array of the radio remnant of SN 1987A. These new images, taken in 2008 Apr and 2008 Oct, substantially extend the frequency range of an ongo- ing monitoring and imaging program conducted between 1.4 and 20 GHz. Our 36.2 GHz images have a diffraction-limited angular resolution of 0.3-0.4 arcseconds, which covers the gap between high resolution, low dynamic range VLBI images of the remnant and low resolution, high dynamic range images at frequencies between 1 and 20 GHz. The radio morphology of the remnant at 36 GHz is an elliptical ring with enhanced emission on the eastern and western sides, similar to that seen previously at lower frequencies. Model fits to the data in the Fourier domain show that the emitting region is consistent with a thick inclined torus of mean radius 0.85 arcsec, and a 2008 Oct flux density of 27 +/- 6 mJy at 36.2 GHz. The spectral index for the remnant at this epoch, determined between 1.4 GHz and 36.2 GHz, is -0.83. There is tentative evidence for an unresolved central source with flatter spectral index.

We present radio timing measurements of six rotating radio transient (RRAT) sources discovered in the Parkes Multibeam Pulsar Survey. These provide four new phase-connected timing solutions and two updated ones, making a total of seven of the original 11 reported RRATs now with high-precision rotational and astrometric parameters. Three of these seven RRATs have magnetic fields greater than 10^13 G, with spin-down properties similar to those of the magnetars and X-ray detected isolated neutron stars. Another two of these RRATs have long periods and large characteristic ages, and lie near the `death-line' for radio pulsar emission. The remaining two RRATs with timing solutions have properties typical of the bulk of the pulsar population. The new solutions offer insights into what might be responsible for the unusual emission properties. We demonstrate that the RRATs have significantly longer periods and higher magnetic fields than normal radio pulsars, and find no correlation with other spin-down parameters. These solutions also provide precise positions, which will facilitate follow-up studies at high energies, crucial for relating these sources with other neutron star populations.

With a collecting area of 70 000 m^2, the Five hundred metre Aperture Spherical Telescope (FAST) will allow for great advances in pulsar astronomy. We have performed simulations to estimate the number of previously unknown pulsars FAST will find with its 19-beam or possibly 100-beam receivers for different survey strategies. With the 19-beam receiver, a total of 5200 previously unknown pulsars could be discovered in the Galactic plane, including about 460 millisecond pulsars (MSPs). Such a survey would take just over 200 days with eight hours survey time per day. We also estimate that, with about 80 six-hour days, a survey of M31 and M33 could yield 50--100 extra-Galactic pulsars. A 19-beam receiver would produce just under 500 MB of data per second and requires about 9 tera-ops to perform the major part of a real time analysis. We also simulate the logistics of high-precision timing of MSPs with FAST. Timing of the 50 brightest MSPs to a signal-to-noise of 500 would take about 24 hours per epoch.

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.

Modern pulsar surveys produce many millions of candidate pulsars, far more than can be individually inspected. Traditional methods for filtering these candidates, based upon the signal-to-noise ratio of the detection, cannot easily distinguish between interference signals and pulsars. We have developed a new method of scoring candidates using a series of heuristics which test for pulsar-like properties of the signal. This significantly increases the sensitivity to weak pulsars and pulsars with periods close to interference signals. By applying this and other techniques for ranking candidates from a previous processing of the Parkes Multi-beam Pulsar Survey, 28 previously unknown pulsars have been discovered. These include an eccentric binary system and a young pulsar which is spatially coincident with a known supernova remnant.

Analysis of pulsar timing data-sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the {\em Fermi Gamma-ray Space Telescope} (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from PSR B1937+21 recorded at Nan\c{c}ay, and using X-ray data on PSR J0218+4232 from XMM-Newton. We demonstrate accuracy of the pulsar phase calculations at the microsecond level. Data Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ .

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.

We report the discovery of PSR J1753-2240 in the Parkes Multibeam Pulsar Survey database. This 95-ms pulsar is in an eccentric binary system with a 13.6-day orbital period. Period derivative measurements imply a characteristic age in excess of 1 Gyr, suggesting that the pulsar has undergone an episode of accretion-induced spin-up. The eccentricity and spin period are indicative of the companion being a second neutron star, so that the system is similar to that of PSR J1811-1736, although other companion types cannot be ruled out at this time. The companion mass is constrained by geometry to lie above 0.48 solar masses, although long-term timing observations will give additional constraints. If the companion is a white dwarf or main sequence star, optical observations may yield a direct detection of the companion. If the system is indeed one of the few known double neutron star systems, it would lie significantly far from the recently proposed spin-period/eccentricity relationship.

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

AGILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to high-energy astrophysics launched in 2007 April. Its 1 microsecond absolute time tagging capability coupled with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the CGRO/EGRET heritage. In this paper we present the first AGILE timing results on the known gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a large number of gamma-ray photons from these pulsars (about 10,000 pulsed counts for Vela) in only few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves.

We present detailed Fourier modeling of the radio remnant of Supernova 1987A, using high-resolution 9 GHz and 18 GHz data taken with the Australia Telescope Compact Array over the period 1992 to 2008. We develop a parameterized three-dimensional torus model for the expanding radio shell, in which the emission is confined to an inclined equatorial belt; our model also incorporates both a correction for light travel-time effects and an overall east-west gradient in the radio emissivity. By deriving an analytic expression for the two-dimensional Fourier transform of the projected three-dimensional brightness distribution, we can fit our spatial model directly to the interferometric visibility data. This provides robust estimates to the radio morphology at each epoch. The best-fit results suggest a constant remnant expansion at 4000 +/- 400 km/s over the 16-year period covered by the observations. The model fits also indicate substantial mid-latitude emission, extending to 40 degree on either side of the equatorial plane. This likely corresponds to the extra-planar structure seen in H$\alpha$ and Ly$\alpha$ emission from the supernova reverse shock, and broadly supports hydrodynamic models in which the complex circumstellar environment was produced by a progression of interacting winds from the progenitor. Our model quantifies the clear asymmetry seen in the radio images: we find that the eastern half of the radio remnant is consistently ~40 brighter than the western half at all epochs, which may result from an asymmetry in the ejecta distribution between these two hemispheres.

Quasi-continuous observations of PSR B03239+54 over 20 days using the Nanshan 25-m telescope at 1540 MHz have been used to study the effects of refractive scintillation on the pulsar flux density and diffractive scintillation properties. Dynamic spectra were obtained from datasets of 90 min duration and diffractive parameters derived from a two-dimensional auto-correlation analysis. Secondary spectra were also computed but these showed no significant evidence for arc structure. Cross correlations between variations in the derived parameters were much lower than predicted by thin screen models and in one case was of opposite sign to the prediction. Observed modulation indices were larger than predicted by thin screen models with a Kolmogorov fluctuation spectrum. Structure functions were computed for the flux density, diffractive timescale and decorrelation bandwidth. These indicated a refractive timescale of $8\pm 2$ h, much shorter than predicted by the thin screen model. The measured structure-function slope of $0.4\pm 0.2$ is also inconsistent with scattering by a single thin screen for which a slope of 2.0 is expected. All observations are consistent with scattering by an extended medium having a Kolmogorov fluctuation spectrum which is concentrated towards the pulsar. This interpretation is also consistent with recent observations of multiple diffuse scintillation arcs for this pulsar.

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 double pulsar system PSR J0737-3039A/B is a highly relativistic double neutron star (DNS) binary, with a 2.4-hour orbital period. The low mass of the second-formed NS, as well the low system eccentricity and proper motion, point to a different evolutionary scenario compared to other known DNS systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to have a low misalignment between its spin and orbital angular momentum axes, with a 68.3% upper limit of 6.1 degrees, assuming emission from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

Six glitches have been recently observed in the rotational frequency of the young pulsar PSR B1737-30 (J1740-3015) using the 25-m Nanshan telescope of Urumqi Observatory. With a total of 20 glitches in 20 years, it is one of the most frequently glitching pulsars of the about 1750 known pulsars. Glitch amplitudes are very variable with fractional increases in rotation rate ranging from 10^{-9} to 10^{-6}. Inter-glitch intervals are also very variable, but no relationship is observed between interval and the size of the preceding glitch. There is a persistent increase in |\dot\nu|, opposite in sign to that expected from slowdown with a positive braking index, which may result from changes in the effective magnetic dipole moment of the star during the glitch.

Detection and study of gravitational waves from astrophysical sources is a major goal of current astrophysics. Ground-based laser-interferometer systems such as LIGO and VIRGO are sensitive to gravitational waves with frequencies of order 100 Hz, whereas space-based systems such as LISA are sensitive in the millihertz regime. Precise timing observations of a sample of millisecond pulsars widely distributed on the sky have the potential to detect gravitational waves at nanohertz frequencies. Potential sources of such waves include binary super-massive black holes in the cores of galaxies, relic radiation from the inflationary era and oscillations of cosmic strings. The Parkes Pulsar Timing Array (PPTA) is an implementation of such a system in which 20 millisecond pulsars have been observed using the Parkes radio telescope at three frequencies at intervals of two -- three weeks for more than two years. Analysis of these data has been used to limit the gravitational wave background in our Galaxy and to constrain some models for its generation. The data have also been used to investigate fluctuations in the interstellar and Solar-wind electron density and have the potential to investigate the stability of terrestrial time standards and the accuracy of solar-system ephemerides.

Variations in the solar wind density introduce variable delays into pulsar timing observations. Current pulsar timing analysis programs only implement simple models of the solar wind, which not only limit the timing accuracy, but can also affect measurements of pulsar rotational, astrometric and orbital parameters. We describe a new model of the solar wind electron density content which uses observations from the Wilcox Solar Observatory of the solar magnetic field. We have implemented this model into the tempo2 pulsar timing package. We show that this model is more accurate than previous models and that these corrections are necessary for high precision pulsar timing applications.

SN 1987A offered a unique opportunity to detect a pulsar at the very beginning of its life and to study its early evolution. Despite many searches at radio and optical wavelengths, no pulsar has yet been detected. Details of a recent search using the Parkes radio telescope are given. Limits on the X-ray, optical and radio luminosity of a point source at the centre of SN 1987A place limits on the properties of a central neutron star. However, neither these nor the pulsar limits preclude the presence of a relatively slowly rotating neutron star (P >~ 100 ms) with a moderate surface dipole magnetic field in SN 1987A. Galactic studies suggest that a significant fraction of pulsars are born with parameters in this range. In view of this, continued searches for a pulsar in SN 1987A are certainly justified.

We investigate the age constraints that can be placed on the double pulsar system using models for the spin-down of the first-born 22.7-ms pulsar A and the 2.77-s pulsar B with characteristic ages of 210 and 50 Myr respectively. Standard models assuming dipolar spin-down of both pulsars suggest that the time since the formation of B is ~50 Myr, i.e. close to B's characteristic age. However, adopting models which account for the impact of A's relativistic wind on B's spin-down we find that the formation of B took place either 80 or 180 Myr ago, depending the interaction mechanism. Formation 80 Myr ago, closer to B's characteristic age, would result in the contribution from J0737-3039 to the inferred coalescence rates for double neutron star binaries increasing by 40%. The 180 Myr age is closer to A's characteristic age and would be consistent with the most recent estimates of the coalescence rate. The new age constraints do not significantly impact recent estimates of the kick velocity, tilt angle between pre and post-supernova orbital planes or pre-supernova mass of B's progenitor.

Supernova 1987A in the Large Magellanic Cloud provides a spectacularly detailed view of the aftermath of a core-collapse explosion. The supernova ejecta initially coasted outward at more than 10% of the speed of light, but in 1990 were observed to decelerate rapidly as they began to encounter dense circumstellar material expelled by the progenitor star. The resulting shock has subsequently produced steadily brightening radio synchrotron emission, which is resolved by the Australia Telescope Compact Array (ATCA) into an expanding limb-brightened shell. Here we present 15 years of ATCA imaging of Supernova 1987A, at an effective angular resolution of 0.4 arcsec. We find that the radio remnant has accelerated in its expansion over this period, from approx 3600 km/s in 1992 to approx 5200 km/s at the end of 2006. The published diameters of the evolving X-ray shell have been ~15% smaller than the corresponding radio values, but a simultaneous Fourier analysis of both radio and X-ray data eliminates this discrepancy, and yields a current diameter for the shell in both wave-bands of approx 1.7 arcsec. An asymmetric brightness distribution is seen in radio images at all ATCA epochs: the eastern and western rims have higher fluxes than the northern and southern regions, indicating that most of the radio emission comes from the equatorial plane of the system, where the progenitor star's circumstellar wind is thought to be densest. The eastern lobe is brighter than and further from the supernova site than the western lobe, suggesting an additional asymmetry in the initial distribution of supernova ejecta.

With its relatively long observation time per pointing, the Parkes multibeam survey was effective in detecting nulling pulsars. We have made 2-hour observations of 23 pulsars which showed evidence for pulse nulling or mode changing in the survey data. Because of the low flux density of these pulsars, in most cases averaging times of between 10 s and 60 s were necessary and so this analysis is insensitive to very short nulls. Seven of the pulsars had null fractions of more than 40% with the largest having a lower limit of 95%. Mode changes were observed in six pulsars with clear relationships between nulling and mode changing in some cases. Combined with earlier results, the data suggest that large null fractions are more related to large characteristic age than to long pulse period. The observations suggest that nulling and mode changing are different manifestations of the same phenomenon.

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.

Pulsar braking indices offer insight into the physics that underlies pulsar spin-down. Only five braking indices have been measured via phase-coherent timing; all measured values are less than 3, the value expected from magnetic dipole radiation. Here we present new measurements for three of the five pulsar braking indices, obtained with phase-coherent timing for PSRs J1846-0258 (n=2.65+/-0.01), B1509-58 (n=2.839+/-0.001) and B0540-69 (n=2.140+/-0.009). We discuss the implications of these results and possible physical explanations for them.

The Parkes Multibeam Survey led to the identification of a number of long-period radio pulsars with magnetic field well above the 'quantum critical field' of ~ 4.4x10^13 G (HBRPs). The HBRPs have similar spin parameters to magnetars, but their emission properties are different, and contradict those theories that predict that radio emission should be suppressed above this critical field. Our observations support the suggestion that initial neutron star spin periods depend on their magnetic fields; in particular, there is a tendency for high-field systems to be born as slow rotators. The aim of this project is to understand the emission properties of HBRPs, using multiple radio frequencies and high time resolution data. One specific objective is to identify HBRPs radio emission characteristics that are different from those of normal pulsars.

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

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

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

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

Tempo2 is a new software package for the analysis of pulsar pulse times of arrival. In this paper we describe in detail the timing model used by tempo2, and discuss limitations on the attainable precision. In addition to the intrinsic slow-down behaviour of the pulsar, tempo2 accounts for the effects of a binary orbital motion, the secular motion of the pulsar or binary system, interstellar, Solar system and ionospheric dispersion, observatory motion (including Earth rotation, precession, nutation, polar motion and orbital motion), tropospheric propagation delay, and gravitational time dilation due to binary companions and Solar system bodies. We believe the timing model is accurate in its description of predictable systematic timing effects to better than one nanosecond, except in the case of relativistic binary systems where further theoretical development is needed. The largest remaining sources of potential error are measurement error, interstellar scattering, Solar system ephemeris errors, atomic clock instability and gravitational waves.

[ABRIDGED] We present the discovery and follow-up observations of 142 pulsars found in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. These new discoveries bring the total number of pulsars found by the survey to 742. In addition to tabulating spin and astrometric parameters, along with pulse width and flux density information, we present orbital characteristics for 13 binary pulsars which form part of the new sample. Combining these results from another recent Parkes multibeam survey at high Galactic latitudes, we have a sample of 1008 normal pulsars which we use to carry out a determination of their Galactic distribution and birth rate. We infer a total Galactic population of 30000 +/- 1100 potentially detectable pulsars (i.e. those beaming towards us) having 1.4-GHz luminosities above 0.1 mJy kpc squared. Using a pulsar current analysis, we derive the birth rate of this population to be 1.4 +/- 0.2 pulsars per century. An important conclusion from our work is that the inferred radial density function of pulsars depends strongly on the assumed distribution of free electrons in the Galaxy. As a result, any analyses using the most recent electron model of Cordes & Lazio predict a dearth of pulsars in the inner Galaxy. We show that this model can also bias the inferred pulsar scale height with respect to the Galactic plane. Combining our results with other Parkes multibeam surveys we find that the population is best described by an exponential distribution with a scale height of 330 pc.

We report the discovery of two highly dispersed pulsars in the direction of the Galactic Centre made during a survey at 3.1 GHz with the Parkes radio telescope. Both PSRs J1745-2912 and J1746-2856 have an angular separation from the Galactic Centre of less than 0.3 degrees and dispersion measures in excess of 1100 cm-3pc, placing them in the top 10 pulsars when ranked on this value. The frequency dependence of the scatter-broadening in PSR J1746-2856 is much shallower than expected from simple theory. We believe it likely that the pulsars are located between 150 and 500 pc from the Galactic Centre on the near side, and are part of an excess population of neutron stars associated with the Centre itself. A second survey made at 8.4 GHz did not detect any pulsars. This implies either that there are not many bright, long-period pulsars at the Galactic Centre or that the scattering is more severe at high frequencies than current models would suggest.

A systematic survey of the Large and Small Magellanic Clouds for radio pulsars using the Parkes radio telescope and the 20-cm multibeam receiver has resulted in the discovery of 14 pulsars and the redetection of five of the eight previously known spin-powered pulsars believed to lie in the Magellanic Clouds. Of the 14 new discoveries, 12 are believed to lie within Clouds, three in the Small Cloud and nine in the Large Cloud, bringing the total number of known spin-powered pulsars in the Clouds to 20. Averaged over all positions within the survey area, the survey had a limiting flux density of about 0.12 mJy. Observed dispersion measures suggest that the mean free electron density in the Magellanic Clouds is similar to that in the disk of our Galaxy. The observed radio luminosities have little or no dependence on pulsar period or characteristic age and the differential luminosity function is consistent with a power-law slope of -1 as is observed for Galactic pulsars.

The recently started Arecibo L-band Feed Array (ALFA) pulsar survey aims to find ~1000 new pulsars. Due to its high time and frequency resolution the survey is especially sensitive to millisecond pulsars, which have the potential to test gravitational theories, detect gravitational waves and probe the neutron-star equation of state. Here we report the results of our preliminary analysis: in the first months we have discovered 21 new pulsars. One of these, PSR J1906+0746, is a young 144-ms pulsar in a highly relativistic 3.98-hr low-eccentricity orbit. The 2.61 +- 0.02 solar-mass system is expected to coalesce in ~300 Myr and contributes significantly to the computed cosmic inspiral rate of compact binary systems.

The Parkes Multibeam Survey led to the identification ofa number of long-period radio pulsars with magnetic field well above the quantum critical field of 4.4 times 10 to 13 G (defined as high magnetic field radio pulsars - HBRPs). Traditional pulsar emission theories postulate that radio emission is suppressed above this critical field. The aim of thisproject is to understand emission properties of HBRPs.

Given sufficient sensitivity, pulsar timing observations can make a direct detection of gravitational waves passing over the Earth. Pulsar timing is most sensitive to gravitational waves with frequencies in the nanoHertz region, with the most likely astronomical sources being binary super-massive black holes in galaxy cores. The Parkes Pulsar Timing Array project uses the Parkes 64-m radio telescope to make precision timing observations of a sample of about 20 millisecond pulsars with a principal goal of making a direct detection of gravitational waves. Observations commenced about one year ago and so far sub-microsecond timing residuals have been achieved for more than half of these pulsars. New receiver and software systems are being developed with the aim of reducing these residuals to the level believed necessary for a positive detection of gravitational waves.

Contemporary pulsar timing experiments have reached a sensitivity level where systematic errors introduced by existing analysis procedures are limiting the achievable science. We have developed tempo2, a new pulsar timing package that contains propagation and other relevant effects implemented at the 1ns level of precision (a factor of ~100 more precise than previously obtainable). In contrast with earlier timing packages, tempo2 is compliant with the general relativistic framework of the IAU 1991 and 2000 resolutions and hence uses the International Celestial Reference System, Barycentric Coordinate Time and up-to-date precession, nutation and polar motion models. Tempo2 provides a generic and extensible set of tools to aid in the analysis and visualisation of pulsar timing data. We provide an overview of the timing model, its accuracy and differences relative to earlier work. We also present a new scheme for predictive use of the timing model that removes existing processing artifacts by properly modelling the frequency dependence of pulse phase.

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.

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

The `radio sky' is relatively unexplored for transient signals, although the potential of radio-transient searches is high, as demonstrated recently by the discovery of a previously unknown type of source which varies on timescales of minutes to hours. Here we report a new large-scale search for radio sources varying on much shorter timescales. This has revealed 11 objects characterized by single, dispersed bursts having durations between 2 and 30 ms. The average time intervals between bursts range from 4 minutes to 3 hours, with radio emission typically detectable for < 1 s per day. From an analysis of the burst arrival times, we have identified periodicities in the range 0.4 - 7 s for ten of the 11 sources, suggesting a rotating neutron star origin. Despite the small number of sources presently detected, their ephemeral nature implies a total Galactic population which significantly exceeds that of the regularly pulsing radio pulsars. Five of the ten sources have periods greater than 4 s, and period derivatives have been measured for three of the sources, with one having a very high inferred magnetic field of 5e13 G, suggesting that this new population is related to other classes of isolated neutron stars observed at X-ray and gamma-ray wavelengths.

We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4-GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, Green Bank, and Parkes, we have measured the spin-down and binary parameters of the pulsar and its basic spectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance 7.57+/-0.03 deg/yr, we infer a total system mass of 2.61+/-0.02 Msun. While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birth rate of binaries similar to PSR J1906+0746 is 60/Myr. This implies that PSR J1906+0746 will make a significant contribution to the computed cosmic inspiral rate of compact binary systems.

We have discovered the pulsar associated with the supernova remnant G21.5-0.9. PSR J1833-1034, with spin period P = 61.8 ms and dispersion measure 169 pc/cc, is very faint, with pulse-averaged flux density of approximately 70 microJy at a frequency of 1.4 GHz, and was first detected in a deep search with the Parkes telescope. Subsequent observations with Parkes and the Green Bank Telescope have confirmed this detection, and yield a period derivative dP/dt = 2.02e-13. These spin parameters imply a characteristic age tau_c = 4.8 kyr and a spin-down luminosity dE/dt = 3.3e37 ergs/s, the latter value exceeded only by the Crab pulsar among the rotation-powered pulsars known in our Galaxy. The pulsar has an unusually steep radio spectrum in the 0.8-2.0 GHz range, with power law index ~3.0, and a narrow single-peaked pulse profile with full-width at half maximum of 0.04P. We have analyzed 350 ks of archival Chandra X-ray Observatory High Resolution Camera (HRC) data, and find a point-like source of luminosity approximately 3e-5dE/dt, offset from the center of an elliptical region of size ~7"x5" and luminosity approximately 1e-3dE/dt within which likely lies the pulsar wind termination shock. We have searched for X-ray pulsations in a 30 ks HRC observation without success, deriving a pulsed fraction upper limit for a sinusoidal pulse shape of about 70% of the pulsar flux. We revisit the distance to G21.5-0.9 based on HI and CO observations, arguing that it is 4.7+-0.4 kpc. We use existing X-ray and radio observations of the pulsar wind nebula, along with the measured properties of its engine and a recent detection of the supernova remnant shell, to argue that G21.5-0.9 and PSR J1833-1034 are much younger than tau_c, and likely their true age is ~<1000 yr. In that case, the initial spin period of the pulsar was >~55 ms.

We have measured the positions of 74 pulsars from regular timing observations using the Nanshan radio telescope at Urumqi Observatory between 2000 January and 2004 August (MJD 51500 -- 53240). Proper motions were determined for these pulsars by comparing their current positions with positions given in pulsar catalogues. We compare our results to earlier measurements in the literature and show that, in general, the values agree. New or improved proper motions are obtained for 16 pulsars. The effect of period fluctuations and other timing noise on the determination of pulsar positions is investigated. For our sample, the mean and rms transverse velocities are 443 and 224 km/s respectively, agreeing with previous work even though we determine distances using the new NE2001 electron density model.

We report the discovery of three binary millisecond pulsars during the Parkes Multibeam Pulsar Survey of the Galactic Plane. The objects are highly recycled and are in orbits of many tens of days about low-mass white-dwarf companions. The eccentricity of one object, PSR J1853+1303, is more than an order of magnitude lower than predicted by the theory of convective fluctuations during tidal circularization. We demonstrate that, under the assumption that the systems are randomly oriented, current theoretical models of the core-mass--orbital-period relation for the progenitors of these systems likely overestimate the white-dwarf masses, strengthening previous concerns about the match of these models to the data. The new objects allow us to update the limits on violation of relativistic equivalence principles to 95% confidence upper limits of 5.6 x 10^-3 for the Strong Equivalence Principle parameter Delta and 4.0 x 10^-20 for the Lorentz-invariance/momentum-conservation parameter alpha_3.

Observations of the radio remnant of Supernova 1987A using the Australia Telescope Compact Array in the 12-mm band on 2003 July 31 (day 6002.7 after the explosion) give the first fully resolved radio image of the supernova remnant. The diffraction-limited image has a resolution of about 0.45 arcsec, a factor of two better than that of the 3-cm images previously obtained. There is excellent agreement between the 12-mm image and a contemporaneous super-resolved 3-cm image. Super-resolution of the 12-mm image gives a further factor of two improvement in resolution, to 0.25 arcsec, albeit with limited dynamic range. While the spatial distributions of the radio and X-ray emission are broadly similar, there are significant differences in detail with no correspondence in the regions of brightest emission. The 12-mm image is well modelled by a thick equatorial ring inclined at 43 deg to the line of sight. This, together with the common east-west asymmetry and the relatively steady increase in the radio flux density, suggests that the reverse shock is the main site for generation of the radio emission.

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

The direct detection of gravitational waves is a major goal of current astrophysics. We provide details of a new method for detecting a stochastic background of gravitational waves using pulsar timing data. Our results show that regular timing observations of 40 pulsars each with a timing accuracy of 100 nano-seconds will be able to make a direct detection of the predicted stochastic background from coalescing black holes within five years. With an improved pre-whitening algorithm, or if the background is at the upper end of the predicted range, a significant detection should be possible with only 20 pulsars.

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

We present polarisation profiles for 48 southern pulsars observed with the new 10-cm receiver at the Parkes telescope. We have exploited the low system temperature and high bandwidth of the receiver to obtain profiles which have good signal to noise for most of our sample at this relatively high frequency. Although, as expected, a number of profiles are less linearly polarised at 3.1 GHz than at lower frequencies, we identify some pulsars and particular components of profiles in other pulsars which have increased linear polarisation at this frequency. We discuss the dependence of linear polarisation with frequency in the context of a model in which emission consists of the superposition of two, orthogonally polarised modes. We show that a simple model, in which the orthogonal modes have different spectral indices, can explain many of the observed properties of the frequency evolution of both the linear polarisation and the total power, such as the high degree of linear polarisation seen at all frequencies in some high spin-down, young pulsars. Nearly all the position angle profiles show deviations from the rotating vector model; this appears to be a general feature of high-frequency polarisation observations.

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

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

We report here on extensive radio observations of the binary system containing PSR B1259-63 and the Be star SS 2883, made around the time of the 2004 periastron. This is the fourth periastron to have been observed in detail. As in previous observations, changes in the pulsar's dispersion and rotation measures are detected over a period spanning 200 days. We show that the eclipse of the pulsed emission lasts from 16 days prior to periastron to 15 days after periastron and is consistent from one periastron to the next. We demonstrate that the timing solution proposed by Wang et al. (2004) provides a good fit through the 2004 periastron. The light curve of the transient unpulsed radio emission is broadly similar from one periastron to the next. For this periastron, however, the light curve is strongly peaked post-periastron with rather little enhancement prior to periastron, in contrast to the 2000 periastron where the peak flux densities were more equal. These observations remain consistent with the interpretation that the pulsar passes through the dense circumstellar disk of the Be star just before and just after periastron. The observed differences from one periastron to the next can be ascribed to variations in the local disk density and magnetic field structure at the time the pulsar enters the disk.

From 2001 January to 2002 June, we monitored PSRs B0329+54, B0823+26, B1929+10, B2020+28 and B2021+51 using the Nanshan 25-m radio telescope of Urumqi Observatory to study their diffractive interstellar scintillation (DISS). The average interval between observations was about 9 days and the observation duration ranged between 2 and 6 hours depending on the pulsar. Wide variations in the DISS parameters were observed over the 18-month data span. Despite this, the average scintillation velocities are in excellent agreement with the proper motion velocities. The average two-dimensional autocorrelation function for PSR B0329+54 is well described by a thin-screen Kolmogorov model, at least along the time and frequency axes. Observed modulation indices for the DISS time scale and bandwidth and the pulsar flux density are greater than values predicted for a Kolmogorov spectrum of electron density fluctuations. Correlated variations over times long compared to the nominal refractive scintillation time are observed, suggesting that larger-scale density fluctuations are important. For these pulsars, the scintillation bandwidth as a function of frequency has a power-law index (~3.6) much less than expected for Kolmogorov turbulence (~4.4). Sloping fringes are commonly observed in the dynamic spectra, especially for PSR B0329+54. The detected range of fringe slopes are limited by our observing resolution. Our observations are sensitive to larger-scale fringes and hence smaller refractive angles, corresponding to the central part of the scattering disk.

We have compiled a new and complete catalog of the main properties of the 1509 pulsars for which published information currently exists. The catalog includes all spin-powered pulsars as well as anomalous X-ray pulsars and soft gamma-ray repeaters showing coherent pulsed emission, but excludes accretion-powered systems. References are given for all data listed. We have also developed a new web interface for accessing and displaying either tabular or plotted data with the option of selecting pulsars to be displayed via logical conditions on parameter expressions. The web interface has an ``expert'' mode giving access to a wider range of parameters and allowing the use of custom databases. For users with locally installed software and database on unix or linux systems, the catalog may be accessed from a command-line interface. C-language functions to access specified parameters are also available. The catalog is updated from time to time to include new information.

Pulsar detection and timing experiments are applications where adaptive filters seem eminently suitable tools for radio-frequency-interference (RFI) mitigation. We describe a novel variant which works well in field trials of pulsar observations centred on an observing frequency of 675 MHz, a bandwidth of 64 MHz and with 2-bit sampling. Adaptive filters have generally received bad press for RFI mitigation in radio astronomical observations with their most serious drawback being a spectral echo of the RFI embedded in the filtered signals. Pulsar observations are intrinsically less sensitive to this as they operate in the (pulsar period) time domain. The field trials have allowed us to identify those issues which limit the effectiveness of the adaptive filter. We conclude that adaptive filters can significantly improve pulsar observations in the presence of RFI.

We report the discovery during the Parkes Multibeam Pulsar Survey of PSR J1756-2251, a 28.5 ms pulsar in a relativistic binary system. Subsequent timing observations showed the pulsar to have an orbital period of 7.67 hrs and an eccentricity of 0.18. They also revealed a significant advance of periastron, 2.585+/-0.002 deg./yr. Assuming this is entirely due to general relativity implies a total system mass (pulsar plus companion) of 2.574+/-0.003 solar mass. This mass and the significant orbital eccentricity suggest that this is a double neutron star system. Measurement of the gravitational redshift, gamma, and an evaluation of the Shapiro delay shape, s, indicate a low companion mass of <1.25 solar mass. The expected coalescence time due to emission of gravitational waves is only ~1.7 Gyr substantially less than a Hubble time. We note an apparent correlation between spin period and eccentricity for normally evolving double neutron star systems.

We have investigated the eclipse of the 23-ms pulsar J0737-3039A by its 2.8-s companion PSR J0737-3039B in the recently discovered double pulsar system using data taken with the Green Bank Telescope at 820 MHz. We find that the pulsed flux density at eclipse is strongly modulated with half the periodicity of the 2.8-s pulsar. The eclipse occurs earlier and is deeper at those rotational phases of B when its magnetic axis is aligned with the line of sight than at phases when its magnetic axis is at right angles to the line of sight. This is consistent with the eclipse of A being due to synchrotron absorption by the shock-heated plasma surrounding B, the asymmetry arising from the higher plasma densities expected in the B magnetosphere's polar cusps.

We present an updated timing solution for the young, energetic pulsar PSR B1509-58 based on 21.3 years of radio timing data and 7.6 years of X-ray timing data. No glitches have occurred in this time span, in contrast to other well-studied young pulsars, which show frequent glitches. We report a measurement of the third frequency derivative of (-1.28+/-0.21)x10^(-31) s^(-4). This value is 1.65 standard deviations from, i.e. consistent with, that predicted by the simple constant magnetic dipole model of pulsar spin-down. We measured the braking index to be n=2.839+/-0.003 and show that it varies by 1.5% over 21.3 yr due to contamination from timing noise. Results of a low-resolution power spectral analysis of the significant noise apparent in the data yield a spectral index of alpha=-4.6+/-1.0 for the red noise component.

We present high-quality observations of PSR B0826-34 at 1374 MHz. The emission from this pulsar exhibits strong bursts of pulses followed by long periods of `null' pulses. When it is strong, the radiation extends through the whole pulse period. We show for the first time that there is weak emission during the `null' phases, which should therefore be considered to be a different mode rather than a null. During this weak mode the profile is similar to that observed in the strong mode at low radio frequency. Using a phase-tracking method, the pattern of drifting subpulses during the strong mode is seen to be coherent across the whole profile. The drift rate is variable and includes positive and negative values. Thirteen subpulse bands have been directly observed, covering the whole longitude range. The subpulses and their spacings ($P_2$) are wider in one half of the profile than those in the other half. This difference, and the variation of observed $P_2$ within the two regions, can be accounted for if the magnetic pole is inclined to the rotation axis by about 0.5 degrees. These two regions appear to represent radiation from outer and inner cones. The intensity modulation of subpulses in all longitude ranges is related to the magnitude of the drift rate.

We have analyzed single pulses from PSR J0737-3039B, the 2.8-s pulsar in the recently discovered double pulsar system, using data taken with the Green Bank Telescope at 820 and 1400 MHz. We report the detection of features similar to drifting subpulses, detectable over only a fraction of the pulse window, with a fluctuation frequency of 0.196 cycles/period. This is exactly the beat frequency between the periods of the two pulsars. In addition, the drifting features have a separation within a given pulse of 23 ms, equal to the pulse period of A. These features are therefore due to the direct influence of PSR J0737-3039A's 44-Hz electromagnetic radiation on PSR J0737-3039B's magnetosphere. We only detect them over a small range of orbital phases, when the radiation from the recycled pulsar PSR J0737-3039A meets our line of sight to PSR J0737-3039B from the side.

The Parkes Multibeam Pulsar Survey is the most successful survey of the Galactic plane ever performed, finding over 600 pulsars in the initial processing. We report on reprocessing of all 40,000 beams with a number of algorithms, including conventional frequency-domain searches and an acceleration search for fast binary pulsars. The very large volume of results coupled with the need to distinguish new candidates from known pulsars and their many harmonics, often with multiple detections from different search algorithms, necessitated the development of a new graphical selection tool tightly linked to a web-based results database. We discuss and demonstrate the benefits of these software systems which are specifically designed for large survey projects. The results of this processing have been encouraging; we have discovered 128 new pulsars including 11 binary and 15 millisecond pulsars, in addition to those previously found in the survey, we have thus far discovered 737 pulsars. In this paper we discuss the discoveries of PSR J1744-3922, a 172 ms mildly recycled pulsar in a 4.6 hr orbit that exhibits nulling behaviour, not previously observed in recycled or binary objects; PSR J1802-2124, an intermediate mass binary pulsar (IMBP) and PSR J1801-1417, a solitary millisecond pulsar.

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

We report on a 19 ks observation of the pulsar J1016-5857 with the Chandra X-ray Observatory. This "Vela-like" pulsar has rotation period 107 ms, characteristic age 21 kyr, and spin-down power 2.6e36 ergs/s. A relatively bright centrally peaked source of radius ~25" around the radio pulsar position has a spectrum that is well fitted by an absorbed power law with photon index 1.32+-0.25. We regard this as a newly identified pulsar wind nebula that we designate PWN G284.0-1.8. We do not detect the pulsar either as a point X-ray source or as a pulsed source in a 55 ks observation with the Rossi X-ray Timing Explorer. The isotropic PWN luminosity is 3e32 ergs/s in the 2-10 keV range, for a distance of 3 kpc that is consistent with the measured neutral hydrogen column density. The unpulsed flux from the pulsar is less than 30% of the measured PWN flux. The brightest component of the PWN, near the pulsar, shows extended emission of size ~2" that may indicate, by analogy with other young pulsars, the wind termination shock. In the Chandra image we also detect a very faint extended structure ~1' by 2' in size that is highly asymmetric about the pulsar position. This structure is a good match in width and position angle to the tip of a "finger" of radio emission that appears to connect to the nearby supernova remnant G284.3-1.8, but we cannot characterize it further with the available data. There is a variable X-ray point source only 1.5 arcmin from the pulsar, which we identify using optical spectroscopy as an accreting binary.

We present an observation with the Chandra X-ray Observatory of the unusual radio source G359.23-0.82 ("the Mouse"), along with updated radio timing data from the Parkes radio telescope on the coincident young pulsar J1747-2958. We find that G359.23-0.82 is a very luminous X-ray source (L_X [0.5-8.0 keV] = 5e34 ergs/s for a distance of 5 kpc), whose morphology consists of a bright head coincident with PSR J1747-2958, plus a 45"-long narrow tail whose power-law spectrum steepens with distance from the pulsar. We thus confirm that G359.23-0.82 is a bow-shock pulsar wind nebula powered by PSR J1747-2958; the nebular stand-off distance implies that the pulsar is moving with a Mach number of ~60, suggesting a space velocity ~600 km/s through gas of density ~0.3 cm^-3. We combine the theory of ion-dominated pulsar winds with hydrodynamic simulations of pulsar bow shocks to show that a bright elongated X-ray and radio feature extending 10'' behind the pulsar represents the surface of the wind termination shock. The X-ray and radio "trails" seen in other pulsar bow shocks may similarly represent the surface of the termination shock, rather than particles in the postshock flow as is usually argued. The tail of the Mouse contains two components: a relatively broad region seen only at radio wavelengths, and a narrow region seen in both radio and X-rays. We propose that the former represents material flowing from the wind shock ahead of the pulsar's motion, while the latter corresponds to more weakly magnetized material streaming from the backward termination shock. This study represents the first consistent attempt to apply our understanding of "Crab-like" nebulae to the growing group of bow shocks around high-velocity pulsars.

In this paper we report unusual glitches in two young pulsars, PSR J1825-0935 (B1822-09) and PSR J1835-1106. For PSR J1825-0935, a slow glitch characterised by a temporary decrease in the slowdown rate occurred between 2000 December 31 to 2001 December 6. This event resulted in a permanent increase in frequency with fractional size $\Delta\nu/\nu\sim31.2(2)\times10^{-9}$, however little effect remained in slowdown rate. The glitch in PSR J1835-1106 occurred abruptly in November 2001 (MJD 52220\pm3) with $\Delta\nu/\nu\sim14.6(4)\times10^{-9}$ and little or no change in the slow-down rate. A significant change in $\ddot\nu$ apparently occurred at the glitch with $\ddot\nu$ having opposite sign for the pre- and post-glitch data.

We report the discovery of a radio pulsar with a characteristic age of 7300 years, making it one of the 10 apparently youngest Galactic pulsars known. PSR J1357-6429, with a spin period of P = 166 ms and spin-down luminosity of 3.1e36 ergs/s, was detected during the Parkes multibeam survey of the Galactic plane. We have measured a large rotational glitch in this pulsar, with Delta P/P = -2.4e-6, similar in magnitude to those experienced occasionally by the Vela pulsar. At a nominal distance of only ~ 2.5 kpc, based on the measured free electron column density of 127 pc/cc and the electron distribution model of Cordes & Lazio, this may be, after the Crab, the nearest very young pulsar known. The pulsar is located near the radio supernova remnant candidate G309.8-2.6.

The Parkes multibeam pulsar survey has led to the discovery of more than 700 pulsars. In this paper, we provide timing solutions, flux densities and pulse profiles for 180 of these new discoveries. Two pulsars, PSRs J1736-2843 and J1847-0130 have rotational periods P > 6s and are therefore among the slowest rotating radio pulsars known. Conversely, with P = 1.8ms, PSR J1843-1113 has the third shortest period of pulsars currently known. This pulsar and PSR J1905+0400 (P = 3.8ms) are both solitary. We also provide orbital parameters for a new binary system, PSR J1420-5625, which has P = 34ms, an orbital period of 40 days and a minimum companion mass of 0.4 solar masses. The 10 degree-wide strip along the Galactic plane that was surveyed is known to contain 264 radio pulsars that were discovered prior to the multibeam pulsar survey. We have redetected almost all of these pulsars and provide new dispersion measure values and flux densities at 20cm for the redetected pulsars.

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

The first ever double pulsar, discovered by our team a few months ago, consists of two pulsars, one with period of 22 ms and the other with a period of 2.7 s. This binary system with a period of only 2.4-hr provides a truly unique laboratory for relativistic gravitational physics. In this contribution we summarize the published results and look at the prospects of future observations.

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.

As discussed elsewhere in these proceedings, the double pulsar system is a magnificent laboratory for gravitational physics and for studying pulsar magnetospheres. Here we consider the uses of having two clocks in the system in the context of its origin and evolution. We find that the ``standard'' evolutionary scenario involving spin-up of the first-born neutron star in an X-ray binary phase is consistent with the observed parameters. Equality of the spin-down ages of the two pulsars requires that the post-accretion spin period of A most likely lies in the range $16 \lapp P_{\rm 0,A} \lapp 21$ ms. The likely age of the system is 30-70 Myr.

Interstellar magnetic fields exist over a broad range of spatial scales, extending from the large Galactic scales ($\sim 10$ kpc) down to the very small dissipative scales ($\ll 1$ pc). In this paper, we use a set of 490 pulsars distributed over roughly one third of the Galactic disk out to a radius $R \simeq 10$ kpc (assuming $R_\odot = 8.5$ kpc) and combine their observed rotation and dispersion measures with their estimated distances to derive the spatial energy spectrum of the Galactic interstellar magnetic field over the scale range $0.5 - 15$ kpc. We obtain a nearly flat spectrum, with a 1D power-law index $\alpha=-0.37\pm0.10$ for $E_{\rm B}(k)=C k^{\alpha}$ and an rms field strength of approximately $6 \mu$G over the relevant scales. Our study complements the derivation of the magnetic energy spectrum over the scale range $0.03 - 100$ pc by \citet{ms96b}, showing that the magnetic spectrum becomes flatter at larger scales. This observational result is discussed in the framework of current theoretical and numerical models.

This paper summarizes the results of 13 years of timing observations of a unique binary pulsar, PSR B1259$-$63, which has a massive B2e star companion. The data span encompasses four complete orbits and includes the periastron passages in 1990, 1994, 1997 and 2000. Changes in dispersion measure occurring around the 1994, 1997 and 2000 periastrons are measured and accounted for in the timing analysis. There is good evidence for a small glitch in the pulsar period in 1997 August, not long after the 1997 periastron, and a significant frequency second derivative indicating timing noise. We find that spin-orbit coupling with secular changes in periastron longitude and projected semi-major axis ($x$) cannot account for the observed period variations over the whole data set. While fitting the data fairly well, changes in pulsar period parameters at each periastron seem ruled out both by X-ray observations and by the large apparent changes in pulsar frequency derivative. Essentially all of the systematic period variations are accounted for by a model consisting of the 1997 August glitch and step changes in $x$ at each periastron. These changes must be due to changes in the orbit inclination, but we can find no plausible mechanism to account for them. It is possible that timing noise may mask the actual changes in orbital parameters at each periastron, but the good fit to the data of the $x$ step-change model suggests that short-term timing noise is not significant.

Mean pulse profiles and polarization properties are presented for nine southern pulsars. The observations were made using the Parkes radio telescope at frequencies near 1330 MHz; three of the nine pulsars were also observed at 660 MHz. A very high degree of circular polarization was observed in PSR J1603-7202. Complex position angle variations which are not well described by the rotating-vector model were observed in PSRs J2124-3358 and J2145-0750, both of which have very wide profiles. Rotation measures were obtained for all nine pulsars, with two implying relatively strong interstellar magnetic fields.

We report on a 10 ksec observation of the double pulsar system J0737-3039 obtained with the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer. We detect the system as a point source with coordinates (J2000) RA = 07:37:51.23, Decl. = -30:39:40.9, making this the first double neutron star system to be detected at X-ray energies. Only 77 source counts are detected. The data are represented satisfactorily by a power-law spectrum with photon index = 2.9 +/- 0.4. The inferred 0.2-10 keV luminosity is approximately 2e30 erg/s for a distance of 0.5 kpc. We do not detect any significant variability with orbital phase in these data, which span one full orbit. These results are consistent with the X-ray emission originating solely from the magnetosphere of the energetic A pulsar, with an efficiency for conversion of rotational energy to X-rays in this band of roughly 0.04%, although we cannot exclude other possibilities.

We present an updated phase-coherent timing solution for the young, energetic pulsar PSR B1509-58 for twenty years of data. Using a partially phase-coherent timing analysis, we show that the second frequency derivative is changing in time implying a third frequency derivative of (-9+/-1)x10^-32 s^(-4). This value is consistent with the simple power law model of pulsar rotation.

The clock-like properties of pulsars moving in the gravitational fields of their unseen neutron-star companions have allowed unique tests of general relativity and provided evidence for gravitational radiation. We report here the detection of the 2.8-sec pulsar J0737-3039B as the companion to the 23-ms pulsar J0737-3039A in a highly-relativistic double-neutron-star system, allowing unprecedented tests of fundamental gravitational physics. We observe a short eclipse of J0737-3039A by J0737-3039B and orbital modulation of the flux density and pulse shape of J0737-3039B, probably due to the influence of J0737-3039A's energy flux upon its magnetosphere. These effects will allow us to probe magneto-ionic properties of a pulsar magnetosphere.

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

Pulsars are among the most highly polarized sources in the universe. The NVSS has catalogued 2 million radio sources with linear polarization measurements, from which we have selected 253 sources, with polarization percentage greater than 25%, as targets for pulsar searches. We believe that such a sample is not biased by selection effects against ultra-short spin or orbit periods. Using the Parkes 64m telescope, we conducted searches with sample intervals of 0.05 ms and 0.08 ms, sensitive to submillisecond pulsars. Unfortunately we did not find any new pulsars.

PSRs J1847-0130 and J1718-37184 have inferred surface dipole magnetic fields greater than those of any other known pulsars and well above the ``quantum critical field'' above which some models predict radio emission should not occur.These fields are similar to those of the anomalous X-ray pulsars (AXPs), which growing evidence suggests are ``magnetars''. The lack of AXP-like X-ray emission from these radio pulsars (and the non-detection of radio emission from the AXPs) creates new challenges for understanding pulsar emission physics and the relationship between these classes of apparently young neutron stars.

We report on calculations of the pulsar birthrate based on the results of the Parkes multibeam survey. From the observed sample of more than 800 pulsars, we compute the pulsar current, accounting as accurately as possible for all known selection effects. The main goal of this work is to understand the pulsar birthrate as a function of the surface dipole magnetic field strengths. We show that pulsars with magnetic fields greater than 10**12.5 G account for about half of the total birthrate.

The number of known pulsars has significantly increased over the previous few years. We have searched the literature to find papers announcing the discovery of pulsars or giving improved parameters for them. Data from the papers have been entered into a new pulsar catalogue that can be accessed via a web interface or from the command line (on Solaris or Linux machines). The user may request over 120 different parameters, select pulsars of interest, generate custom variables and choose between different ways of displaying or tabulating the data. Full bibliographic references are available on all observed parameters.

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

We report the discovery of PSR J1847-0130, a radio pulsar with a 6.7-s spin period, in the Parkes multibeam survey of the Galactic plane. The slowdown rate for the pulsar, 1.3x10^{-12} s/s, is high and implies a surface dipole magnetic field strength of 9.4x10^{13} G. This inferred dipolar magnetic field strength is the highest by far among all known radio pulsars and over twice the ``quantum critical field'' above which some models predict radio emission should not occur. The inferred dipolar magnetic field strength and period of this pulsar are in the same range as those of the anomalous X-ray pulsars, which have been identified as being "magnetars" whose luminous X-ray emission is powered by their large magnetic fields. We have examined archival ASCA data and place an upper limit on the X-ray luminosity of J1847-0130 which is lower than the luminosities of all but one AXP. The properties of this pulsar prove that inferred dipolar magnetic field strength and period cannot alone be responsible for the unusual high-energy properties of the magnetars and creates new challenges for understanding the possible relationship between these two manifestations of young neutron stars.

The Parkes multibeam pulsar survey has unlocked vast areas of the Galactic plane which were previously invisible to earlier low-frequency and less-sensitive surveys. The survey has discovered more than 600 new pulsars so far, including many that are young and exotic. In this paper we report the discovery of 200 pulsars for which we present positional and spin-down parameters, dispersion measures, flux densities and pulse profiles. A large number of these new pulsars are young and energetic, and we review possible associations of $\gamma$-ray sources with the sample of about 1300 pulsars for which timing solutions are known. Based on a statistical analysis, we estimate that about $19\pm6$ associations are genuine. The survey has also discovered 12 pulsars with spin properties similar to those of the Vela pulsar, nearly doubling the known population of such neutron stars. Studying the properties of all known `Vela-like' pulsars, we find their radio luminosities to be similar to normal pulsars, implying that they are very inefficient radio sources. Finally, we review the use of the newly discovered pulsars as Galactic probes and discuss the implications of the new NE2001 Galactic electron density model for the determination of pulsar distances and luminosities.

We have carried out a deep search at 1.4 GHz for radio pulsed emission from six soft X-ray transient sources observed during their X-ray quiescent phase. The commonly accepted model for the formation of the millisecond radio pulsars predicts the presence of a rapidly rotating, weakly magnetized neutron star in the core of these systems. The sudden drop in accretion rate associated with the end of an X-ray outburst causes the Alfv\`en surface to move outside the light cylinder, allowing the pulsar emission process to operate. No pulsed signal was detected from the sources in our sample. We discuss several mechanisms that could hamper the detection and suggest that free-free absorption from material ejected from the system by the pulsar radiation pressure could explain our null result.

The Parkes Multibeam (PM) survey has found nearly 50% of the >1400 known pulsars and detected 20% of the pulsars in binary systems (75% within the survey area), few of these are millisecond pulsars with orbits of less than one day. It is known that the normal search techniques used in the PM Survey have selected against finding binary pulsars, due to the long integration times used. The relative performance of the search techniques used and the improvements made prior to the complete reprocessing of the data currently being undertaken are discussed. The processing is being done on COBRA - the new supercomputer at Jodrell Bank. The results of the reprocessing, as well as adding more solitary pulsars, will give more confidence on the population and types of binaries in the Galactic plane. However, it can be shown that for millisecond pulsars in binary orbits of between 20min and a few hours, the search techniques are still relatively insensitive.

Pulsars provide unique probes for the {\it large-scale} interstellar magnetic field in the Galactic disk. The Parkes multibeam pulsar survey has discovered many distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. The fields are found to be {\it coherent in direction} over a linear scale of $\sim 10$ kpc between the Carina-Sagittarius and Crux-Scutum arms from $l\sim45\degr$ to $l\sim305\degr$. The coherent spiral structures and field direction reversals, including the newly determined counterclockwise field near the Norma arm, are consistent with bi-symmetric spiral model for the disk field. However, the antisymmetric rotation measure sky from the Galactic halo and the dipole field in the Galactic center suggest that the A0 dynamo is operating there.

The young pulsar J1740-3052 is in an 8-month orbit with a companion of at least 11 solar masses. We present multifrequency GBT and Parkes timing observations, and discuss implications for the nature of the companion.

We review the properties of the 22 millisecond radio pulsars currently known in the globular cluster 47 Tucanae, and their implications for the mass and gas content of the cluster. Further details can be found in the publications from this project to date (Camilo et al. 2000; Freire 2000; Freire et al. 2001a,b; Freire et al. 2003). Throughout the review, we look ahead to the future results anticipated from this fascinating cluster.

We report the discovery of PSR J1747-2958, a radio pulsar with period P = 98 ms and dispersion measure DM = 101 pc/cc, in a deep observation with the Parkes telescope of the axially-symmetric "Mouse" radio nebula (G359.23-0.82). Timing measurements of the newly discovered pulsar reveal a characteristic age Pdt/2dP = 25 kyr and spin-down luminosity dE/dt = 2.5e36 erg/s. The pulsar (timing) position is consistent with that of the Mouse's "head". The distance derived from the DM, ~2 kpc, is consistent with the Mouse's distance limit from HI absorption, < 5.5 kpc. Also, the X-ray energetics of the Mouse are compatible with being powered by the pulsar. Therefore we argue that PSR J1747-2958, moving at supersonic speed through the local interstellar medium, powers this unusual non-thermal nebula. The pulsar is a weak radio source, with period-averaged flux density at 1374 MHz of 0.25 mJy and luminosity ~1 mJy kpc^2.

We present the results of first deep optical observations of the field of the old ($\sim 10^8$ yr), nearby, isolated pulsar J0108-1431, in an attempt to detect its optical counterpart. The observations were performed using the FORS1 instrument at the focus of the European Southern Observatory Antu Telescope of the VLT. Observations with the Australia Telescope Compact Array (ATCA) were made to determine an accurate position for the radio pulsar at the current epoch. The imaging data, obtained in the V, B, and U passbands reveal no counterpart at the revised radio position down to $V \simeq 28$, $B\simeq 28.6$ and $U \simeq 26.4$. For a distance of 130 pc, estimated from the pulsar's dispersion measure, our constraints on the optical flux put an upper limit of $T=4.5\times 10^4$ K for the surface temperature of the neutron star, assuming a stellar radius $R_\infty=13$ km. Our new radio position allows us to place an upper limit on the pulsar proper motion of 82 mas yr$^{-1}$ which, for d=130 pc, implies a transverse velocity $\la 50$ km sec$^{-1}$.

We report here on a sequence of 28 observations of the binary pulsar system PSR B1259-63/SS2883 at four radio frequencies made with the Australia Telescope Compact Array around the time of the 2000 periastron passage. Observations made on 2000 Sep 1 show that the pulsar's apparent rotation measure (RM) reached a maximum of $-14800 \pm 1800$ rad m$^{-2}$, some 700 times the value measured away from periastron, and is the largest astrophysical RM measured. This value, combined with the dispersion measure implies a magnetic field in the Be star's wind of 6 mG. We find that the light curve of the unpulsed emission is similar to that obtained during the 1997 periastron but that differences in detail imply that the emission disc of the Be star is thicker and/or of higher density. The behaviour of the light curve at late times is best modelled by the adiabatic expansion of a synchrotron bubble formed in the pulsar/disc interaction. The expansion rate of the bubble $\sim 12$ km s$^{-1}$ is surprisingly low but the derived magnetic field of 1.6 G close to that expected.

PSR J1806-2125 is a pulsar discovered in the Parkes multibeam pulsar survey with a rotational period of 0.4s and a characteristic age of 65kyr. Between MJDs 51462 and 51894 this pulsar underwent an increase in rotational frequency of $\Delta \nu/\nu \approx 16 \times 10^{-6}$. The magnitude of this glitch is approximately 2.5 times greater than any previously observed in any pulsar and 16 times greater than the mean glitch size. This letter gives the parameters of the glitch and compares its properties to previously observed events. The existence of such large and rare glitches offers new hope for attempts to observe thermal X-ray emission from the internal heat released following a glitch, and suggests that pulsars which previously have not been observed to glitch may do so on long timescales.

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

The Parkes multibeam pulsar survey is a sensitive survey of a strip of the Galactic plane with $|b|<5\degr$ and $260\degr < l < 50\degr$ at 1374 MHz. Here we report the discovery of 120 new pulsars and subsequent timing observations, primarily using the 76-m Lovell radio telescope at Jodrell Bank. The main features of the sample of 370 published pulsars discovered during the multibeam survey are described. Furthermore, we highlight two pulsars: PSR J1734$-$3333, a young pulsar with the second highest surface magnetic field strength among the known radio pulsars, $B_s = 5.4\times10^{13}$ G, and PSR J1830$-$1135, the second slowest radio pulsar known, with a 6-s period.

Pulsars provide unique probes of the large-scale interstellar magnetic field in the Galactic disk. Up to now, the limited Galactic distribution of the known pulsar population has restricted these investigations to within a few kiloparsec of the Sun. The Parkes multibeam pulsar survey has discovered many more-distant pulsars which enables us for the first time to explore the magnetic field in most of the nearby half of the Galactic disk. Here we report the detection of counterclockwise magnetic fields in the Norma spiral arm using pulsar rotation measures. The fields are coherent in direction over a linear scale of $\sim 5$ kpc along the arm and have a strength of $-4.4\pm0.9$ $\mu$G. The magnetic field between the Carina-Sagittarius and Crux-Scutum arms is confirmed to be coherent from $l\sim45\degr$ to $l\sim305\degr$ over a length of $\sim 10$ kpc. These results strengthen arguments for a bisymmetric spiral model for the field configuration in the Galactic disk.

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

We report on spectral and timing observations of the nearest millisecond pulsar J0437-4715 with Chandra. The pulsar spectrum, detected up to 7 keV, cannot be described by a simple one-component model. We suggest that it consists of two components, a nonthermal power-law spectrum generated in the pulsar magnetosphere, with a photon index about 2, and a thermal spectrum emitted by heated polar caps, with a temperature decreasing outwards from 2 MK to 0.5 MK. The lack of spectral features in the thermal component suggests that the neutron star surface is covered by a hydrogen (or helium) atmosphere. The timing analysis shows one X-ray pulse per period, with a pulsed fraction of about 40% and the peak at the same pulse phase as the radio peak. No synchrotron pulsar-wind nebula is seen in the X-rays.

We report the results of the spectral and timing analysis of observations of the Vela pulsar with the Chandra X-ray Observatory. The spectrum shows no statistically significant spectral lines in the observed 0.25--8.0 keV band. It consists of two distinct continuum components. The softer component can be modeled as either a magnetic hydrogen atmosphere spectrum with kT = 59 +- 3 eV, R = 15.5 +- 1.5 km, or a standard blackbody with kT = 129 +- 4 eV, R = 2.1 +- 0.2 km (the radii are for a distance of 250 pc). The harder component, modeled as a power-law spectrum, gives photon indices depending on the model adopted for the soft component: gamma = 1.5 +- 0.3 for the magnetic atmosphere soft component, and gamma = 2.7 +- 0.4 for the blackbody soft component. Timing analysis shows three peaks in the pulse profile, separated by about 0.3 in phase. Energy-resolved timing provides evidence for pulse profile variation with energy. The higher energy (E > 1.8 keV) profile shows significantly higher pulsed fraction.

The Parkes multibeam pulsar survey is covering a 10deg-wide strip of the southern Galactic plane from l=260deg to l=50deg. It utilizes a 13-beam receiver operating in the 20-cm band on the Parkes 64-m radio telescope and is much more sensitive than any previous large-scale survey. Most of the 608 pulsars discovered so far are relatively distant and many are young, with 37 having a characteristic age of less than 10^5 years. At least one of these is associated with a supernova remnant and four other probable associations are suggested. Several multibeam pulsars have high values of the parameter Edot/d^2 and are within the position error contours of unidentified EGRET gamma-ray sources. These possible associations will be tested with the advent of new gamma-ray telescopes.

We describe a sensitive targeted search campaign at Parkes to find faint young radio pulsars using the center beam of the 20-cm multibeam receiver. The high frequency of the receiver mitigates scattering effects which can dominate profiles of distant pulsars at low frequencies and prevent detection. The targets included several composite supernova remnants (SNRs), two fast X-ray pulsars, long-period anomalous X-ray pulsars, and a soft gamma repeater. The sensitivity limits in this search are significantly better than those previously reported for these objects. We have discovered three new radio pulsars, two of which are not associated with the target objects. The third pulsar, PSR J1713-3949, is possibly associated with SNR G347.3-0.5; the estimated ages and distances of the two objects are consistent, but a more significant positional coincidence must be established with the central X-ray source in the SNR through pulsar timing to confirm or refute the association.

We report on Australia Telescope Compact Array observations in the direction of the young high-magnetic-field radio pulsar PSR J1119-6127. In the resulting images we identify a non-thermal radio shell of diameter 15 arcminutes, which we classify as a previously uncataloged young supernova remnant (SNR), G292.2-0.5. This SNR is positionally coincident with PSR J1119-6127, and we conclude that the two objects are physically associated. No radio emission is detected from any pulsar wind nebula (PWN) associated with the pulsar; our observed upper limits are consistent with the expectation that high-magnetic-field pulsars produce radio nebulae which fade rapidly. This system suggests a possible explanation for the lack of associated radio pulsars and/or PWNe in many SNRs.

The development of the radio remnant of SN 1987A has been followed using the Australia Telescope Compact Array since its first detection in 1990 August. The remnant has been observed at four frequencies, 1.4, 2.4, 4.8 and 8.6 GHz, at intervals of 4 - 6 weeks since the first detection. These data are combined with the 843 MHz data set of Ball et al. (2001) obtained at Molonglo Observatory to study the spectral and temporal variations of the emission. These observations show that the remnant continues to increase in brightness, with a larger rate of increase at recent times. They also show that the radio spectrum is becoming flatter, with the spectral index changing from -0.97 to -0.88 over the 11 years. In addition, at roughly yearly intervals since 1992, the remnant has been imaged at 9 GHz using super-resolution techniques to obtain an effective synthesised beamwidth of about 0".5. The imaging observations confirm the shell-morphology of the radio remnant and show that it continues to expand at ~3000 km/s. The bright regions of radio emission seen on the limb of the shell do not appear to be related to the optical hotspots which have subsequently appeared in surrounding circumstellar material.

A pulsar timing system has been operating in the 18-cm band at the Urumqi Astronomical Observatory 25-m telescope since mid-1999. Frequency resolution allowing dedispersion of the pulsar signals is provided by a 2X128X2.5 MHz filterbank/digitiser system. Observations of 74 pulsars over more than 12 months have resulted in updated pulsar periods and period derivatives, as well as improved positions. Comparison with previous measurements showed that the changes in period and period derivative tend to have the same sign and to be correlated in amplitude. A model based on unseen glitches gives a good explanation of the observed changes, suggesting that long-term fluctuations in period and period derivatives are dominated by glitches. In 2000 July, we detected a glitch of relative amplitude Delta_nu/nu~24X1e-9 in the Crab pulsar. The post-glitch decay appears similar to other large Crab glitches.

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

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

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 report the detection of ionized intracluster gas in the globular cluster 47 Tucanae. Pulsars in this cluster with a negative period derivative, which must lie in the distant half of the cluster, have significantly higher measured integrated electron column densities than the pulsars with a positive period derivative. We derive the plasma density within the central few pc of the cluster using two different methods which yield consistent values. Our best estimate of n_e = (0.067+-0.015)/cm^3 is about 100 times the free electron density of the ISM in the vicinity of 47 Tucanae, and the ionized gas is probably the dominant component of the intracluster medium.

We report the discovery of a young and energetic pulsar in the Parkes multibeam survey of the Galactic plane. PSR J1016-5857 has a rotation period of 107 ms and period derivative of 8e-14, implying a characteristic age of 21 kyr and spin-down luminosity of 2.6e36 erg/s. The pulsar is located just outside, and possibly interacting with, the shell supernova remnant G284.3-1.8. Archival X-ray data show a source near the pulsar position which is consistent with emission from a pulsar wind nebula. The pulsar is also located inside the error box of the unidentified EGRET source 3EG J1013-5915, for which it represents a plausible counterpart.

We report on radio polarimetry observations for nine southern pulsars. Six of the nine in the sample are young, with characteristic ages under 100 kyr and high spin-down luminosities. All six show a significant degree of linear polarization. We also confirm a previously noticed trend in which the degree of linear polarization increases with spin-down luminosity. Where possible we have used the rotating-vector model of the pulsar emission geometry to fit the observed position angle data. Our fit for PSR J1513-5908 (B1509-58) in particular is useful for directly testing the magnetospheric model of Melatos (1997) in combination with further timing observations. For this pulsar we find that a magnetic inclination angle greater than or equal to 60 degrees is excluded at the 3-sigma level, and that the geometry suggested by the morphology of an apparent bipolar X-ray outflow is marginally inconsistent with the Melatos model. We also report on polarimetry of three older pulsars: PSR J0045-7319, PSR J1627-4845, and PSR J1316-6232 (whose discovery we also report). Of these, only PSR J0045-7319 shows significant polarization.

The Parkes multibeam pulsar survey is a sensitive survey of a strip along the Galactic plane with |b|<5 deg and l=260 deg to l=50 deg. It uses a 13-beam receiver on the 64-m Parkes radio telescope, receiving two polarisations per beam over a 288 MHz bandwidth centred on 1374 MHz. Receiver and data acquisition systems are described in some detail. For pulsar periods in the range 0.1 - 2 s and dispersion measures of less than 300 cm^{-3} pc, the nominal limiting flux density of the survey is about 0.2 mJy. At shorter or longer periods or higher dispersions, the sensitivity is reduced. Timing observations are carried out for pulsars discovered in the survey for 12 - 18 months after confirmation to obtain accurate positions, spin parameters, dispersion measures, pulse shapes and mean flux densities. The survey is proving to be extremely successful, with more than 600 pulsars discovered so far. We expect that, when complete, this one survey will come close to finding as many pulsars as all previous pulsar surveys put together. The newly discovered pulsars tend to be young, distant and of high radio luminosity. They will form a valuable sample for studies of pulsar emission properties, the Galactic distribution and evolution of pulsars, and as probes of interstellar medium properties. This paper reports the timing and pulse shape parameters for the first 100 pulsars timed at Parkes, including three pulsars with periods of less than 100 ms which are members of binary systems. These results are briefly compared with the parameters of the previously known population.

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

In the last 10 years 20 millisecond pulsars have been discovered in the globular cluster 47 Tucanae. Hitherto, only 3 of these had published timing solutions. Here we improve upon these 3 and present 12 new solutions. These measurements can be used to determine a variety of physical properties of the pulsars and of the cluster. The 15 pulsars have positions determined with typical uncertianties of only a few milliarcsec and they are all located within 1.2 arcmin of the cluster centre. We have also measured the proper motions of 5 of the pulsars, which are consistent with the proper motion of 47 Tuc based on Hipparcos data. The period derivatives measured for many of the pulsars are dominated by the dynamical effects of the cluster gravitational field, and are used to constrain the surface mass density of the cluster. All pulsars have characteristic ages T > 170 Myr and magnetic fields B < 2.4e9 Gauss, and the average T > 1 Gyr. We have measured the rate of advance of periastron for the binary pulsar J0024-7204H, implying a total system mass 1.4+-0.8 solar masses.

We report on the discovery of a binary pulsar, PSR J1740-3052, during the Parkes multibeam survey. Timing observations of the 570-ms pulsar at Jodrell Bank and Parkes show that it is young, with a characteristic age of 350 kyr, and is in a 231-day, highly eccentric orbit with a companion whose mass exceeds 11 M_sun. An accurate position for the pulsar was obtained using the Australia Telescope Compact Array. Near-infrared 2.2-um observations make with the telescopes at the Siding Spring observatory reveal a late-type star coincident with the pulsar position. However, we do not believe that this star is the pulsar's companion, because a typical star of this spectral type and required mass would extend beyond the pulsar's orbit. Furthermore, the measured advance of periastron of the pulsar suggests a more compact companion, for example, a main-sequence star with radius only a few times that of the sun. Such a companion is also more consistent with the small dispersion measure variations seen near periastron. Although we cannot conclusively rule out a black-hole companion, we believe the companion is probably an early B star, making the system similar to the binary PSR J0045-7319.

We report on Australia Telescope Compact Array observations in the direction of the young high magnetic-field pulsar J1119-6127. In the resulting images we identify a non-thermal radio shell of diameter 15', which we classify as a previously uncatalogued young supernova remnant, G292.2-0.5. This supernova remnant is positionally coincident with PSR J1119-6127, and we conclude that the two objects are physically associated. No radio emission is detected from any pulsar wind nebula (PWN) associated with the pulsar; our observed upper limits are consistent with the expectation that high magnetic-field pulsars produce radio nebulae which fade rapidly. This system suggests a possible explanation for the lack of an associated radio pulsar and/or PWN in many supernova remnants.

We report on five binary pulsars discovered in the Parkes multibeam Galactic plane survey. All of the pulsars are old, with characteristic ages 1-11 Gyr, and have relatively small inferred magnetic fields, 5-90e8 G. The orbital periods range from 1.3 to 15 days. As a group these objects differ from the usual low-mass binary pulsars (LMBPs): their spin periods of 9-88 ms are relatively long; their companion masses, 0.2-1.1 Msun, are, in at least some cases, suggestive of CO or more massive white dwarfs; and some of the orbital eccentricities, 1e-5 < e < 0.002, are unexpectedly large. We argue that these observed characteristics reflect binary evolution that is significantly different from that of LMBPs. We also note that intermediate-mass binary pulsars apparently have a smaller scale-height than LMBPs.

We report the results of a survey of the Small Magellanic Cloud (SMC) for radio pulsars conducted with the 20-cm multibeam receiver of the Parkes 64-meter telescope. This survey targeted a more complete region of the SMC than a previous pulsar search and had an improvement in sensitivity by a factor of about two for most pulsar periods. This survey is much more sensitive to fast young pulsars (with P < 100 ms) and is the first survey of the SMC with any sensitivity to millisecond pulsars. Two new pulsars were discovered in the survey, one of which is located within the SMC. The number of pulsars found in the survey is consistent with the expected number derived using several methods. We also report the serendipitous discovery of a new pulsar in the 30 Doradus region of the Large Magellanic Cloud (LMC). These discoveries bring the total number of rotation-powered pulsars currently known in the Magellanic Clouds to eight. We have also made refined timing measurements for the new discoveries as well as for three previously known LMC pulsars. The age distribution of luminous Magellanic Cloud pulsars supports the conjecture that pulsars younger than about 5 Myr are more luminous on average than older pulsars.

Using all available multi-component radio pulse profiles for pulsars with medium to long periods and good polarisation data, we have constructed a two-dimensional image of the mean radio beamshape. This shows a peak near the centre of the beam but is otherwise relatively uniform with only mild enhancements in a few regions. This result supports the patchy-beam model for emission beams in which the mean beam shape represents the properties of the emission mechanism and observed pulse components result from emission sources distributed randomly across the beam.

We report the discovery using the Parkes radio telescope of binary millisecond pulsars in four clusters for which no associated pulsars were previously known. The four pulsars have pulse periods lying between 3 and 6 ms. All are in circular orbits with low-mass companions and have orbital periods of a few days or less. One is in a 1.7-hour orbit with a companion of planetary mass. Another is eclipsed by a wind from its companion for 40% of the binary period despite being in a relatively wide orbit. These discoveries result from the use of improved technologies and prove that many millisecond pulsars remain to be found in globular clusters.

We report the detection of a new 3.06 ms binary pulsar in the globular cluster NGC6544 using a Fourier-domain ``acceleration'' search. With an implied companion mass of ~0.01 solar masses and an orbital period of only P_b~1.7 hours, it displays very similar orbital properties to many pulsars which are eclipsed by their companion winds. The orbital period is the second shortest of known binary pulsars after 47 Tuc R. The measured flux density of 1.3 +/- 0.4 mJy at 1332 MHz indicates that the pulsar is almost certainly the known steep-spectrum point source near the core of NGC6544.

There are many reasons why it is important to increase the number of known pulsars. Not only do pulsar searches continue to improve statistical estimates of, for example, pulsar birthrates, lifetimes and the Galactic distribution, but they continue to turn up interesting and, in some cases, unique individual pulsars. In the early days of pulsar astronomy, the Molonglo radio telescope led the world as a pulsar detection instrument. However, the Parkes radio telescope, with its frequency versatility and greater tracking ablility, combined with sensitive receivers and powerful computer detection algorithms, is now the world's most successful telescope at finding pulsars. The Parkes multibeam survey, begun in 1997, by itself will come close to doubling the number of known pulsars. Parkes has also been very successful at finding millisecond pulsars, especially in globular clusters. One third of the known millisecond pulsars have been found in just one cluster, 47 Tucanae.

We report here on multifrequency radio observations of the pulsed emission from PSR B1259-63 around the time of the closest approach to its B2e companion star. There was a general increase in the pulsar's dispersion measure and scatter broadening, and a decrease in the flux density towards periastron although changes in these parameters were seen on timescales as short as minutes. The pulsed emission disappeared 16 days prior to periastron and remained undetectable until 16 days after periastron. The observations are used to determine the parameters of the wind from the Be star. We show that a simple model, in which the wind density varies with radius as r^-2, provides a good fit to the data. The wind is clumpy with size scales <= 10^10 cm, densities of 10^6 cm^-3 and a velocity of 2000 km/s at a distance of 20 - 50 stellar radii. We find a correlation between dispersion measure variations and the pulse scattering times, suggesting that the same electrons are responsible for both effects.

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

We report on the discovery of PSR J1141-6545, a radio pulsar in an eccentric, relativistic 5-hr binary orbit. The pulsar shows no evidence for being recycled, having pulse period P = 394 ms, characteristic age tau_c = 1.4 x 10^6 yr, and inferred surface magnetic dipole field strength B = 1.3 x 10^12 G. From the mass function and measured rate of periastron advance, we determine the total mass in the system to be (2.300 +/- 0.012) solar masses, assuming that the periastron advance is purely relativistic. Under the same assumption, we constrain the pulsar's mass to be M_p < 1.348 solar masses and the companion's mass to be M_c > 0.968 solar masses (both 99% confidence). Given the total system mass and the distribution of measured neutron star masses, the companion is probably a massive white dwarf which formed prior to the birth of the pulsar. Optical observations can test this hypothesis.

We report the discovery of two young isolated radio pulsars with very high inferred magnetic fields. PSR J1119-6127 has period P = 0.407 s, and the largest period derivative known among radio pulsars, Pdot = 4.0e-12. Under standard assumptions these parameters imply a characteristic spin-down age of only tau = 1.6 kyr and a surface dipole magnetic field strength of B = 4.1e13 G. We have measured a stationary period-second-derivative for this pulsar, resulting in a braking index of n = 2.91+-0.05. We have also observed a glitch in the rotation of the pulsar, with fractional period change Delta_P/P = -4.4e-9. Archival radio imaging data suggest the presence of a previously uncataloged supernova remnant centered on the pulsar. The second pulsar, PSR J1814-1744, has P = 3.975 s and Pdot = 7.4e-13. These parameters imply tau = 85 kyr, and B = 5.5e13 G, the largest of any known radio pulsar. Both PSR J1119-6127 and PSR J1814-1744 show apparently normal radio emission in a regime of magnetic field strength where some models predict that no emission should occur. Also, PSR J1814-1744 has spin parameters similar to the anomalous X-ray pulsar (AXP) 1E 2259+586, but shows no discernible X-ray emission. If AXPs are isolated, high magnetic field neutron stars (``magnetars''), these results suggest that their unusual attributes are unlikely to be merely a consequence of their very high inferred magnetic fields.

We report on searches of the globular cluster Terzan 5 for low luminosity and accelerated radio pulsars using the 64-m Parkes radio telescope. One new millisecond pulsar, designated PSR J1748-2446C, was discovered, having a period of 8.44ms. Timing measurements using the 76-m Lovell radio telescope at Jodrell Bank show that it is a solitary pulsar and lies close to the core of the cluster. We also present the results of timing measurements which show that the longer-period pulsar PSR J1748-2444 (formerly known as PSR B1744-24B) lies 10 arcmin from the core of the cluster and is unlikely to be associated with the cluster. We conclude that there are further pulsars to be detected in the cluster.

We report the discovery of two isolated radio pulsars having the largest inferred surface dipole magnetic fields yet seen in the population: 4.1e10^13 G and 5.5e10^13 G. These pulsars show apparently normal radio emission in a regime of magnetic field strength where some models predict no emission should occur. They have spin parameters and magnetic fields similar to those of some magnetar candidates, but exhibit very different radiative properties. This demonstrates that if the putative magnetars are indeed isolated neutron stars, their unusual attributes cannot be only a consequence of their large inferred magnetic fields.

A high-frequency survey of the Galactic plane for radio pulsars is in progress, using the multibeam receiver on the 64-m Parkes radiotelescope. We describe the survey motivations, the observing plan and the inital results. The survey is discovering many pulsars, more than 500 so far. Eight of the new pulsars are binary, one with a massive companion. At least eight are young, with characteristic ages of less than 100 kyr. Two of these (Kaspi et al, this Conference) have surface dipole magnetic field strengths greater than any other known radio pulsar.

Recent observations of the globular cluster 47 Tuc, made with the Parkes telescope at a wavelength of 20 cm, have resulted in the discovery of nine new millisecond pulsars, all in binary systems. The number of timing solutions available has risen from two to 14. These results will make possible a more detailed study of the cluster dynamics.

Measurement of accurate positions, pulse periods and period derivatives is an essential follow-up to any pulsar survey. The procedures being used to obtain timing parameters for the pulsars discovered in the Parkes multibeam pulsar survey are described. Completed solutions have been obtained so far for about 80 pulsars. They show that the survey is preferentially finding pulsars with higher than average surface dipole magnetic fields. Eight pulsars have been shown to be members of binary systems and some of the more interesting results relating to these are presented.

The Parkes multibeam pulsar survey began in 1997 and is now about 50% complete. It has discovered more than 400 new pulsars so far, including a number of young, high magnetic field, and relativistic binary pulsars. Early results, descriptions of the survey and follow up timing programs can be found in papers by Lyne et al. (1999 MNRAS in press, astro-ph/9911313), Camilo et al. (astro-ph/9911185), and Manchester et al. (astro-ph/9911319). This paper describes the data release policy and how you can gain access to the raw data and details on the pulsars discovered.

We are undertaking a high-frequency survey of the Galactic plane for radio pulsars, using the 13-element multibeam receiver on the 64-m Parkes radio telescope. We describe briefly the survey system and some of the initial results. PSR J1811-1736, one of the first pulsars discovered with this system, has a rotation period of 104 ms. Subsequent timing observations using the 76-m radio telescope at Jodrell Bank show that it is in an 18.8-day, highly-eccentric binary orbit. We have measured the rate of advance of periastron which indicates a total system mass of 2.6 +- 0.9 Msun, and the minimum companion mass is about 0.7 Msun. This, the high orbital eccentricity and the recycled nature of the pulsar suggests that this system is composed of two neutron stars, only the fourth or fifth such system known in the disk of the Galaxy.

Polarization profiles of several hundred pulsars have been published recently (eg. Gould & Lyne 1998, GL98; Weisberg et al. 1999). In this report, we summarize the characteristics of circular and high linear polarizations of pulsar profiles, based on all previously published data.

We have used a new observing system on the Parkes radio telescope to carry out a series of pulsar observations of the globular cluster 47 Tucanae at 20-cm wavelength. We detected all 11 previously known pulsars, and have discovered nine others, all of which are millisecond pulsars in binary systems. We have searched the data for relatively short orbital period systems, and found one pulsar with an orbital period of 96 min, the shortest of any known radio pulsar. The increased rate of detections with the new system resulted in improved estimates of the flux density of the previously known pulsars, determination of the orbital parameters of one of them, and a coherent timing solution for another one. Five of the pulsars now known in 47 Tucanae have orbital periods of a few hours and implied companion masses of only ~ 0.03 Msun. Two of these are eclipsed at some orbital phases, while three are seen at all phases at 20 cm but not always at lower frequencies. Four and possibly six of the other binary systems have longer orbital periods and companion masses ~ 0.2 Msun, with at least two of them having relatively large orbital eccentricities. All 20 pulsars have rotation periods in the range 2-8 ms.

We present radio polarimetry results for nine Southern pulsars. Six of the nine are young, with characteristic ages less than 100 kyr and high spin-down luminosities. All six show significant linear polarization, and we confirm a previously noticed trend in which the degree of linear polarization increases with spin-down luminosity. We have used the rotating vector model to fit the observed position angle data for PSR J1513-5908 (B1509-58). We find that a magnetic inclination angle alpha > 60 degrees is excluded at the 3-sigma level in the fit, and that the geometry suggested by the morphology of an apparent bipolar X-ray outflow is marginally inconsistent with a recent model of the pulsar magnetosphere.

The Parkes multibeam pulsar survey uses a 13-element receiver operating at a wavelength of 20 cm to survey the inner Galactic plane with remarkable sensitivity. To date we have collected and analyzed data from 45% of the survey region (|b| < 5 deg.; 260 deg. < l < 50 deg.), and have discovered 440 pulsars, in addition to re-detecting 190 previously known ones. Most of the newly discovered pulsars are at great distances, as inferred from a median dispersion measure (DM) of 400 pc/cm^3.

The interaction between the ejecta from Supernova 1987A and surrounding material is producing steadily brightening radio and X-ray emission. The new-born supernova remnant has been significantly decelerated by this interaction, while its morphology reflects the axisymmetric nature of the progenitor wind.

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

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

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

The Parkes multibeam pulsar survey is a high-frequency, fast-sampled survey of the Galactic Plane, expected to discover at least 500 new pulsars. To date, over 200 pulsars have been found, including several young pulsars and at least one with a very high magnetic field. Seven of the new stars are in binary systems; this number includes one probable double-neutron-star system, and one pulsar with an extremely massive companion.

We have obtained 63 rotation measures (RMs) from polarization observations of southern pulsars, of which 54 are new measurements and 3 are varied from previous values. The new pulsar RM data at high Galactic latitudes are mostly consistent with the antisymmetric RM distribution found previously. For the Galactic disc, evidence for a field reversal near the Perseus arm, and possibly another beyond it, is presented. Inside the Solar Circle, in addition to the two known field reversals in or near the Carina-Sagittartus arm and the Crux-Scutum arm, a further reversal in the Norma arm is tentatively identified. These reversals, together with the pitch angle derived from pulsar RM and stellar polarization distributions, are consistent with bisymmetric spiral (BSS) models for the large-scale magnetic field structure in the disc of our Galaxy. However, discrimination between models is complicated by the presence of smaller-scale irregularities in the magnetic field, as well as uncertainties in the theoretical modelling.

We describe radio continuum and spectral line observations of the source G320.4-01.2 (MSH 15-52) and the coincident young pulsar B1509-58 made with the Australia Telescope Compact Array. Based on a comparison between X-ray and radio observations, we argue that the two main radio components of G320.4-01.2 are a single supernova remnant (SNR), which HI absorption indicates is at a distance of 5.2+/-1.4 kpc. A high-resolution correspondence between radio and X-rays argues that the pulsar is interacting with the SNR via an opposed pair of collimated outflows. The outflow itself is seen as an elongated X-ray feature surrounded by a highly polarized radio sheath, while the interaction with the SNR manifests itself as a ring of radio/X-ray knots within the optical nebula RCW 89. We reject the hypothesis that the pulsar outflow powers the entire RCW 89 region. SNR G320.4-01.2 and PSR B1509-58 agree in distance and in rotation measure, and appear to be interacting. We thus conclude that the two objects are associated and have an age of <~1700 yr. We propose that the SNR resulted from a high-energy or low-mass supernova which occurred near the edge of an elongated cavity. Such a model can account for the SNR's bilateral appearance, its large apparent age, the significant offset of the pulsar from the SNR's centre and the faintness of the pulsar-powered nebula at radio wavelengths.

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

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

G320.4-01.2 is a complex radio and X-ray source, coinciding on the sky with the young energetic pulsar B1509-58. A young pulsar embedded in a SNR would seem to accord with expectations, but previous observations suggest that all may not be what it seems. Controversy persists over whether the pulsar and the SNR are associated, and as to what causes the remnant's unusual appearance. To answer these questions, we have undertaken a set of high-resolution radio observations of the system. We present the results of this study, which provide new evidence that PSR B1509-58 is associated with and is interacting with G320.4-01.2.

Using the Parkes radio telescope, we have searched for pulsed radio emission from the recently discovered X-ray pulsars AX J1811.5-1926 and PSR J0537-6910 in the supernova remnants G11.2-0.3 and N157B, respectively. We detected no significant pulsed radio emission from these pulsars and have set an upper limit of 0.07 mJy on the 1374 MHz flux from AX J1811.5-1926 and upper limits of 0.18 mJy and 0.06 mJy for the flux from PSR J0537-6910 at 660 MHz and 1374 MHz, respectively. Assuming a power law radio spectral index of 2, these flux limits correspond to luminosity limits at 400 MHz of 20 mJy kpc^2 for AX J1811.5-1926 and 1100 mJy kpc^2 for PSR J0537-6910. Our luminosity limit for AX J1811.5-1926 is lower than the observed luminosities of other young radio pulsars. The upper limit on the luminosity for PSR J0537-6910 in N157B is not significantly constraining due to the large distance to the pulsar. We have also searched for giant radio pulses from both pulsars and have found no convincing candidates.

We present a systematic study of the circular polarization in pulsar integrated profiles, based on published polarization data. For core components, we find no significant correlation between the sense-change of circular polarization and the sense of linear position angle variation. Circular polarization is not restricted to core components and, in some cases, reversals of circular polarization sense are observed across the conal emission. In conal double profiles, the sense of circular polarization is found to be correlated with the sense of position-angle variation. Pulsars with a high degree of linear polarization often have one hand of circular polarization across the whole profile. For most pulsars, the sign of circular polarization is the same at 50-cm and 20-cm wavelength, and the degree of polarization is similar, albeit with a wide scatter. However, at least two cases of frequency-dependent sign reversals are known. This diverse behaviour may require more than one mechanism to generate circular polarization.

We report the detection of the radio pulsar counterpart to the 69 ms X-ray pulsar discovered near the supernova remnant RCW 103 (G332.4-0.4). Our detection confirms that the pulsations arise from a rotation-powered neutron star, which we name PSR J1617-5055. The observed barycentric period derivative confirms that the pulsar has a characteristic age of only 8 kyr, the sixth smallest of all known pulsars. The unusual apparent youth of the pulsar and its proximity to a young remnant requires that an association be considered. Although the respective ages and distances are consistent within substantial uncertainties, the large inferred pulsar transverse velocity is difficult to explain given the observed pulsar velocity distribution, the absence of evidence for a pulsar wind nebula, and the symmetry of the remnant. Rather, we argue that the objects are likely superposed on the sky; this is reasonable given the complex area. Without an association, the question of where is the supernova remnant left behind following the birth of PSR J1617-5055 remains open. We also discuss a possible association between PSR J1617-5055 and the gamma-ray source 2CG 333+01. Though an association is energetically plausible, it is unlikely given that EGRET did not detect 2CG 333+01.

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

We present Australia Telescope Compact Array observations of the supernova remnant (SNR) G309.2-00.6. In a 1.3-GHz continuum image the remnant appears as a near-circular shell, but with two brightened and distorted arcs of emission on opposite sides. HI absorption against the SNR yields a distance in the range 5.4 to 14.1 kpc, corresponding to an age (1-20)x10^3 yr. On the basis of the SNR's morphology we argue that it is a younger analogue of the W50/SS433 system, and that its unusual appearance is a result of opposed jets or outflows from a central source. A jet-like feature and breaks in the shell can both be seen along the axis of proposed outflow, providing further support for this interpretation; the central source itself is not detected. The SNR may be interacting with the adjacent HII region RCW 80 through an extension of the proposed outflow beyond its shell. This would put the SNR at the lower limit of its distance range and would imply an age <4000 yr. We consider other SNRs similar to G309.2-00.6, and propose remnants whose shells are affected by jets as one of several classes of SNR from which the presence of a central source can be inferred.

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

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

Observations with the Rossi X-ray Timing Explorer and the Jodrell Bank, Parkes, and Green Bank telescopes have enabled us to determine the time delay between radio and X-ray pulses in the two isolated pulsars B1821-24 and B1509-58. For the former we find that the narrow X-ray and radio pulse components are close to being coincident in time, with the radio peak leading by 0.02 period (60 +/- 20 microsec), while the wide X-ray pulse component lags the last of the two wider radio components by about 0.08 period. For the latter pulsar we find, using the standard value for the dispersion measure, that the X-ray pulse lags the radio by about 0.27 period, with no evidence for any energy-dependence in the range 2-100 keV. However, uncertainties in the history of the dispersion measure for this pulsar make a comparison with previous results difficult. It is clear that there are no perceptable variations in either the lag or the dispersion measure at time scales of a year or less.

We present Australia Telescope Compact Array (ATCA) observations of the supernova remnant (SNR) G296.8-00.3. A 1.3-GHz continuum image shows the remnant to have a complex multi-shelled appearance, with an unusual rectangular strip running through its centre. HI absorption yields a kinematic distance to the remnant of 9.6+/-0.6 kpc, from which we estimate an age in the range (2-10) x 10^3 yr. The ATCA's continuum mode allows a measurement of the Faraday rotation across the band, from which we derive a mean rotation measure towards the SNR of 430 rad m^(-2). We consider possible explanations for the morphology of G296.8-00.3, and conclude that it either has a biannular structure, as might be produced through interaction with an asymmetric progenitor wind, or that its appearance is caused by the effects of the surrounding interstellar medium. We argue that the adjacent pulsar J1157-6224 is at a similar distance as the SNR, but that a physical association is highly unlikely. The pulsar is the only detectable source in the field in circular polarization, suggesting a method for finding pulsars during aperture synthesis.

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

We present seven years of radio observations of SN 1987A made with the Australia Telescope Compact Array. At 1.4, 2.4, 4.8 and 8.6 GHz, the flux density of the radio remnant has increased monotonically since emission was redetected 1200 days after the explosion. On day 3200, the remnant was expanding at 2800 +/- 400 km/s, which we interpret as indicating significant deceleration of the fastest moving ejecta. Since day 1787 the spectral index has remained constant at alpha = -0.95 +/- 0.04. These observations are all consistent with the shock having encountered a denser, shocked, component of the progenitor's stellar wind. At the current rate of expansion, the shock is expected to encounter the inner optical ring in the year 2006 +/- 3. Using super-resolution, we have also obtained 9 GHz images of the remnant (resolution approx 0".5) at four epochs. The emission is distributed around the rim of a near-circular shell, but has become increasingly asymmetric with time. There are two "hotspots" to the east and west, aligned along the major axis of the optical ring. This morphology is most likely indicative of an axisymmetric circumstellar medium into which the shock is expanding, consistent with present understanding of the progenitor star and its environment.

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