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

The gravitational strong equivalence principle (SEP) is a cornerstone of the general theory of relativity (GR). The extreme difference in binding energy between neutron stars and white dwarfs allows for precision tests of the SEP via the technique of pulsar timing. To date, the best limit on the validity of SEP under strong-field conditions was obtained with a unique pulsar in a triple stellar system, PSR J0337+1715. We report here on an improvement of this test using an independent data set acquired over 6 years with the Nan\c cay radio telescope (NRT). The improvements arise from a uniformly sampled data set, a theoretical analysis, and a treatment that fixes some short-comings in the previously published results, leading to better precision and reliability of the test. In contrast to the previously published test, we use a different long-term timing data set, developed a new timing model and an independent numerical integration of the motion of the system, and determined the masses and orbital parameters with a different methodology that treats the parameter $\Delta$, describing a possible strong-field SEP violation, identically to all other parameters. We obtain a violation parameter $\Delta = (+0.5 \pm 1.8) \times 10^{-6}$ at 95\% confidence level, which is compatible with and improves upon the previous study by 30\%. This result is statistics-limited and avoids limitation by systematics as previously encountered. We find evidence for red noise in the pulsar spin frequency, which is responsible for up to 10\% of the reported uncertainty. We use the improved limit on SEP violation to place constraints on a class of well-studied scalar-tensor theories, in particular we find $\omega_{\rm BD} > 140\,000$ for the Brans-Dicke parameter. The conservative limits presented here fully take into account current uncertainties in the equation for state of neutron-star matter.

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

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

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

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.

Binary pulsars are affected by general relativity (GR), causing the spin axes of each pulsar to precess. We present polarimetric radio observations of PSR J1906+0746 that demonstrate the validity of the geometrical model of pulsar polarisation. We reconstruct the (sky-projected) polarisation emission map over the pulsar's magnetic pole and predict the disappearance of the detectable emission by 2028. Two additional tests of GR are performed in this system, including the spin-precession for strongly self-gravitating bodies. We constrain the relativistic treatment of the pulsar polarisation model and measure the pulsar beaming fraction, with implications for the population of neutron stars and the expected rate of neutron star mergers.

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

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

We report on the timing observations of the millisecond pulsar PSR J2055+3829 originally discovered as part of the SPAN512 survey conducted with the Nan\c{c}ay Radio Telescope. The pulsar has a rotational period of 2.089 ms, and is in a tight 3.1 hr orbit around a very low mass ($0.023 \leq m_c \lesssim 0.053$ M$_\odot$, 90\% c.l.) companion. Our 1.4 GHz observations reveal the presence of eclipses of the pulsar's radio signal caused by the outflow of material from the companion, for a few minutes around superior conjunction of the pulsar. The very low companion mass, the observation of radio eclipses, and the detection of time variations of the orbital period establish PSR J2055+3829 as a `black widow' (BW) pulsar. Inspection of the radio signal from the pulsar during ingress and egress phases shows that the eclipses in PSR J2055+3829 are asymmetric and variable, as is commonly observed in other similar systems. More generally, the orbital properties of the new pulsar are found to be very similar to those of other known eclipsing BW pulsars. No gamma-ray source is detected at the location of the pulsar in recent \textit{Fermi}-LAT source catalogs. We used the timing ephemeris to search ten years of \textit{Fermi} Large Area Telescope (LAT) data for gamma-ray pulsations, but were unable to detect any, possibly because of the pulsar's large distance. We finally compared the mass functions of eclipsing and non-eclipsing BW pulsars and confirmed previous findings that eclipsing BWs have higher mass functions than their non-eclipsing counterparts. Larger inclinations could explain the higher mass functions of eclipsing BWs. On the other hand, the mass function distributions of Galactic disk and globular cluster BWs appear to be consistent, suggesting, despite the very different environments, the existence of common mechanisms taking place in the last stages of evolution of BWs.

