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

Pulsar radio emission and its polarization are observed to evolve with frequency. This frequency dependence is key to the emission mechanism and the structure of the radio beam. With the new Ultra-Wideband receiver (UWL) on the Parkes radio telescope we are able, for the first time, to observe how pulsar profiles evolve over a broad continuous bandwidth of 700-4000 MHz. We describe here a technique for processing broadband polarimetric observations to establish a meaningful alignment and visualize the data across the band. We apply this to observations of PSRs J1056-6258 and J1359-6038, chosen due to previously unresolved questions about the frequency evolution of their emission. Application of our technique reveals that it is possible to align the polarization position angle (PA) across a broad frequency range when constrained to applying only corrections for dispersion and Faraday rotation to do so. However, this does not correspond to aligned intensity profiles for these two sources. We find that it is possible to convert these misalignments into emission height range estimates that are consistent with published and simulated values, suggesting that they can be attributed to relativistic effects in the magnetosphere. We discuss this work in the context of the radio beam structure and prepare the ground for a wider study of pulsar emission using broadband polarimetric data.

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

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

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 report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

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

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

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

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

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

The globular cluster Omega Centauri is the most massive and luminous cluster in the Galaxy. The $\gamma$-ray source FL8Y J1326.7$-$4729 is coincident with the core of the cluster, leading to speculation that hitherto unknown radio pulsars or annihilating dark matter may be present in the cluster core. Here we report on the discovery of five millisecond pulsars in Omega Centauri following observations with the Parkes radio telescope. Four of these pulsars are isolated with spin periods of 4.1, 4.2, 4.6 and 6.8 ms. The fifth has a spin period of 4.8 ms and is in an eclipsing binary system with an orbital period of 2.1 hours. Deep radio continuum images of the cluster centre with the Australian Telescope Compact Array reveal a small population of compact radio sources making it likely that other pulsars await discovery. We consider it highly likely that the millisecond pulsars are the source of the $\gamma$-ray emission. The long-term timing of these pulsars opens up opportunities to explore the dynamics and interstellar medium of the cluster.

We establish that for PSR B0031-07 the orthogonal polarization modes switch at a single pulse level synchronously with the periodic drifting subpulses seen in total intensity. There are only four other pulsars known for which this phenomenon is observed. PSR B0031-07 is unique as it is the only source in this group which has multiple stable drift modes. For both drift modes visible at our observing frequency centered at 1369 MHz, the modulation of polarization modes is synchronous with the drifting subpulses. In one of the drift modes, a discontinuity in the modulation pattern of polarization properties occurs halfway through the pulse, coinciding with a slight change in the slope of the intensity drift band. In contrast to what has been suggested for this pulsar in the past, this, plus other differences in the polarization of the modulated emission observed for the two drift modes, suggests that a drift mode change is more than a change in the underlying carousel radius and that magnetospheric propagation effects play an important role. The ellipticity evolves asymmetrically in time during the modulation cycle, which in the framework of a carousel model implies that the polarized sub-beams are asymmetric with respect to the sense of circulation, something which is not observed for other pulsars. Birefringence in the magnetosphere, resulting in the orthogonal polarization modes to spatially separate, is not enough to explain these results. It is argued that more complex magnetospheric processes, which possibly allow conversion between orthogonal polarization modes, play a role.

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.

In the canonical picture of pulsars, radio emission arises from a narrow cone centered on the star's magnetic axis but many basic details remain unclear. We use high-quality polarization data taken with the Parkes radio telescope to constrain the geometry and emission heights of pulsars showing interpulse emission, and include the possibility that emission heights in the main and interpulse may be different. We show that emission heights are low in the centre of the beam, typically less than 3\% of the light cylinder radius. The emission beams are under-filled in longitude, with an average profile width only 60% of the maximal beam width and there is a strong preference for the visible emission to be located on the trailing part of the beam. We show substantial evidence that the emission heights are larger at the beam edges than in the beam centre. There is some indication that a fan-like emission beam explains the data better than conal structures. Finally, there is a strong correlation between handedness of circular polarization in the main and interpulse profiles which implies that the hand of circular polarization is determined by the hemisphere of the visible emission.

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.

It is now known that the emission from radio pulsars can vary over a wide range of timescales, from fractions of seconds to decades. However, it is not yet known if long- and short-term emission variability are caused by the same physical processes. It has been observed that long-term emission variability is often correlated with rotational changes in the pulsar. We do not yet know if the same is true of short-term emission variability, as the rotational changes involved cannot be directly measured over such short timescales. To remedy this, we propose a continuous pulsar monitoring technique that permits the statistical detection of any rotational changes in nulling and mode-changing pulsars with certain properties. Using a simulation, we explore the range of pulsar properties over which such an experiment would be possible.

The frequency widening of pulsar profiles is commonly attributed to lower frequencies being produced at greater heights above the surface of the pulsar; so-called radius-to-frequency mapping. The observer's view of pulsar emission is a 1D cut through a 3D magnetosphere: we can only see that emission which points along our line of sight. However, by comparing the frequency evolution of many single pulses positioned at different phases, we can build up an understanding of the shape of the active emission region. We use single pulses observed with the Giant Metrewave Radio Telescope to investigate the emission region of PSR J1136+1551 and test radius-to-frequency mapping. Assuming that emission is produced tangential to the magnetic field lines and that each emission frequency corresponds to a single height, we simulate the single pulse profile evolution resulting from the canonical conal beam model and a fan beam model. Comparing the results of these simulations with the observations, we conclude that the emission region of PSR J1136+1551 is better described by the fan beam model. The diversity of profile widening behaviour observed for the single pulses can be explained by orthogonally polarized modes propagating along differing frequency-dependent paths in the magnetosphere.

Radio pulsars are often used as clocks in a wide variety of experiments. Imperfections in the clock, known as timing noise, have the potential to reduce the significance of, or even thwart e.g. the attempt to find a stochastic gravitational wave (GW) background. We measure the timing noise in a group of 129 mostly middle-aged pulsars (i.e. characterstic ages near 1~Myr) observed with the Parkes radio telescope on a monthly basis since 2014. We examine four different metrics for timing noise, but it remains unclear which, if any, provides the best determination. In spite of this, it is evident that these pulsars have significantly less timing noise than their younger counterparts, but significantly more than the (much older) millisecond pulsars (MSPs). As with previous authors, we find a strong correlation between timing noise and the pulsar spin-down rate, $\dot{\nu}$. However, for a given $\dot{\nu}$ there is a spread of about a factor 30 in the strength of the timing noise likely indicating that nuclear conditions in the interior of the stars differs between objects. We briefly comment on the implications for GW detection through pulsar timing arrays as the level of timing noise in MSPs may be less than predicted.

PKS B1322-110 is a radio quasar that is located only 8.5' in angular separation from the bright B star Spica. It exhibits intra-day variability in its flux density at GHz frequencies attributed to scintillations from plasma inhomogeneities. We have tracked the rate of scintillation of this source for over a year with the Australia Telescope Compact Array, recording a strong annual cycle that includes a near-standstill in August and another in December. The cycle is consistent with scattering by highly anisotropic plasma microstructure, and we fit our data to that model in order to determine the kinematic parameters of the plasma. Because of the low ecliptic latitude of PKS B1322-110, the orientation of the plasma microstructure is poorly constrained. Nonetheless at each possible orientation our data single out a narrow range of the corresponding velocity component, leading to a one-dimensional constraint in a two-dimensional parameter space. The constrained region is consistent with a published model in which the scattering material is associated with Spica and consists of filaments that are radially oriented around the star. This result has a 1% probability of arising by chance.

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

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

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

In the standard picture of radio pulsars, the radio emission arises from a set of open magnetic field lines, the extent of which is primarily determined by the pulsar's spin period, P, and the emission height. We have used a database of parameters from 600 pulsars to show that the observed profile width, W, follows W proportional to P^-0.3 albeit with a large scatter, emission occurs from heights below 400 km and that the beam is underfilled. Furthermore, the prevalence in the data for long period pulsars to have relatively wide profiles can only be explained if the angle between the magnetic and rotation axis decays with time.

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.

We have conducted the largest investigation to date into the origin of phase resolved, apparent RM variations in the polarized signals of radio pulsars. From a sample of 98 pulsars based on observations at 1.4 GHz with the Parkes radio telescope, we carefully quantified systematic and statistical errors on the measured RMs. A total of 42 pulsars showed significant phase resolved RM variations. We show that both magnetospheric and scattering effects can cause these apparent variations. There is a clear correlation between complex profiles and the degree of RM variability, in addition to deviations from the Faraday law. Therefore, we conclude that scattering cannot be the only cause of RM variations, and show clear examples where magnetospheric effects dominate. It is likely that, given sufficient signal-to-noise, such effects will be present in all radio pulsars. These signatures provide a tool to probe the propagation of the radio emission through the magnetosphere.

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

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

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year.

We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies.

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

We report on the measurement and investigation of pulsed high-energy $\gamma$-ray emission from the Vela pulsar, PSR B0833$-$45, based on 40.3 hours of observations with the largest telescope of H.E.S.S., CT5, in monoscopic mode, and on 8 years of data obtained with the Fermi-LAT. A dedicated very-low-threshold event reconstruction and analysis pipeline was developed and, together with the CT5 telescope response model, was validated using the Fermi-LAT data as reference. A pulsed $\gamma$-ray signal at a significance level of more than $15\sigma$ is detected from the P2 peak of the Vela pulsar light curve. Of a total of 15835 events, more than 6000 lie at an energy below 20 GeV, implying a significant overlap between H.E.S.S. II-CT5 and the Fermi-LAT. While the investigation of the pulsar light curve with the LAT confirms characteristics previously known up to 20 GeV, in the tens of GeV energy range, CT5 data show a change in the pulse morphology of P2, i.e., an extreme sharpening of its trailing edge, together with the possible onset of a new component at 3.4$\sigma$ significance level. Assuming a power-law model for the P2 spectrum, an excellent agreement is found for the photon indices ($\Gamma \simeq$ 4.1) obtained with the two telescopes above 10 GeV and an upper bound of 8% is derived on the relative offset between their energy scales. Using both instruments data, it is however shown that the spectrum of P2 in the 10-100 GeV has a pronounced curvature, i.e. a confirmation of the sub-exponential cutoff form found at lower energies with the LAT. This is further supported by the weak evidence for an emission above 100 GeV obtained with CT5. In contrast, converging indications are found from both CT5 and LAT data for the emergence of a hard component above 50 GeV in the leading wing (LW2) of P2, which possibly extends beyond 100 GeV.

The pulsar/massive star binary system PSR B1259-63 / LS 2883 is one of the best-studied gamma-ray binaries, a class of systems whose bright gamma-ray flaring can provide important insights into high-energy physics. Using the Australian Long Baseline Array we have conducted very long baseline interferometric observations of PSR B1259-63 over 4.4 years, fully sampling the 3.4-year orbital period. From our measured parallax of $0.38\pm0.05$ mas we use a Bayesian approach to infer a distance of $2.6^{+0.4}_{-0.3}$ kpc. We find that the binary orbit is viewed at an angle of $154\pm3$ degrees to the line of sight, implying that the pulsar moves clockwise around its orbit as viewed on the sky. Taking our findings together with previous results from pulsar timing observations, all seven orbital elements for the system are now fully determined. We use our measurement of the inclination angle to constrain the mass of the stellar companion to lie in the range 15-31$M_{\odot}$. Our measured distance and proper motion are consistent with the system having originated in the Cen OB1 association and receiving a modest natal kick, causing it to have moved $\sim$8 pc from its birthplace over the past $\sim3\times10^5$ years. The orientation of the orbit on the plane of the sky matches the direction of motion of the X-ray synchrotron-emitting knot observed by the Chandra X-ray Observatory to be moving away from the system.

The star KIC 8462852 (Boyajian's Star) displays both fast dips of up to 20% on time scales of days, plus long-term secular fading by up to 19% on time scales from a year to a century. We report on CCD photometry of KIC 8462852 from 2015.75 to 2018.18, with 19,176 images making for 1,866 nightly magnitudes in BVRI. Our light curves show a continuing secular decline (by 0.023 +- 0.003 mags in the B-band) with three superposed dips with duration 120-180 days. This demonstrates that there is a continuum of dip durations from a day to a century, so that the secular fading is seen to be by the same physical mechanism as the short-duration Kepler dips. The BVRI light curves all have the same shape, with the slopes and amplitudes for VRI being systematically smaller than in the B-band by factors of 0.77 +- 0.05, 0.50 +- 0.05, and 0.31 +- 0.05. We rule out any hypothesis involving occultation of the primary star by any star, planet, solid body, or optically thick cloud. But these ratios are the same as that expected for ordinary extinction by dust clouds. This chromatic extinction implies dust particle sizes going down to ~0.1 micron, suggesting that this dust will be rapidly blown away by stellar radiation pressure, so the dust clouds must have formed within months. The modern infrared observations were taken at a time when there was at least 12.4% +- 1.3% dust coverage (as part of the secular dimming), and this is consistent with dimming originating in circumstellar dust.

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

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

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

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

The Murchison Widefield Array (MWA), and its recently-developed Voltage Capture System (VCS), facilitates extending the low-frequency range of pulsar observations at high-time and -frequency resolution in the Southern Hemisphere, providing further information about pulsars and the ISM. We present the results of an initial time-resolved census of known pulsars using the MWA. To significantly reduce the processing load, we incoherently sum the detected powers from the 128 MWA tiles, which yields ~10% of the attainable sensitivity of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185 MHz), allowing multiple pulsars to be observed simultaneously. We developed a WIde-field Pulsar Pipeline (WIPP) that processes the data from each observation and automatically folds every known pulsar located within the beam. We have detected 50 pulsars to date, 6 of which are millisecond pulsars. This is consistent with our expectation, given the telescope sensitivity and the sky coverage of the processed data (~17,000 deg2). For ten pulsars, we present the lowest-frequency detections published. For a subset of the pulsars, we present multi-frequency pulse profiles by combining our data with published profiles from other telescopes. Since the MWA is a low-frequency precursor to the Square Kilometre Array (SKA), we use our census results to forecast that a survey using Phase 1 of SKA-Low (SKA1-Low) can potentially detect around 9400 pulsars.

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

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

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

