We report on the detection of seven bursts from the periodically active, repeating fast radio burst (FRB) source FRB 180916.J0158+65 in the 300-400-MHz frequency range with the Green Bank Telescope (GBT). Emission in multiple bursts is visible down to the bottom of the GBT band, suggesting that the cutoff frequency (if it exists) for FRB emission is lower than 300 MHz. Observations were conducted during predicted periods of activity of the source, and had simultaneous coverage with the Low Frequency Array (LOFAR) and the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. We find that one of the GBT-detected bursts has potentially associated emission in the CHIME band (400-800 MHz) but we detect no bursts in the LOFAR band (110-190 MHz), placing a limit of $\alpha > -1.0$ on the spectral index of broadband emission from the source. We also find that emission from the source is severely band-limited with burst bandwidths as low as $\sim$40 MHz. In addition, we place the strictest constraint on observable scattering of the source, $<$ 1.7 ms, at 350 MHz, suggesting that the circumburst environment does not have strong scattering properties. Additionally, knowing that the circumburst environment is optically thin to free-free absorption at 300 MHz, we find evidence against the association of a hyper-compact HII region or a young supernova remnant (age $<$ 50 yr) with the source.

The discovery that at least some Fast Radio Bursts (FRBs) repeat has ruled out cataclysmic events as the progenitors of these particular bursts. FRB~121102 is the most well-studied repeating FRB but despite extensive monitoring of the source, no underlying pattern in the repetition has previously been identified. Here, we present the results from a radio monitoring campaign of FRB~121102 using the 76-m Lovell telescope. Using the pulses detected in the Lovell data along with pulses from the literature, we report a detection of periodic behaviour of the source over the span of five years of data. We predict that the source is currently `off' and that it should turn `on' for the approximate MJD range $59002-59089$ (2020-06-02 to 2020-08-28). This result, along with the recent detection of periodicity from another repeating FRB, highlights the need for long-term monitoring of repeating FRBs at a high cadence. Using simulations, we show that one needs at least 100 hours of telescope time to follow-up repeating FRBs at a cadence of 0.5--3 days to detect periodicities in the range of 10--150 days. If the period is real, it shows that repeating FRBs can have a large range in their activity periods that might be difficult to reconcile with neutron star precession models.

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 bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while many other low frequency surveys in the past have failed to find any FRBs. Here, we present results from a FRB survey recently conducted at the Jodrell Bank Observatory at 332 MHz with the 76-m Lovell telescope for a total of 58 days. We did not detect any FRBs in the survey and report a 90$\%$ upper limit of 5500 FRBs per day per sky for a Euclidean Universe above a fluence threshold of 46 Jy ms. We discuss the possibility of absorption as the main cause of non-detections in low frequency (< 800 MHz) searches and invoke different absorption models to explain the same. We find that Induced Compton Scattering alone cannot account for absorption of radio emission and that our simulations favour a combination of Induced Compton Scattering and Free-Free Absorption to explain the non-detections. For a free-free absorption scenario, our constraints on the electron density are consistent with those expected in the post-shock region of the ionized ejecta in Super-Luminous SuperNovae (SLSNe).

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

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.

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

We are using the LOw-Frequency ARray (LOFAR) to perform the LOFAR Tied-Array All-Sky (LOTAAS) survey for pulsars and fast transients. Here we present the astrometric and rotational parameters of 20 pulsars discovered as part of LOTAAS. These pulsars have regularly been observed with LOFAR at 149 MHz and the Lovell telescope at 1532 MHz, supplemented by some observations with the Lovell telescope at 334 MHz and the Nancay Radio Telescope at 1484 MHz. Timing models are calculated for the 20 pulsars, some of which are among the slowest-spinning pulsars known. PSR J1236-0159 rotates with a period P ~ 3.6 s, while 5 additional pulsars show P > 2 s. Also, the spin-down rates Pdot are, on average, low, with PSR J0815+4611 showing Pdot ~ 4E-18. Some of the pulse profiles, generically single-peaked, present complex shapes evolving with frequency. Multi-frequency flux measurements show that these pulsars have generically relatively steep spectra but exceptions are present, with values ranging between ~ -4 and -1. Among the pulsar sample, a large fraction shows large single-pulse variability, with 4 pulsars being undetectable more than 15% of the time and one tentatively classified as a Rotating Radio Transient. Two single-peaked pulsars show drifting sub-pulses.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

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.

We investigate the impact of noise processes on high-precision pulsar timing. Our analysis focuses on the measurability of the second spin frequency derivative $\ddot{\nu}$. This $\ddot{\nu}$ can be induced by several factors including the radial velocity of a pulsar. We use Bayesian methods to model the pulsar times-of-arrival in the presence of red timing noise and dispersion measure variations, modelling the noise processes as power laws. Using simulated times-of-arrival that both include red noise, dispersion measure variations and non-zero $\ddot{\nu}$ values, we find that we are able to recover the injected $\ddot{\nu}$, even when the noise model used to inject and recover the input parameters are different. Using simulations, we show that the measurement uncertainty on $\ddot{\nu}$ decreases with the timing baseline $T$ as $T^\gamma$, where $\gamma=-7/2+\alpha/2$ for power law noise models with shallow power law indices $\alpha$ ($0<\alpha<4$). For steep power law indices ($\alpha>8$), the measurement uncertainty reduces with $T^{-1/2}$. We applied this method to times-of-arrival from the European Pulsar Timing Array and the Parkes Pulsar Timing Array and determined $\ddot{\nu}$ probability density functions for 49 millisecond pulsars. We find a statistically significant $\ddot{\nu}$ value for PSR\,B1937+21 and consider possible options for its origin. Significant (95 per cent C.L.) values for $\ddot{\nu}$ are also measured for PSRs\,J0621+1002 and J1022+1001, thus future studies should consider including it in their ephemerides. For binary pulsars with small orbital eccentricities, like PSR\,J1909$-$3744, extended ELL1 models should be used to overcome computational issues. The impacts of our results on the detection of gravitational waves are also discussed.

We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array.

We report on Westerbork Synthesis Radio Telescope observations of the Crab pulsar at 350 MHz from 2012 November 24 until 2015 June 21. During this period we consistently observe variations in the pulse profile of the Crab. Both variations in the scattering width of the pulse profile as well as delayed copies, also known as echoes, are seen regularly. These observations support the classification of two types of echoes: those that follow the truncated exponential shape expected for the thin-screen scattering approximation, and echoes that show a smoother, more Gaussian shape. During a sequence of high-cadence observations in 2015, we find that these non-exponential echoes evolve in time by approaching the main pulse and interpulse in phase, overlapping the main pulse and interpulse, and later receding. We find a pulse scatter-broadening time scale, $\tau$, scaling with frequency as $\nu^{\alpha}$, with $\alpha=-3.9\pm0.5$, which is consistent with expected values for thin-screen scattering models

FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination of a dispersion measure (DM); we argue that it is appropriate to use a DM metric that maximizes frequency-averaged pulse structure, as opposed to peak signal-to-noise, and find DM = 560.57 +/- 0.07 pc/cc at MJD 57644. After correcting for dispersive delay, we find that the sub-bursts have characteristic frequencies that typically drift lower at later times in the total burst envelope. In the 1.1-1.7 GHz band, the ~ 0.5-1-ms sub-bursts have typical bandwidths ranging from 100-400 MHz, and a characteristic drift rate of ~ 200 MHz/ms towards lower frequencies. At higher radio frequencies, the sub-burst bandwidths and drift rate are larger, on average. While these features could be intrinsic to the burst emission mechanism, they could also be imparted by propagation effects in the medium local to the source. Comparison of the burst DMs with previous values in the literature suggests an increase of Delta(DM) ~ 1-3 pc/cc in 4 years, though this could be a stochastic variation as opposed to a secular trend. This implies changes in the local medium or an additional source of frequency-dependent delay. Overall, the results are consistent with previously proposed scenarios in which FRB 121102 is embedded in a dense 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 present the discovery of PSR J0250+5854, a radio pulsar with a spin period of 23.5 s. This is the slowest-spinning radio pulsar known. PSR J0250+5854 was discovered by the LOFAR Tied-Array All-Sky Survey (LOTAAS), an all-Northern-sky survey for pulsars and fast transients at a central observing frequency of 135 MHz. We subsequently detected pulsations from the pulsar in the interferometric images of the LOFAR Two-metre Sky Survey, allowing for sub-arcsecond localization. This, along with a pre-discovery detection 2 years prior, allowed us to measure the spin-period derivative to be $\dot{P}=2.7 \times 10^{-14}$ s s$^{-1}$. The observed spin period derivative of PSR J0250+5854 indicates a surface magnetic field strength, characteristic age and spin-down luminosity of $2.6 \times 10^{13}$G, $13.7$ Myr and $8.2 \times 10^{28}$ erg s$^{-1}$ respectively, for a dipolar magnetic field configuration. This also places the pulsar beyond the conventional pulsar death line, where radio emission is expected to cease. The spin period of PSR J0250+5854 is similar to those of the high-energy-emitting magnetars and X-ray dim isolated neutron stars (XDINSs). However, the pulsar was not detected by the Swift/XRT in the energy band of 0.3-10 keV, placing a bolometric luminosity limit of $1.5 \times 10^{32}$ erg s$^{-1}$ for an assumed $N_{\rm H}=1.35\times10^{21}$ cm$^{-2}$ and a temperature of 85 eV (typical of XDINSs). We discuss the implications of the discovery for models of the pulsar death line as well as the prospect of finding more similarly long-period pulsars, including the advantages provided by LOTAAS for this.

Searches for millisecond-duration, dispersed single pulses have become a standard tool used during radio pulsar surveys in the last decade. They have enabled the discovery of two new classes of sources: rotating radio transients and fast radio bursts. However, we are now in a regime where the sensitivity to single pulses in radio surveys is often limited more by the strong background of radio frequency interference (RFI, which can greatly increase the false-positive rate) than by the sensitivity of the telescope itself. To mitigate this problem, we introduce the Single-pulse Searcher (SpS). This is a new machine-learning classifier designed to identify astrophysical signals in a strong RFI environment, and optimized to process the large data volumes produced by the new generation of aperture array telescopes. It has been specifically developed for the LOFAR Tied-Array All-Sky Survey (LOTAAS), an ongoing survey for pulsars and fast radio transients in the northern hemisphere. During its development, SpS discovered 7 new pulsars and blindly identified ~80 known sources. The modular design of the software offers the possibility to easily adapt it to other studies with different instruments and characteristics. Indeed, SpS has already been used in other projects, e.g. to identify pulses from the fast radio burst source FRB 121102. The software development is complete and SpS is now being used to re-process all LOTAAS data collected to date.

We report on a time-domain search for pulsars in 44 steep spectrum radio sources originally identified from recent imaging surveys. The time-domain search was conducted at 327 MHz using the Ooty radio telescope, and utilized a semi-coherent dedispersion scheme retaining the sensitivity even for sub-millisecond periods up to reasonably high dispersion measures. No new pulsars were found. We discuss the nature of these steep spectrum sources and argue that majority of the sources in our sample should either be pulsars or a new category of Galactic sources. Several possibilities that could hinder detection of these sources as pulsars, including anomalously high scattering or alignment of the rotation and magnetic axes, are discussed in detail, and we suggest unconventional search methods to further probe these possibilities.