Recycled pulsars are old ($\gtrsim10^{8}$ yr) neutron stars that are descendants from close, interacting stellar systems. In order to understand their evolution and population, we must find and study the largest number possible of recycled pulsars in a way that is as unbiased as possible. In this work, we present the discovery and timing solutions of five recycled pulsars in binary systems (PSRs J0509$+$0856, J0709$+$0458, J0732$+$2314, J0824$+$0028, J2204$+$2700) and one isolated millisecond pulsar (PSR J0154$+$1833). These were found in data from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327). All these pulsars have a low dispersion measure (DM) ($\lesssim 45 \, \rm{pc}\, cm^{-3}$), and have a DM-determined distance of $\lesssim$ 3 kpc. Their timing solutions, have data spans ranging from 1 to $\sim$ 7 years, include precise estimates of their spin and astrometric parameters, and for the binaries, precise estimates of their Keplerian binary parameters. Their orbital periods range from about 4 to 815 days and the minimum companion masses (assuming a pulsar mass of 1.4 $\rm{M_{\odot}}$) range from $\sim$ 0.06--1.11 $\rm{M_{\odot}}$. For two of the binaries we detect post-Keplerian parameters; in the case of PSR~J0709$+$0458 we measure the component masses but with a low precision, in the not too distant future the measurement of the rate of advance of periastron and the Shapiro delay will allow very precise mass measurements for this system. Like several other systems found in the AO327 data, PSRs J0509$+$0854, J0709$+$0458 and J0732$+$2314 are now part of the NANOGrav timing array for gravitational wave detection.

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

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

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.

PSR~J1618$-$3921 is an $11.99$-ms pulsar in a $22.7$-d orbit around a likely low-mass He white dwarf companion, discovered in a survey of the intermediate Galactic latitudes at 1400 MHz conducted with the Parkes radio telescope in the late 1990s. Although PSR~J1618$-$3921 was discovered more than 15 years ago, only limited information has been published about this pulsar which has a surprisingly large orbital eccentricity ($e \simeq 0.027$), considering its high spin frequency and likely small companion mass. The focus of this work is a precise measurement of the spin, astrometric and orbital characteristics of PSR J1618$-$3921. This was done with timing observations made at the Nan\c{c}ay Radio Telescope, from 2009 to 2017. We analyzed the timing data recorded at the Nan\c{c}ay Radio Telescope over several years to characterize the properties of PSR~J1618$-$3921. A rotation ephemeris for this pulsar was obtained by analyzing the times of arrival of the radio pulses at the telescope. We confirm the unusual eccentricity of PSR~J1618$-$3921 and discuss several hypotheses regarding its formation, in the context of other discoveries of recycled pulsars in eccentric orbits.

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

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.

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

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.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

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

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

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

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

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.

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

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

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

Millisecond pulsars (MSPs) represent nearly half of the more than 160 currently known $\gamma$-ray pulsars detected by the Large Area Telescope on the \textit{Fermi} satellite, and a third of all known MSPs are seen in $\gamma$ rays. The least energetic $\gamma$-ray MSPs enable us to probe the so-called deathline for high-energy emission, i.e., the spin-down luminosity limit under which pulsars (PSRs) cease to produce detectable high-energy radiation. Characterizing the MSP luminosity distribution helps to determine their contribution to the Galactic diffuse $\gamma$-ray emission. We made use of the high-quality pulsar timing data recorded at the Nan\c{c}ay Radio Telescope over several years to characterize the properties of a selection of MSPs. For one of the pulsars, the dataset was complemented with Westerbork Synthesis Radio Telescope observations. The rotation ephemerides derived from this analysis were also used to search the LAT data for new $\gamma$-ray MSPs. For the MSPs considered in this study, we obtained new transverse proper motion measurements or updated the existing ones, and placed new distance constraints for some of them, with four new timing parallax measurements. We discovered significant GeV $\gamma$-ray signals from four MSPs, i.e., PSRs J0740+6620, J0931$-$1902, J1455$-$3330, and J1730$-$2304. The latter is now the least energetic $\gamma$-ray pulsar found to date. Despite the improved $\dot E$ and $L_\gamma$ estimates, the relationship between these two quantities remains unclear, especially at low $\dot E$ values.

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.

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

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

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

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

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.

PSR J1906+0746 is a young pulsar in the relativistic binary with the second-shortest known orbital period, of 3.98 hours. We here present a timing study based on five years of observations, conducted with the 5 largest radio telescopes in the world, aimed at determining the companion nature. Through the measurement of three post-Keplerian orbital parameters we find the pulsar mass to be 1.291(11) M_sol, and the companion mass 1.322(11) M_sol respectively. These masses fit well in the observed collection of double neutron stars, but are also compatible with other white dwarfs around young pulsars such as J1906+0746. Neither radio pulsations nor dispersion-inducing outflows that could have further established the companion nature were detected. We derive an HI-absorption distance, which indicates that an optical confirmation of a white dwarf companion is very challenging. The pulsar is fading fast due to geodetic precession, limiting future timing improvements. We conclude that young pulsar J1906+0746 is likely part of a double neutron star, or is otherwise orbited by an older white dwarf, in an exotic system formed through two stages of mass transfer.

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

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.