We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 au, and strong flux density enhancement before and after the ESE may arise from refractive focusing. Long observations reveal scintillation arcs characteristic of interference between rays scattered at large angles, and the clearest arcs appear during the ESE. The arcs suggest scattering by a screen 100-200 pc from the earth, perhaps ionized filamentary structure associated with the boundary of the local bubble(s). Towards PSR J1740-3015 we observed a "double dip" in the measured flux density similar to ESEs observed towards compact extragalactic radio sources. The observed shape is consistent with that produced by a many-au scale diverging plasma lens with electron density ~500cm$^{-3}$ . The continuing ESE is at least 1500 d long, making it the longest detected event to date. These detections, with materially different observational signatures, indicate that well-calibrated pulsar monitoring is a keen tool for ESE detection and ISM diagnostics. They illustrate the strong r\^ole au-scale non-Kolmogorov density fluctuations and the local ISM structure play in such events and are key to understanding both their intrinsic physics and their impact on other phenomena, particularly fast radio bursts.

Over the past 13 years, the Parkes radio telescope has observed a large number of pulsars using digital filterbank backends with high time and frequency resolution and the capability for Stokes recording. Here we use archival data to present polarimetry data at an observing frequency of 1.4 GHz for 600 pulsars with spin-periods ranging from 0.036 to 8.5 s. We comment briefly on some of the statistical implications from the data and highlight the differences between pulsars with high and low spin-down energy. The dataset, images and table of properties for all 600 pulsars are made available in a public data archive maintained by the CSIRO.

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

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

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

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

In Dai et al. (2016) we developed a formalism for computing variance images from standard, interferometric radio images containing time and frequency information. Variance imaging with future radio continuum surveys allows us to identify radio pulsars and serves as a complement to conventional pulsar searches which are most sensitive to strictly periodic signals. Here, we carry out simulations to predict the number of pulsars that we can uncover with variance imaging on future continuum surveys. We show that the Australian SKA Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) survey can find $\sim30$ normal pulsars and $\sim40$ millisecond pulsars (MSPs) over and above the number known today, and similarly an all-sky continuum survey with SKA-MID can discover $\sim140$ normal pulsars and $\sim110$ MSPs with this technique. Variance imaging with EMU and SKA-MID will detect pulsars with large duty cycles and is therefore a potential tool for finding MSPs and pulsars in relativistic binary systems. Compared with current pulsar surveys at high Galactic latitudes in the southern hemisphere, variance imaging with EMU and SKA-MID will be more sensitive, and will enable detection of pulsars with dispersion measures between $\sim10$ and 100 cm$^{-3}$ pc.

Radio observations of the region surrounding PSR J1301-6305 at 5.5 GHz and 7.5 GHz were conducted with ATCA on September 5th, 2013. They were dedicated to the search of the radio counterpart of the evolved pulsar wind nebula HESS J1303-631, detected in X-rays and GeV-TeV gamma-rays. The collected data do not reveal any significant extended emission associated with PSR J1301-6305. In addition, archival 1.384 GHz and 2.368 GHz data do not show any evidence for a radio counterpart of HESS J1303-631. Archival 1.384 GHz observations reveal a detection of an extended structure centred at an angular distance of 190 from the pulsar. This extended structure might be a Supernova remnant (SNR) and a potential birth place of PSR J1301-6305. The implications of the lack of radio counterpart of HESS J1303-631 on the understanding of the nature of the PWN are discussed.

We report the discovery of a synchrotron nebula, G283.1-0.59, associated with PSR J1015-5719. Radio observations using the Molonglo Observatory Synthesis Telescope and the Australia Telescope Compact Array at 36, 16, 6, and 3 cm reveal a complex morphology. The pulsar is embedded in the "head" of the nebula with fan-shaped diffuse emission. This is connected to a circular bubble of 20" radius and a collimated tail extending over 1'. Polarization measurements show a highly ordered magnetic field in the nebula. It wraps around the edge of the head and shows an azimuthal configuration near the pulsar, then switches direction quasi-periodically near the bubble and in the tail. Together with the flat radio spectrum observed, we suggest that this system is most plausibly a pulsar wind nebula (PWN), with the head as a bow shock that has a low Mach number and the bubble as a shell expanding in a dense environment. The bubble could act as a magnetic bottle trapping the relativistic particles. A comparison with other bow-shock PWNe with higher Mach numbers shows similar structure and B-field geometry, implying that pulsar velocity may not be the most critical factor in determining the properties of these systems. We also derive analytic expressions for the projected standoff distance and shape of an inclined bow shock. It is found that the projected distance is always larger than the true distance in three dimensions. On the other hand, the projected shape is not sensitive to the inclination after rescaling with the projected standoff distance.

During February 2016, CSIRO Astronomy and Space Science and the Max-Planck-Institute for Radio Astronomy installed, commissioned and carried out science observations with a phased array feed (PAF) receiver system on the 64m diameter Parkes radio telescope. Here we demonstrate that the PAF can be used for pulsar observations and we highlight some unique capabilities. We demonstrate that the pulse profiles obtained using the PAF can be calibrated and that multiple pulsars can be simultaneously observed. Significantly, we find that an intrinsic polarisation leakage of -31dB can be achieved with a PAF beam offset from the centre of the field of view. We discuss the possibilities for using a PAF for future pulsar observations and for searching for fast radio bursts with the Parkes and Effelsberg telescopes.

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

We report radio observations of the quasar PMN J1106-3647. Our data, taken with the Australia Telescope Compact Array, show large variations in the amplitude and shape of its spectrum, on a short time-scale. A great variety of spectral features is evident, including: sharp kinks; broad spectral peaks; and wiggles. No two spectra are alike. We interpret the variations as interstellar scintillation of a radio source that is compact, but not point-like. Under this interpretation, complex spectral structure can arise purely refractively, under high magnification conditions, or from interference between waves that have been scattered by small-scale density fluctuations (diffractive scintillation). Both effects may be playing a role in J1106-3647, and we tentatively identify kinks with the former, and wiggles with the latter. Diffractive scintillation of AGN is uncommon, as the fringe visibility is low for all but the most compact radio sources. Refractive interpretation of the kink implies that the source has a sharp, concave boundary. Our data are consistent with a mildly boosted synchrotron source, provided the scattering material is at a distance ~50 pc from us.

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

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

We use data on extreme radio scintillation to demonstrate that this phenomenon is associated with hot stars in the solar neighbourhood. The ionized gas responsible for the scattering is found at distances up to 1.75pc from the host star, and on average must comprise 1.E5 distinct structures per star. We detect azimuthal velocities of the plasma, relative to the host star, up to 9.7 km/s, consistent with warm gas expanding at the sound speed. The circumstellar plasma structures that we infer are similar in several respects to the cometary knots seen in the Helix, and in other planetary nebulae. There the ionized gas appears as a skin around tiny molecular clumps. Our analysis suggests that molecular clumps are ubiquitous circumstellar features, unrelated to the evolutionary state of the star. The total mass in such clumps is comparable to the stellar mass.

When a massive star in a binary system explodes as a supernova its companion star may be polluted with heavy elements from the supernova ejecta. Such a pollution had been detected in a handful of post-supernova binaries (Gonzalez Hernandez et al. 2011), but none of them is associated with a supernova remnant. We report the discovery of a solar-type star in a close, eccentric binary system with a neutron star within the young Galactic supernova remnant RCW 86. Our discovery implies that the supernova progenitor was a moving star, which exploded near the edge of its wind bubble and lost most of its initial mass due to common-envelope evolution shortly before core collapse. We find that the solar-type star is strongly polluted with calcium and other elements, which places the explosion within the class of calcium-rich supernovae -- faint and fast transients (Filippenko et al. 2003; Kasliwal et al. 2012}, whose origin is strongly debated (Kawabata et al. 2010; Waldman et al. 2011), and provides the first observational evidence that supernovae of this type can arise from core-collapse explosions (Kawabata et al. 2010; Moriya et al. 2010).

We report on a timing and spectral analysis of the young, high magnetic field rotation-powered pulsar (RPP) B1509-58 using Chandra continuous-clocking mode observation. The pulsar's X-ray light curve can be fit by the two Gaussian components and the pulsed fraction shows moderate energy dependence over the Chandra band. The pulsed X-ray spectrum is well described by a power law with a photon index 1.16(4), which is harder than the values measured with RXTE PCA and NuSTAR. This result supports the log-parabolic model for the broadband X-ray spectrum. With the unprecedented angular resolution of Chandra, we clearly identified off-pulse X-ray emission from the pulsar, and its spectrum is best fit by a power law plus blackbody model. The latter component has a temperature of ~0.14 keV with a bolometric luminosity comparable to the luminosities of other young and high magnetic field RPPs, and it lies between the temperature of magnetars and typical RPPs. In addition, we found that the nonthermal X-ray emission of PSR B1509-58 is significantly softer in the off-pulse phase than in the pulsed phase, with the photon index varying between 1.0 and 1.8 and anticorrelated with the flux. This is similar to the behavior of three other young pulsars. We interpreted it as different contributions of pair-creation processes at different altitudes from the neutron star surface according to the outer-gap model.

We present low-frequency spectral energy distributions of 60 known radio pulsars observed with the Murchison Widefield Array (MWA) telescope. We searched the GaLactic and Extragalactic All-sky MWA (GLEAM) survey images for 200-MHz continuum radio emission at the position of all pulsars in the ATNF pulsar catalogue. For the 60 confirmed detections we have measured flux densities in 20 x 8 MHz bands between 72 and 231 MHz. We compare our results to existing measurements and show that the MWA flux densities are in good agreement.

The location of radio pulsars in the period-period derivative (P-Pdot) plane has been a key diagnostic tool since the early days of pulsar astronomy. Of particular importance is how pulsars evolve through the P-Pdot diagram with time. Here we show that the decay of the inclination angle (alpha-dot) between the magnetic and rotation axes plays a critical role. In particular, alpha-dot strongly impacts on the braking torque, an effect which has been largely ignored in previous work. We carry out simulations which include a negative alpha-dot term, and show that it is possible to reproduce the observational P-Pdot diagram without the need for either pulsars with long birth periods or magnetic field decay. Our best model indicates a birth rate of 1 radio pulsar per century and a total Galactic population of ~20000 pulsars beaming towards Earth.

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

We present observations made with the Australia Telescope Compact Array (ATCA), the Jansky Very Large Array (JVLA) and the Giant Metre-Wave Telescope of the radio source within the galaxy WISE~J071634.59-190039.2, claimed to be host of FRB~150418 by Keane et al. (2016). We have established a common flux density scale between the ATCA and JVLA observations, the main result of which is to increase the flux densities obtained by Keane et al. At a frequency of 5.5 GHz, the source has a mean flux density of 140uJy and is variable on short timescales with a modulation index of 0.36. Statistical analysis of the flux densities shows that the variations seen are consistent with refractive interstellar scintillation of the weak active galactic nucleus at the centre of the galaxy. It may therefore be the case that the FRB and the galaxy are not associated. However, taking into account the rarity of highly variable sources in the radio sky, and our lack of knowledge of the progenitors of FRBs as a class, the association between WISE~J071634.59-190039.2 and FRB~150418 remains a possibility.

We compare radio profile widths of young, energetic gamma-ray-detected and non-gamma-ray-detected pulsars. We find that the latter typically have wider radio profiles, with the boundary between the two samples exhibiting a dependence on the rate of rotational energy loss. We also find that within the sample of gamma-ray-detected pulsars, radio profile width is correlated with both the separation of the main gamma-ray peaks and the presence of narrow gamma-ray components. These findings lead us to propose that these pulsars form a single population where the main factors determining gamma ray detectability are the rate of rotational energy loss and the proximity of the line of sight to the rotation axis. The expected magnetic inclination angle distribution will be different for radio pulsars with and without detectable gamma rays, naturally leading to the observed differences. Our results also suggest that the geometry of existing radio and outer-magnetosphere gamma-ray emission models are at least qualitatively realistic, implying that information about the viewing geometry can be extracted from profile properties of pulsars.

A fading radio source, coincident in time and position with the fast radio burst FRB150418, has been associated with the galaxy WISE J071634.59-190039.2. Subsequent observations of this galaxy have revealed that it contains a persistent, but variable, radio source. We present e-MERLIN, VLBA, and ATCA radio observations and Subaru optical observations of WISE J071634.59-190039.2 and find that the persistent radio source is unresolved and must be compact (<0.01 kpc), and that its location is consistent with the optical centre of the galaxy. We conclude that it is likely that WISE J071634.59-190039.2 contains a weak radio AGN.

Pulsars are the only cosmic radio sources known to be sufficiently compact to show diffractive interstellar scintillations. Images of the variance of radio signals in both time and frequency can be used to detect pulsars in large-scale continuum surveys using the next generation of synthesis radio telescopes. This technique allows a search over the full field of view while avoiding the need for expensive pixel-by-pixel high time resolution searches. We investigate the sensitivity of detecting pulsars in variance images. We show that variance images are most sensitive to pulsars whose scintillation time-scales and bandwidths are close to the subintegration time and channel bandwidth. Therefore, in order to maximise the detection of pulsars for a given radio continuum survey, it is essential to retain a high time and frequency resolution, allowing us to make variance images sensitive to pulsars with different scintillation properties. We demonstrate the technique with Murchision Widefield Array data and show that variance images can indeed lead to the detection of pulsars by distinguishing them from other radio sources.

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.

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.

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

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

We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field.

Extreme scattering events (ESEs) are distinctive fluctuations in the brightness of astronomical radio sources caused by occulting plasma lenses in the interstellar medium. The inferred plasma pressures of the lenses are $\sim 10^3$ times the ambient pressure, challenging our understanding of gas conditions in the Milky Way. Using a new survey technique, we have discovered an ESE while it was in progress. We report radio and optical follow-up observations. Modelling of the radio data demonstrates that the lensing structure is a density enhancement and that the lens is diverging, ruling out one of two competing physical models. Our technique will uncover many more ESEs, addressing a long-standing mystery of the small-scale gas structure of the Galaxy.

The Boolardy Engineering Test Array is a 6 x 12 m dish interferometer and the prototype of the Australian Square Kilometre Array Pathfinder (ASKAP), equipped with the first generation of ASKAP's phased array feed (PAF) receivers. These facilitate rapid wide-area imaging via the deployment of simultaneous multiple beams within a 30 square degree field of view. By cycling the array through 12 interleaved pointing positions and using 9 digitally formed beams we effectively mimic a traditional 1 hour x 108 pointing survey, covering 150 square degrees over 711 - 1015 MHz in 12 hours of observing time. Three such observations were executed over the course of a week. We verify the full bandwidth continuum imaging performance and stability of the system via self-consistency checks and comparisons to existing radio data. The combined three epoch image has arcminute resolution and a 1-sigma thermal noise level of 375 micro-Jy per beam, although the effective noise is a factor 3 higher due to residual sidelobe confusion. From this we derive a catalogue of 3,722 discrete radio components, using the 35 percent fractional bandwidth to measure in-band spectral indices for 1,037 of them. A search for transient events reveals one significantly variable source within the survey area. The survey covers approximately two-thirds of the Spitzer South Pole Telescope Deep Field. This pilot project demonstrates the viability and potential of using PAFs to rapidly and accurately survey the sky at radio wavelengths.

