• 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.
  • In 8 years of operation, the Large Area Telescope (LAT) on the Fermi satellite has impacted our understanding of gamma-ray pulsars dramatically. The LAT now sees over two hundred pulsars: the largest class of GeV sources in the Milky Way. They are diverse -- radio loud versus quiet, young versus millisecond, in evolving binary systems versus isolated, and so on. Relatively few of the GeV pulsars have also been seen in soft gamma rays. After an overview, we present 10 new radio pulsars, six young and four recycled, for which we detect gamma-ray pulsations.
  • 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 Square Kilometre Array (SKA) is an integral part of the next-generation observatories that will survey the Universe across the electromagnetic spectrum, and beyond, revolutionizing our view of fundamental physics, astrophysics and cosmology. Owing to their extreme nature and clock-like properties, pulsars discovered and monitored by SKA will enable a broad range of scientific endeavour and play a key role in this quest. This chapter summarizes the pulsar-related science goals that will be reached with coordinated efforts among SKA and other next-generation astronomical facilities.
  • We report on our second-year campaign of X-ray follow-up observations of unidentified Fermi-LAT \gamma-ray sources at high Galactic latitudes (|b|>10 degree) using the X-ray Imaging Spectrometer onboard the Suzaku X-ray Observatory. In this second year of the project, seven new targets were selected from the First Fermi-LAT Catalog, and studied with 20-40 ks effective Suzaku exposures. We detected an X-ray point source coincident with the position of the recently discovered millisecond pulsar PSR J2302+4442 within the 95% confidence error circle of 1FGL J2302.8+4443. The X-ray spectrum of the detected counterpart was well fit by a blackbody model with temperature of kT ~0.3 keV, consistent with an origin of the observed X-ray photons from the surface of a rotating magnetized neutron star. For four other targets which were also recently identified with a normal pulsar (1FGL J0106.7+4853) and millisecond pulsars (1FGL J1312.6+0048, J1902.0-5110, and J2043.2+1709), only upper limits in the 0.5-10 keV band were obtained at the flux levels of ~10^{-14} erg cm^{-2} s^{-1}. A weak X-ray source was found in the field of 1FGL J1739.4+8717, but its association with the variable \gamma-ray emitter could not be confirmed with the available Suzaku data alone. For the remaining Fermi-LAT object 1FGL J1743.8-7620 no X-ray source was detected within the LAT 95% error ellipse. We briefly discuss the general properties of the observed high Galactic-latitude Fermi-LAT objects by comparing their multiwavelength properties with those of known blazars and millisecond pulsars.
  • Radio timing observations of millisecond pulsars (MSPs) in support of Fermi LAT observations of the gamma-ray sky enhance the sensitivity of high-energy pulsation searches. With contemporaneous ephemerides we have detected gamma-ray pulsations from PSR B1937+21, the first MSP ever discovered, and B1957+20, the first known black-widow system. The two MSPs share a number of properties: they are energetic and distant compared to other gamma-ray MSPs, and both of them exhibit aligned radio and gamma-ray emission peaks, indicating co-located emission regions in the outer magnetosphere of the pulsars. However, radio observations are also crucial for revealing MSPs in Fermi unassociated sources. In a search for radio pulsations at the position of such unassociated sources, the Nan\c{c}ay Radio Telescope discovered two MSPs, PSRs J2017+0603 and J2302+4442, increasing the sample of known Galactic disk MSPs. Subsequent radio timing observations led to the detection of gamma-ray pulsations from these two MSPs as well. We describe multiwavelength timing and spectral analysis of these four pulsars, and the modeling of their gamma-ray light curves in the context of theoretical models.
  • The Fermi observatory was launched on June 11, 2008. It hosts the \emph{Large Area Telescope} (LAT), sensitive to $\gamma$-ray photons from 20 MeV to over 300 GeV. When the LAT began its activity, nine young and energetic pulsars were known in $\gamma$ rays. At least several tens of pulsar detections by the LAT were predicted before launch. The LAT also allowed the study of millisecond pulsars (MSPs), never firmly detected in $\gamma$ rays before Fermi. This thesis first presents the pulsar timing campaign for the LAT, in collaboration with large radiotelescopes and X-ray telescopes, allowing for high sensitivity pulsed searches. Furthermore, it lead to quasi-homogeneous coverage of the galactic MSPs, so that the search for pulsations in LAT data for this population of stars was not affected by an \emph{a priori} bias. We present a search for pulsations from these objects in LAT data. For the first time, eight galactic MSPs have been detected as sources of pulsed $\gamma$-ray emission over 100 MeV. In addition, a couple of good candidates for future detection are seen. A similar search for globular cluster MSPs has not succeeded so far. Comparison of the phase-aligned $\gamma$-ray and radio light curves, as well as the spectral shapes, leads to the conclusion that their $\gamma$-ray emission is similar to that of normal pulsars, and is probably produced in the outer-magnetosphere. This discovery suggests that many unresolved $\gamma$-ray sources are unknown MSPs.