In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys using image pattern recognition with deep neural nets---the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interferences by looking for patterns from candidate. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of up to thousands pixel of image data. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ~9000 neurons. Different from other pulsar selection programs which use pre-designed patterns, the PICS AI teaches itself the salient features of different pulsars from a set of human-labeled candidates through machine learning. The deep neural networks in this AI system grant it superior ability in recognizing various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates different from the training set. In this completely independent test, PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars, to the top 961 (1%) of 90008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

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.

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

We have made extensive observations of 35 distant slow (non-recycled) pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Timing observations of these pulsars over several years at Arecibo Observatory and Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation properties. Despite being a relatively distant population, these pulsars have properties that mirror those of the previously known pulsar population. Many of the sources exhibit timing noise, and one underwent a small glitch. We have used multifrequency data to measure the interstellar scattering properties of these pulsars. We find scattering to be higher than predicted along some lines of sight, particularly in the Cygnus region. Lastly, we present XMM-Newton and Chandra observations of the youngest and most energetic of the pulsars, J1856+0245, which has previously been associated with the GeV-TeV pulsar wind nebula HESS J1857+026.

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

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

Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.

Shortly after the discovery of PSR J1906+0746, some hints of profile variations were already interpreted as first signs of relativistic spin-precession occuring. Using observations from the Nan\c{c}ay, Arecibo and Green Bank Radio Observatories, we report here the measurement of pulse profile and polarimetric variations. Using the Rotating Vector Model, we show that PSR J1906+0746 is likely to be an orthogonal rotator ($\alpha \simeq 80^\circ$). Fitting our polarimetric data to a precession model, we determined the geometry of the pulsar and found a wide misalignment angle ($\delta = 89_{-44}^{+85}$ deg, 95% C.L.), although the uncertainty is large. Assuming this geometry, we constructed the beam maps of both magnetic poles.

We present an ongoing survey with the Nan\c{c}ay Radio Telescope at L-band. The targeted area is $74^\circ \lesssim l <150^\circ$ and $3.5^\circ < |b| < 5^\circ$. This survey is characterized by a long integration time (18 min), large bandwidth (512 MHz) and high time and frequency resolution (64 $\mu$s and 0.5 MHz) giving a nominal sensitivity limit of 0.055 mJy for long period pulsars. This is about 2 times better than the mid-latitude HTRU survey, and is designed to be complementary with current large scale surveys. This survey will be more sensitive to transients (RRATs, intermittent pulsars), distant and faint millisecond pulsars as well as scintillating sources (or any other kind of radio faint sources) than all previous short-integration surveys.

The on-going PALFA survey is searching the Galactic plane (|b| < 5 deg., 32 < l < 77 deg. and 168 < l < 214 deg.) for radio pulsars at 1.4 GHz using ALFA, the 7-beam receiver installed at the Arecibo Observatory. By the end of August 2012, the PALFA survey has discovered 100 pulsars, including 17 millisecond pulsars (P < 30 ms). Many of these discoveries are among the pulsars with the largest DM/P ratios, proving that the PALFA survey is capable of probing the Galactic plane for millisecond pulsars to a much greater depth than any previous survey. This is due to the survey's high sensitivity, relatively high observing frequency, and its high time and frequency resolution. Recently the rate of discoveries has increased, due to a new more sensitive spectrometer, two updated complementary search pipelines, the development of online collaborative tools, and access to new computing resources. Looking forward, focus has shifted to the application of artificial intelligence systems to identify pulsar-like candidates, and the development of an improved full-resolution pipeline incorporating more sophisticated radio interference rejection. The new pipeline will be used in a complete second analysis of data already taken, and will be applied to future survey observations. An overview of recent developments, and highlights of exciting discoveries will be presented.