Compact radio sources sometimes exhibit intervals of large, rapid changes in their flux-density, due to lensing by interstellar plasma crossing the line-of-sight. A novel survey program has made it possible to discover these "Extreme Scattering Events" (ESEs) in real time, resulting in a high-quality dynamic spectrum of an ESE observed in PKS 1939-315. Here we present a method for determining the column-density profile of a plasma lens, given only the dynamic radio spectrum of the lensed source, under the assumption that the lens is either axisymmetric or totally anisotropic. Our technique relies on the known, strong frequency dependence of the plasma refractive index in order to determine how points in the dynamic spectrum map to positions on the lens. We apply our method to high-frequency (4.2-10.8 GHz) data from the Australia Telescope Compact Array of the PKS 1939-315 ESE. The derived electron column-density profiles are very similar for the two geometries we consider, and both yield a good visual match to the data. However, the fit residuals are substantially above the noise level, and deficiencies are evident when we compare the predictions of our model to lower-frequency (1.6-3.1 GHz) data on the same ESE, thus motivating future development of more sophisticated inversion techniques.

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

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

Pulsars are known to display short-term variability. Recently, examples of longer-term emission variability have emerged that are often correlated with changes in the rotational properties of the pulsar. To further illuminate this relationship, we have developed techniques to identify emission and rotation variability in pulsar data, and determine correlations between the two. Individual observations may be too noisy to identify subtle changes in the pulse profile. We use Gaussian process (GP) regression to model noisy observations and produce a continuous map of pulse profile variability. Generally, multiple observing epochs are required to obtain the pulsar spin frequency derivative ($\dot{\nu}$). GP regression is, therefore, also used to obtain $\dot{\nu}$, under the hypothesis that pulsar timing noise is primarily caused by unmodelled changes in $\dot{\nu}$. Our techniques distinguish between two types of variability: changes in the total flux density versus changes in the pulse shape. We have applied these techniques to 168 pulsars observed by the Parkes radio telescope, and see that although variations in flux density are ubiquitous, substantial changes in the shape of the pulse profile are rare. We reproduce previously published results and present examples of profile shape changing in seven pulsars; in particular, a clear new example of correlated changes in profile shape and rotation is found in PSR~J1602$-$5100. In the shape changing pulsars, a more complex picture than the previously proposed two state model emerges. We conclude that our simple assumption that all timing noise can be interpreted as $\dot{\nu}$ variability is insufficient to explain our dataset.

PSR J1357$-$6429 is a young and energetic radio pulsar detected in X-rays and $\gamma$-rays. It powers a compact pulsar wind nebula with a jet visible in X-rays and a large scale plerion detected in X-ray and TeV ranges. Previous multiwavelength studies suggested that the pulsar has a significant proper motion of about 180 mas yr$^{-1}$ implying an extremely high transverse velocity of about 2000 km s$^{-1}$. In order to verify that, we performed radio-interferometric observations of PSR J1357$-$6429 with the the Australia Telescope Compact Array (ATCA) in the 2.1 GHz band. We detected the pulsar with a mean flux density of $212\pm5$ $\mu$Jy and obtained the most accurate pulsar position, RA = 13:57:02.525(14) and Dec = $-$64:29:29.89(15). Using the new and archival ATCA data, we did not find any proper motion and estimated its 90 per cent upper limit $\mu < 106$ mas yr$^{-1}$. The pulsar shows a highly polarised single pulse, as it was earlier observed at 1.4 GHz. Spectral analysis revealed a shallow spectral index $\alpha_{\nu}$ = $0.5 \pm 0.1$. Based on our new radio position of the pulsar, we disclaim its optical counterpart candidate reported before.

In a search for periodic variation in the arrival times of pulses from 151 young, energetic pulsars, we have identified seven cases of modulation consistent with one or two harmonics of a single fundamental with time-scale 0.5-1.5 yr. We use simulations to show that these modulations are statistically significant and of high quality (sinusoidal) even when contaminated by the strong stochastic timing noise common to young pulsars. Although planetary companions could induce such modulation, the large implied masses and 2:1 mean motion resonances challenge such an explanation. Instead, the modulation is likely to be intrinsic to the pulsar, arising from quasi-periodic switching between stable magnetospheric states, and we propose that precession of the neutron star may regulate this switching.

We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-\alpha}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058.

Based on its energy-dependent morphology the initially unidentified very high energy (VHE; E>100 GeV) gamma-ray source HESS J1303-631 was recently associated with the pulsar PSR J1301-6305. Subsequent detection of X-ray and GeV counterparts also supports the identification of the H.E.S.S. source as evolved pulsar wind nebula (PWN). We report here on recent radio observations of the PSR J1301-6305 field of view (FOV) with ATCA dedicated to search for the radio counterpart of this evolved PWN. Observations at 5.5 GHz and 7.5 GHz do not reveal any extended emission associated with the pulsar. The analysis of the archival 1.384 GHz and 2.368 GHz data also does not show any significant emission. The 1.384 GHz data reveal a hint of an extended shell-like emission in the PSR J1301-6305 FOV which might be a supernova remnant. We discuss the implications of the non-detection at radio wavelengths on the nature and evolution of the PWN as well as the possibility of the SNR candidate being the birth place of PSR J1301-6305.

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

We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$.

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

We have searched a sample of 151 young, energetic pulsars for periodic variation in pulse time-of-arrival arising from the influence of planetary companions. We are sensitive to objects with masses two orders of magnitude lower than those detectable with optical transit timing, but we find no compelling evidence for pulsar planets. For the older pulsars most likely to host planets, we can rule out Mercury analogues in one third of our sample and planets with masses $>0.4 M_{\oplus}$ and periods $P_b<1$ yr in all but 5% of such systems. If pulsar planets form primarily from supernova fallback disks, these limits imply that such disks do not form, are confined to $<0.1$ AU radii, are disrupted, or form planets more slowly ($>2$ Myr) than their protoplanetary counterparts.

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

We present HI imaging of the galaxy group IC 1459 carried out with six antennas of the Australian SKA Pathfinder equipped with phased-array feeds. We detect and resolve HI in eleven galaxies down to a column density of $\sim10^{20}$ cm$^{-2}$ inside a ~6 deg$^2$ field and with a resolution of ~1 arcmin on the sky and ~8 km/s in velocity. We present HI images, velocity fields and integrated spectra of all detections, and highlight the discovery of three HI clouds -- two in the proximity of the galaxy IC 5270 and one close to NGC 7418. Each cloud has an HI mass of $10^9$ M$_\odot$ and accounts for ~15% of the HI associated with its host galaxy. Available images at ultraviolet, optical and infrared wavelengths do not reveal any clear stellar counterpart of any of the clouds, suggesting that they are not gas-rich dwarf neighbours of IC 5270 and NGC 7418. Using Parkes data we find evidence of additional extended, low-column-density HI emission around IC 5270, indicating that the clouds are the tip of the iceberg of a larger system of gas surrounding this galaxy. This result adds to the body of evidence on the presence of intra-group gas within the IC 1459 group. Altogether, the HI found outside galaxies in this group amounts to several times $10^9$ M$_\odot$, at least 10% of the HI contained inside galaxies. This suggests a substantial flow of gas in and out of galaxies during the several billion years of the group's evolution.

We examine the effect of Galactic diffractive interstellar scintillation as a means of explaining the reported deficit of Fast Radio Burst (FRB) detections at low Galactic latitude. We model the unknown underlying FRB flux density distribution as a power law with a rate scaling as $S_\nu^{-5/2+\delta}$ and account for the fact that the FRBs are detected at unknown positions within the telescope beam. We find that the event rate of FRBs located off the Galactic plane may be enhanced by a factor ~30-300% relative to objects near the Galactic plane without necessarily affecting the slope of the distribution. For FRBs whose flux densities are subject to relatively weak diffractive scintillation, as is typical for events detected at high Galactic latitudes, we demonstrate that an effect associated with Eddington bias is responsible for the enhancement. The magnitude of the enhancement increases with the steepness of the underlying flux density distribution, so that existing limits on the disparity in event rates between high and low Galactic latitudes suggest that the FRB population has a steep differential flux density distribution, scaling as $S_\nu^{-3.5}$ or steeper. Existing estimates of the event rate in the flux density range probed by the High Time Resolution Universe (HTRU) survey overestimate the true rate by a factor of ~3.

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

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

We have detected sporadic, bright, short-duration radio pulses from PSR J0901$-$4624. These pulses are emitted simultaneously with persistent, periodic emission that dominates the flux density when averaging over many periods of the pulsar. The bright pulses have energies that are consistent with a power-law distribution. The integrated profile of PSR J0901$-$4624 is highly polarized and shows four distinct components. The bright pulses appear to originate near the magnetic pole of the pulsar and have polarization properties unlike that of the underlying emission at the same pulse phase. We conclude that the bright pulses represent a secondary giant-micropulse emission process, possibly from a different region in the pulsar magnetosphere.

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.

Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.

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

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

Several studies have shown the distribution of pulsars' magnetic inclination angles to be skewed towards low values compared with the distribution expected if the rotation and magnetic axes are placed randomly on the star. Here we focus on a sample of 28 $\gamma$-ray-detected pulsars using data taken as part of the Parkes telescope's \emph{FERMI} timing program. In doing so we find a preference in the sample for low magnetic inclination angles, $\alpha$, in stark contrast to both the expectation that the magnetic and rotation axes are orientated randomly at the birth of the pulsar and to $\gamma$-ray-emission-model-based expected biases. In this paper, after exploring potential explanations, we conclude that there are two possible causes of this preference, namely that low $\alpha$ values are intrinsic to the sample, or that the emission regions extend outside what is traditionally thought to be the open-field-line region in a way which is dependent on the magnetic inclination. Each possibility is expected to have important consequences, ranging from supernova physics to population studies of pulsars and considerations of the radio beaming fraction. We also present a simple conversion scheme between the observed and intrinsic magnetic inclinations which is valid under the assumption that the observed skew is not intrinsic and which can be applied to all existing measurements. We argue that extending the active field-line region will help to resolve the existing tension between emission geometries derived from radio polarisation measurements and those required to model $\gamma$-ray light curves.

We present radio intensity and polarisation profiles of 28 $\gamma$-ray-detected pulsars with the aim of putting constraints on their viewing geometries using data from the Parkes telescope. Constraints are formed both from the goodness-of-fit of the position angles to the Rotating Vector Model and from the beam opening angle considering aberration and retardation effects. Uncertainties on the relevant parameters are systematically taken into account in order to produce a more robust constraint, using a new approach. Surprisingly, we find that the distribution of the magnetic inclination angle ($\alpha$) in this subset of pulsars peaks at low values, contrary to the predictions of $\gamma$-ray models. We find a lack of correlation between these and a set of $\alpha$ values which were derived using $\gamma$-ray light curves, suggesting a problem in the interpretation of the data in one or both of these domains. Finally, we also show that the $\alpha$ distribution of pulsars with multiple radio components is no different to that of single-component pulsars.

This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interferometer, fitted with prototype signal processing hardware capable of forming at least 9 dual-polarisation beams simultaneously, allowing several square degrees to be imaged in a single pointed observation. The main purpose of the test array is to develop beamforming and wide-field calibration methods for use with the full telescope, but it will also be capable of limited early science demonstrations.

PSR J1717-4054 is one of the small class of pulsars which null on intermediate (~hour) timescales. Such pulsars represent an important link between classical nullers, whose emission vanishes for a few rotations, and the intermittent pulsars which null for months and years. Using the Parkes radio telescope and the Australia Telescope Compact Array, we have studied the emission from PSR J1717-4054 over intervals from single pulses to years. We have identified and characterised nulling at three discrete timescales: the pulsar emits during "active states" separated by nulls lasting many thousands of rotations, while active states themselves are interrupted by nulls with a bimodal distribution of durations; one to two rotations, or tens of rotations. We detect no subpulse structure or flux density variations during active states, and we place stringent limits (<0.1% of the mean active flux density) on pulsed emission during inactive states. Finally, our high-quality data have also allowed us to measure for the first time many important properties of PSR J1717-4054, including its position, spindown rate, spectrum, polarization characteristics, and pulse broadening in the interstellar medium.

Context: GeV gamma-ray pulsations from over 140 pulsars have been characterized using the Fermi Large Area Telescope, enabling improved understanding of the emission regions within the neutron star magnetospheres, and the contributions of pulsars to high energy electrons and diffuse gamma rays in the Milky Way. The first gamma-ray pulsars to be detected were the most intense and/or those with narrow pulses. Aims: As the Fermi mission progresses, progressively fainter objects can be studied. In addition to more distant pulsars (thus probing a larger volume of the Galaxy), or ones in high background regions (thus improving the sampling uniformity across the Galactic plane), we detect pulsars with broader pulses or lower luminosity. Adding pulsars to our catalog with inclination angles that are rare in the observed sample, and/or with lower spindown power, will reduce the bias in the currently known gamma-ray pulsar population. Methods: We use rotation ephemerides derived from radio observations to phase-fold gamma rays recorded by the Fermi Large Area Telescope, to then determine the pulse profile properties. Spectral analysis provides the luminosities and, when the signal-to-noise ratio allows, the cutoff energies. We constrain the pulsar distances by different means in order to minimize the luminosity uncertainties. Results: We present six new gamma-ray pulsars with an eclectic mix of properties. Three are young, and three are recycled. They include the farthest, the lowest power, two of the highest duty-cycle pulsars seen, and only the fourth young gamma-ray pulsar with a radio interpulse. We discuss the biases existing in the current gamma-ray pulsar catalog, and steps to be taken to mitigate the bias.

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

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

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$\alpha$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$\alpha$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$\gamma$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$\alpha$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

Debris disks and asteroid belts are expected to form around young pulsars due to fallback material from their original supernova explosions. Disk material may migrate inwards and interact with a pulsar's magnetosphere, causing changes in torque and emission. Long term monitoring of PSR J0738-4042 reveals both effects. The pulse shape changes multiple times between 1988 and 2012. The torque, inferred via the derivative of the rotational period, changes abruptly from September 2005. This change is accompanied by an emergent radio component that drifts with respect to the rest of the pulse. No known intrinsic pulsar processes can explain these timing and radio emission signatures. The data lead us to postulate that we are witnessing an encounter with an asteroid or in-falling debris from a disk.

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.

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

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

