The investigation of pulsars between millimetre and optical wavelengths is challenging due to the faintness of the pulsar signals and the relative low sensitivity of the available facilities compared to 100-m class telescopes operating in the centimetre band. The Kinetic Inductance Detector (KID) technology offers large instantaneous bandwidths and a high sensitivity that can help to substantially increase the ability of existing observatories at short wavelengths to detect pulsars and transient emission. To investigate the feasibility of detecting pulsars with KIDs, we observed the anomalous X-ray pulsar XTE J1810-197 with the New IRAM KIDs Array-2 (NIKA2) camera installed at the IRAM 30-m Telescope in Spain. We detected the pulsations from the pulsar with NIKA2 at its two operating frequency bands, 150 and 260 GHz ($\lambda$=2.0 and 1.15 mm, respectively). This is the first time that a pulsar is detected with a receiver based on KID technology in the millimetre band. In addition, this is the first report of short millimetre emission from XTE J1810-197 after its reactivation in December 2018, and it is the first time that the source is detected at 260 GHz, which gives us new insights into the radio emission process of the star.

Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

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

We report on simultaneous radio and X-ray observations of the repeating fast radio burst source FRB 180916.J0158+65 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME), Effelsberg, and Deep Space Network (DSS-14 and DSS-63) radio telescopes and the Chandra X-ray Observatory. During 33 ks of Chandra observations, we detect no radio bursts in overlapping Effelsberg or Deep Space Network observations and a single radio burst during CHIME/FRB source transits. We detect no X-ray events in excess of the background during the Chandra observations. These non-detections imply a 5-$\sigma$ limit of $<5\times10^{-10}$ erg cm$^{-2}$ for the 0.5--10 keV fluence of prompt emission at the time of the radio burst and $1.3\times10^{-9}$ erg cm$^{-2}$ at any time during the Chandra observations at the position of FRB 180916.J0158+65. Given the host-galaxy redshift of FRB 180916.J0158+65 ($z\sim0.034$), these correspond to energy limits of $<1.6\times10^{45}$ erg and $<4\times10^{45}$ erg, respectively. We also place a 5-$\sigma$ limit of $<8\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ on the 0.5--10\,keV absorbed flux of a persistent source at the location of FRB 180916.J0158+65. This corresponds to a luminosity limit of $<2\times10^{40}$ erg s$^{-1}$. Using Fermi/GBM data we search for prompt gamma-ray emission at the time of radio bursts from FRB 180916.J0158+65 and find no significant bursts, placing a limit of $4\times10^{-9}$ erg cm$^{-2}$ on the 10--100 keV fluence. We also search Fermi/LAT data for periodic modulation of the gamma-ray brightness at the 16.35-day period of radio-burst activity and detect no significant modulation. We compare these deep limits to the predictions of various fast radio burst models, but conclude that similar X-ray constraints on a closer fast radio burst source would be needed to strongly constrain theory.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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.

We study the timing stability of three black widow pulsars, both in terms of their long-term spin evolution and their shorter-term orbital stability. The erratic timing behaviour and radio eclipses of the first two black widow pulsar systems discovered (PSRs B1957+20 and J2051$-$0827) was assumed to be representative for this class of pulsars. With several new black widow systems added to this population in the last decade, there are now several systems known that do not show these typical orbital variations or radio eclipses. We present timing solutions using 7$-$8 yrs of observations from four of the European Pulsar Timing Array telescopes for PSRs\,J0023+0923, J2214+3000 and J2234+0944, and confirm that two of these systems do not show any significant orbital variability over our observing time span, both in terms of secular or orbital parameters. The third pulsar PSR J0023+0923 shows orbital variability and we discuss the implications for the timing solution. Our results from the long-term timing of these pulsars provide several new or improved parameters compared to earlier works. We discuss our results regarding the stability of these pulsars, and the stability of the class of black widow pulsars in general, in the context of the binary parameters, and discuss the potential of the Roche-lobe filling factor of the companion star being an indicator for stability of these systems.

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

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