We present two millisecond pulsar discoveries from the PALFA survey of the Galactic plane with the Arecibo telescope. PSR J1955+2527 is an isolated pulsar with a period of 4.87 ms, and PSR J1949+3106 has a period of 13.14 ms and is in a 1.9-day binary system with a massive companion. Their timing solutions, based on 4 years of timing measurements with the Arecibo, Green Bank, Nan\c{c}ay and Jodrell Bank telescopes, allow precise determination of spin and astrometric parameters, including precise determinations of their proper motions. For PSR J1949+3106, we can clearly detect the Shapiro delay. From this we measure the pulsar mass to be 1.47(+0.43/-0.31) solar masses, the companion mass to be 0.85(+0.14/-0.11) solar masses and the orbital inclination to be i = 79.9(+1.6/-1.9) degrees, where uncertainties correspond to +/- 1-\sigma\ confidence levels. With continued timing, we expect to also be able to detect the advance of periastron for the J1949+3106 system. This effect, combined with the Shapiro delay, will eventually provide very precise mass measurements for this system and a test of general relativity.

We present the discovery and phase-coherent timing of four highly dispersed millisecond pulsars (MSPs) from the Arecibo PALFA Galactic plane survey: PSRs J1844+0115, J1850+0124, J1900+0308, and J1944+2236. Three of the four pulsars are in binary systems with low-mass companions, which are most likely white dwarfs, and which have orbital periods on the order of days. The fourth pulsar is isolated. All four pulsars have large dispersion measures (DM > 100 pc cm-3), are distant (> 3.4 kpc), faint at 1.4 GHz (< 0.2 mJy), and are fully recycled (with spin periods P between 3.5 and 4.9 ms). The three binaries also have very small orbital eccentricities, as expected for tidally circularized, fully recycled systems with low-mass companions. These four pulsars have DM/P ratios that are among the highest values for field MSPs in the Galaxy. These discoveries bring the total number of confirmed MSPs from the PALFA survey to fifteen. The discovery of these MSPs illustrates the power of PALFA for finding weak, distant MSPs at low-Galactic latitudes. This is important for accurate estimates of the Galactic MSP population and for the number of MSPs that the Square Kilometer Array can be expected to detect.

We report the discovery of the millisecond pulsar PSR J2043+1711 in a search of a Fermi Large Area Telescope (LAT) source with no known associations, with the Nancay Radio Telescope. The new pulsar, confirmed with the Green Bank Telescope, has a spin period of 2.38 ms, is relatively nearby (d <~ 2 kpc), and is in a 1.48 day orbit around a low mass companion, probably a He-type white dwarf. Pulsed gamma-ray emission was detected in the data recorded by the Fermi LAT. The gamma-ray light curve and spectral properties are typical of other gamma-ray millisecond pulsars seen with Fermi. X-ray observations of the pulsar with Suzaku and the Swift/XRT yielded no detection. At 1.4 GHz we observe strong flux density variations because of interstellar diffractive scintillation, however a sharp peak can be observed at this frequency during bright scintillation states. At 327 MHz the pulsar is detected with a much higher signal-to-noise ratio and its flux density is far more steady. However, at that frequency the Arecibo instrumentation cannot yet fully resolve the pulse profile. Despite that, our pulse time-of-arrival measurements have a post-fit residual rms of 2 \mus. This and the expected stability of this system has made PSR J2043+1711 one of the first new Fermi-selected millisecond pulsars to be added to pulsar gravitational wave timing arrays. It has also allowed a significant measurement of relativistic delays in the times of arrival of the pulses due to the curvature of space-time near the companion, but not yet with enough precision to derive useful masses for the pulsar and the companion. A mass for the pulsar between 1.7 and 2.0 solar masses can be derived if a standard millisecond pulsar formation model is assumed. In this article we also present a comprehensive summary of pulsar searches in Fermi LAT sources with the Nancay Radio Telescope to date.

A new spectrometer for the Green Bank Telescope (GBT) is being built jointly by the NRAO and the CASPER, University of California, Berkeley. The spectrometer uses 8 bit ADCs and will be capable of processing up to 1.25 GHz bandwidth from 8 dual polarized beams. This mode will be used to process data from focal plane arrays. The spectrometer supports observing mode with 8 tunable digital sub-bands within the 1.25 GHz bandwidth. The spectrometer can also be configured to process a bandwidth of up to 10 GHz with 64 tunable sub-bands from a dual polarized beam. The vastly enhanced backend capabilities will support several new science projects with the GBT.