We have carried out observations of the newly-discovered magnetar in the direction of Sagittarius A* using the Australia Telescope Compact Array in four frequency bands from 4.5 to 20 GHz. Radio pulsations are clearly detected at all frequencies. We measure the pulsar's dispersion measure to be 1650 +/- 50 pc/cm^3, the highest of any of the known pulsars. Once Faraday rotation has been taken into account, the pulse profile is more than 80% linearly polarized at all frequencies and has a small degree (5%) of circular polarization. The rotation measure of -67000 +/- 500$ rad/m^2 is the largest (in magnitude) ever measured for a pulsar but still a factor 8 smaller than Sgr A* itself. The combination of the dispersion and rotation measures implies an integrated magnetic field strength of -50uG along the line of sight. The flux density appears to have increased by about a factor of two between observations made 30 days apart. This object therefore joins the small class of radio emitting magnetars.

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

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

Recent discoveries of highly dispersed millisecond radio bursts by Thornton et al. in a survey with the Parkes radio telescope at 1.4 GHz point towards an emerging population of sources at cosmological distances whose origin is currently unclear. Here we demonstrate that the scattering effects at lower radio frequencies are less than previously thought, and that the bursts could be detectable at redshifts out to about $z=0.5$ in surveys below 1 GHz. Using a source model in which the bursts are standard candles with bolometric luminosities $\sim 8 \times 10^{44}$ ergs/s uniformly distributed per unit comoving volume, we derive an expression for the observed peak flux density as a function of redshift and use this, together with the rate estimates found by Thornton et al. to find an empirical relationship between event rate and redshift probed by a given survey. The non-detection of any such events in Arecibo 1.4 GHz survey data by Deneva et al., and the Allen Telescope Array survey by Simeon et al. is consistent with our model. Ongoing surveys in the 1--2 GHz band should result in further discoveries. At lower frequencies, assuming a typical radio spectral index $\alpha=-1.4$, the predicted peak flux densities are 10s of Jy. As a result, surveys of such a population with current facilities would not necessarily be sensitivity limited and could be carried out with small arrays to maximize the sky coverage. We predict that sources may already be present in 350-MHz surveys with the Green Bank Telescope. Surveys at 150 MHz with 30 deg$^2$ fields of view could detect one source per hour above 30 Jy.

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

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep observations to search for radio pulsations in the directions of unidentified Fermi Large Area Telescope (LAT) gamma-ray sources. We report the discovery of an eclipsing black-widow millisecond pulsar, PSR J1544+4937, identified with the un-cataloged gamma-ray source Fermi J1544.2+4941. This 2.16 ms pulsar is in a 2.9 hours compact circular orbit with a very low-mass companion (Mc > 0.017 Msun). At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Variations in the eclipse ingress phase are observed, indicating a clumpy and variable eclipsing medium. Moreover, additional short-duration absorption events are observed around the eclipse boundaries. Using the radio timing ephemeris we were able to detect gamma-ray pulsations from this pulsar, confirming it as the source powering the gamma-ray emission.

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

PSR J0742-2822 exhibits two distinct emission states that are identified by discrete changes in the observed pulse profile. These changes have previously been shown to correlate with changes in the derivative of the pulse frequency. In this paper, we use observations with the Parkes radio telescope at a centre frequency of 1369 and 3100 MHz to produce high phase resolution polarisation profiles for the two modes and perform a detailed study of the correlation between observed pulse shape and spin-down rate. We find no correlation for at least 200 days prior to a glitch in the pulsar at MJD 55022, following which the correlation becomes strong. This suggests a link between the emission state switching phenomenon and glitch events. We discuss the possibility that emission state switching is driven by the interaction between the magnetosphere and the interior of the neutron star.

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.

This paper presents a detailed investigation of the dependence of pulsar spin-velocity alignment, which has been observed for a sample of 58 pulsars, on pulsar age. At first, our study considers only pulsar characteristic ages, resulting in no change in the degree of correlation as a function of age, up to at least 100 Myr. Subsequently, we consider a more reliable estimate of pulsar age, the kinematic age, assuming that pulsars are born near the Galactic plane. We derive kinematic ages for 52 pulsars, based on the measured pulsar proper motions and positions, by modelling the trajectory of the pulsars in a Galactic potential. The sample of 52 pulsar kinematic ages constitutes the largest number of independently estimated pulsar ages to date. Using only the 33 most reliable kinematic ages from our simulations, we revisit the evolution of spin- velocity alignment, this time as a function of kinematic age. We find that the strong correlation seen in young pulsars is completely smeared out for pulsars with kinematic ages above 10 Myr, a length of time beyond which we expect the gravitational pull of the Galaxy to have a significant effect on the directions of pulsar velocities. In the discussion, we investigate the impact of large distance uncertainties on the reliability of the calculated kinematic ages. Furthermore, we present a detailed investigation of the implications of our revised pulsar ages for the braking-index and birth-period distributions. Finally, we discuss the predictions of various SN-kick mechanisms and their compatibility with our results.

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have performed a targeted radio pulsar survey of 289 unassociated gamma-ray sources discovered by the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL catalogue (Abdo et al., 2010). This survey resulted in the discovery of millisecond pulsar J1745+1017, which resides in a short-period binary system with a low-mass companion, Mmin ~ 0.0137 Msun, indicative of `Black Widow' type systems. A two-year timing campaign has produced a refined radio ephemeris, accurate enough to allow for phase-folding of the LAT photons, resulting in the detection of a dual-peaked gamma-ray light-curve, proving that PSR J1745+1017 is the source responsible for the gamma-ray emission seen in 1FGL J1745.5 + 1018 (2FGL J1745.6+1015; Nolan et al., 2012). We find the gamma-ray spectrum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to contain a further 10 newly discovered pulsars which were undetected in this survey. Our radio observations of these sources are discussed and in all cases limiting flux densities are calculated. The reasons behind the seemingly low yield of discoveries are also discussed.

We report the detection of radio emission from PSR J1311-3430, the first millisecond pulsar discovered in a blind search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio pulsations at 2 GHz, visible for <10% of ~4.5-hrs of observations using the Green Bank Telescope (GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nancay, and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint detection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferometric imaging observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311-3430, is not radio quiet and provides additional evidence that the radio beaming fraction of millisecond pulsars is very large. The radio detection yields a distance estimate of 1.4 kpc for the system, yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay in the radio pulsar as the pulsar appears from eclipse and we speculate on possible mechanisms for the non-detections of the pulse at other orbital phases and observing frequencies.

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.

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

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

The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered by young and energetic pulsars. HESS J1857+026 is a spatially extended gamma-ray source detected by H.E.S.S. and classified as a possible pulsar wind nebula candidate powered by PSR J1856+0245. We search for gamma-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1.5 GHz, we phase-fold 36 months of gamma-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete gamma-ray spectral and morphological analysis. No gamma-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident with the TeV source HESS J1857+026. The gamma-ray spectrum is well described by a simple power-law with a spectral index of 1.53 \pm 0.11_{\rm stat} \pm 0.55_{\rm syst}$ and an energy flux of $G(0.1$--100 GeV$)=(2.71 \pm 0.52_{\rm stat} \pm 1.51_{\rm syst}) \times 10^{-11}$ ergs cm$^{-2}$ s$^{-1}$. The $\gamma$-ray luminosity is $L_{PWN}^{\gamma} (0.1$--100 GeV$)=(2.5 \pm 0.5_{stat} \pm 1.5_{syst}) \times 10^{35} (\frac{d}{9 kpc})^2$ ergs s$^{-1}$, assuming a distance of 9 kpc. This implies a $\gamma-$ray efficiency of $\sim$ 5% for $\dot{E}=4.6 \times 10^{36}$ erg $s^{-1}$, in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new constraints on its pulsar wind nebula nature.

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

The reported alignment between the projected spin-axes and proper motion directions of pulsars is revisited in the light of new data from Jodrell Bank and Effelsberg. The present investigation uses 54 pulsars, the largest to date sample of pulsars with proper-motion and absolute polarisation, to study this effect. Our study has found strong evidence for pulsar spin-velocity alignment, excluding that those two vectors are completely uncorrelated, with >99% confidence. Although we cannot exclude the possibility of orthogonal spin-velocity configurations, comparison of the data with simulations shows that the scenario of aligned vectors is more likely than that of the orthogonal case. Moreover, we have determined the spread of velocities that a spin-aligned and spin-orthogonal distribution of kicks must have to produce the observed distribution of spin-velocity angle offsets. If the observed distribution of spin-velocity offset angles is the result of spin-aligned kicks, then we find that the distribution of kick-velocity directions must be broad with {\sigma}_v~30\degree if the orthogonal-kick scenario is assumed, then the velocity distribution is much narrower with {\sigma}_v<10\degree. Finally, in contrast to previous studies, we have performed robustness tests on our data, in order to determine whether our conclusions are the result of a statistical and/or systematic bias. The conclusion of a correlation between the spin and velocity vectors is independent of a bias introduced by subsets in the total sample. Moreover, we estimate that the observed alignment is robust to within 10% systematic uncertainties on the determination of the spin-axis direction from polarisation data.

We present a detailed single pulse study of PSR B1055-52 based on observations at the Parkes radio telescope. The radio emission is found to have a complex modulation dominated by a periodicity of ~20 times its rotational period P (0.197s), whose phase and strength depends on pulse longitude. This periodicity exhibits a phase-locked delay of about 2.5P between the main pulse (MP) and interpulse (IP), presumed to be the opposite poles of the pulsar. This delay corresponds to a light travel distance of many times the light cylinder radius. More complex modulations are found within the MP on timescales down to about 9P, and both these and the principal modulation vary strongly across the (at least) 7 components which the MP and IP exhibit. The nature of the single pulse emission, which ranges from smooth and longitudinally extended to `spiky', is also component-dependent. Despite these disparities, the total pulse intensity distributions at the MP and IP are virtually identical in shape, suggesting a common emission mechanism. In an attempt to account for the complex modulations we examine a number of physical models, both intrinsic (which presuppose the pulsar to be an isolated neutron star) and extrinsic (appealing to the presence of circumstellar material to modulate the emission). Significant objections can be made to each model, so this pulsar's behaviour patterns remain a crucial challenge to theorists.

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

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

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

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

In a search for radio pulsations from the magnetar 1E 1841-045, we have discovered the unrelated pulsar J1841-0500, with rotation period P=0.9 s and characteristic age 0.4 Myr. One year after discovery with the Parkes telescope at 3 GHz, radio emission ceased from this bright pulsar. After 580 days, emission resumed as before. The P-dot during both on states is 250% of the average in the off state. PSR J1841-0500 is a second example of an extremely intermittent pulsar, although with a much longer off period and larger ratio of spin-down rates than PSR B1931+24. The new pulsar is hugely scattered by the ISM, with a fitted timescale referenced to 1 GHz of tau_1=2 s. Based on polarimetric observations at 5 GHz with the Green Bank Telescope, the intrinsic pulse profile has not obviously changed between the two on states observed so far, although relatively small variations cannot be excluded. The magnitude of its rotation measure is the largest known, RM=-3000 rad/m^2, and with a dispersion measure DM=532 pc/cc implies a large electron-weighted average magnetic field strength along the line of sight, 7 microG.

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi-LAT sources, we have discovered the middle-aged pulsar J2030+3641, associated with 1FGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the entire Fermi mission. With a rotation period of 0.2 s, spin-down luminosity of 3e34 erg/s, and characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM=246 pc/cc. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive -- PSR J2030+3641 would have been found blindly in gamma rays if only >0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.

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

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

PSR B1259-63 is a 48 ms pulsar in a highly eccentric 3.4 year orbit around the young massive star LS 2883. During the periastron passage the system displays transient non-thermal unpulsed emission from radio to very high energy gamma rays. It is one of the three galactic binary systems clearly detected at TeV energies, together with LS 5039 and LS I +61 303. We observed PSR B1259-63 after the 2007 periastron passage with the Australian Long Baseline Array at 2.3 GHz to trace the milliarcsecond (mas) structure of the source at three different epochs. We have discovered extended and variable radio structure. The peak of the radio emission is detected outside the binary system near periastron, at projected distances of 10-20 mas (25-45 AU assuming a distance of 2.3 kpc). The total extent of the emission is ~50 mas (~120 AU). This is the first observational evidence that non-accreting pulsars orbiting massive stars can produce variable extended radio emission at AU scales. Similar structures are also seen in LS 5039 and LS I +61 303, in which the nature of the compact object is unknown. The discovery presented here for the young non-accreting pulsar PSR B1259-63 reinforces the link with these two sources and supports the presence of pulsars in these systems as well. A simple kinematical model considering only a spherical stellar wind can approximately trace the extended structures if the binary system orbit has a longitude of the ascending node of omega ~ -40 deg and a magnetization parameter of sigma~0.005.

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.