After spending almost a decade in a radio-quiet state, the Anomalous X-ray Pulsar XTE J1810-197 turned back on in early December 2018. We have observed this radio magnetar at 1.5 GHz with ~daily cadence since the first detection of radio re-activation on 8 December 2018. In this paper, we report on the current timing properties of XTE J1810-197 and find that the magnitude of the spin frequency derivative has increased by a factor of 2.6 over our 48-day data set. We compare our results with the spin-down evolution reported during its previous active phase in the radio band. We also present total intensity pulse profiles at five different observing frequencies between 1.5 and 8.4 GHz, collected with the Lovell and the Effelsberg telescopes. The profile evolution in our data set is less erratic than what was reported during the previous active phase, and can be seen varying smoothly between observations. Profiles observed immediately after the outburst show the presence of at least five cycles of a very stable ~50-ms periodicity in the main pulse component that lasts for at least tens of days. This remarkable structure is seen across the full range of observing frequencies.

Ofek (2017) identified FIRST J141918.9+394036 (hereafter FIRST J1419+3940) as a radio source sharing similar properties and host galaxy type to the compact, persistent radio source associated with the first known repeating fast radio burst, FRB 121102. Law et al. (2018) showed that FIRST J1419+3940 is a transient source decaying in brightness over the last few decades. One possible interpretation is that FIRST J1419+3940 is a nearby analogue to FRB 121102 and that the radio emission represents a young magnetar nebula (as several scenarios assume for FRB 121102). Another interpretation is that FIRST J1419+3940 is the afterglow of an `orphan' long gamma-ray burst (GRB). The environment is similar to where most such events are produced. To distinguish between these hypotheses, we conducted radio observations using the European VLBI Network at 1.6 GHz to spatially resolve the emission and to search for millisecond-duration radio bursts. We detect FIRST J1419+3940 as a compact radio source with a flux density of $620 \pm 20\ \mathrm{\mu Jy}$ (on 2018 September 18) and a source size of $3.9 \pm 0.7\ \mathrm{mas}$ (i.e. $1.6 \pm 0.3\ \mathrm{pc}$ given the angular diameter distance of $83\ \mathrm{Mpc}$). These results confirm that the radio emission is non-thermal and imply an average expansion velocity of $(0.10 \pm 0.02)c$. Contemporaneous high-time-resolution observations using the 100-m Effelsberg telescope detected no millisecond-duration bursts of astrophysical origin. The source properties and lack of short-duration bursts are consistent with a GRB jet expansion, whereas they disfavor a magnetar birth nebula.

We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $\alpha = -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $\mu$s, with power law indices $\alpha_{\text{low}} = -3.48 \pm 0.04$ and $\alpha_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs.

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

We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.

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.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

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

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

The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nan\c{c}ay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data.

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

Polarised radio emission from PSR J1745-2900 has already been used to investigate the strength of the magnetic field in the Galactic Centre, close to Sagittarius A*. Here we report how persistent radio emission from this magnetar, for over four years since its discovery, has revealed large changes in the observed Faraday rotation measure, by up to 3500 rad m$^{-2}$ (a five per cent fractional change). From simultaneous analysis of the dispersion measure, we determine that these fluctuations are dominated by variations in the projected magnetic field, rather than the integrated free electron density, along the changing line of sight to the rapidly moving magnetar. From a structure function analysis of rotation measure variations, and a recent epoch of rapid change of rotation measure, we determine a minimum scale of magnetic fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR J1745-2900 is just ~ 0.1 pc behind an additional scattering screen.

We report the discovery and the results of follow-up timing observations of PSR J2045+3633 and PSR J2053+4650, two binary pulsars found in the Northern High Time Resolution Universe pulsar survey being carried out with the Effelsberg radio telescope. Having spin periods of 31.7 ms and 12.6 ms respectively, and both with massive white dwarf companions, $M_{c}\, > \, 0.8\, M_{\odot}$, the pulsars can be classified as mildly recycled. PSR J2045+3633 is remarkable due to its orbital period (32.3 days) and eccentricity $e\, = \, 0.01721244(5)$ which is among the largest ever measured for this class. After almost two years of timing the large eccentricity has allowed the measurement of the rate of advance of periastron at the 5-$\sigma$ level, 0.0010(2)$^\circ~\rm yr^{-1}$. Combining this with a detection of the orthometric amplitude of the Shapiro delay, we obtained the following constraints on the component masses (within general relativity): $M_{p}\, = \, 1.33^{+0.30}_{-0.28}\, M_{\odot}$, and $M_{c}\, = \, 0.94^{+0.14}_{-0.13}\, M_{\odot}$. PSR J2053+4650 has a 2.45-day circular orbit inclined to the plane of the sky at an angle $i\, = \, 85.0^{+0.8}_{-0.9}\,{^\circ}$. In this nearly edge-on case the masses can be obtained from the Shapiro delay alone. Our timing observations resulted in a significant detection of this effect giving: $M_{p}\, = \, 1.40^{+0.21}_{-0.18}\, M_{\odot}$, and $M_{c}\, = \, 0.86^{+0.07}_{-0.06}\, M_{\odot}$.

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.

The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.

Recent studies have shown possible connections between highly magnetized neutron stars ("magnetars"), whose X-ray emission is too bright to be powered by rotational energy, and ordinary radio pulsars. In addition to the magnetar SGR J1745-2900, one of the radio pulsars in the Galactic centre (GC) region, PSR J1746-2850, had timing properties implying a large magnetic field strength and young age, as well as a flat spectrum. All characteristics are similar to those of rare, transient, radio-loud magnetars. Using several deep non-detections from the literature and two new detections, we show that this pulsar is also transient in the radio. Both the flat spectrum and large amplitude variability are inconsistent with the light curves and spectral indices of 3 radio pulsars with high magnetic field strengths. We further use frequent, deep archival imaging observations of the GC in the past 15 years to rule out a possible X-ray outburst with a luminosity exceeding the rotational spin down rate. This source, either a transient magnetar without any detected X-ray counterpart or a young, strongly magnetized radio pulsar producing magnetar-like radio emission, further blurs the line between the two categories. We discuss the implications of this object for the radio emission mechanism in magnetars and for star and compact object formation in the GC.

In Torne et al. (2015), we showed detections of SGR J1745-2900 up to 225 GHz (1.33 mm); at that time the highest radio frequency detection of pulsar emission. In this work, we present the results of new observations of the same magnetar with detections up to 291 GHz (1.03 mm), together with evidence of linear polarization in its millimetre emission. SGR J1745-2900 continues to show variability and is, on average, a factor $\sim$4 brighter in the millimetre band than in our observations of July 2014. The new measured spectrum is slightly inverted, with $\left<\alpha\right> = +0.4\pm0.2$ (for $S_{\nu} \propto \nu^{\alpha})$. However, the spectrum does not seem to be well described by a single power law, which might be due to the intrinsic variability of the source, or perhaps a turn-up somewhere between 8.35 and 87 GHz. These results may help us to improve our still incomplete model of pulsar emission and, in addition, they further support the search for and study of pulsars located at the Galactic Centre using millimetre wavelengths.

We have resolved the scatter-broadened image of PSR B0329+54 and detected substructure within it. These results are not influenced by any extended structure of a source but instead are directly attributed to the interstellar medium. We obtained these results at 324 MHz with the ground-space interferometer RadioAstron which included the space radio telescope (SRT), ground-based Westerbork Synthesis Radio Telescope and 64-m Kalyazin Radio Telescope on baseline projections up to 330,000 km in 2013 November 22 and 2014 January 1 to 2. At short 15,000 to 35,000 km ground-space baseline projections the visibility amplitude decreases with baseline length providing a direct measurement of the size of the scattering disk of 4.8$\pm$0.8 mas. At longer baselines no visibility detections from the scattering disk would be expected. However, significant detections were obtained with visibility amplitudes of 3 to 5% of the maximum scattered around a mean and approximately constant up to 330,000 km. These visibilities reflect substructure from scattering in the interstellar medium and offer a new probe of ionized interstellar material. The size of the diffraction spot near Earth is 17,000$\pm$3,000 km. With the assumption of turbulent irregularities in the plasma of the interstellar medium, we estimate that the effective scattering screen is located 0.6$\pm$0.1 of the distance from Earth toward the pulsar.

Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of 7 and 3 pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.

Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.

A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.

Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.

We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) \times 10^{-19}\,\text{Hz} \ \text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${\Delta \nu/\nu=2.5(1) \times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $\sim 50$% that another glitch will be observed in a timing array pulsar within 10 years.

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely-separated antennas. This instrument included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14x25-m Westerbork Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function show two exponential scales, with characteristic delays of $\tau_1=4.1\pm 0.3\ \mu{\rm s}$ and $\tau_2=23\pm 3\ \mu{\rm s}$, indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles.

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.

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

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

The PSRIX backend is the primary pulsar timing instrument of the Effelsberg 100-m radio telescope since early 2011. This new ROACH-based system enables bandwidths up to 500 MHz to be recorded, significantly more than what was possible with its predecessor, the Effelsberg-Berkeley Pulsar Processor (EBPP). We review the first four years of PSRIX timing data for 33 pulsars collected as part of the monthly European Pulsar Timing Array (EPTA) observations. We describe the automated data analysis pipeline, CoastGuard, that we developed to reduce these observations. We also introduce TOASTER, the EPTA timing database used to store timing results, processing information and observation metadata. Using these new tools, we measure the phase-averaged flux densities at 1.4 GHz of all 33 pulsars. For 7 of these pulsars, our flux density measurements are the first values ever reported. For the other 26 pulsars, we compare our flux density measurements with previously published values. By comparing PSRIX data with EBPP data, we find an improvement of ~2-5 times in signal-to-noise ratio achievable, which translates to an increase of ~2-5 times in pulse time-of-arrival (TOA) precision. We show that such an improvement in TOA precision will improve the sensitivity to the stochastic gravitational wave background. Finally, we showcase the flexibility of the new PSRIX backend by observing several millisecond-period pulsars (MSPs) at 5 and 9 GHz. Motivated by our detections, we discuss the potential for complementing existing pulsar timing array data sets with MSP monitoring campaigns at these frequencies.

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.

The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are: the Effelsberg telescope, the Lovell telescope at Jodrell Bank, the Nan\c cay radio telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time-series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centered at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data, and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time-series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $\Omega^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

We report on simultaneous observations of the magnetar SGR J1745-2900 at frequencies $\nu = 2.54$ to $225\,\rm{GHz}$ using the Nancay 94-m equivalent, Effelsberg 100-m, and IRAM 30-m radio telescopes. We detect SGR J1745-2900 up to 225 GHz, the highest radio frequency detection of pulsed emission from a neutron star to date. Strong single pulses are also observed from 4.85 up to 154 GHz. At the millimetre band we see significant flux density and spectral index variabilities on time scales of tens of minutes, plus variability between days at all frequencies. Additionally, SGR J1745-2900 was observed at a different epoch at frequencies 296 to 472 GHz using the APEX 12-m radio telescope, with no detections. Over the period MJD 56859.83-56862.93 the fitted spectrum yields a spectral index of $\left<\alpha\right> = -0.4 \pm 0.1$ for a reference flux density $\left< S_{154} \right> = 1.1 \pm 0.2\rm{\,mJy}$ (with $S_{\nu} \propto {\nu}^{\alpha})$, a flat spectrum alike those of the other radio-loud magnetars. These results show that strongly magnetized neutron stars can be effective radio emitters at frequencies notably higher to what was previously known and that pulsar searches in the Galactic Centre are possible in the millimetre band.

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.

Millisecond pulsars (MSPs) are known as highly stable celestial clocks. Nevertheless, recent studies have revealed the unstable nature of their integrated pulse profiles, which may limit the achievable pulsar timing precision. In this paper, we present a case study on the pulse profile variability of PSR J1022+1001. We have detected approximately 14,000 sub-pulses (components of single pulses) in 35-hr long observations, mostly located at the trailing component of the integrated profile. Their flux densities and fractional polarisation suggest that they represent the bright end of the energy distribution in ordinary emission mode and are not giant pulses. The occurrence of sub-pulses from the leading and trailing components of the integrated profile is shown to be correlated. For sub-pulses from the latter, a preferred pulse width of approximately 0.25 ms has been found. Using simultaneous observations from the Effelsberg 100-m telescope and the Westerbork Synthesis Radio Telescope, we have found that the integrated profile varies on a timescale of a few tens of minutes. We show that improper polarisation calibration and diffractive scintillation cannot be the sole reason for the observed instability. In addition, we demonstrate that timing residuals generated from averages of the detected sub-pulses are dominated by phase jitter, and place an upper limit of ~700 ns for jitter noise based on continuous 1-min integrations.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

Simultaneous observations at four different frequencies viz. 313, 607, 1380 and 4850 MHz, for three pulsars, PSRs B0031-07, B0809+74 and B2319+60, are reported in this paper. Identified null and burst pulses are highly concurrent across more than decade of frequency. Small fraction of non-concurrent pulses (less than or equal to 3%) are observed, most of which occur at the transition instances. We report, with very high significance for the first time, full broadband nature of the nulling phenomenon in these three pulsars. These results suggest that nulling invokes changes on the global magnetospheric scale.

In recent years, instrumentation enabling pulsar observations with unprecedentedly high fractional bandwidth has been under development which can be used to substantially improve the precision of pulsar timing experiments. The traditional template-matching method used to calculate pulse times-of-arrival (ToAs), may not function effectively on these broadband data due to a variety of effects such as diffractive scintillation in the interstellar medium, profile variation as a function of frequency, dispersion measure (DM) evolution and so forth. In this paper, we describe the channelised Discrete Fourier Transform method that can greatly mitigate the influence of the aforementioned effects when measuring ToAs from broadband timing data. The method is tested on simulated data, and its potential in improving timing precision is shown. We further apply the method to PSR J1909$-$3744 data collected at the Nan\c{c}ay Radio Telescope with the Nan\c{c}ay Ultimate Pulsar Processing Instrument. We demonstrate a removal of systematics due to the scintillation effect as well as improvement on ToA measurement uncertainties. Our method also determines temporal variations in dispersion measure, which are consistent with multi-channel timing approaches used earlier.

In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

We report the discovery of four gamma-ray pulsars, detected in computing-intensive blind searches of data from the Fermi Large Area Telescope (LAT). The pulsars were found using a novel search approach, combining volunteer distributed computing via Einstein@Home and methods originally developed in gravitational-wave astronomy. The pulsars PSRs J0554+3107, J1422-6138, J1522-5735, and J1932+1916 are young and energetic, with characteristic ages between 35 and 56 kyr and spin-down powers in the range $6\times10^{34}$ - $10^{36}$ erg s$^{-1}$. They are located in the Galactic plane and have rotation rates of less than 10 Hz, among which the 2.1 Hz spin frequency of PSR J0554+3107 is the slowest of any known gamma-ray pulsar. For two of the new pulsars, we find supernova remnants coincident on the sky and discuss the plausibility of such associations. Deep radio follow-up observations found no pulsations, suggesting that all four pulsars are radio-quiet as viewed from Earth. These discoveries, the first gamma-ray pulsars found by volunteer computing, motivate continued blind pulsar searches of the many other unidentified LAT gamma-ray sources.

We present temporal scattering measurements of single pulses and average profiles of PSR J1745--2900, a magnetar recently discovered only 3 arcsec away from Sagittarius A* (Sgr A*), from 1.2 - 18.95 GHz using the Effelsberg 100-m Radio Telescope, the Nan\c{c}ay Decimetric Radio Telescope, and the Jodrell Bank Lovell Telescope. Single pulse analysis shows that the integrated pulse profile above 2 GHz is dominated by pulse jitter, while below 2 GHz the pulse profile shape is dominated by scattering. The high dispersion measure and rotation measure of the magnetar suggest that it is close to Sgr A* (within ~0.1 pc). This is the first object in the GC with both pulse broadening and angular broadening measurements. We measure a pulse broadening spectral index of alpha = -3.8 +/- 0.2 and a pulse broadening time scale at 1 GHz of tau_GHz = 1.3 +/- 0.2 s, which is several orders of magnitude lower than the scattering predicted by the NE2001 model (Cordes and Lazio 2002). If this scattering timescale is representative of the GC as a whole, then previous surveys should have detected many pulsars. The lack of detections implies either our understanding of scattering in the GC is incomplete or there are fewer pulsars in the GC than previously predicted. Given that magnetars are a rare class of radio pulsar, we believe that there are many canonical and millisecond pulsars in the GC, and not surprisingly, scattering regions in the GC have complex spatial structures.

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

We report on the setup and initial discoveries of the Northern High Time Resolution Universe survey for pulsars and fast transients, the first major pulsar survey conducted with the 100-m Effelsberg radio telescope and the first in 20 years to observe the whole northern sky at high radio frequencies. Using a newly developed 7-beam receiver system combined with a state-of-the-art polyphase filterbank, we record an effective bandwidth of 240 MHz in 410 channels centred on 1.36 GHz with a time resolution of 54 $\mu$s. Such fine time and frequency resolution increases our sensitivity to millisecond pulsars and fast transients, especially deep inside the Galaxy, where previous surveys have been limited due to intra-channel dispersive smearing. To optimise observing time, the survey is split into three integration regimes dependent on Galactic latitude, with 1500-s, 180-s and 90-s integrations for latitude ranges $|b|<3.5^{\circ}$, $|b|<15^{\circ}$ and $|b|>15^{\circ}$, respectively. The survey has so far resulted in the discovery of 15 radio pulsars, including a pulsar with a characteristic age of $\sim18$ kyr, {PSR J2004+3429}, and a highly eccentric, binary millisecond pulsar, {PSR J1946+3417}. All newly discovered pulsars are timed using the 76-m Lovell radio telescope at the Jodrell Bank Observatory and the Effelsberg radio telescope. We present timing solutions for all newly discovered pulsars and discuss potential supernova remnant associations for {PSR J2004+3429}.

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

Radio pulsars in relativistic binary systems are unique tools to study the curved space-time around massive compact objects. The discovery of a pulsar closely orbiting the super-massive black hole at the centre of our Galaxy, Sgr A*, would provide a superb test-bed for gravitational physics. To date, the absence of any radio pulsar discoveries within a few arc minutes of Sgr A* has been explained by one principal factor: extreme scattering of radio waves caused by inhomogeneities in the ionized component of the interstellar medium in the central 100 pc around Sgr A*. Scattering, which causes temporal broadening of pulses, can only be mitigated by observing at higher frequencies. Here we describe recent searches of the Galactic centre region performed at a frequency of 18.95 GHz with the Effelsberg radio telescope.

In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of observing time constant. The optimized schedule depends on the timing noise characteristics of each individual pulsar as well as the performance of instrumentation. Several examples are given to illustrate the effects of different types of noise. A method to select the most suitable pulsars to be included in a pulsar timing array project is also presented.

Abbreviated: We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax. We find that, at low redshift, a PTA is able to detect the nano-Hertz GW from super massive black hole binary systems with masses of $\sim10^8 - 10^{10}\,M_{\sun}$ less than $\sim10^5$\,years before the final merger, and those with less than $\sim10^3 - 10^4$ years before merger may allow us to detect the evolution of binaries. We derive an analytical expression to describe the accuracy of a pulsar distance measurement via timing parallax. We consider five years of bi-weekly observations at a precision of 15\,ns for close-by ($\sim 0.5 - 1$\,kpc) pulsars. Timing twenty pulsars would allow us to detect a GW source with an amplitude larger than $5\times 10^{-17}$. We calculate the corresponding GW and binary orbital parameters and their measurement precision. The accuracy of measuring the binary orbital inclination angle, the sky position, and the GW frequency are calculated as functions of the GW amplitude. We note that the "pulsar term", which is commonly regarded as noise, is essential for obtaining an accurate measurement for the GW source location. We also show that utilizing the information encoded in the GW signal passing the Earth also increases the accuracy of pulsar distance measurements. If the gravitational wave is strong enough, one can achieve sub-parsec distance measurements for nearby pulsars with distance less than $\sim 0.5 - 1$\,kpc.

The low frequency (110--180 MHz) capabilities of the Westerbork Synthesis Radio Telescope (WSRT) are used to characterise a large collection of single pulses from three low magnetic field pulsars. Using the Pulsar Machine II (PuMa-II) to acquire and coherently dedisperse the pulsar signals, we examine whether the bright pulses observed in these pulsars are related to the classical giant pulse emission. Giant pulses are reported from PSR B1112$+$50 and bright pulses from the PSRs B1133$+$16 and B0031$-$07. These pulsars also exhibit large intensity modulations observed as rapid changes in the single pulse intensity. Evidence of global magnetospheric effects is provided by our detection of bright double pulses in PSRs B0031$-$07 and B1133$+$16. Using the multi-frequency observations, we accurately determine the dispersion measures (4.844$\pm$0.002 for B1133$+$16 and 9.1750$\pm$0.0001 for B1112$+$50), derive the radio emission height in PSR B1133$+$16 and report on the properties of subpulse drift modes in these pulsars. We also find that these pulsars show a much larger intensity modulation at low sky frequencies resulting in narrow and bright emissions.

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

The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.

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

With LOFAR beginning operation in 2008 there is huge potential for studying pulsars with high signal to noise at low frequencies. We present results of observations made with the Westerbork Synthesis Radio Telescope to revisit, with modern technology, this frequency range. Coherently dedispersed profiles of millisecond pulsars obtained simultaneously between 115-175 MHz are presented. We consider the detections and non-detections of 14 MSPs in light of previous observations and the fluxes, dispersion measures and spectral indices of these pulsars. The excellent prospects for LOFAR finding new MSPs and studying the existing systems are then discussed in light of these results.