We have measured variations in scattering time scales in the Crab Pulsar over a 30-year period, using observations made at 610 MHz with the 42-ft telescope at Jodrell Bank Observatory. Over more recent years, where regular Lovell Telescope observations at frequencies around 1400 MHz were available, we have also determined the dispersion measure variations, after disentangling the scattering delay from the dispersive delay. We demonstrate a relationship between scattering and dispersion measure variations, with a correlation coefficient of $0.56\pm0.01$. The short time scales over which these quantities vary, the size of the variations, and the close correlation between scattering and dispersion measure all suggest that the effects are due to discrete structures within the Crab Nebula, with size scales of $\sim6$ AU (corresponding to an angular size of $\sim2$ mas at an assumed distance of 2200 pc). We mitigate the effects of scattering on the observed pulse shape by using the measured scattering information to modify the template used for generating the pulse arrival times, thus improving the precision to which the pulsar can be timed. We test this on timing data taken during periods of high scattering, and obtain a factor of two improvement in the root mean square of the timing residuals.

We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $\Delta \nu = 1.530 \times 10^{-5}$ Hz, corresponding to a fractional increase of $\Delta \nu / \nu = 0.516 \times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time observations of the Crab pulsar we are able to partially resolve a fraction of the total spin-up of the star. This delayed spin-up occurred over a timescale of $\sim$1.7 days and is similar to the behaviour seen in the 1989 and 1996 large Crab pulsar glitches. The spin-down rate also increased at the glitch epoch by $\Delta \dot{\nu} / \dot{\nu} = 7 \times 10^{-3}$. In addition to being the largest such event observed in the Crab, the glitch occurred after the longest period of glitch inactivity since at least 1984 and we discuss a possible relationship between glitch size and waiting time. No changes to the shape of the pulse profile were observed near the glitch epoch at 610 MHz or 1520 MHz, nor did we identify any changes in the X-ray flux from the pulsar. The long-term recovery from the glitch continues to progress as $\dot{\nu}$ slowly rises towards pre-glitch values. In line with other large Crab glitches, we expect there to be a persistent change to $\dot{\nu}$. We continue to monitor the long-term recovery with frequent, high quality observations.

We investigate the impact of intrinsic, kinematic and gravitational effects on high precision pulsar timing. We present an analytical derivation and a numerical computation of the impact of these effects on the first and second derivative of the pulsar spin frequency. In addition, in the presence of white noise, we derive an expression to determine the expected measurement uncertainty of a second derivative of the spin frequency for a given timing precision, observing cadence and timing baseline and find that it strongly depends on the latter ($\propto t^{-7/2}$). We show that for pulsars with significant proper motion, the spin frequency second derivative is dominated by a term dependent on the radial velocity of the pulsar. Considering the data sets from three Pulsar Timing Arrays, we find that for PSR J0437$-$4715 a detectable spin frequency second derivative will be present if the absolute value of the radial velocity exceeds 33 km s$^{-1}$. Similarly, at the current timing precision and cadence, continued timing observations of PSR J1909$-$3744 for about another eleven years, will allow the measurement of its frequency second derivative and determine the radial velocity with an accuracy better than 14 km s$^{-1}$. With the ever increasing timing precision and observing baselines, the impact of the, largely unknown, radial velocities of pulsars on high precision pulsar timing can not be neglected.

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.

We present the longest-term timing study so far of three Rotating Radio Transients (RRATs) - J1819-1458, J1840-1419 and J1913+1330 - performed using the Lovell, Parkes and Green Bank telescopes over the past decade. We study long-term and short- term variations of the pulse emission rate from these RRATs and report a marginal indication of a long-term increase in pulse detection rate over time for PSR J1819-1458 and J1913+1330. For PSR J1913+1330, we also observe a two orders of magnitude variation in the observed pulse detection rates across individual epochs, which may constrain the models explaining the origin of RRAT pulses. PSR J1913+1330 is also observed to exhibit a weak persistent emission mode. We investigate the post-glitch timing properties of J1819-1458 (the only RRAT for which glitches are observed) and discuss the implications for possible glitch models. Its post-glitch over-recovery of the frequency derivative is magnetar-like and similar behaviour is only observed for two other pulsars, both of which have relatively high magnetic field strengths. Following the over-recovery we also observe that some fraction of the pre-glitch frequency derivative is gradually recovered.