Most pulsars observed by the Fermi LAT have gamma-ray luminosities scaling with spindown power Edot as L_gamma (Edot x 10^33 erg/s)^{1/2}. However, there exist one detection and several upper limits an order of magnitude or more fainter than this trend. We describe these `sub-luminous' gamma-ray pulsars, and discuss the case for this being an orientation effect. Of the 12 known young radio pulsars with Edot>10^34 erg/s and d<2kpc several are substantially sub-luminous. The limited available geometrical constraints favor aligned geometries for these pulsars, although no one case for alignment is compelling. In this scenario GeV emission detected from such sub-luminous pulsars can be due to a lower altitude, lower-power accelerator gap.

We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind nebula HESS J1825-137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of sigma = 0.56{\deg} $\pm$ 0.07{\deg} for an assumed Gaussian model. The 1 - 100 GeV LAT spectrum of this source is well described by a power-law with a spectral index of 1.38 $\pm$ 0.12 $\pm$ 0.16 and an integral flux above 1 GeV of (6.50 $\pm$ 0.21 $\pm$ 3.90) x 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV-TeV emission in a cooling-limited nebula powered by the pulsar PSR J1826-1334.

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

Here we describe the Compact Array Broadband Backend (CABB) and present first results obtained with the upgraded Australia Telescope Compact Array (ATCA). The 16-fold increase in observing bandwidth, from 2 x 128 MHz to 2 x 2048 MHz, high bit sampling, and addition of 16 zoom windows (each divided into a further 2048 channels) provide major improvements for all ATCA observations. The benefits of the new system are: (1) hugely increased radio continuum and polarization sensitivity as well as image fidelity, (2) substantially improved capability to search for and map emission and absorption lines over large velocity ranges, (3) simultaneous multi-line and continuum observations, (4) increased sensitivity, survey speed and dynamic range due to high-bit sampling, and (5) high velocity resolution, while maintaining full polarization output. The new CABB system encourages all observers to make use of both spectral line and continuum data to achieve their full potential. Given the dramatic increase of the ATCA capabilities in all bands (ranging from 1.1 to 105 GHz) CABB enables scientific projects that were not feasible before the upgrade, such as simultaneous observations of multiple spectral lines, on-the-fly mapping, fast follow-up of radio transients (e.g., the radio afterglow of new supernovae) and maser observations at high velocity resolution and full polarization. The first science results presented here include wide-band spectra, high dynamic-range images, and polarization measurements, highlighting the increased capability and discovery potential of the ATCA.

One of the tenets of the radio pulsar observational picture is that the integrated pulse profiles are constant with time. This assumption underpins much of the fantastic science made possible via pulsar timing. Over the past few years, however, this assumption has come under question with a number of pulsars showing pulse shape changes on a range of timescales. Here, we show the dramatic appearance of a bright component in the pulse profile of PSR J0738-4042 (B0736-40). The component arises on the leading edge of the profile. It was not present in 2004 but strongly present in 2006 and all observations thereafter. A subsequent search through the literature shows the additional component varies in flux density over timescales of decades. We show that the polarization properties of the transient component are consistent with the picture of competing orthogonal polarization modes. Faced with the general problem of identifying and characterising average profile changes, we outline and apply a statistical technique based on a Hidden Markov Model. The value of this technique is established through simulations, and is shown to work successfully in the case of low signal-to-noise profiles.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

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

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

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

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

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 searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

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

We demonstrate that the high-magnetic field pulsar J1119-6127 exhibits three different types of behaviour in the radio band. Trailing the "normal" profile peak there is an "intermittent" peak and these components are flanked by two additional components showing very erratic "RRAT-like" emission. Both the intermittent and RRAT-like events are extremely rare and are preceded by a large amplitude glitch in the spin-down parameters. The post-glitch spin-down rate is smaller than the pre-glitch rate. This type of relaxation is very unusual for the pulsar population as a whole, but is observed in the glitch recovery of a RRAT. The abnormal emission behaviour in PSR J1119-6127 was observed up to three months after the epoch of the large glitch, suggestive of changes in the magnetospheric conditions during the fast part of the recovery process. We argue that both the anomalous recoveries and the emission changes could be related to reconfigurations of the magnetic field. Apart from the glitches, the spin-down of PSR J1119-6127 is relatively stable, allowing us to refine the measurement of the braking index (n=2.684\pm0.002) using more than 12 years of timing data. The properties of this pulsar are discussed in light of the growing evidence that RRATs do not form a distinct class of pulsar, but rather are a combination of different extreme emission types seen in other neutron stars. Different sub-classes of the RRATs can potentially be separated by calculating the lower limit on the modulation index of their emission. We speculate that if the abnormal behaviour in PSR J1119-6127 is indeed glitch induced then there might exist a population of neutron stars which only become visible in the radio band for a short duration in the immediate aftermath of glitch activity. These neutron stars will be visible in the radio band as sources that only emit some clustered pulses every so many years.

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

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

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

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

We report radio polarization observations of G319.9-0.7 (MSC 319.9-0.7) at 3 and 6 cm obtained with the Australia Telescope Compact Array. The source shows a highly elongated morphology with the energetic pulsar J1509-5850 located at the tip. We found a flat radio spectrum of index \alpha=-0.26 +/- 0.04 and a high degree of linear polarization. These results confirm G319.9-0.7 as a bow-shock pulsar wind nebula. The polarization maps suggest a helical magnetic field trailing the pulsar, with the symmetry axis parallel to the system's inferred direction of motion. This is the first time such a field geometry has been seen in a bow-shock nebula, and it may be the result of an alignment between the pulsar spin axis and its space velocity. Compared to other bow-shock examples, G319.9-0.7 exhibits very different properties in the field structure and surface brightness distribution, illustrating the large diversity of the population.

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

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.

We present polarisation observations of five pulsars whose profiles exhibit two distinct emission regions separated by close to 180 degrees of longitude. We fitted the position angle of the linear polarisation using the rotating vector model and convincingly show that all the pulsars have the angle between their magnetic and rotation axes close to 90 degrees. The simplest interpretation of the results is that we see `main pulse' emission from one pole and `interpulse' emission from the opposite pole. We have attempted to produce emission maps of the magnetosphere above the polar caps for each pulsar and find that the maps support the view that the emission region in pulsars is complex, even when the profile appears simple. For three pulsars, we can derive emission heights and polar maps which are consistent with emission regions located symmetrically about the magnetic axis and confined to the open field lines. For two pulsars, we find that either the emission arises from `closed' field lines or that the profiles are highly asymmetric with respect to the magnetic axis.

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.

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

We have detected significant Rotation Measure variations for 9 bright pulsars, as a function of pulse longitude. An additional sample of 10 pulsars showed a rather constant RM with phase, yet a small degree of RM fluctuation is visible in at least 3 of those cases. In all cases, we have found that the rotation of the polarization position angle across our 1.4 GHz observing band is consistent with the wavelength-squared law of interstellar Faraday Rotation. We provide for the first time convincing evidence that RM variations across the pulse are largely due to interstellar scattering, although we cannot exclude that magnetospheric Faraday Rotation may still have a minor contribution; alternative explanations of this phenomenon, like erroneous de-dispersion and the presence of non-orthogonal polarization modes, are excluded. If the observed, phase-resolved RM variations are common amongst pulsars, then many of the previously measured pulsar RMs may be in error by as much as a few tens of rad m-2.

Using gamma-ray data collected by the Astrorivelatore Gamma ad Immagini LEggero (AGILE) satellite over a period of almost one year (from 2007 July to 2008 June), we searched for pulsed signals from 35 potentially interesting radio pulsars, ordered according to $F_{\gamma}\propto \sqrt{\dot{E}} d^{-2}$ and for which contemporary or recent radio data were available. AGILE detected three new top-ranking nearby and Vela-like pulsars with good confidence both through timing and spatial analysis. Among the newcomers we find pulsars with very high rotational energy losses, such as the remarkable PSR B1509-58 with a magnetic field in excess of 10^13 Gauss, and PSR J2229+6114 providing a reliable identification for the previously unidentified EGRET source 3EG 2227+6122. Moreover, the powerful millisecond pulsar B1821-24, in the globular cluster M28, is detected during a fraction of the observations. Four other promising gamma-ray pulsar candidates, among which is the notable J2043+2740 with an age in excess of 1 million years, show a possible detection in the timing analysis only and deserve confirmation.

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 that will be 50 times more sensitive than any existing radio facility. Most of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from a few hundred MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is a technology demonstrator aimed in the mid-frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phased-array feed systems on parabolic reflectors. The large field-of-view makes ASKAP an unprecedented synoptic telescope that will make substantial advances in SKA key science. 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 two sites selected by the international community as a potential location for the SKA. In this paper, we outline the ASKAP project and summarise its headline science goals as defined by the community at large.

We have simulated a population of young spin-powered pulsars and computed the beaming pattern and lightcurves for the three main geometrical models: polar cap emission, two-pole caustic ("slot gap") emission and outer magnetosphere emission. The light curve shapes depend sensitively on the magnetic inclination alpha and viewing angle zeta. We present the results as maps of observables such as peak multiplicity and gamma-ray peak separation in the (alpha, zeta) plane. These diagrams can be used to locate allowed regions for radio-loud and radio-quiet pulsars and to convert observed fluxes to true all-sky emission.

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

Multi-wavelength studies at radio, infrared, optical, X-ray, and TeV wavelengths have discovered probable counterparts to many Galactic sources of GeV emission detected by EGRET. These include pulsar wind nebulae, high mass X-ray binaries, and mixed morphology supernova remnants. Here we provide an overview of the observational properties of Galactic sources which emit across 19 orders of magnitude in energy. We also present new observations of several sources.

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

In this paper we compare the characteristics of pulsars with a high spin-down energy loss rate (Edot) against those with a low Edot. We show that the differences in the total intensity pulse morphology between the two classes are in general rather subtle. A much more significant difference is the fractional polarization which is very high for high Edot pulsars and low for low Edot pulsars. The Edot at the transition is very similar to the death line predicted for curvature radiation. This suggests a possible link between high energy and radio emission in pulsars and could imply that gamma-ray efficiency is correlated with the degree of linear polarization in the radio band. The degree of circular polarization is in general higher in the second component of doubles, which is possibly caused by the effect of co-rotation on the curvature of the field lines in the inertial observer frame. The most direct link between the high energy emission and the radio emission could be the sub-group of pulsars which we call the energetic wide beam pulsars. These young pulsars have very wide profiles with steep edges and are likely to be emitted from a single magnetic pole. The similarities with the high energy profiles suggest that both types of emission are produced at the same extended height range in the magnetosphere. Alternatively, the beams of the energetic wide beam pulsars could be magnified by propagation effects in the magnetosphere. This would naturally lead to decoupling of the wave modes, which could explain the high degree of linear polarization. As part of this study, we have discovered three previous unknown interpulse pulsars (and we detected one for the first time at 20 cm). We also obtained rotation measures for 18 pulsars whose values had not previously been measured.

The observed fraction of pulsars with interpulses, their period distribution and the observed pulse width versus pulse period correlation is shown to be inconsistent with a model in which the angle alpha between the magnetic axis and the rotation axis is random. This conclusion appears to be unavoidable, even when non-circular beams are considered. Allowing the magnetic axis to align from a random distribution at birth with a timescale of 70 Myrs years can, however, explain those observations well. The timescale derived is consistent with that obtained via independent methods. The probability that a pulsar beam intersects the line of sight is a function of the angle alpha and therefore beam evolution has important consequences for evolutionary models and for estimations of the total number of neutron stars. The validity of the standard formula for the spin-down rate, which is independent of alpha appears to be questionable.

High time resolution observations of PSR B0906-49 (or PSR J0908-4913) over a wide range of frequencies have enabled us to determine the geometry and beam shape of the pulsar. We have used the position angle traverse to determine highly-constrained solutions to the rotating vector model which show conclusively that PSR B0906-49 is an orthogonal rotator. The accuracy obtained in measuring the geometry is unprecedented. This may allow tests of high-energy emission models, should the pulsar be detected with GLAST. Although the impact parameter, beta, appears to be frequency dependent, we have shown that this is due to the effect of interstellar scattering. As a result, this pulsar provides some of the strongest evidence yet that the position angle swing is indeed related to a geometrical origin, at least for non-recycled pulsars. We show that the beam structures of the main pulse and interpulse in PSR B0906-49 are remarkably similar. The emission comes from a height of ~230 km and is consistent with originating in a patchy cone located about half way to the last open field lines. The rotation axis and direction of motion of the pulsar appear to be aligned.

We present a new survey for pulsars in the error boxes of the low-latitude EGRET sources 3EG J1027-5817, 3EG J1800-2338 and 3EG J1810-1032. Although all of these sources have been covered by previous pulsar surveys, the recent discovery of the young, energetic pulsar PSR J1410-6132 at 6.7 GHz has shown that pulsars of this type can be hidden from low frequency surveys. Using an observing frequency of 3.1 GHz we discovered a 91-ms pulsar, PSR J1028-5819, which observations made at the Parkes telescope and the Australia Telescope Compact Array have shown to be young and energetic. We believe this pulsar is likely to be powering the unidentified EGRET source 3EG J1027-5817. Like other energetic pulsars, PSR J1028-5819 is highly linearly polarised, but astonishingly has a pulse duty cycle of only 0.4%, one of the smallest in the entire pulsar catalogue.

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

We have observed a total of 67 pulsars at five frequencies ranging from 243 to 3100 MHz. Observations at the lower frequencies were made at the Giant Metre Wave Telescope in India and those at higher frequencies at the Parkes Telescope in Australia. We present profiles from 34 of the sample with the best signal to noise ratio and the least scattering. The general `rules' of pulsar profiles are seen in the data; profiles get narrower, the polarization fraction declines and outer components become more prominent as the frequency increases. Many counterexamples to these rules are also observed, and pulsars with complex profiles are especially prone to rule breaking. We hypothesise that the location of pulsar emission within the magnetosphere evolves with time as the the pulsar spins down. In highly energetic pulsars, the emission comes from a confined range of high altitudes, in the middle range of spin down energies the emission occurs over a wide range of altitudes whereas in pulsars with low spin-down energies it is confined to low down in the magnetosphere.

We measured a sample of 150 pulsar Rotation Measures (RMs) using the 20-cm receiver of the Parkes 64-m radio telescope. 46 of the pulsars in our sample have not had their RM values previously published, whereas 104 pulsar RMs have been revised. We used a novel quadratic fitting algorithm to obtain an accurate RM from the calibrated polarisation profiles recorded across 256 MHz of receiver bandwidth. The new data are used in conjunction with previously known Dispersion Measures (DMs) and the NE2001 electron-density model to study models of the direction and magnitude of the Galactic magnetic field.

We have investigated the radio emission from the anomalous X-ray pulsar 1E 1547.0-5408 (PSR J1550-5418) using the Parkes telescope and the Australia Telescope Compact Array. The flux density of the pulsar is roughly the same between 1.4 and 45 GHz, but shows time variability. The radiation is nearly 100% linearly polarized between frequencies of 45 and 3.2 GHz, but from 2.3 to 1.4 GHz it becomes increasingly more depolarized. The rotation measure of -1860 rad/m^2 is the largest for any known pulsar, and implies an average magnetic field strength along the line of sight of 2.7 microG. The pulse profiles are circularly polarized at all frequencies observed, more so at lower frequencies, at the ~15% level. The observed swing of the position angle of linear polarization as a function of pulse phase suggests that in this neutron star the rotation and magnetic axes are nearly aligned, and that its radio emission is only detectable within a small solid angle. Timing measurements indicate that the period derivative of this 2 s pulsar has increased by nearly 40% in a 6-month period. The flat spectrum and variability in flux density and pulse profiles are reminiscent of the properties of XTE J1810-197, the only other known radio-emitting magnetar, and are anomalous by comparison with those of ordinary radio pulsars.

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 that will be 50 times more sensitive than any existing radio facility. Most of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from a few hundred MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is a technology demonstrator aimed in the mid-frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phased-array feed systems on parabolic reflectors. The large field-of-view makes ASKAP an unprecedented synoptic telescope that will make substantial advances in SKA key science. 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 two sites selected by the international community as a potential location for the SKA. In this paper, we outline an ambitious science program for ASKAP, examining key science such as understanding the evolution, formation and population of galaxies including our own, understanding the magnetic Universe, revealing the transient radio sky and searching for gravitational waves.