Compared to the previous INPOP versions, the INPOP10a planetary and lunar ephemeris has several improvements. For the planets of our solar system, no big change was brought in the dynamics but improvements were implemented in the fitting process, the data sets used in the fit and in the selection of fitted parameters. We report here the main characteristics of the planetary part of INPOP10a like the fit of the product of the Solar mass with the gravitational constant (GM$_{\odot}$) instead of the astronomical unit. Determinations of PPN parameters as well as adjustments of the Sun J2 and of asteroid masses are also presented. New advances of nodes and perihelia of planets were also estimated and are given here. As for INPOP08, INPOP10a provides to the user, positions and velocities of the planets, the moon, the rotation angles of the Earth and the Moon as well as TT-TDB chebychev polynomials at http://www.imcce.fr/inpop

We report the discovery of the 20.7 ms binary pulsar J1952+2630, made using the distributed computing project Einstein@Home in Pulsar ALFA survey observations with the Arecibo telescope. Follow-up observations with the Arecibo telescope confirm the binary nature of the system. We obtain a circular orbital solution with an orbital period of 9.4 hr, a projected orbital radius of 2.8 lt-s, and a mass function of f = 0.15 solar masses by analysis of spin period measurements. No evidence of orbital eccentricity is apparent; we set a 2-sigma upper limit e < 1.7e-3. The orbital parameters suggest a massive white dwarf companion with a minimum mass of 0.95 solar masses, assuming a pulsar mass of 1.4 solar masses. Most likely, this pulsar belongs to the rare class of intermediate mass binary pulsars. Future timing observations will aim to determine the parameters of this system further, measure relativistic effects, and elucidate the nature of the companion star.

Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $\alpha$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $\alpha=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.

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

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

Pulsars PSR J0248+6021 (rotation period P=217 ms and spin-down power Edot = 2.13E35 erg/s) and PSR J2240+5832 (P=140 ms, Edot = 2.12E35 erg/s) were discovered in 1997 with the Nancay radio telescope during a northern Galactic plane survey, using the Navy-Berkeley Pulsar Processor (NBPP) filter bank. GeV gamma-ray pulsations from both were discovered using the Fermi Large Area Telescope. Twelve years of radio and polarization data allow detailed investigations. The two pulsars resemble each other both in radio and in gamma-ray data. Both are rare in having a single gamma-ray pulse offset far from the radio peak. The high dispersion measure for PSR J0248+6021 (DM = 370 pc cm^-3) is most likely due to its being within the dense, giant HII region W5 in the Perseus arm at a distance of 2 kpc, not beyond the edge of the Galaxy as obtained from models of average electron distributions. Its high transverse velocity and the low magnetic field along the line-of-sight favor this small distance. Neither gamma-ray, X-ray, nor optical data yield evidence for a pulsar wind nebula surrounding PSR J0248+6021. The gamma-ray luminosity for PSR J0248+6021 is L_ gamma = (1.4 \pm 0.3)\times 10^34 erg/s. For PSR J2240+5832, we find either L_gamma = (7.9 \pm 5.2) \times 10^34 erg/s if the pulsar is in the Outer arm, or L_gamma = (2.2 \pm 1.7) \times 10^34 erg/s for the Perseus arm. These luminosities are consistent with an L_gamma ~ sqrt(Edot) rule. Comparison of the gamma-ray pulse profiles with model predictions, including the constraints obtained from radio polarization data, favor emission in the far magnetosphere. These two pulsars differ mainly in their inclination angles and acceleration gap widths, which in turn explains the observed differences in the gamma-ray peak widths.

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.

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

We report on the high precision timing analysis of the pulsar-white dwarf binary PSR J1012+5307. Using 15 years of multi-telescope data from the European Pulsar Timing Array (EPTA) network, a significant measurement of the variation of the orbital period is obtained. Using this ideal strong-field gravity laboratory we derive theory independent limits for both the dipole radiation and the variation of the gravitational constant.

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

Nancay radiotelescope is involved in high precision timing since 20 years. Since 2004, a coherent dedispersion instrumentation enables numerous routine observations on more than 200 pulsars using half of the time if this 100-meters class radiotelescope. Two main programs are currently conducted. A large set of young and old pulsars is timed for a multi-wavelength approach, complementary to the very successful high energy observations of pulsars done by FERMI. A set of highly stable millisecond pulsars is monitored as our contribution to the European Pulsar Timing Array in order to probe the cosmological Gravitational Wave Background.

We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nan\c{c}ay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of $\pi = 1.2(3)$ mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of $e<8.4\times 10^{-7}$ is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of $\dot{P}_{b} = 5(1)\times 10^{-14}$ is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant $\dot{G}$.