We present SMA and NOEMA observations of the host galaxy of FRB 121102 in the CO 3-2 and 1-0 transitions, respectively. We do not detect emission from either transition. We set $3\sigma$ upper limits to the CO luminosity $L_{CO} < 2.5 \times 10^7\,{\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 3-2 and $L_{CO} < 2.3 \times 10^9\, {\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 1-0. For Milky-Way-like star formation properties, we set a $3\sigma$ upper limit on the $H_2$ mass of $2.5 \times 10^8 \rm\ M_{\odot}$, slightly less than the predictions for the $H_2$ mass based on the star formation rate. The true constraint on the $H_2$ mass may be significantly higher, however, because of the reduction in CO luminosity that is common forlow-metallicity dwarf galaxies like the FRB host galaxy. These results demonstrate the challenge of identifying the nature of FRB progenitors through study of the host galaxy molecular gas. We also place a limit of 42 $\mu$Jy ($3\sigma$) on the continuum flux density of the persistent radio source at 97 GHz, consistent with a power-law extrapolation of the low frequency spectrum, which may arise from an AGN or other nonthermal source.

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 have observed and analysed the eclipses of the black widow pulsar J1810+1744 at low radio frequencies. Using LOw-Frequency ARray (LOFAR) and Westerbork Synthesis Radio Telescope observations between 2011--2015 we have measured variations in flux density, dispersion measure and scattering around eclipses. High-time-resolution, simultaneous beamformed and interferometric imaging LOFAR observations show concurrent disappearance of pulsations and total flux from the source during the eclipses, with a $3\sigma$ upper limit of 36 mJy ($<10\%$ of the pulsar's averaged out-of-eclipse flux density). The dispersion measure variations are highly asymmetric, suggesting a tail of material swept back due to orbital motion. The egress deviations are variable on timescales shorter than the 3.6 hr orbital period and are indicative of a clumpy medium. Additional pulse broadening detected during egress is typically $<20\%$ of the pulsar's spin period, showing no evidence of scattering the pulses beyond detectability in the beamformed data. The eclipses, lasting $\sim13\%$ of the orbit at 149 MHz, are shown to be frequency-dependent with total duration scaling as $\propto\nu^{-0.41\pm0.03}$. The results are discussed in the context of the physical parameters of the system, and an examination of eclipse mechanisms reveals cyclotron-synchrotron absorption as the most likely primary cause, although non-linear scattering mechanisms cannot be quantitatively ruled out. The inferred mass loss rate is a similar order-of-magnitude to the mean rate required to fully evaporate the companion in a Hubble time.

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.

Fast radio bursts (FRBs) are millisecond-duration, extragalactic radio flashes of unknown physical origin. FRB 121102, the only known repeating FRB source, has been localized to a star-forming region in a dwarf galaxy at redshift z = 0.193, and is spatially coincident with a compact, persistent radio source. The origin of the bursts, the nature of the persistent source, and the properties of the local environment are still debated. Here we present bursts that show ~100% linearly polarized emission at a very high and variable Faraday rotation measure in the source frame: RM_src = +1.46 x 10^5 rad m^-2 and +1.33 x 10^5 rad m^-2 at epochs separated by 7 months, in addition to narrow (< 30 mus) temporal structure. The large and variable rotation measure demonstrates that FRB 121102 is in an extreme and dynamic magneto-ionic environment, while the short burst durations argue for a neutron star origin. Such large rotation measures have, until now, only been observed in the vicinities of massive black holes (M_BH > 10^4 MSun). Indeed, the properties of the persistent radio source are compatible with those of a low-luminosity, accreting massive black hole. The bursts may thus come from a neutron star in such an environment. However, the observed properties may also be explainable in other models, such as a highly magnetized wind nebula or supernova remnant surrounding a young neutron star.

We have used LOFAR to perform targeted millisecond pulsar surveys of Fermi gamma-ray sources. Operating at a center frequency of 135 MHz, the surveys use a novel semi-coherent dedispersion approach where coherently dedispersed trials at coarsely separated dispersion measures are incoherently dedispersed at finer steps. Three millisecond pulsars have been discovered as part of these surveys. We describe the LOFAR surveys and the properties of the newly discovered pulsars.

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.

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

One-opposition near-Earth asteroids (NEAs) are growing in number, and they must be recovered to prevent loss and mismatch risk, and to improve their orbits, as they are likely to be too faint for detection in shallow surveys at future apparitions. We aimed to recover more than half of the one-opposition NEAs recommended for observations by the Minor Planet Center (MPC) using the Isaac Newton Telescope (INT) in soft-override mode and some fractions of available D-nights. During about 130 hours in total between 2013 and 2016, we targeted 368 NEAs, among which 56 potentially hazardous asteroids (PHAs), observing 437 INT Wide Field Camera (WFC) fields and recovering 280 NEAs (76% of all targets). Engaging a core team of about ten students and amateurs, we used the THELI, Astrometrica, and the Find_Orb software to identify all moving objects using the blink and track-and-stack method for the faintest targets and plotting the positional uncertainty ellipse from NEODyS. Most targets and recovered objects had apparent magnitudes centered around V~22.8 mag, with some becoming as faint as V~24 mag. One hundred and three objects (representing 28% of all targets) were recovered by EURONEAR alone by Aug 2017. Orbital arcs were prolonged typically from a few weeks to a few years; our oldest recoveries reach 16 years. The O-C residuals for our 1,854 NEA astrometric positions show that most measurements cluster closely around the origin. In addition to the recovered NEAs, 22,000 positions of about 3,500 known minor planets and another 10,000 observations of about 1,500 unknown objects (mostly main-belt objects) were promptly reported to the MPC by our team. Four new NEAs were discovered serendipitously in the analyzed fields, increasing the counting to nine NEAs discovered by the EURONEAR in 2014 and 2015.

We present coordinated Chandra X-ray Observatory and Karl G. Jansky Very Large Array observations of the transitional millisecond pulsar PSR J1023+0038 in its low-luminosity accreting state. The unprecedented five hours of strictly simultaneous X-ray and radio continuum coverage for the first time unambiguously show a highly reproducible, anti-correlated variability pattern. The characteristic switches from the X-ray high mode into a low mode are always accompanied by a radio brightening with duration that closely matches the X-ray low mode interval. This behavior cannot be explained by a canonical inflow/outflow accretion model where the radiated emission and the jet luminosity are powered by, and positively correlated with, the available accretion energy. We interpret this phenomenology as alternating episodes of low-level accretion onto the neutron star during the X-ray high mode that are interrupted by rapid ejections of plasma by the active rotation-powered pulsar, possibly initiated by a reconfiguration of the pulsar magnetosphere, that cause a transition to a less luminous X-ray mode. The observed anti-correlation between radio and X-ray luminosity has an additional consequence: transitional MSPs can make excursions into a region of the radio/X-ray luminosity plane previously thought to be occupied solely by black hole X-ray binary sources. This complicates the use of this luminosity relation to identify candidate black holes, suggesting the need for additional discriminants when attempting to establish the true nature of the accretor.

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5$\sigma$ upper limit of $3\times10^{-11}$ erg cm$^{-2}$ on the 0.5--10 keV fluence for X-ray bursts at the time of radio bursts for durations $<700$ ms, which corresponds to a burst energy of $4\times10^{45}$ erg at the measured distance of FRB 121102. We also place limits on the 0.5--10 keV fluence of $5\times10^{-10}$ erg cm$^{-2}$ and $1\times10^{-9}$ erg cm$^{-2}$ for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5$\sigma$ upper limit on the 10--100 keV fluence of $4\times10^{-9}$ erg cm$^{-2}$ ($5\times10^{47}$ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5$\sigma$ upper limit on the 0.5--10 keV flux of $4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($3\times10^{41}$ erg~s$^{-1}$ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general.

Using LOFAR, we have performed a very-low-frequency (115-155 MHz) radio survey for millisecond pulsars (MSPs). The survey targeted 52 unidentified Fermi $\gamma$-ray sources. Employing a combination of coherent and incoherent dedispersion, we have mitigated the dispersive effects of the interstellar medium while maintaining sensitivity to fast-spinning pulsars. Toward 3FGL J1553.1+5437 we have found PSR J1552+5437, the first MSP to be discovered (through its pulsations) at a radio frequency $<$ 200 MHz. PSR J1552+5437 is an isolated MSP with a 2.43 ms spin period and a dispersion measure of 22.9 pc cm$^{-3}$. The pulsar has a very steep radio spectral index ($\alpha < -$2.8 $\pm$ 0.4). We obtain a phase-connected timing solution combining the 0.74 years of radio observations with $\gamma$-ray photon arrival times covering 7.5 years of Fermi observations. We find that the radio and $\gamma$-ray pulse profiles of PSR J1552+5437 appear to be nearly aligned. The very steep spectrum of PSR J1552+5437, along with other recent discoveries, hints at a population of radio MSPs that have been missed in surveys using higher observing frequencies. Detecting such steep spectrum sources is important for mapping the population of MSPs down to the shortest spin periods, understanding their emission in comparison to slow pulsars, and quantifying the prospects for future surveys with low-frequency radio telescopes like SKA-Low and its precursors.

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.

We report on the results of a 4-year timing campaign of PSR~J2222$-0137$, a 2.44-day binary pulsar with a massive white dwarf (WD) companion, with the Nan\c{c}ay, Effelsberg and Lovell radio telescopes. Using the Shapiro delay for this system, we find a pulsar mass $m_{p}=1.76,\pm\,0.06,M_\odot$ and a WD mass $m_{c}\,=\,1.293\,\pm\,0.025\, M_\odot$. We also measure the rate of advance of periastron for this system, which is marginally consistent with the GR prediction for these masses. The short lifetime of the massive WD progenitor star led to a rapid X-ray binary phase with little ($< \, 10^{-2} \, M_\odot$) mass accretion onto the neutron star (NS); hence, the current pulsar mass is, within uncertainties, its birth mass; the largest measured to date. We discuss the discrepancy with previous mass measurements for this system; we conclude that the measurements presented here are likely to be more accurate. Finally, we highlight the usefulness of this system for testing alternative theories of gravity by tightly constraining the presence of dipolar radiation. This is of particular importance for certain aspects of strong-field gravity, like spontaneous scalarization, since the mass of PSR~J2222$-0137$ puts that system into a poorly tested parameter range.

We present optical, near- and mid-infrared imaging of the host galaxy of FRB 121102 with the Gemini North telescope, the Hubble Space Telescope and the Spitzer Space Telescope. The FRB 121102 host galaxy is resolved, revealing a bright star forming region located in the outskirts of the irregular, low-metallicity dwarf galaxy. The star forming region has a half-light radius of 0.68 kpc (0.20 arcsec), encompassing the projected location of the compact (<0.7 pc), persistent radio source that is associated with FRB 121102. The half-light diameter of the dwarf galaxy is 5 to 7 kpc, and broadband spectral energy distribution fitting indicates that it has a total stellar mass of M*~10^8 Msun. The metallicity of the host galaxy is low, 12+log10 ([O/H])=8.0+-0.1. The properties of the host galaxy of FRB 121102 are comparable to those of extreme emission line galaxies, also known to host hydrogen-poor superluminous supernovae and long-duration gamma-ray bursts. The projected location of FRB 121102 within the star forming region supports the proposed connection of FRBs with newly born neutron stars or magnetars.

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate that the pulsar is in either a low-eccentricity ($0.33\lesssim e\lesssim0.4$) smaller orbit with a low mass companion (such as a main sequence star, white dwarf, neutron star, or stellar mass black hole) or a high-eccentricity ($e\gtrsim0.9$) larger orbit with a massive companion. The cluster mass properties and the observed properties of 4U 1820$-$30 and the other pulsars in the cluster argue against the low-eccentricity possibility. The high-eccentricity solution reveals that the pulsar is most likely orbiting around an intermediate-mass black hole (IMBH) of mass $> 7,500$~M$_\odot$ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central black hole (BH), predicts an acceleration that is commensurate with that measured for the pulsar. It further predicts that the model-dependent minimum mass of the IMBH is $\sim60,000$~M$_\odot$. Accounting for the associated contribution to the observed period derivative indicates that the $\gamma$-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central black holes and they may be revealed by the study of new pulsars found sufficiently close to their centres. Note that we found an erratum in Section 5 and thus, the $\sim$60,000~M$_\odot$ mass mentioned above has to be replaced by the correct model-dependent mass limit of $\sim$20,000~M$_\odot$. See the erratum appended.

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.

The precise localization of the repeating fast radio burst (FRB 121102) has provided the first unambiguous association (chance coincidence probability $p\lesssim3\times10^{-4}$) of an FRB with an optical and persistent radio counterpart. We report on optical imaging and spectroscopy of the counterpart and find that it is an extended ($0.6^{\prime\prime}-0.8^{\prime\prime}$) object displaying prominent Balmer and [OIII] emission lines. Based on the spectrum and emission line ratios, we classify the counterpart as a low-metallicity, star-forming, $m_{r^\prime} = 25.1$ AB mag dwarf galaxy at a redshift of $z=0.19273(8)$, corresponding to a luminosity distance of 972 Mpc. From the angular size, the redshift, and luminosity, we estimate the host galaxy to have a diameter $\lesssim4$ kpc and a stellar mass of $M_*\sim4-7\times 10^{7}\,M_\odot$, assuming a mass-to-light ratio between 2 to 3$\,M_\odot\,L_\odot^{-1}$. Based on the H$\alpha$ flux, we estimate the star formation rate of the host to be $0.4\,M_\odot\,\mathrm{yr^{-1}}$ and a substantial host dispersion measure depth $\lesssim 324\,\mathrm{pc\,cm^{-3}}$. The net dispersion measure contribution of the host galaxy to FRB 121102 is likely to be lower than this value depending on geometrical factors. We show that the persistent radio source at FRB 121102's location reported by Marcote et al (2017) is offset from the galaxy's center of light by $\sim$200 mas and the host galaxy does not show optical signatures for AGN activity. If FRB 121102 is typical of the wider FRB population and if future interferometric localizations preferentially find them in dwarf galaxies with low metallicities and prominent emission lines, they would share such a preference with long gamma ray bursts and superluminous supernovae.

Fast radio bursts are astronomical radio flashes of unknown physical nature with durations of milliseconds. Their dispersive arrival times suggest an extragalactic origin and imply radio luminosities orders of magnitude larger than any other kind of known short-duration radio transient. Thus far, all FRBs have been detected with large single-dish telescopes with arcminute localizations, and attempts to identify their counterparts (source or host galaxy) have relied on contemporaneous variability of field sources or the presence of peculiar field stars or galaxies. These attempts have not resulted in an unambiguous association with a host or multi-wavelength counterpart. Here we report the sub-arcsecond localization of FRB 121102, the only known repeating burst source, using high-time-resolution radio interferometric observations that directly image the bursts themselves. Our precise localization reveals that FRB 121102 originates within 100 mas of a faint 180 uJy persistent radio source with a continuum spectrum that is consistent with non-thermal emission, and a faint (25th magnitude) optical counterpart. The flux density of the persistent radio source varies by tens of percent on day timescales, and very long baseline radio interferometry yields an angular size less than 1.7 mas. Our observations are inconsistent with the fast radio burst having a Galactic origin or its source being located within a prominent star-forming galaxy. Instead, the source appears to be co-located with a low-luminosity active galactic nucleus or a previously unknown type of extragalactic source. [Truncated] If other fast radio bursts have similarly faint radio and optical counterparts, our findings imply that direct sub-arcsecond localizations of FRBs may be the only way to provide reliable associations.

The millisecond-duration radio flashes known as Fast Radio Bursts (FRBs) represent an enigmatic astrophysical phenomenon. Recently, the sub-arcsecond localization (~ 100mas precision) of FRB121102 using the VLA has led to its unambiguous association with persistent radio and optical counterparts, and to the identification of its host galaxy. However, an even more precise localization is needed in order to probe the direct physical relationship between the millisecond bursts themselves and the associated persistent emission. Here we report very-long-baseline radio interferometric observations using the European VLBI Network and the 305-m Arecibo telescope, which simultaneously detect both the bursts and the persistent radio emission at milliarcsecond angular scales and show that they are co-located to within a projected linear separation of < 40pc (< 12mas angular separation, at 95% confidence). We detect consistent angular broadening of the bursts and persistent radio source (~ 2-4mas at 1.7GHz), which are both similar to the expected Milky Way scattering contribution. The persistent radio source has a projected size constrained to be < 0.7pc (< 0.2mas angular extent at 5.0GHz) and a lower limit for the brightness temperature of T_b > 5 x 10^7K. Together, these observations provide strong evidence for a direct physical link between FRB121102 and the compact persistent radio source. We argue that a burst source associated with a low-luminosity active galactic nucleus or a young neutron star energizing a supernova remnant are the two scenarios for FRB121102 that best match the observed data.

The recent discovery of a population of eccentric (e ~ 0.1) millisecond pulsar (MSP) binaries with low-mass white dwarf companions in the Galactic field represents a challenge to evolutionary models that explain MSP formation as recycling: all such models predict that the orbits become highly circularised during a long period of accretion. The members of this new population exhibit remarkably similar properties (orbital periods, eccentricities, companion masses, spin periods) and several models have been put forward that suggest a common formation channel. In this work we present the results of an extensive timing campaign focusing on one member of this new population, PSR J1946+3417. Through measurement of the both the advance of periastron and Shapiro delay for this system, we determine the mass of the pulsar, companion and the inclination of the orbit to be 1.828(22) Msun, 0.2656(19) Msun and 76.4(6) , under the assumption that general relativity is the true description of gravity. Notably, this is the third highest mass measured for any pulsar. Using these masses and the astrometric properties of PSR J1946+3417 we examine three proposed formation channels for eccentric MSP binaries. While our results are consistent with eccentricity growth driven by a circumbinary disk or neutron star to strange star phase transition, we rule out rotationally delayed accretion-induced collapse as the mechanism responsible for the configuration of the PSR J1946+3417 system.

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 report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that is three times the predicted maximum Galactic value. The 2-GHz bursts have highly variable spectra like those at 1.4 GHz, indicating that the frequency structure seen across the individual 1.4 and 2-GHz bandpasses is part of a wideband process. X-ray observations of the FRB 121102 field with the Swift and Chandra observatories show at least one possible counterpart; however, the probability of chance superposition is high. A radio imaging observation of the field with the Jansky Very Large Array at 1.6 GHz yields a 5$\sigma$ upper limit of 0.3 mJy on any point-source continuum emission. This upper limit, combined with archival WISE 22-$\mu$m and IPHAS H$\alpha$ surveys, rules out the presence of an intervening Galactic HII region. We update our estimate of the FRB detection rate in the PALFA survey to be 1.1$^{+3.7}_{-1.0} \times 10^4$ FRBs sky$^{-1}$ day$^{-1}$ (95% confidence) for peak flux density at 1.4 GHz above 300 mJy. We find that the intrinsic widths of the 12 FRB 121102 bursts from Arecibo are, on average, significantly longer than the intrinsic widths of the 13 single-component FRBs detected with the Parkes telescope.

Transitional millisecond pulsars (tMSPs) switch, on roughly multi-year timescales, between rotation-powered radio millisecond pulsar (RMSP) and accretion-powered low-mass X-ray binary (LMXB) states. The tMSPs have raised several questions related to the nature of accretion flow in their LMXB state and the mechanism that causes the state switch. The discovery of coherent X-ray pulsations from PSR J1023+0038 (while in the LMXB state) provides us with the first opportunity to perform timing observations and to compare the neutron star's spin variation during this state to the measured spin-down in the RMSP state. Whereas the X-ray pulsations in the LMXB state likely indicate that some material is accreting onto the neutron star's magnetic polar caps, radio continuum observations indicate the presence of an outflow. The fraction of the inflowing material being ejected is not clear, but it may be much larger than that reaching the neutron star's surface. Timing observations can measure the total torque on the neutron star. We have phase-connected nine XMM-Newton observations of PSR J1023+0038 over the last 2.5 years of the LMXB state to establish a precise measurement of spin evolution. We find that the average spin-down rate as an LMXB is 26.8+/-0.4% faster than the rate (-2.39x10^-15 Hz s-1) determined during the RMSP state. This shows that negative angular momentum contributions (dipolar magnetic braking and outflow) exceed positive ones (accreted material), and suggests that the pulsar wind continues to operate at a largely unmodified level. We discuss implications of this tight observational constraint in the context of possible accretion models.

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.

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.

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

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 an implementation of the coherent dedispersion algorithm capable of dedispersing high-time-resolution radio observations to many different dispersion measures (DMs). This approach allows the removal of the dispersive effects of the interstellar medium and enables searches for pulsed emission from pulsars and other millisecond-duration transients at low observing frequencies and/or high DMs where time broadening of the signal due to dispersive smearing would otherwise severely reduce the sensitivity. The implementation, called 'cdmt', for Coherent Dispersion Measure Trials, exploits the parallel processing capability of general-purpose graphics processing units to accelerate the computations. We describe the coherent dedispersion algorithm and detail how cdmt implements the algorithm to efficiently compute many coherent DM trials. We present the concept of a semi-coherent dedispersion search, where coherently dedispersed trials at coarsely separated DMs are subsequently incoherently dedispersed at finer steps in DM. The software is used in an ongoing LOFAR pilot survey to test the feasibility of performing semi-coherent dedispersion searches for millisecond pulsars at 135MHz. This pilot survey has led to the discovery of a radio millisecond pulsar -- the first at these low frequencies. This is the first time that such a broad and comprehensive search in DM-space has been done using coherent dedispersion, and we argue that future low-frequency pulsar searches using this approach are both scientifically compelling and feasible. Finally, we compare the performance of cdmt with other available alternatives.

As part of the Chandra Galactic Bulge Survey (GBS), we present a catalogue of optical sources in the GBS footprint. This consists of two regions centered at Galactic latitude b = 1.5 degrees above and below the Galactic Centre, spanning (l x b) = (6x1) degrees. The catalogue consists of 2 or more epochs of observations for each line of sight in r', i' and H{\alpha} filters. It is complete down to r' = 20.2 and i' = 19.2 mag; the mean 5{\sigma} depth is r' = 22.5 and i' = 21.1 mag. The mean root-mean-square residuals of the astrometric solutions is 0.04 arcsec. We cross-correlate this optical catalogue with the 1640 unique X-ray sources detected in Chandra observations of the GBS area, and find candidate optical counterparts to 1480 X-ray sources. We use a false alarm probability analysis to estimate the contamination by interlopers, and expect ~ 10 per cent of optical counterparts to be chance alignments. To determine the most likely counterpart for each X-ray source, we compute the likelihood ratio for all optical sources within the 4{\sigma} X-ray error circle. This analysis yields 1480 potential counterparts (~ 90 per cent of the sample). 584 counterparts have saturated photometry (r'<17, i'<16), indicating these objects are likely foreground sources and the real counterparts. 171 candidate counterparts are detected only in the i'-band. These sources are good qLMXB and CV candidates as they are X-ray bright and likely located in the Bulge.

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

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

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

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

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 report on photometric and spectroscopic observations of white dwarf companions to four binary radio millisecond pulsars, leading to the discovery of companions to PSRs J0614-3329, J1231-1411 and J2017+0603. We place limits on the brightness of the companion to PSR J0613-0200. Optical spectroscopy of the companion to PSR J0614-3329 identifies it as a DA type white dwarf with a temperature of Teff=6460+-80 K, a surface gravity log g=7.0+-0.2 cgs and a mass of Mwd=0.24+-0.04 Msun. We find that the distance to PSR J0614-3329 is smaller than previously estimated, removing the need for the pulsar to have an unrealistically high gamma-ray efficiency. Comparing the photometry with predictions from white dwarf cooling models allows us to estimate temperatures and cooling ages of the companions to PSRs J0613-0200, J1231-1411 and J2017+0603. We find that the white dwarfs in these systems are cool Teff<4000 K and old >5 Gyr. Thin Hydrogen envelopes are required for these white dwarfs to cool to the observed temperatures, and we suggest that besides Hydrogen shell flashes, irradiation driven mass loss by the pulsar may have been important.

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.

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.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Accreting millisecond X-ray pulsars are an important subset of low-mass X-ray binaries in which coherent X-ray pulsations can be observed during occasional, bright outbursts (X-ray luminosity $L_X\sim 10^{36}$ erg s$^{-1}$). These pulsations show that matter is being channeled onto the neutron star's magnetic poles. However, such sources spend most of their time in a low-luminosity, quiescent state ($L_X\lesssim 10^{34}$ erg s$^{-1}$), where the nature of the accretion flow onto the neutron star (if any) is not well understood. Here we report that the millisecond pulsar/low-mass X-ray binary transition object PSR J1023+0038 intermittently shows coherent X-ray pulsations at luminosities nearly 100 times fainter than observed in any other accreting millisecond X-ray pulsar. We conclude that in spite of its low luminosity PSR J1023+0038 experiences episodes of channeled accretion, a discovery that challenges existing models for accretion onto magnetized neutron stars.

The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary. We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent low-mass X-ray binary state. Two long XMM-Newton observations reveal that the system spends ~70% of the time in a $\approx$$3\times10^{33}$ erg/s X-ray luminosity mode, which, as shown in Archibald et al. (2014), exhibits coherent X-ray pulsations. This emission is interspersed with frequent lower flux mode intervals with $\approx$$5\times 10^{32}$ erg/s and sporadic flares reaching up to $\approx$$10^{34}$ erg/s, with neither mode showing significant X-ray pulsations. The switches between the three flux modes occur on timescales of order 10 s. In the UV and optical, we observe occasional intense flares coincident with those observed in X-rays. Our radio timing observations reveal no pulsations at the pulsar period during any of the three X-ray modes, presumably due to complete quenching of the radio emission mechanism by the accretion flow. Radio imaging detects highly variable, flat-spectrum continuum emission from PSR J1023+0038, consistent with an origin in a weak jet-like outflow. Our concurrent X-ray and radio continuum data sets do not exhibit any correlated behavior. The observational evidence we present bears qualitative resemblance to the behavior predicted by some existing propeller and trapped disk accretion models although none can fully account for all aspects of the rich phenomenology of this system.

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

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering.

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.

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.

We study the spin-down properties of PSR B0919+06 based on almost 30 years of radio observations. We confirm that the time derivative of the rotational frequency $\dot \nu$ is modulated quasi-periodically and show that it exhibits a repeating double-peaked structure throughout the entire observation span. We model the $\dot \nu$ variation of the pulsar assuming two spin-down rates with sudden switches between them in time. Our results show that the double-peak structure in $\dot \nu$ has a repetition time of about 630 days until MJD 52000 (April 2001) and 550 days since then. During this cycle, the pulsar spin varies from the lower spin-down rate to the upper spin-down rate twice with different amounts of time spent in each state, resulting in a further quasi-stable secondary modulation of the two-state switching. This particular spin-down state switching is broadly consistent with free precession of the pulsar, however, a strong evidence linked with this mechanism is not clearly established. We also confirm that the pulsar occasionally emits groups of pulses which appear early in pulse phase, so-called "flares", and these events significantly contribute to the pulse profile shape. We find the $\dot \nu$ modulation and the pulse shape variations are correlated throughout the observations. However, the flare-state is not entirely responsible for this correlation. In addition to the flare-state, we detect flare-like events from the pulsar in single pulse observations. During these events, the shift in pulse phase is small compared to that of the main flare-state and clearly visible only in single pulse observations.

We present the results of a search for gravitational waves (GWs) from individual sources using high cadence observations of PSR B1937+21. The data were acquired from an intensive observation campaign with the Lovell telescope at Jodrell Bank, between June 2011 and May 2013. The almost daily cadence achieved, allowed us to be sensitive to GWs with frequencies up to $4.98\times10^{-6}\,\rm {Hz}$, extending the upper bound of the typical frequency range probed by Pulsar Timing Arrays. We used observations taken at three different radio frequencies with the Westerbork Synthesis Radio Telescope in order to correct for dispersion measure effects and scattering variances. The corrected timing residuals exhibited an unmodeled periodic noise with an amplitude $~150\,\rm {ns}$ and a frequency of $3.4\rm {yr}^{-1}$. As the signal is not present in the entire data set, we attributed it to the rotational behaviour of the pulsar, ruling out the possibilities of being either due to a GW or an asteroid as the cause. After removing this noise component, we placed limits on the GW strain of individual sources equaling to $h_{\rm s}=1.53\times10^{-11}$ and $h_{\rm s}=4.99\times10^{-14}$ at $10^{-7}\,\rm {Hz}$ for random and optimal sources locations respectively.

We present optical lightcurves of variable stars consistent with the positions of X-ray sources identified with the Chandra X-ray Observatory for the Chandra Galactic Bulge Survey. Using data from the Mosaic-II instrument on the Blanco 4m Telescope at CTIO, we gathered time-resolved photometric data on timescales from $\sim2$ hr to 8 days over the $\frac{3}{4}$ of the X-ray survey containing sources from the initial GBS catalog. Among the lightcurve morphologies we identify are flickering in interacting binaries, eclipsing sources, dwarf nova outbursts, ellipsoidal variations, long period variables, spotted stars, and flare stars. $87\%$ of X-ray sources have at least one potential optical counterpart. $24\%$ of these candidate counterparts are detectably variable; a much greater fraction than expected for randomly selected field stars, which suggests that most of these variables are real counterparts. We discuss individual sources of interest, provide variability information on candidate counterparts, and discuss the characteristics of the variable population.

We report NuSTAR observations of the millisecond pulsar - low mass X-ray binary (LMXB) transition system PSR J1023+0038 from June and October 2013, before and after the formation of an accretion disk around the neutron star. Between June 10-12, a few days to two weeks before the radio disappearance of the pulsar, the 3-79 keV X-ray spectrum was well fit by a simple power law with a photon index of Gamma=1.17 +/-0.08 (at 90% confidence) with a 3-79 keV luminosity of 7.4+/-0.4 x 10^32 erg/s. Significant orbital modulation was observed with a modulation fraction of 36+/-10%. During the October 19-21 observation, the spectrum is described by a softer power law (Gamma=1.66+/-0.06) with an average luminosity of 5.8+/-0.2 x 10^33 erg/s and a peak luminosity of ~1.2 x 10^34 erg/s observed during a flare. No significant orbital modulation was detected. The spectral observations are consistent with previous and current multi-wavelength observations and show the hard X-ray power law extending to 79 keV without a spectral break. Sharp edged, flat bottomed `dips' are observed with widths between 30-1000 s and ingress and egress time-scales of 30-60 s. No change in hardness ratio was observed during the dips. Consecutive dip separations are log-normal in distribution with a typical separation of approximately 400 s. These dips are distinct from dipping activity observed in LMXBs. We compare and contrast these dips to observations of dips and state changes in the similar transition systems PSR J1824-2452I and XSS J1227.0-4859 and discuss possible interpretations based on the transitions in the inner disk.

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

We present XMM-Newton and Chandra observations of the low-mass X-ray binary XSS J12270--4859, which experienced a dramatic decline in optical/X-ray brightness at the end of 2012, indicative of the disappearance of its accretion disk. In this new state, the system exhibits previously absent orbital-phase-dependent, large-amplitude X-ray modulations with a decline in flux at superior conjunction. The X-ray emission remains predominantly non-thermal but with an order of magnitude lower mean luminosity and significantly harder spectrum relative to the previous high flux state. This phenomenology is identical to the behavior of the radio millisecond pulsar binary PSR J1023+0038 in the absence of an accretion disk, where the X-ray emission is produced in an intra-binary shock driven by the pulsar wind. This further demonstrates that XSS J12270-4859 no longer has an accretion disk and has transformed to a full-fledged eclipsing "redback" system that hosts an active rotation-powered millisecond pulsar. There is no evidence for diffuse X-ray emission associated with the binary that may arise due to outflows or a wind nebula. An extended source situated 1.5' from XSS J12270--4859 is unlikely to be associated, and is probably a previously uncatalogued galaxy cluster.

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.

Millisecond radio pulsars acquire their rapid rotation rates through mass and angular momentum transfer in a low-mass X-ray binary system. Recent studies of PSR J1824-2452I and PSR J1023+0038 have observationally demonstrated this link, and they have also shown that such systems can repeatedly transition back-and-forth between the radio millisecond pulsar and low-mass X-ray binary states. This also suggests that a fraction of such systems are not newly born radio millisecond pulsars but are rather suspended in a back-and-forth state switching phase, perhaps for giga-years. XSS J12270-4859 has been previously suggested to be a low-mass X-ray binary, and until recently the only such system to be seen at MeV-GeV energies. We present radio, optical and X-ray observations that offer compelling evidence that XSS J12270-4859 is a low-mass X-ray binary which transitioned to a radio millisecond pulsar state between 2012 November 14 and 2012 December 21. Though radio pulsations remain to be detected, we use optical and X-ray photometry/spectroscopy to show that the system has undergone a sudden dimming and no longer shows evidence for an accretion disk. The optical observations constrain the orbital period to 6.913+-0.002 hr.

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 present radio, X-ray, and gamma-ray observations which reveal that the binary millisecond pulsar / low-mass X-ray binary transition system PSR J1023+0038 has undergone a transformation in state. Whereas until recently the system harbored a bright millisecond radio pulsar, the radio pulsations at frequencies between 300 to 5000 MHz have now become undetectable. Concurrent with this radio disappearance, the gamma-ray flux of the system has quintupled. We conclude that, though the radio pulsar is currently not detectable, the pulsar mechanism is still active and the pulsar wind, as well as a newly formed accretion disk, are together providing the necessary conditions to create the gamma-ray increase. The system is the first example of a transient, compact, low-mass gamma-ray binary and will continue to provide an exceptional test bed for better understanding the formation of millisecond pulsars as well as accretion onto neutron stars in general.

PSR J1023+0038 is an exceptional system for understanding how slowly rotating neutron stars are spun up to millisecond rotational periods through accretion from a companion star. Observed as a radio pulsar from 2007-2013, optical data showed that the system had an accretion disk in 2000/2001. Starting at the end of 2013 June, the radio pulsar has become undetectable, suggesting a return to the previous accretion-disk state, where the system more closely resembles an X-ray binary. In this Letter we report the first targeted X-ray observations ever performed of the active phase and complement them with UV/Optical and radio observations collected in 2013 October. We find strong evidence that indeed an accretion disk has recently formed in the system and we report the detection of fast X-ray changes spanning about two orders of magnitude in luminosity. No radio pulsations are seen during low flux states in the X-ray light-curve or at any other times.

Pulsars are highly-magnetised rotating neutron stars and are well-known for the stability of their signature pulse shapes, allowing high-precision studies of their rotation. However, during the past 22 years, the radio pulse profile of the Crab pulsar has shown a steady increase in the separation of the main pulse and interpulse components at 0.62$^{\rm o}\pm$0.03$^{\rm o}$ per century. There are also secular changes in the relative strengths of several components of the profile. The changing component separation indicates that the axis of the dipolar magnetic field, embedded in the neutron star, is moving towards the stellar equator. This evolution of the magnetic field could explain why the pulsar does not spin down as expected from simple braking by a rotating dipolar magnetic field.

(Abridged:) We present the identification of optical counterparts to 23 Galactic Bulge Survey X-ray sources. We report their accurate coordinates and optical spectra acquired at the VLT and Magellan. All sources are classified as accreting binaries according to their emission line characteristics. To distinguish accreting binaries from chromospherically active objects we develop and explain criteria based on Halpha and HeI 5786,6678 emission line properties available in the literature. The spectroscopic properties and photometric variability of all the objects are discussed and a classification of the source is given where possible. Among the 23 systems, at least 9 of them show an accretion-dominated optical spectrum (CX28, CX63, CX70, CX128, CX142, CX207, CX522, CX794, CX1011) and another 6 show photospheric lines from a late-type donor star in addition to accretion disc emission (CX44, CX93, CX137, CX154, CX377 and CX1004) indicating that they are probably accreting binaries in quiescence or in a low accretion rate state. Two sources are confirmed to be eclipsing: CX207 and CX794. CX207 shows a broad asymmetric Halpha profile blue-shifted by >300 km/s. Such line profile characteristics are consistent with a magnetic (Polar) CV. CX794 is an eclipsing nova-like CV in the period gap. Time-resolved photometry and the large broadening of the Halpha emission lines in CX446 (2100 km/s FWHM) suggest that this is also an eclipsing or high-inclination accreting binary. Finally, the low-accretion rate source CX1004 shows a double-peaked Halpha profile with a FWHM of 2100 km/s. This supports a high inclination or even eclipsing system. Whether the compact object is a white dwarf in an eclipsing CV or a black hole primary in a high-inclination LMXB remains to be established.

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

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.

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

We present upper limits on the X-ray emission for three neutron stars. For PSR J1840$-$1419, with a characteristic age of 16.5 Myr, we calculate a blackbody temperature upper limit (at 99% confidence) of $kT_{\mathrm{bb}}^{\infty}<24^{+17}_{-10}$ eV, making this one of the coolest neutron stars known. PSRs J1814$-$1744 and J1847$-$0130 are both high magnetic field pulsars, with inferred surface dipole magnetic field strengths of $5.5\times10^{13}$ and $9.4\times10^{13}$ G, respectively. Our temperature upper limits for these stars are $kT_{\mathrm{bb}}^{\infty}<123^{+20}_{-33}$ eV and $kT_{\mathrm{bb}}^{\infty}<115^{+16}_{-33}$ eV, showing that these high magnetic field pulsars are not significantly hotter than those with lower magnetic fields. Finally, we put these limits into context by summarizing all temperature measurements and limits for rotation-driven neutron stars.

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.

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest orbital period, making it especially interesting for a range of astrophysical and gravity related questions. We present a spectroscopic and photometric study of the white dwarf companion and infer its radial velocity curve, effective temperature, surface gravity and luminosity. We find that the white dwarf has properties consistent with those of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding mass-radius relation to infer its mass; M_WD = 0.181 +/- +0.007/-0.005 solar masses. Combined with the mass ratio q=8.1 +/- 0.2 inferred from the radial velocities and the precise pulsar timing ephemeris, the neutron star mass is constrained to M_PSR = 1.47 +/- +0.07/-0.06 solar masses. Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. We find that, even if the birth mass of the neutron star was only 1.20 solar masses, more than 60% of the matter that left the surface of the white dwarf progenitor escaped the system. The accurate determination of the component masses transforms this system in a laboratory for fundamental physics by constraining the orbital decay predicted by general relativity. Currently, the agreement is within 1 sigma of the observed decay. Further radio timing observations will allow precise tests of white dwarf models, assuming the validity of general relativity.

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.

We present phase-resolved optical spectroscopy of the counterpart to the neutron star low mass X-ray binary EXO 0748-676, almost one year after it turned into quiescence. The spectra display prominent Hbeta and Hgamma and weak Fe II lines in emission. An average of all the spectra (corrected for the orbital motion) also exhibits a very weak line from Mg I. Tomographic reconstructions show that the accretion disc is not contributing to the optical line emission, which is instead dominated by the irradiated hemisphere of the companion star facing the neutron star. We could not detect absorption features from the mass donor star in the spectra. The emission lines appear broad, with an intrinsic FWHM of 255+-22 km/s. Under the assumption that the width of the Fe II emission lines is dominated by rotational broadening, we obtain a lower limit on the compact object mass which is inconsistent with a NS accretor. We discuss this incongruity and conclude that either the lines are blends of unresolved features (although this requires some fine tuning) or they are broadened by additional effects such as bulk gas motion in an outflow. The fact that the Fe II lines slightly lag in phase with respect to the companion star can be understood as outflowing gas consistent with a Black-Widow like scenario. Nevertheless, we can not rule out the possibility that blends of various emission lines cause the apparent phase lag of the Fe II emission lines as well as their large width.

We study the electron-proton temperature equilibration behind several shocks of the RCW 86 supernova remnant. To measure the proton temperature, we use published and new optical spectra, all from different locations on the remnant. For each location, we determine the electron temperature from X-ray spectra, and correct for temperature equilibration between the shock front and the location of the X-ray spectrum. We confirm the result of previous studies that the electron and proton temperatures behind shock fronts are consistent with equilibration for slow shocks and deviate for faster shocks. However, we can not confirm the previously reported trend of the electron temperature to proton temperature ratio of 1/v^2.

We introduce the Galactic Bulge Survey (GBS) and we provide the Chandra source list for the region that has been observed to date. Among the goals of the GBS are constraining the neutron star equation of state and the black hole mass distribution via the identification of eclipsing neutron star and black hole low-mass X-ray binaries. The latter goal will, in addition, be obtained by significantly enlarging the number of black hole systems for which a black hole mass can be derived. Further goals include constraining X-ray binary formation scenarios, in particular the common envelope phase and the occurrence of kicks, via source-type number counts and an investigation of the spatial distribution of X-ray binaries, respectively. The GBS targets two strips of 6x1 degrees (12 square degrees in total), one above (1<b<2 degrees) and one below (-2<b<-1 degrees) the Galactic plane in the direction of the Galactic Center at both X-ray and optical wavelengths. By avoiding the Galactic plane (-1<b<1 degrees) we limit the influence of extinction on the X-ray and optical emission but still sample relatively large number densities of sources. The survey is designed such that a large fraction of the X-ray sources can be identified from their optical spectra. The X-ray survey, by design, covers a large area on the sky while the depth is shallow using 2 ks per Chandra pointing. In this way we maximize the predicted number ratio of (quiescent) low-mass X-ray binaries to Cataclysmic Variables. The survey is approximately homogeneous in depth to an 0.5-10 keV flux of 7.7x10^-14 erg cm-2 s-1. So far, we have covered about two-thirds (8.3 square degrees) of the projected survey area with Chandra providing over 1200 unique X-ray sources. We discuss the characteristics and the variability of the brightest of these sources.

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.

We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.

We report on the identification of a near-infrared counterpart to the massive (>11 Msun) binary companion of pulsar J1740-3052. An accurate celestial position of PSR J1740-3052 is determined from interferometric radio observations. Adaptive optics corrected near-infrared imaging observations show a counterpart at the interferometric position of the pulsar. The counterpart has Ks=15.87+-0.10 and J-Ks>0.83. Based on distance and absorption estimates from models of the Galactic electron and dust distributions these observed magnitudes are consistent with those of a main-sequence star as the binary companion. We argue that this counterpart is the binary companion to PSR J1740-3052 and thus rule out a stellar mass black hole as the pulsar companion.

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

Highly sensitive and precise X-ray imaging from Chandra, combined with the superb spatial resolution of HST optical images, dramatically enhances our empirical understanding of compact binaries such as cataclysmic variables and low mass X-ray binaries, their progeny, and other stellar X-ray source populations deep into the cores of globular clusters. Our Chandra X-ray images of the globular cluster NGC 362 reveal 100 X-ray sources, the bulk of which are likely cluster members. Using HST color-magnitude and color-color diagrams, we quantitatively consider the optical content of the NGC 362 Chandra X-ray error circles, especially to assess and identify the compact binary population in this condensed-core globular cluster. Despite residual significant crowding in both X-rays and optical, we identify an excess population of H{\alpha}-emitting objects that is statistically associated with the Chandra X-ray sources. The X-ray and optical characteristics suggest that these are mainly cataclysmic variables, but we also identify a candidate quiescent low mass X-ray binary. A potentially interesting and largely unanticipated use of observations such as these may be to help constrain the macroscopic dynamic state of globular clusters.

Pulsars with compact companions in close eccentric orbits are unique laboratories for testing general relativity and alternative theories of gravity. Moreover, they are excellent targets for future gravitational wave experiments like LISA and they are also highly important for understanding the equation of state of super-dense matter and the evolution of massive binaries. Here we report on optical observations of the 1.02 M_sun companion to the pulsar PSR J1141-6545. We detect an optical counterpart with apparent magnitudes V=25.08(11) and R=24.38(14), consistent with the timing position of the pulsar. We demonstrate that our results are in agreement with a white dwarf companion. However the latter is redder than expected and the inferred values are not consistent with the theoretical cooling tracks, preventing us from deriving the exact age. Our results confirm the importance of the PSR J1141-6545 system for gravitational experiments.

We investigate a sample of eleven Galactic X-ray sources recently discovered with INTEGRAL or RXTE with the goal of identifying their optical and/or near-infrared (NIR) counterpart. For this purpose new Chandra positions of nine objects are presented together with follow-up observations of all the targets in the optical and NIR. For the four sources IGR J16194-2810, IGRJ 16479-4514, IGR J16500-3307 and IGR J19308+530, the Chandra position confirms an existing association with an optical/NIR object, while for two sources (XTE J1716-389 and 18490-0000) it rules out previously proposed counterparts indicating new ones. In the case of IGR J17597-220, a counterpart is selected out of the several possibilities proposed in the literature and we present the first association with an optical/NIR source for J16293-4603 and XTE J1743-363. Moreover, optical/NIR observations are reported for XTE J1710-281 and IGR J17254-3257: we investigate the counterpart to the X-ray sources based on their XMM-Newton positions. We discuss the nature of each system considering its optical/NIR and X-ray properties.

In August 2008, the accreting milli-second X-ray pulsar (AMXP), IGR J00291+5934, underwent an outburst lasting ~ 100 days, the first since its discovery in 2004. We present data from the double-peaked outburst from Faulkes Telescope North, the INT, the Keck Telescope, PAIRITEL, the Westerbork Synthesis Radio Telescope and the Swift, XMM-Newton and RXTE X-ray missions. We study the outburst's evolution at various wavelengths. We study the light curve morphology, presenting the first radio-X-ray Spectral Energy Distributions (SEDs) for this source and the most detailed UV-IR SEDs for any outbursting AMXP. We show simple models that attempt to identify the emission mechanisms responsible. We analyse short-timescale optical variability, and compare a medium resolution optical spectrum with those from 2004. The outburst morphology is unusual for an AMXP, comprising two peaks, the second containing a 'plateau' of ~ 10 days at maximum brightness within 30 days of the initial activity. This has implications on duty cycles of short-period X-ray transients. The X-ray spectrum can be fitted by a single, hard power-law. We detect optical variability of ~ 0.05 magnitudes, on timescales of minutes, but find no periodic modulation. In the optical, the SEDs contain a blue component, indicative of an irradiated disc, and a transient near-infrared (NIR) excess. This excess is consistent with a simple model of an optically thick synchrotron jet (as seen in other outbursting AMXPs). The optical spectrum shows a double-peaked H alpha profile, a diagnostic of an accretion disc, but we do not clearly see other lines (e.g. He I, II) reported in 2004. Optical/IR observations of AMXPs are excellent for studying the evolution of both the outer accretion disc and the inner jet, and may eventually provide us with tight constraints to model disc-jet coupling in accreting neutron stars.

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.

We have timed four millisecond pulses, PSRs J1721-2457, J1745-0952, J1810-2005, and J1918-0642, for up to a total of 10.5 years each using multiple telescopes in the European Pulsar Timing Array network: the Westerbork Synthesis Radio Telescope in The Netherlands, the Nancay Radio Telescope in France and the Lovell telescope at Jodrell Bank in the UK. The long time span has enabled us to measure the proper motions of J1745-0952 and J1918-0642, indicating that they have transverse velocities of 200(50) and 54(7) km/s respectively. We have obtained upper limits on the proper motion of J1721-2457 and J1810-2005, which imply that they have transverse velocities less than 140 and 400 km/s respectively. In all cases, the velocities lie in the range typical of millisecond pulsars. We present pulse profiles for each pulsar taken from observations at multiple frequencies in the range of 350 to 2600 MHz, and show that J1810-2005 shows significant profile evolution in this range. Using our multi-frequency observations, we measured the spectral indices for all four pulsars, and for J1810-2005 it appears to be very flat. The flux density of J1918-0642 shows extensive modulation which we attribute to the combined effects of refractive and diffractive scintillation. We discuss the possible use of including J1721-2457 or J1918-0642 in a pulsar timing array, and find that J1918-0642 will be useful to include when the timing precision of this pulsar is improved over the next few years. We have searched archival optical observations to detect companions of the binary pulsars, but none were detected. However, we provide lower limits on the masses of the white dwarf companions of PSRs J1745-0952 and J1918-0642.

We report on the Chandra X-ray Observatory and Hubble Space Telescope observation of two low core density globular clusters, NGC6144 and E3. By comparing the number of X-ray sources inside the half-mass radius to those outside, we found 6 X-ray sources within the half-mass radius of NGC6144, among which 4 are expected to be background sources; 3 X-ray sources are also found within the half-mass radius of E3, of which 3 is expected to be background source. Therefore, we cannot exclude that all our sources are background sources. However, combining the results from X-ray and optical observations, we found that 1-2 sources in NGC6144 and 1 source in E3 are likely to be cataclysmic variables and that 1 source in NGC6144 is an active binary, based on the X-ray and optical properties. The number of faint X-ray sources in NGC6144 and E3 found with Chandra and HST is higher than a prediction based on collision frequency, but is closer to that based on mass. Our observations strongly suggest that the compact binary systems in NGC6144 and E3 are primordial in origin.

We study a Chandra X-ray Observatory ACIS-S observation of the Galactic globular clusterM12. With a 26 ks exposure time, we detect 6 X-ray sources inside the half-mass radius (2.16 arcminutes) of which two are inside the core radius (0.72 arcminutes) of the cluster. If we assume these sources are all associated with globular cluster M12, the luminosity L_{X} among these sources between 0.3-7.0 keV varies roughly from 10^{30} to 10^{32} ergs s^{-1}. For identification, we also analyzed the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field and Planetary Camera 2 (WFPC2) data and identified the optical counterparts to five X-ray sources inside the HST ACS field of view. According to the X-ray and optical features, we found 2-5 candidate active binaries (ABs) or cataclysmic variables (CVs) and 0-3 background galaxies within the HST ACS field of view. Based on the assumption that the number of X-ray sources scales with the encounter rate and the mass of the globular cluster, we expect 2 X-ray source inside M12, and the expectation is consistent with our observational results. Therefore, the existence of identified X-ray sources (possible CVs or ABs) in M12 suggests the primordial origin of X-ray sources in globular clusters which is in agreement with previous studies.

An X-ray monitoring campaign with Chandra and Swift confirms that 1RXS J171824.2-402934 is accreting at the lowest rate among the known persistently accreting low-mass X-ray binaries. A thermonuclear X-ray burst was detected with the all-sky monitor on RXTE. This is only the second such burst seen in 1RXS J171824.2-402934 in more than 20 Ms of observations done over 19 years. The low burst recurrence rate is in line with the low accretion rate. The persistent nature and low accretion rate can be reconciled within accretion disk theory if the binary system is ultracompact. An unprecedentedly short orbital period of less than approximately 7 minutes would be implied. An ultracompact nature, together with the properties of the type I X-ray burst, suggests, in turn, that helium-rich material is accreted. Optical follow-up of the Chandra error region does not reveal an unambiguous counterpart.

We present phase resolved optical spectroscopy and X-ray timing of the neutron star X-ray binary EXO 0748-676 after the source returned to quiescence in the fall of 2008. The X-ray light curve displays eclipses consistent in orbital period, orbital phase and duration with the predictions and measurements before the return to quiescence. Halpha and He I emission lines are present in the optical spectra and show the signature of the orbit of the binary companion, placing a lower limit on the radial velocity semi-amplitude of K2>405 km/s. Both the flux in the continuum and the emission lines show orbital modulations, indicating that we observe the hemisphere of the binary companion that is being irradiated by the neutron star. Effects due to this irradiation preclude a direct measurement of the radial velocity semi-amplitude of the binary companion; in fact no stellar absorption lines are seen in the spectrum. Nevertheless, our observations place a stringent lower limit on the neutron star mass of M1>1.27 Msun. For the canonical neutron star mass of M1=1.4 Msun, the mass ratio is constrained to 0.075<q<0.105.

Cosmic rays are the most energetic particles arriving at earth. Although most of them are thought to be accelerated by supernova remnants, the details of the acceleration process and its efficiency are not well determined. Here we show that the pressure induced by cosmic rays exceeds the thermal pressure behind the northeast shock of the supernova remnant RCW 86, where the X-ray emission is dominated by synchrotron radiation from ultra-relativistic electrons. We determined the cosmic-ray content from the thermal Doppler broadening measured with optical spectroscopy, combined with a proper-motion study in X- rays. The measured post-shock proton temperature in combination with the shock velocity does not agree with standard shock heating, implying that >50% of the post-shock pressure is produced by cosmic rays.

We report on optical imaging of the X-ray binary SAX J1808.4-3658 with the 8-m Gemini South Telescope. The binary, containing an accretion-powered millisecond pulsar, appears to have a large periodic modulation in its quiescent optical emission. In order to clarify the origin of this modulation, we obtained three time-resolved $r'$-band light curves (LCs) of the source in five days. The LCs can be described by a sinusoid, and the long time-span between them allows us to determine optical period P=7251.9 s and phase 0.671 at MJD 54599.0 (TDB; phase 0.0 corresponds to the ascending node of the pulsar orbit), with uncertainties of 2.8 s and 0.008 (90 % confidence), respectively. This periodicity is highly consistent with the X-ray orbital ephemeris. By considering this consistency and the sinusoidal shape of the LCs, we rule out the possibility of the modulation arising from the accretion disk. Our study supports the previous suggestion that the X-ray pulsar becomes rotationally powered in quiescence, with its energy output irradiating the companion star, causing the optical modulation. While it has also been suggested that the accretion disk would be evaporated by the pulsar, we argue that the disk exists and gives rise to the persistent optical emission. The existence of the disk can be verified by long-term, multi-wavelength optical monitoring of the source in quiescence, as an increasing flux and spectral changes from the source would be expected based on the standard disk instability model.

We present new Chandra X-ray observations and detailed astrometry of the field of the type Ia supernova 2007on, for which the detection of a likely progenitor in archival Chandra data was recently reported. No source is detected in the new Chandra images, taken six weeks after optical maximum. We calculate a 90-99% probability that any X-ray source near the position of the supernova (SN) is fainter than in the pre-outburst images, depending on the choice of aperture, which supports the identification of the archival X-ray source with the SN. Detailed astrometry of the X-ray and new optical images, however, gives an offset between the supernova and the measured X-ray source position of 1.15+/-0.27". Extensive simulations show that the probability of finding an offset of this magnitude is ~1%, equal to the (trial-corrected) probability of a chance alignment with any X-ray source in the field. This casts doubt on the identification of the X-ray source with the progenitor, although the scenario in which at least some of the observed X-rays are connected to the supernova may be the least unlikely based on all available data. After a brief review of the auxiliary evidence, we conclude that only future X-ray observations can shed further light on the proposed connection between the X-ray source and the progenitor of SN 2007on, and thus whether an accreting white dwarf scenario is truly favoured for this SN Ia.

PSR B1931+24 is the first intermittent radio pulsar discovered to date, characterized by a 0.8 s pulsation which turns on and off quasi-periodically every ~35 days, with a duty cycle of ~10%. We present here X-ray and optical observations of PSR B1931+24 performed with the Chandra X-ray Observatory and Isaac Newton Telescope, respectively. Simultaneous monitoring from the Jodrell Bank Observatory showed that this intermittent pulsar was in the radio-on phase during our observations. We do not find any X-ray or optical counterpart to PSR B1931+24 translating into an upper limit of ~2x10^{31} erg/s on the X-ray luminosity, and of g' > 22.6 on the optical magnitude. If the pulsar is isolated, these limits cannot constrain the dim X-ray and optical emission expected for a pulsar of that age (~1.6 Myr). We discuss the possibility that the quasi-periodic intermittent behavior of PSR B1931+24 is due to the presence of a low mass companion star or gaseous planet, tight with the pulsar in an eccentric orbit. In order to constrain the parameters of this putative binary system we re-analysed the pulsar radio timing residuals and we found that (if indeed hosted in a binary system), PSR B1931+24 should have a very low mass companion and an orbit of low inclination.

We present Chandra X-ray Observatory ACIS-S3 X-ray imaging observations and VLT/FORS2 and Hubble Space Telescope optical observations of two low-density Galactic globular clusters; NGC 6366 and M55. We detect 16 X-ray sources with 0.5-6.0 keV luminosities above Lx=4E30 erg/s within the half-mass radius of M55, of which 8 or 9 are expected to be background sources, and 5 within the half-mass radius of NGC 6366, of which 4 are expected to be background sources. Optical counterparts are identified for several X-ray sources in both clusters and from these we conclude that 3 of the X-ray sources in M55 and 2 or 3 of the X-ray sources in NGC 6366 are probably related to the cluster. Combining these results with those for other clusters, we find the best fit for a predicted number of X-ray sources in a globular cluster Nc=1.2 Gamma+1.1 Mh, where Gamma is the collision number and Mh is (half of) the cluster mass, both normalized to the values for the globular cluster M4. Some sources tentatively classified as magnetically active binaries are more luminous in X-rays than the upper limit of Lx~0.001Lbol of such binaries in the solar neighbourhood. Comparison with XMM and ROSAT observations lead us to conclude that the brightest X-ray source in M55, a dwarf nova, becomes fainter in X-rays during the optical outburst, in accordance with other dwarf novae. The brightest X-ray source in NGC 6366 is a point source surrounded by a slightly offset extended source. The absence of galaxies and Halpha emission in our optical observations argues against a cluster of galaxies and against a planetary nebula, and we suggest that the source may be an old nova.

Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.

We present HST/WFPC2, GALEX and Chandra observations of the position of the type Ia supernova SN2007sr in the Antennae galaxies, taken before the explosion. No source is found in any of the observations, allowing us to put interesting constraints on the progenitor luminosity. In total there is about 450 kilosecond of Chandra data, spread over 7 different observations. Limiting magnitudes of FUV (23.7 AB mag), NUV (23.8 AB mag), F555W (26.5 Vega mag) and F814W (24.5-25 Vega mag) are derived. The distance to the Antennae is surprisingly poorly known, with almost a factor of 2 difference between the latest distance based on the tip of the red giant branch (13.3 Mpc) and the distance derived from the SN2007sr light curve (25 Mpc). Using these distances we derive limits on absolute optical and UV magnitudes of any progenitor but these are still above the brightest (symbiotic) proposed progenitors. From the Chandra data a 3 sigma upper limit to the X-ray luminosity of 0.5 -- 8.0 10^37 erg/s in the 0.3-1 keV range is found. This is below the X-ray luminosity of the potential progenitor of the type Ia supernova 2007on that we recently discovered and for which we report a corrected X-ray luminosity. If that progenitor is confirmed it suggests the two supernovae have different progenitors. The X-ray limit is comparable to the brightest supersoft X-ray sources in the Galaxy, the LMC and the SMC and significantly below the luminosities of the brightest supersoft and quasi-soft X-ray sources found in nearby galaxies, ruling out such sources as progenitors of this type Ia supernova.

Results are discussed of an X-ray and optical observation campaign of the low-mass X-ray binary A1246-58 performed with instruments on Satellite per Astronomia X ('BeppoSAX'), the Rossi X-ray Timing Explorer (RXTE), the X-ray Multi-mirror Mission ('XMM-Newton'), the Swift mission, and the Very Large Telescope. Spectra and flux time histories are studied. The most important results are the lack of hydrogen spectral features in the optical spectrum, supporting the proposition that this is an ultracompact X-ray binary (UCXB), the determination of a 4.3 kpc distance from time-resolved spectroscopy of thermonuclear X-ray bursts, and the detection of intermediately long thermonuclear bursts as seen in a number of other UCXBs. There is evidence for a Ne/O abundance ratio in the line of sight that is higher than solar and variable. This may be due to different changes in the ionization degrees of Ne and O, which may be related to the variable irradiating flux. We discuss the spectral variability and the peculiarities of the long-term light curve.

We report on optical and infrared observations of the anomalous X-ray pulsar (AXP) 1E 1048.1-5937, made during its ongoing X-ray flare which started in 2007 March. We detected the source in the optical I and near-infrared Ks bands in two ground-based observations and obtained deep flux upper limits from four observations, including one with the Spitzer Space Telescope at 4.5 and 8.0 microns. The detections indicate that the source was approximately 1.3--1.6 magnitudes brighter than in 2003--2006, when it was at the tail of a previous similar X-ray flare. Similar related flux variations have been seen in two other AXPs during their X-ray outbursts, suggesting common behavior for large X-ray flux variation events in AXPs. The Spitzer flux 1E 1048.1-5937 limits are sufficiently deep that we can exclude mid-infrared emission similar to that from the AXP 4U 0142+61, which has been interpreted as arising from a dust disk around the AXP. The optical/near-infrared emission from probably has a magnetospheric origin. The similarity in the flux spectra of 4U 0142+61 and 1E 1048.1-5937 challenges the dust disk model proposed for the latter.

We present the results of X-ray and near-IR observations of the anomalous X-ray pulsar 1E 1048.1-5937, believed to be a magnetar. This AXP underwent a period of extreme variability during 2001-2004, but subsequently entered an extended and unexpected quiescence in 2004-2006, during which we monitored it with RXTE, CXO, and HST. Its timing properties were stable for >3 years throughout the quiescent period. 1E 1048.1-5937 again went into outburst in March 2007, which saw a factor of >7 total X-ray flux increase which was anti-correlated with a pulsed fraction decrease, and correlated with spectral hardening, among other effects. The near-IR counterpart also brightened following the 2007 event. We discuss our findings in the context of the magnetar and other models.

(Abridged) We report on new and archival X-ray and near-infrared observations of the anomalous X-ray pulsar 1E 1048.1-5937 performed between 2001-2007 with RXTE, CXO, Swift, HST, and VLT. During its ~2001-2004 active period, 1E 1048.-5937 exhibited two large, long-term X-ray pulsed-flux flares as well as short bursts, and large (>10x) torque changes. Monitoring with RXTE revealed that the source entered a phase of timing stability in 2004; at the same time, a series of four simultaneous observations with CXO and HST in 2006 showed that its X-ray flux and spectrum and near-IR flux, all variable prior to 2005, stabilized. The near-IR flux, when detected by HST (H~22.7 mag) and VLT (K_S~21.0 mag), was considerably fainter than previously measured. Recently, in 2007 March, this newfound quiescence was interrupted by a sudden flux enhancement, X-ray spectral changes and a pulse morphology change, simultaneous with a large spin-up glitch and near-IR enhancement. Our RXTE observations revealed a sudden pulsed flux increase by a factor of ~3 in the 2-10 keV band. In observations with CXO and Swift, we found that the total X-ray flux increased much more than the pulsed flux, reaching a peak value of >7 times the quiescent value (2-10 keV). With these recent data, we find a strong anti-correlation between X-ray flux and pulsed fraction, and a correlation between X-ray spectral hardness and flux. Simultaneously with the radiative and timing changes, we observed a significant X-ray pulse morphology change such that the profile went from nearly sinusoidal to having multiple peaks. We compare these remarkable events with other AXP outbursts and discuss implications in the context of the magnetar model and other models of AXP emission.

Low-mass X-ray binaries, recycled pulsars, cataclysmic variables and magnetically active binaries are observed as X-ray sources in globular clusters. We discuss the classification of these systems, and find that some presumed active binaries are brighter than expected. We discuss a new statistical method to determine from observations how the formation of X-ray sources depends on the number of stellar encounters and/or on the cluster mass. We show that cluster mass is not a proxy for the encounter number, and that optical identifications are essential in proving the presence of primordial binaries among the low-luminosity X-ray sources.

We report on our analysis of a 20 ksec Chandra X-ray observation of the quasi-persistent neutron star soft X-ray transient (SXT) 1M1716-315 in quiescence. Only one source was detected in the HEAO-I error region. Its luminosity is 1.6E32-1.3E33 erg s-1. In this the range is dominated by the uncertainty in the source distance. The source spectrum is well described by an absorbed soft spectrum, e.g. a neutron star atmosphere or black body model. No optical or near-infrared counterpart is present at the location of the X-ray source, down to a magnitude limit of I> 23.5 and K_s> 19.5. The positional evidence, the soft X-ray spectrum together with the optical and near-infrared non-detections provide strong evidence that this source is the quiescent neutron star SXT. The source is 10-100 times too bright in X-rays in order to be explained by stellar coronal X-ray emission. Together with the interstellar extinction measured in outburst and estimates for the source distance, the reported optical and near-infrared limit give an upper limit on the absolute magnitude of the counterpart of I>8.6 and K_s>5.1. This implies that the system is either an ultra-compact X-ray binary having P_orb<1 hr or the companion star is an M-dwarf. We reconstructed the long term X-ray lightcurve of the source. 1M1716-315 has been active for more than 12 years before returning to quiescence, the reported Chandra observation started 16.9+-4.1 years after the outburst ended.

In this Paper we report on phase resolved I-band optical spectroscopic and photometric observations of CirX-1 obtained with the Very Large Telescope. The spectra are dominated by Paschen absorption lines at nearly all orbital phases except near phase zero (coinciding with the X-ray dip) when the absorption lines are filled-in by broad Paschen emission lines. The radial velocity curve of the absorption lines corresponds to an eccentric orbit (e=0.45) whose period and time of periastron passage are consistent with the period and phase predicted by the most recent X-ray dip ephemeris. We found that the I-band magnitude decreases from 17.6 to ~16.8 near phase 0.9-1.0, this brightening coincides in phase with the X-ray dip. Even though it is likely that the absorption line spectrum is associated with the companion star of CirX-1, we cannot exclude the possibility that the spectrum originates in the accretion disc. However, if the spectrum belongs to the companion star, it must be a supergiant of spectral type B5-A0. If we assume that the compact object does not move through the companion star at periastron, the companion star mass is constrained to ~<10 Msun for a 1.4 Msun neutron star, whereas the inclination has to be ~> 13.7 degrees. Alternatively, the measured absorption lines and their radial velocity curve can be associated with the accretion disc surrounding a 1.4 Msun neutron star and its motion around the centre of mass. An absorption line spectrum from an accretion disc is typically found when our line-of-sight passes through the accretion disc rim implying a high inclination. However, from radio observations it was found that the angle between the line-of-sight and the jet axis is smaller than 5 degrees implying that the jet ploughs through the accretion disc in this scenario.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

We present spectroscopic and photometric observations of the optical counterpart to PSR J1911-5958A, a millisecond pulsar located towards the globular cluster NGC 6752. We measure radial velocities from the spectra and determine the systemic radial velocity of the binary and the radial-velocity amplitude of the white-dwarf orbit. Combined with the pulsar orbit obtained from radio timing, we infer a mass ratio of Mpsr/Mwd=7.36+-0.25. The spectrum of the counterpart is that of a hydrogen atmosphere, showing Balmer absorption lines upto H12, and we identify the counterpart as a helium-core white dwarf of spectral type DA5. Comparison of the spectra with hydrogen atmosphere models yield a temperature Teff=10090+-150 K and a surface gravity log g=6.44+-0.20 cm s^-2. Using mass-radius relations appropriate for low-mass helium-core white dwarfs, we infer the white-dwarf mass Mwd=0.18+-0.02 Msun and radius Rwd=0.043+-0.009 Rsun. Combined with the mass ratio, this constrains the pulsar mass to Mpsr=1.40^+0.16_-0.10 Msun. If we instead use the white-dwarf spectrum and the distance of NGC 6752 to determine the white-dwarf radius, we find Rwd=0.058+-0.004 Rsun. For the observed temperature, the mass-radius relations predict a white-dwarf mass of Mwd=0.175+-0.010 Msun, constraining the pulsar mass to Mpsr=1.34+-0.08 Msun. We find that the white-dwarf radius determined from the spectrum and the systemic radial velocity of the binary are only marginally consistent with the values that are expected if PSR J1911-5958A is associated with NGC 6752. We discuss possible causes to explain this inconsistency, but conclude that our observations do not conclusively confirm nor disprove the assocation of the pulsar binary with the globular cluster.

We report on the Chandra X-ray Observatory observations of the globular cluster NGC 288. We detect four X-ray sources within the core radius and seven additional sources within the half-mass radius down to a limiting luminosity of Lx=7e30 erg/s (assuming cluster membership) in the 0.3-7 keV band. We also observed the cluster with the Hubble Space Telescope Advanced Camera for Surveys and identify optical counterparts to seven X-ray sources out of the nine sources within the HST field-of-view. Based on the X-ray and optical properties, we find 2-5 candidates of cataclysmic variables (CVs) or chromospherically active binaries, and 2-5 background galaxies inside the half-mass radius. Since the core density of NGC 288 is very low, the faint X-ray sources of NGC 288 found in the Chandra and HST observations is higher than the prediction on the basis of the collision frequency. We suggest that the CVs and chromospherically active binaries are primordial in origin, in agreement with theoretical expectation.

In this Paper we report on our analysis of a ~25 ksec. Chandra X-ray observation of the neutron star soft X-ray transient (SXT) 1H1905+000 in quiescence. Furthermore, we discuss our findings of the analysis of optical photometric observations which we obtained using the Magellan telescope and photometric and spectroscopic observations which we obtained using the Very Large Telescope at Paranal. The X-ray counterpart of 1H1905+000 was not detected in our Chandra data, with a 95 per cent confidence limit to the source count rate of 1.2x10^-4 counts s^-1. For different spectral models this yields an upper limit on the luminosity of 1.8x10^31 erg s^-1 (for an upper limit on the distance of 10 kpc.) This luminosity limit makes 1H1905+000 the faintest neutron star SXT in quiescence observed to date. The neutron star luminosity is so low that it is similar to the lowest luminosities derived for black hole SXTs in quiescence. This low luminosity for a neutron star SXT challanges the hypothesis presented in the literature that black hole SXTs in quiescence have lower luminosities than neutron star SXTs as a result of the presence of a black hole event horizon. Furthermore, the limit on the neutron star luminosity obtained less than 20 years after the outburst has ceased, constrains the thermal conductivity of the neutron star crust. Finally, the neutron star core must be so cold that unless the time averaged mass accretion rate is lower than 2x10^-12 M_sun yr^-1, core cooling has to proceed via enhanced neutrino emission processes. We derive a limit on the absolute I-band magnitude of the quiescent counterpart of M_I>7.8 assuming the source is at 10 kpc. This is in line with 1H1905+000 being an ultra-compact X-ray binary, as has been proposed based on the low outburst V-band absolute magnitude.

We present optical and near-infrared observations with Keck of the binary millisecond pulsar PSR J0751+1807. We detect a faint, red object - with R=25.08+-0.07, B-R=2.5+-0.3, and R-I=0.90+-0.10 - at the celestial position of the pulsar and argue that it is the white dwarf companion of the pulsar. The colours are the reddest among all known white dwarfs, and indicate a very low temperature, Teff~4000 K. This implies that the white dwarf cannot have the relatively thick hydrogen envelope that is expected on evolutionary grounds. Our observations pose two puzzles. First, while the atmosphere was expected to be pure hydrogen, the colours are inconsistent with this composition. Second, given the low temperature, irradiation by the pulsar should be important, but we see no evidence for it. We discuss possible solutions to these puzzles.

We present the identification of the optical counterparts to the low-mass X-ray binaries 1A 1246-588 and 4U 1812-12. We determine the X-ray position of 1A 1246-588 from ROSAT/PSPC observations and find within the error circle a blue star with V=19.45, B-V=0.22 and R-I=0.22 which we identify as the counterpart. Within the Chandra error circle of 4U 1812-12, a single star is present which appears blue with respect to the stars in the vicinity. It has R=22.15, R-I=1.53. Distance estimates for both systems indicate that the optical counterparts are intrinsically faint, suggesting that they are ultra-compact X-ray binaries. These identifications would increase the number of candidate ultra-compact X-ray binaries from 2 to 4, whereas orbital periods are measured for only 7 systems in the Galactic disk.

Optical observations of the companions of pulsars can help determine the properties of the binaries, as well as those of their components, and give clues to the preceding evolution. In this review, we first describe the different classes of binary pulsars, and present a table with a summary what is known about their optical counterparts. Next, we focus on the class of pulsars that have low-mass, helium-core white dwarf companions. We discuss attempts to determine the masses of both components using optical spectroscopy, and compare the pulsar spin-down ages with cooling ages of the white dwarfs. We confirm that for a given age, the lowest-mass white dwarfs are much hotter than the more massive ones, consistent with recent evolutionary models, although with one glaring exception. We discuss the case of PSR B0820+02, where the cooling age indicates a braking index less than 3, and conclude by describing how cooling ages can be used to test formation scenarios for PSR J1911-5958A, a pulsar binary in the outskirts of NGC 6752.

We present optical astrometry of archival ground and space based imaging of the companion to PSR J1740-5340. The optical position of the companion is significantly offset from the timing position of the pulsar. We briefly investigate the effects of this inconsistency on other timing parameters and compare our position with an improved position of PSR J1740-5340 from recent, preliminary, timing results.

We report on the Chandra X-ray Observatory ACIS-S3 imaging observation of the Galactic globular cluster M4 (NGC 6121). We detect 12 X-ray sources inside the core and 19 more within the cluster half-mass radius. The limiting luminosity of this observation is Lx~10e29 erg/sec for sources associated with the cluster, the deepest X-ray observation of a globular cluster to date. We identify 6 X-ray sources with known objects and use ROSAT observations to show that the brightest X-ray source is variable. Archival data from the Hubble Space Telescope allow us to identify optical counterparts to 16 X-ray sources. Based on the X-ray and optical properties of the identifications and the information from the literature, we classify two (possibly three) sources as cataclysmic variables, one X-ray source as a millisecond pulsar and 12 sources as chromospherically active binaries. Comparison of M4 with 47 Tuc and NGC 6397 suggests a scaling of the number of active binaries in these clusters with the cluster (core) mass.

We report the discovery of two binary millisecond pulsars in the core-collapsed globular cluster M30 using the Green Bank Telescope (GBT) at 20 cm. PSR J2140-2310A (M30A) is an eclipsing 11-ms pulsar in a 4-hr circular orbit and PSR J2140-23B (M30B) is a 13-ms pulsar in an as yet undetermined but most likely highly eccentric (e>0.5) and relativistic orbit. Timing observations of M30A with a 20-month baseline have provided precise determinations of the pulsar's position (within 4" of the optical centroid of the cluster), and spin and orbital parameters, which constrain the mass of the companion star to be m_2 >~ 0.1Msun. The position of M30A is coincident with a possible thermal X-ray point source found in archival Chandra data which is most likely due to emission from hot polar caps on the neutron star. In addition, there is a faint (V_555 ~ 23.8) star visible in archival HST F555W data that may be the companion to the pulsar. Eclipses of the pulsed radio emission from M30A by the ionized wind from the compact companion star show a frequency dependent duration (\propto\nu^{-\alpha} with \alpha ~ 0.4-0.5) and delay the pulse arrival times near eclipse ingress and egress by up to 2-3 ms. Future observations of M30 may allow both the measurement of post-Keplerian orbital parameters from M30B and the detection of new pulsars due to the effects of strong diffractive scintillation.

We report on the identification of the optical counterpart of the binary millisecond pulsar PSR J1911-5958A, located in the outskirts of the globular cluster NGC 6752. At the position of the pulsar we find an object with V=22.08, B-V=0.38, U-B=-0.49. The object is blue with respect to the cluster main sequence by 0.8 magnitudes in B-V. We argue that the object is the white dwarf companion of the pulsar. Comparison with white dwarf cooling models shows that this magnitude and colors are consistent with a low-mass white dwarf at the distance of NGC 6752. If associated with NGC 6752, the white dwarf is relatively young, <2 Gyr, which sets constraints on the formation of the binary and its ejection from the core of the globular cluster.

We report on Keck optical BVRI images and spectroscopy of the companion of the binary millisecond pulsar PSR J0218+4232. A faint bluish (V=24.2, B-V=0.25) counterpart is observed at the pulsar location. Spectra of this counterpart reveal Balmer lines which confirm that the companion is a Helium-core white dwarf. We find that the white dwarf has a temperature of Teff=8060+-150 K. Unfortunately, the spectra are of insufficient quality to put a strong constraint on the surface gravity, although the best fit is for low log g and hence low mass (~0.2 Msun), as expected. We compare predicted white dwarf cooling ages with the characteristic age of the pulsar and find similar values for white dwarf masses of 0.19 to 0.3 Msun. These masses would imply a distance of 2.5 to 4 kpc to the system. The spectroscopic observations also enable us to estimate the mass ratio between the white dwarf and the pulsar. We find q=7.5+-2.4, which is consistent with the current knowledge of white dwarf companions to millisecond pulsars.