We have conducted observations of 22 pulsars at frequencies of 0.7, 1.4 and 3.1 GHz and present their polarization profiles. The observations were carried out for two main purposes. First we compare the orientation of the spin and velocity vectors to verify the proposed alignment of these vectors by Johnston et al. (2005). We find, for the 14 pulsars for which we were able to determine both vectors, that 7 are plausibly aligned, a fraction which is lower than, but consistent with, earlier measurements. Secondly, we use profiles obtained simultaneously at widely spaced frequencies to compute the radio emission heights. We find, similar to other workers in the field, that radiation from the centre of the profile originates from lower in the magnetosphere than the radiation from the outer parts of the profile.

Motivated by recent results on the location of the radio emission in pulsar magnetospheres, we have developed a model which can account for the large diversity found in the average profile shapes of pulsars. At the centre of our model lies the idea that radio emission at a particular frequency arises from a wide range of altitudes above the surface of the star and that it is confined to a region close to the last open field lines. We assert that the radial height range over which emission occurs is responsible for the complex average pulse shapes rather than the transverse (longitudinal) range proposed in most current models. By implementing an abrupt change in the height range to discriminate between young, short-period, highly-energetic pulsars and their older counterparts, we obtain the observed transition between the simple and complex average pulse profiles observed in each group respectively. Monte Carlo simulations are used to demonstrate the match of our model to real observations.

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

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

We revisit the phenomenon of pulse nulling using high-quality single-pulse data of PSR B1133+16 from simultaneous multifrequency observations. Observations were made at 325, 610, 1400 and 4850 MHz as part of a joint program between the European Pulsar Network and the Giant Metrewave Radio Telescope. The pulse energy time series are analysed to derive improved statistics of nulling pulses as well as to investigate the frequency dependence of the phenomenon. The pulsar is observed to be in null state for approximately 15% of the time; however, we find that nulling does not always occur simultaneously at all four frequencies of observation. We characterise this "selective nulling'' as a function of frequency, separation in frequency, and combination of frequencies. The most remarkable case is a significantly large number of nulls ($\approx$6%) at lower frequencies, that are marked by the presence of a fairly narrow emission feature at the highest frequency of 4850 MHz. We refer to these as "low frequency (LF) nulls." Our analysis shows that this high frequency emission tends to occur preferentially over a narrow range in longitude and with pulse widths typically of the order of a few milliseconds. We discuss the implications of our results for the pulsar emission mechanism in general and for the broadbandness of nulling phenomenon in particular. Our results signify the presence of an additional process of emission which does not turn off when the pulsar nulls at low frequencies, and becomes more prominent at higher frequencies. Our analysis also hints at a possible outer gap origin for this new population of pulses, and thus a likely connection to some high-energy emission processes that occur in the outer parts of the pulsar magnetosphere.

We present the results of a first campaign of radio continuum observations toward supergiants with spectral types in the range O8 - B3. Three targets out of twelve were detected with the ATCA and the VLA. The measured fluxes and the upper limits yielded values of the stellar mass-loss rates dM/dt and the wind efficiencies. The comparison of predicted and derived values of dM/dt shows a good agreement.

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.

We present polarization data for 32 mainly southern pulsars at 8.4 GHz. The observations show that the polarization fraction is low in most pulsars at this frequency except for the young, energetic pulsars which continue to show polarization fractions in excess of 60 per cent. All the pulsars in the sample show evidence for conal emission with only one third also showing core emission. Many profiles are asymmetric, with either the leading or the trailing part of cone not detectable. Somewhat surprisingly, the asymmetric profiles tend to be more polarized than the symmetrical profiles. Little or no pulse narrowing is seen between 1 and 8.4 GHz. The spectral behaviour of the orthogonal polarization modes and radius to frequency mapping can likely account for much of the observational phenomenology. Highly polarized components may orginate from higher in the magnetosphere than unpolarized components.

We have observed the supernova remnant (SNR) G290.1-0.8 in the 21-cm HI line and the 20-cm radio continuum using the Australia Telescope Compact Array (ATCA). The HI data were combined with data from the Southern Galactic Plane Survey to recover the shortest spatial frequencies. In contrast, HI absorption was analyzed by filtering extended HI emission, with spatial frequencies shorter than 1.1 k-lambda. The low-resolution ATCA radio continuum image of the remnant shows considerable internal structure, resembling a network of filaments across its 13 arcmin diameter. A high-resolution ATCA radio continuum image was also constructed to study the small scale structure in the SNR. It shows that there are no structures smaller than ~17", except perhaps for a bright knot to the south, which is probably an unrelated object. The HI absorption study shows that the gas distribution and kinematics in front of SNR G290.1-0.8 are complex. We estimate that the SNR probably lies in the Carina arm, at a distance 7 (+/- 1) kpc. In addition, we have studied nearby sources in the observed field using archival multiwavelength data to determine their characteristics.

We present polarization profiles at 1.4 and 3.1 GHz for 14 young pulsars with characteristic ages less than 75 kyr. Careful calibration ensures that the absolute position angle of the linearly polarized radiation at the pulsar is obtained. In combination with previously published data we draw three main conclusions about the pulse profiles of young pulsars. (1) Pulse profiles are simple and consist of either one or two prominent components. (2) The linearly polarized fraction is nearly always in excess of 70 per cent. (3) In profiles with two components the trailing component nearly always dominates, only the trailing component shows circular polarization and the position angle swing is generally flat across the leading component and steep across the trailing component. Based on these results we can make the following generalisations about the emission beams of young pulsars. (1) There is a single, relatively wide cone of emission from near the last open field lines. (2) Core emission is absent or rather weak. (3) The height of the emission is between 1 and 10 per cent of the light cylinder radius.

PSR B1259-63 is in a highly eccentric 3.4 year orbit with a Be star and crosses the Be star disc twice per orbit, just prior to and just after periastron. Unpulsed radio, X-ray and gamma-ray emission observed from the binary system is thought to be due to the collision of pulsar wind with the wind of Be star. We present here the results of new XMM-Newton observations of the PSR B1259-63 system during the beginning of 2004 as the pulsar approached the disc of Be star.We combine these results with earlier unpublished X-ray data from BeppoSAX and XMM-Newton as well as with ASCA data. The detailed X-ray lightcurve of the system shows that the pulsar passes (twice per orbit) through a well-defined gaussian-profile disk with the half-opening angle (projected on the pulsar orbit plane) ~18.5 deg. The intersection of the disk middle plane with the pulsar orbital plane is inclined at ~70 deg to the major axis of the pulsar orbit. Comparing the X-ray lightcurve to the TeV lightcurve of the system we find that the increase of the TeV flux some 10--100 days after the periastron passage is unambiguously related to the disk passage. At the moment of entrance to the disk the X-ray photon index hardens from 1.8 up to 1.2 before returning to the steeper value 1.5. Such behaviour is not easily accounted for by the model in which the X-ray emission is synchrotron emission from the shocked pulsar wind. We argue that the observed hardening of the X-ray spectrum is due to the inverse Compton or bremsstrahlung emission from 10-100 MeV electrons responsible for the radio synchrotron emission.

Aim: We observed the center of the supernova remnant Vela Jr in radio continuum in order to search for a counterpart to the compact central X-ray source CXOU J085201.4-461753, possibly a neutron star candidate which could be the remnant of the supernova explosion. Method: Observations were made with the Australia Telescope Compact Array at 13 and 20 cm. Spectral indices were obtained using flux density correlations of the data which were spatially filtered to have the same u-v coverage. A multiwavelength search for counterparts to the compact central X-ray source was made. Results: We compiled a new catalogue of 31 small diameter radio sources, including the previously known source PMN J0853-4620, listing the integrated flux densities at 20 cm and, for half of the sources, the flux densities at 13 cm with the corresponding spectral indices. All sources are unresolved at the present angular resolution except for Source 18, which is clearly elongated and lies strikingly close to CXOU J085201.4-461753. Our observations show no evidence for the existence of a pulsar wind driven nebula associated with the point X-ray source. Furthermore, Source 18 has a thermal spectrum with index +0.8 +/- 0.4, and appears to be the counterpart of the optical source Wray 16-30. In spite of the absence of OIII emission lines as reported in the literature, we find that this object could be explained as a low emission planetary nebula belonging to the ``butterfly'' morphological class. Conclusions: We conclude that if the radio source 18 is actually a planetary nebula, then CXOU J085201.4-461753 is more likely to be related to it rather than to Vela Jr.

We present strong observational evidence for a relationship between the direction of a pulsar's motion and its rotation axis. We show carefully calibrated polarization data for 25 pulsars, 20 of which display linearly polarized emission from the pulse longitude at closest approach to the magnetic pole. Such data allow determination of the position angle of the linear polarisation which in turn reflects the position angle of the rotation axis. Of these 20 pulsars, 10 show an offset between the velocity vector and the polarisation position angle which is either less than 10\degr or more than 80\degr, a fraction which is very unlikely by random chance. We believe that the bimodal nature of the distribution arises from the presence of orthogonal polarisation modes in the pulsar radio emission. In some cases this orthogonal ambiguity is resolved by observations at other wavelengths so that we conclude that the velocity vector and the rotation axis are aligned at birth. Strengthening the case is the fact that 4 of the 5 pulsars with ages less than 3 Myr show this relationship, including the Vela pulsar. We discuss the implications of these findings in the context of the Spruit & Phinney (1998)\nocite{sp98} model of pulsar birth-kicks. We point out that, contrary to claims in the literature, observations of double neutron star systems do not rule out aligned kick models and describe a possible observational test involving the double pulsar system.

The high-energy gamma-ray light curve of the binary system PSR B1259 -63, is computed using the approach that successfully predicted the spectrum at periastron. The simultaneous INTEGRAL and H.E.S.S. spectra taken 16 days after periastron currently permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. In this paper we shown how the long-term light curve may help to lift this degeneracy.

We present here a direct comparison of the polarisation position angle (PA) profiles of 17 pulsars, observed at 1.4 and 3.1 GHz. Absolute PAs are obtained at each frequency, permitting a measurement of the difference in the profiles. By doing this, we obtain more precise rotation measure (RM) values for some of the pulsars in the current catalogue. We find that, apart from RM corrections, there are small, pulse longitude dependent differences in PA with frequency. Such differences go beyond the interpretation of a geometrical origin. We describe in detail the PA evolution between the two frequencies and discuss possible causes, such as orthogonal and non-orthogonal polarisation modes of emission. We also use the PA and total power profiles to estimate the difference in emission height at which the two frequencies originate. In our data sample, there are changes in the relative strengths of different pulse components, especially overlapping linearly polarised components, which coincide with intrinsic changes of the PA profile, resulting in interesting PA differences between the two frequencies.

The high-energy gamma-ray emission discovered using the H.E.S.S. telescopes from the binary system PSR B1259 -63, is modelled using an extension of the approach that successfully predicted it. We find that the simultaneous INTEGRAL and H.E.S.S. data permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. Additional, simultaneous, X-ray and TeV data sets are needed to lift this degeneracy.

The sheer number of pulsars discovered by the SKA, in combination with the exceptional timing precision it can provide, will revolutionize the field of pulsar astrophysics. The SKA will provide a complete census of pulsars in both the Galaxy and in Galactic globular clusters that can be used to provide a detailed map of the electron density and magnetic fields, the dynamics of the systems, and their evolutionary histories. This complete census will provide examples of nearly every possible outcome of the evolution of massive stars, including the discovery of very exotic systems such as pulsar black-hole systems and sub-millisecond pulsars, if they exist. These exotic systems will allow unique tests of the strong field limit of relativistic gravity and the equation of state at extreme densities. Masses of pulsars and their binary companions -- planets, white dwarfs, other neutron stars, and black holes -- will be determined to $\sim 1$% for hundreds of objects. With the SKA we can discover and time highly-stable millisecond pulsars that comprise a pulsar-timing array for the detection of low-frequency gravitational waves. The SKA will also provide partial censuses of nearby galaxies through periodicity and single-pulse detections, yielding important information on the intergalactic medium.

Stimulated emission of radiation has not been directly observed in astrophysical situations up to this time. Here we demonstrate that photons from pulsar B1641-45 stimulate pulses of excess 1720 MHz line emission in an interstellar OH cloud. As this stimulated emission is driven by the pulsar, it varies on a few millisecond timescale, orders of magnitude shorter than the quickest OH maser variations previously detected. Our 1612 MHz spectra are inverted copies of the 1720 MHz spectra. This "conjugate line" phenomenon enables us to constrain the properties of the interstellar OH line-producing gas. We also show that pulsar signals suffer significantly deeper OH absorption than do other background sources; confirming earlier tentative findings that OH clouds are clumpier on small scales than neutral hydrogen clouds.

PSR J0537-6910 is a young, energetic, rotation-powered X-ray pulsar with a spin period of 16 ms located in the Large Magellanic Cloud. We have searched for previously undetected radio pulsations (both giant and standard) from this pulsar in a 12-hour observation taken at 1400 MHz with the Parkes 64-m radio telescope. The very large value of the magnetic field at the light cylinder radius suggests that this pulsar might be emitting giant radio pulses like those seen in other pulsars with similar field strengths. No radio emission of either kind was detected from the pulsar, and we have established an upper limit of ~ 25 mJy kpc^2 for the average 1400-MHz radio luminosity of PSR J0537-6910. The 5-sigma single-pulse detection threshold was ~ 750 mJy for a single 80-us sample. These limits are likely to be the best obtainable until searches with greatly improved sensitivity can be made with next-generation radio instruments.

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 present a high-resolution (24 arcsec) study of the HI interstellar gas distribution around the radio-quiet neutron star Geminga. Based on Very Large Array (VLA) and MPIfR Effelsberg telescope data, we analyzed a 40' x 40' field around Geminga. These observations have revealed the presence of a neutral gas shell, 0.4 pc in radius, with an associated HI mass of 0.8 Msun, which surrounds Geminga at a radial velocity compatible with the kinematical distance of the neutron star. In addition, morphological agreement is observed between the internal face of the HI shell and the brightest structure of Geminga's tail observed in X-rays.We explore the possibility that this morphological agreement is the result of a physical association.

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.

The sensitivity of the SKA enables a number of tests of theories of gravity. A Galactic Census of pulsars will discover most of the active pulsars in the Galaxy beamed toward us. In this census will almost certainly be pulsar-black hole binaries as well as pulsars orbiting the super-massive black hole in the Galactic centre. These systems are unique in their capability to probe the ultra-strong field limit of relativistic gravity. These measurements can be used to test the Cosmic Censorship Conjecture and the No-Hair theorem. The large number of millisecond pulsars discovered with the SKA will also provide a dense array of precision clocks on the sky. These clocks will act as the multiple arms of a huge gravitational wave detector, which can be used to detect and measure the stochastic cosmological gravitational wave background that is expected from a number of sources.

PSR B0540-69 is one of a small handful of pulsars known to emit giant pulses and the only extra-galactic pulsar known to do so. We observed the pulsar for a total of 72 hr over a 6 month interval and detected 141 giant pulses. We have obtained correct phasing between the radio arrival times of the giant pulses and the X-ray pulse profile. The giant pulses occur in two phase windows, located 6.7 ms before and 5.0 ms after the midpoint of the X-ray profile. We have detected the integrated profile of the pulsar at 1.4 GHz with a flux density of only 24 micro-Jy. The statistics of the giant pulses are clearly power-law and it is likely that the giant pulses contribute only a few per cent to the integrated pulse flux density, similar to other pulsars with giant pulse emission. Simultaneous X-ray and radio observations show no significant increase in the X-ray flux at the time of the radio giants. The relative enhancements in the X-ray emission must be at least 5.5 times smaller than in the radio.

In most pulsars, the circularly polarized component, Stokes $V$, is weak in the average pulse profiles. By forming the average profile of $|V|$ from single pulses we can distinguish between pulsars where $V$ is weak in the individual pulses and those where large $V$ of variable handedness is observed from one pulse to the other. We show how $|V|$ profiles depend on the signal-to-noise ratio of $V$ in the single pulses and demonstrate that it is possible to simulate the observed, broad distributions of $V$ by assuming a model where $|V|$ is distributed around a mean value and the handedness of $V$ is permitted to change randomly. The $|V|$ enhanced profiles of 13 pulsars are shown, 5 observed at 1.41 GHz and 8 observed at 4.85 GHz, to complement the set in Karastergiou et al. (2003b). It is argued that the degree of circular polarization in the single pulses is related to the orthogonal polarization mode phenomenon and not to the classification of the pulse components as cone or core.

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.

We report on the first polarimetric observations at 19 GHz made with the upgraded Australia Telescope Compact Array. Observations were made of the Galactic supernova remnant (SNR) G328.8+0.2. We find the SNR has circular morphology with a strong central bar, similar to that seen at lower frequencies. The SNR has high linear polarization throughout, with fractional polarization in the bar up to 50 per cent. The orientation of the magnetic field lines follow the filamentary structure of the SNR. The magnetic field at the edge of the SNR is generally toroidal, interspersed with radial fingers, likely caused by Rayleigh-Taylor instabilities. Although the SNR has been identified as Crab-like, we prefer an interpretation in which the bar is a pulsar powered wind nebula with the rest of the SNR consisting of the shell. The proposed pulsar parameters make the SNR / pulsar system more like SNR G11.2--0.3 than the Crab Nebula.

The integrated profile of PSR B1641-45 can be decomposed into four Gaussian components. The sum of these Gaussian components can be made to replicate the total intensity, linear and circular polarization of the integrated profile, and the discontinuity in the swing of the position angle. Surprisingly, the single pulses from PSR B1641-45 can also be decomposed into these same Gaussian components with only the amplitude of the components as a free parameter. The distributions of the fitted amplitudes are log-normal. Under the assumption that emission consists of 100% polarized orthogonal modes which are emitted simultaneously, I show that the polarization properties of the single pulses can also be replicated in a simple way. This lends support to models involving the superposition of orthogonal modes rather than disjoint emission of the modes.

We have carried out an HI survey towards X-ray central compact objects (CCOs) inside supernova remnants (SNRs) which shows that many of them are placed within local HI minima. The nature of these minima is not clear, but the most likely explanation is that the CCOs have evacuated the neighboring gas. This survey also allowed us to detect a weak, diffuse radio nebula inside the SNR Vela Jr, probably created by the winds of its associated CCO.

We have carried out a study of the neutral hydrogen in the direction of the X-ray source 1E 161348-5055, a compact central object (CCO) located in the interior of the supernova remnant (SNR) RCW 103. The HI 21 cm line observations were carried out using the Australia Telescope Compact Array, complemented with single dish data from the Parkes radio telescope to recover information at all spatial scales. We derive a distance to RCW 103 of 3.3 kpc, in agreement with previous distance measurements. We have also detected a small hole in the HI emission which is positionally and kinematically coincident with the location of the CCO which confirms the association between the SNR and the CCO. This is the third case of a depression in HI emission seemingly associated with CCOs in SNRs. The characteristic parameters of the holes such as their size, eccentricity and evacuated mass are similar in all three cases. We estimate the absorbing HI column density towards 1E 161348-5055 to be ~6 x 10^{21} cm^{-2}, a value compatible with a blackbody solution for the CCO X-ray emission. However, the implied brightness temperature is very high comparedto most neutron stars. Moreover, the strong long-term variability in X-rays favours the hypothesis that 1E 161348-5055 is an accreting binary sourcerather than an isolated, cooling neutron star. An analysis of the continuum image obtained at 1.4 GHz from these observations shows no trace of a pulsar wind nebula around 1E 161348-5055, in spite of it being a young object.

We present a study of single pulses from nine bright northern pulsars to investigate the behaviour of circular polarisation, V. The observations were conducted with the Effelsberg 100-m radio telescope at 1.41 GHz and 4.85 GHz and the Westerbork radio telescope at 352 MHz. For the first time, we present the average profile of the absolute circular polarisation |V| in the single pulses. We demonstrate that the average profile of |V| is the distinguishing feature between pulse components that exhibit low V in the single pulses and components that exhibit high V of either handedness, despite both cases resulting in a low mean. We also show that the |V| average profile remains virtually constant with frequency, which is not generally the case for V, leading us to the conclusion that |V| is a key quantity in the pulsar emission problem.

We have observed the HI distribution in an area 40' x 40' around the neutron star candidate RX J0822-4300, which is located in the supernova remnant Puppis A. The observations of the 21 cm line were obtained with the Australia Telescope Compact Array (ATCA) and were combined with single dish data from the Southern Galactic Plane Survey. The spatial resolution is 90", and the velocity resolution, 1 km/s . A sensitivity of ~ 0.7 K was achieved. The results revealed a double lobed feature of reduced emission at +16 km/s, centered on the central compact object (CCO), and aligned with an HI hole blueshifted by 13 km/s. The HI depressions have probably been created by the sweeping up of ~ 2 solar masses. The alignement between the lobes and the optical expansion centre of Puppis A suggests that the CCO could be ejecting two opposite jets. The velocity at which the two lobes are best defined allowed us to confirm that the distance to Puppis A is 2.2 kpc, based on a systemic velocity of +16 km/s. The hydrogen column density computed using this systemic velocity is consistent with estimates from models for X-ray spectra, thus reinforcing our conclusion that the kinematic distance is 2.2 k pc.

In this paper we demonstrate that a large, unexplored reservoir of information about pulsar emission exists, that is directly linked to the radiating particles and their radiation process: We present a study of flux density measurements of individual pulses simultaneously observed at four different frequencies. Correcting for effects caused by the interstellar medium, we derive intrinsic flux density spectra of individual radio pulses observed at several frequencies for the first time. Pulsar B0329+54 was observed at 238, 626, 1412 and 4850 MHz, while observations of PSR B1133+16 were made at 341, 626, 1412 and 4850 MHz. We derive intrinsic pulse-to-pulse modulation indices which show a minimum around 1 GHz. Correlations between the flux densities of different frequency pairs worsen as the frequency separation widens and also tend to be worse for outer profile components. The single pulse spectra of PSR B0329+54 resemble the spectra of the integrated profile. However, the spectral index distributions for the single pulses of PSR B1133+16 show significant deviations from a Gaussian. This asymmetry is caused by very strong pulses with flux densities exceeding the mean value by more than a factor of ten. These strong pulses occur preferentially at the trailing edge of the leading component and appear to be broadband in most cases. Their properties are similar to those of so-called giant pulses, suggesting that these phenomena are related.

We report the discovery of the first giant pulses from an extragalactic radio pulsar. Observations of PSR B0540-69 in the Large Magellanic Cloud made with the Parkes radio telescope at 1.38 GHz show single pulses with energy more than 5000 times that of the average pulse energy. This is only the second young pulsar, after the Crab, to show giant pulse emission. Similar to the Crab pulsar, the giant pulses occur in two distinct phase ranges and have significant arrival time jitter within these ranges. The location of the giant pulses appears to lag the peak of the (sinusoidal) X-ray profile by 0.37 and 0.64 phase, although absolute timing between the radio and X-ray data is not yet secure. The dispersion measure of the giant pulses is 146.5 pc/cc, in agreement with the detection of the pulsar at 0.64 GHz by Manchester et al. (1993). The giant pulses are scatter broadened at 1.4 GHz with an exponential scattering time of 0.4 ms and have an emission bandwidth of at least 256 MHz. In 8 hr of integration we have failed to detect any integrated flux density from the pulsar to a level of 13 microJy, assuming a duty cycle of 10%. This implies the spectral index between 0.64 and 1.38 GHz is steeper than -4.4.

The variability of the Vela pulsar (PSR B0833-45) corresponds to well-defined field statistics that vary with pulsar phase, ranging from Gaussian intensity statistics off-pulse to approximately power-law statistics in a transition region and then lognormal statistics on-pulse, excluding giant micropulses. These data are analyzed here in terms of two superposed wave populations, using a new calculation for the amplitude statistics of two vectorially-combined transverse fields. Detailed analyses show that the approximately power-law and lognormal distributions observed are fitted well at essentially all on-pulse phases by Gaussian-lognormal and double-lognormal combinations, respectively. These good fits, plus the smooth but significant variations in fit parameters across the source, provide strong evidence that the approximately power-law statistics observed in the transition region are not intrinsic. Instead, the data are consistent with normal pulsar emission having lognormal statistics at all phases. This is consistent with generation in an inhomogeneous source obeying stochastic growth theory (SGT) and with the emission mechanism being purely linear (either direct or indirect). A nonlinear mechanism is viable only if it produces lognormal statistics when suitably ensemble-averaged. Variations in the SGT fit parameters with phase imply that the radiation is relatively more variable near the pulse edges than near the center, as found in earlier work. In contrast, Vela's giant micropulses come from a very restricted phase range and have power-law statistics with indices ($6.7 \pm 0.6$) not inconsistent with nonlinear wave collapse. These results imply that normal pulses have a different source and generation mechanism than giant micropulses, as suggested previously on other grounds.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations do not exist for this variability or for the mechanism producing the radiation. The variability corresponds to the field statistics, particularly the distribution of wave field amplitudes, which are predicted by theories for wave growth in inhomogeneous media. This paper shows that the field statistics of the Vela pulsar (PSR B0833-45) are well-defined and vary as a function of pulse phase, evolving from Gaussian intensity statistics off-pulse to approximately power-law and then lognormal distributions near the pulse peak to approximately power-law and eventually Gaussian statistics off-pulse again. Detailed single-component fits confirm that the variability corresponds to lognormal statistics near the peak of the pulse profile and Gaussian intensity statistics off-pulse. The lognormal field statistics observed are consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. The simplest interpretations are that the pulsar's variability is a direct manifestation of an SGT state and the emission mechanism is linear (either direct or indirect), with no evidence for nonlinear mechanisms like modulational instability and wave collapse which produce power-law field statistics. Stringent constraints are placed on nonlinear mechanisms: they must produce lognormal statistics when suitably ensemble-averaged. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.

We investigate circular polarization in pulsar radio emission through simultaneous observations of PSR B1133+16 at two frequencies. In particular, we investigate the association of the handedness of circular polarization with the orthogonal polarization mode phenomenon at two different frequencies. We find the association to be significant across the pulse for PSR B1133+16, making a strong case for orthogonal polarization modes determining the observed circular polarization. The association however is not perfect and decreases with frequency. Based on these results and assuming emission occurs in superposed orthogonal polarization modes, we present a technique of mode decomposition based on single pulses. Average profiles of the polarization of each mode can then be computed by adding the individual mode-separated single pulses. We show that decomposing single pulses produces different average profiles for the orthogonal polarization modes from decomposing average profiles. Finally, we show sample single pulses and discuss the implications of the frequency dependence of the correlation of the circular polarization with the orthogonal polarization mode phenomenon.

Radio observations of two isolated neutron stars, RXJ0720.4-3125 and RX J0806.4-4132, have been made with the Australia Telescope Compact Array at a frequency of 1.4 GHz. No continuum emission is detected from either object with a 3 sigma upper limit of 0.2 mJy. The data were also folded synchronously with the known rotation periods of 8.4 and 11.4 s respectively. No pulsed emission was detected. If the pulse duty cycle is small, the upper limit on pulsed emission can be reduced still further to 0.04 mJy. The best evidence seems to indicate that the isolated neutron stars detected by ROSAT all-sky survey are relatively young objects, born with a very high magnetic field.

We have measured 21-cm absorption spectra in the direction of three southern pulsars, PSRs B0736-40, B1451-68 and B1641-45, in three separate epochs spread over 2.5 yr. We see no evidence for any changes in the absorption spectra over this time span in spite of good velocity resolution and sensitivity. Towards PSR B1641-45 we place an upper limit of 10^19 per cm^2 on the change in the column density of the cold, neutral gas. In addition, we observed PSR B1557-50 and compared its HI absorption spectrum with spectra taken in 1980 and 1994. A prominent deep absorption feature seen at -110 km/s in 1994 is weaker in 2000 and likely was not present at all in 1980. Using the standard interpretation which links the distance traversed by the pulsar to a physical cloud size, this results in a cloud size of ~1000 AU and density of 2x10^4 per cc, parameters typical of those seen in other observations. These results are also consistent with an alternative model by Deshpande (2000) who expects the largest variations in optical depth to be seen against the longest time intervals.

We have made radio continuum, HI and X-ray observations in the direction of the unidentified EGRET source 3EG J1410-6147, using the Australia Telescope Compact Array and the Chandra X-ray Observatory. The observations encompass the supernova remnant (SNR) G312.4-0.4 and the two young pulsars PSRs J1412-6145 and J1413-6141. We derive a lower distance limit of 6 kpc to the SNR, although interpretation of positive velocity features in the HI spectrum may imply the SNR is more distant than 14 kpc. PSR J1412-6145, with an age of 50 kyr, is the pulsar most likely associated with SNR G312.4-0.4. X-rays are not detected from either pulsar and diffuse X-ray emission near the bright western edge of the SNR is weak. Although there is circumstantial evidence that this western region is a pulsar wind nebula (PWN), the embedded pulsar PSR J1412-6145 is apparently not sufficiently powerful to explain the radio enhancement. The origin of the electron acceleration in this region and of the gamma-rays remain unidentified, unless the distance to PSR J1413-6141 is at least a factor of 3 lower than its dispersion measure distance.

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.

In this paper, we present a study of orthogonal polarization modes in the radio emission of PSR B1133+16, conducted within the frame of simultaneous, multi-frequency, single-pulse observations. Simultaneously observing at two frequencies (1.41 GHz and 4.85 GHz) provides the means to study the bandwidth of polarization features such as the polarization position angle. We find two main results. First, that there is a high degree of correlation between the polarization modes at the two frequencies. Secondly, the modes occur more equally and the fractional linear polarization decreases towards higher frequencies. We discuss this frequency evolution and propose propagation effects in the pulsar magnetosphere as its origin.

Using the formalism of Hamaker et al. (1996), I derive a method for the polarization calibration of observations made with a single radio telescope. This method is particularly appropriate for observations of pulsars, where the sign and magnitude of the circular polarization are useful for understanding the emission processes at work. I apply the method to observations of PSR J1359-6038 made using the multibeam receiver on the Parkes radio telescope.

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

We have carried out a survey for `giant pulses' in 6 young, Vela-like pulsars. In no cases did we find single pulses with flux densities more than 10 times the mean flux density. However, in PSR B1706-44 we have detected giant micro-pulses very similar to those seen in the Vela pulsar. In PSR B1706-44 these giant micro-pulses appear on the trailing edge of the profile and have an intrinsic width of ~1 ms. The cumulative probability distribution of their intensities is best described by a power-law. If the power-law continues to higher intensities, then 3.7x10^6 rotations are required to obtain a pulse with 20 times the mean pulse flux. This number is similar to the giant pulse rate in PSR B1937+21 and PSR B1821-24 but significantly higher than that for the Crab.

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to pulse, from the average profile, and from phase to phase within a pulse. Widely accepted explanations for pulsar radio emission and its time variability do not exist. Here, by analysing data near the peak of the Vela pulsar's average profile, we show that Vela's variability corresponds to lognormal field statistics, consistent with the prediction of stochastic growth theory (SGT) for a purely linear system close to marginal stability. Vela's variability is therefore a direct manifestation of an SGT state and the field statistics constrain the emission mechanism to be linear (either direct or indirect), ruling out nonlinear mechanisms like wave collapse. Field statistics are thus a powerful, potentially widely applicable tool for understanding variability and constraining mechanisms and source characteristics of coherent astrophysical and space emissions.

We present X-ray, radio, and infrared observations of the 68 ms pulsar PSR J1420-6048 and its surrounding nebula, a possible counterpart of the gamma-ray source GeV J1417-6100/3EG J1420-6038. Pulsed X-ray emission at the radio period is marginally detected by ASCA from a source embedded in the hard spectrum X-ray nebula AX J1420.1-6049. At radio wavelengths, the pulsar is found to be strongly linearly and circularly polarized, and the polarization sweep is measured. A comparison of high resolution ATCA radio imaging of the Kookaburra's upper wing (G313.6+0.3), which contains the pulsar and the X-ray nebula, with infrared images suggests the radio emission is partly non-thermal.

We have carried out a survey for `giant pulses' in 5 millisecond pulsars. We detect individual pulses from the high Edot pulsar PSR B1821-24 with energies exceeding 50x the mean pulse energy. These giant pulses are concentrated in a narrow phase window coincident with the power-law non-thermal pulse seen in hard X-rays. This is the third example of the giant pulse phenomenon. It supports the idea that large B fields in the outer magnetosphere are critical to the formation of such pulses and further suggests a direct connection between giant pulses and high energy emission.

We present X-ray and radio studies of sources which are brightabove 1 GeV (F_{>1GeV} > 4e-8 ph/cm^2/s. Only 11 out of ~30 of these gamma-ray sources have been identified with lower energy counterparts: 5 blazars and 6 pulsars. Three of these pulsars are surrounded by radio pulsar wind nebulae (PWN), two of which are also seen as bright, extended X-ray synchrotron nebulae. The ASCA X-ray telescope has observed 28 of the bright GeV sources, revealing an excess of F_{2-10keV} > 10e-12 ergs/cm^2/s sources within the {\it EGRET} error contours of the unidentified sources. Although several supernova remnants are positionally coincident with these sources, we find no X-ray evidence of high energy particle production in SNR shell shocks consistent with the GeV positions. We also present initial results from follow on radio imaging studies of several fields containing unidentified sources. We have discovered new X-ray/radio nebulae in three of these fields which are strong candidates for PWN. These sources, along with a similar nebula in CTA 1 and the PWN around PSR B1853+01 in W44, are all positionally coincident with variable EGRET sources. This suggests a class of variable gamma-ray sources associated with synchrotron emitting regions powered by the winds of young pulsars.

Scintillation observations of the nearby single millisecond pulsar (MSP) PSR J0030+0451 were carried out with the Parkes 64m radiotelescope at three different epochs in 1999. From analysis of the dynamic spectrum we obtained the amplitude of the electron density power spectrum log C_N^2 ~= -3.33$ and a scintillation velocity V_iss <~ 15 km s^-1. This result shows that the Shklovskii effect on the spin-down rate \dot P is negligible. We also performed a correlation analysis between pulsar proper motions (V_pm) and scintillation velocities (V_iss) using updated measurements for a sample of 77 objects, 17 of which are MSPs. The full sample shows a correlation coefficient r_s ~= 80% at an extremely high significance level, while for the MSP sub-sample (excluding 2 outliers) we obtain r_s ~= 90%.

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

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

We have measured 21-cm absorption and emission spectra in the direction of a further 7 southern pulsars with the Parkes telescope to derive their kinematic distances and to study the interstellar medium. For the first time we have successfully obtained HI absorption measurements for PSRs J1602-5100, J1740-3015 and J1745-3040. We have also significantly improved the sensitivity and resolution on PSRs J1600-5044, J1752-2806, and J1825-0935, whose spectra have previously been measured, and have corrected an error in the published distance to PSR J1824-1945. Since the Frail & Weisberg summary of pulsar distances in 1990, a further 23 pulsars now have measured HI distances, mainly through the efforts of the current group. We discuss the Taylor & Cordes electron density model in light of these new measurements and find that, although the model towards the Galactic Centre appears good, the line of sight through the Carina spiral arm is poorly fit by the model.

We report on a sensitive survey for radio pulsar wind nebulae (PWN) towards 27 energetic and/or high velocity pulsars. Observations were carried out at 1.4 GHz using the Very Large Array and the Australia Telescope Compact Array, and utilised pulsar-gating to search for off-pulse emission. These observing parameters resulted in a considerably more sensitive search than previous surveys, and could detect PWN over a much wider range of spatial scales (and hence ambient densities and pulsar velocities). However, no emission clearly corresponding to a PWN was discovered. Based on these non-detections we argue that the young and energetic pulsars in our sample have winds typical of young pulsars, but produce unobservable PWN because they reside in low density (n approx 0.003 cm^-3) regions of the ISM. However, non-detections of PWN around older and less energetic pulsars can only be explained if the radio luminosity of their winds is less than 1e-5 of their spin-down luminosity, implying an efficiency at least an order of magnitude smaller than that seen for young pulsars.

We have measured the properties of the diffractive scintillation toward the Vela pulsar under the extremely strong scattering conditions encountered at 660 MHz. We obtain a decorrelation bandwidth of $\nu_d = 244 \pm 4$ Hz and diffractive decorrelation timescale of $t_{\rm diff} = 3.3\pm 0.3$ s. Our measurement of the modulation indices $m=0.87\pm 0.003\pm 0.05$ and $m=0.93\pm 0.03 \pm 0.05$ (one for each polarization stream), are at variance with the modulation index of the Vela pulsar obtained at 2.3 GHz by Gwinn et al. (1997) {\it if} the deviation from a modulation index of unity is ascribed to a source size effect.

We summarise the proceedings of the SRCfTA workshop on ``Supernova Remnants, Pulsars and the Interstellar Medium'' that was held at the University of Sydney on Mar 18 and 19, 1999.

Using the standard equation for the slowdown of a neutron star, we derive a formula for the braking index via integration rather than the conventional differentiation. The new formula negates the need to measure the second time derivative of the rotation frequency. We show that the method gives similar braking indices for PSR B1509-58 and the Crab pulsar to those already in the literature. We point out that our method is useful for obtaining the braking indices of moderate aged pulsars without the need for long, phase-connected timing solutions. We applied the method to 20 pulsars and discuss the implications of the results. We find that virtually all the derived braking indices are dominated by the effects of (unseen) glitches, the recovery from which corrupts the value of the frequency first derivative. However, any real, large, positive braking index has implications for magnetic field decay and offers support to recent models of pulsar evolution.

Using the Australia Telescope Compact Array (ATCA) we have imaged the fields around five promising pulsar candidates to search for radio pulsar wind nebulae (PWNe). We have used the ATCA in its pulsar gating mode; this enables an image to be formed containing only off-pulse visibilities, thereby dramatically improving the sensitivity to any underlying PWN. Data from the Molonglo Observatory Synthesis Telescope were also used to provide sensitivity on larger spatial scales. This survey found a faint new PWN around PSR B0906-49; here we report on non-detections of PWNe towards PSRs B1046-58, B1055-52, B1610-50 and J1105-6107. Our radio observations of the field around PSR B1055-52 argue against previous claims of an extended X-ray and radio PWNe associated with the pulsar. If these pulsars power unseen, compact radio PWN, upper limits on the radio flux indicate that less than 1e-6 of their spin-down energy is used to power this emission. Alternatively PSR B1046-58 and PSR B1610-50 may have relativistic winds similar to other young pulsars and the unseen PWN is resolved and fainter than our surface brightness sensitivity threshold. We can then determine upper limits on the local ISM density of 2.2e-3 cm^-3 and 1e-2 cm^-3, respectively. Furthermore we constrain the spatial velocities of these pulsars to be less than ~450 km/s and thus rule out the association of PSR B1610-50 with SNR G332.4+00.1 (Kes 32). Strong limits on the ratio of unpulsed to pulsed emission are also determined for three pulsars.

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 discuss the interpretation of transient, unpulsed radio emission detected from the unique pulsar/Be-star binary system PSR B1259-63. Extensive monitoring of the 1994 and 1997 periastron passages has shown that the source flares over a 100-day interval around periastron, varying on time-scales as short as a day and peaking at 60 mJy (~100 times the apastron flux density) at 1.4 GHz. Interpreting the emission as synchrotron radiation, we show that (i) the observed variations in flux density are too large to be caused by the shock interaction between the pulsar wind and an isotropic, radiatively driven, Be-star wind, and (ii) the radio emitting electrons do not originate from the pulsar wind. We argue instead that the radio electrons originate from the circumstellar disk of the Be star and are accelerated at two epochs, one before and one after periastron, when the pulsar passes through the disk. A simple model incorporating two epochs of impulsive acceleration followed by synchrotron cooling reproduces the essential features of the radio light curve and spectrum and is consistent with the system geometry inferred from pulsed radio data.

We have mapped the radio emission in the error ellipse of GeV J1417-6100 (2EGS J1418-6049) at 13cm and 20cm using the Australia Telescope Compact Array. We find a large shell with extended wings, at the edge of which is a non-thermal, polarized structure with a center filled morphology (the `Rabbit'), coincident with an extended, hard X-ray source. We discuss the various sources seen within the ellipse as potential counterparts of the gamma-ray source. We conclude that the most likely scenario is that the Rabbit is a wind nebula surrounding a radio-quiet gamma-ray pulsar.

Previous searches for radio pulsar wind nebulae (PWN) suggest that only the youngest and most energetic pulsars power PWN. However the selection effects associated with such searches are severe, in that emission from a faint compact PWN will be masked by that from its associated pulsar. This has motivated us to search for radio PWN using the technique of pulsar gating, which allows us to ``switch off'' the pulsed emission and see what is hidden underneath. This search has resulted in the detection of the faintest radio PWN yet discovered, while non-detections around other pulsars have implications for previous claims of nebular emission and for associations with nearby supernova remnants.

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.

We report on Australia Telescope Compact Array observations of the ~10^5 yr old pulsar PSR B0906-49. In an image containing only off-pulse emission, we find a weak, slightly extended source coincident with the pulsar's position, which we argue is best interpreted as a pulsar wind nebula (PWN). A trail of emission extending behind the pulsar aligns with the major axis of the PWN, and implies that the pulsar is moving north-west with projected velocity ~60 km/s, consistent with its scintillation speed. The consequent density we infer for the pulsar's environment is >2 cm^(-3), so that the PWN around PSR B0906-49 is confined mainly by the high density of its surroundings rather than by the pulsar's velocity. Other properties of the system such as the PWN's low luminosity and apparent steep spectrum, and the pulsar's large characteristic age, lead us to suggest that this nebula is substantially different from other radio PWNe, and may represent a transition between young pulsars with prominent radio PWNe and older pulsars for which no radio PWN has been detected. We recommend that further searches for radio PWNe should be made as here: at low frequencies and with the pulsed emission subtracted.

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

We collate the evidence for rotation-powered neutron stars that are visible as X-ray sources and not as radio pulsars. Quantitative upper limits exist for the radio fluxes of eight objects, revealing a population at least an order of magnitude less luminous at radio wavelengths than known radio pulsars of similar power or age. These objects could have intrinsically low radio luminosities, but a simple alternative explanation is that they are pulsars in which the radio beams are directed away from Earth. Where hard X-ray or gamma-ray fluxes are seen, the beaming explanation implies different emission sites for the high-energy radiation and the unseen radio beams. We estimate that (a) the radio beaming fraction of young pulsars is roughly 50%, certainly much less than 100%, and (b) one neutron star is born in the Galaxy every ~90 years. We conclude that probably all neutron stars are born as radio pulsars and that most young, nearby pulsars have already been discovered.

PSR B1259-63 is in a highly eccentric 3.4 yr orbit around the Be-star SS2883. The system is located in the direction of the Sagittarius-Carina spiral arm at a distance of $\sim$1.5 kpc. We have made scintillation observations of the pulsar far from periastron at 4.8 GHz and 8.4 GHz, determining the diffractive bandwidth, $\Delta \nu_{\rm d}$, and timescale, $\Delta t_{\rm d}$, at both frequencies. We find no dependency on orbital phase until within 30 days of periastron. The data indicate that the scintillation is caused, not by the circumstellar environment, but by an HII region within the Sagittarius-Carina spiral arm located at least three-quarters of the way to the pulsar. Close to periastron, when the line-of-sight to the pulsar intersects the disk of the Be star, the electron densities within the disk are sufficient to overcome the `lever-arm effect' and produce a reduction in the scintillation bandwidth by six orders of magnitude.