• We apply the algorithm published by Liang et al. (2014) to describe the Double Pulsar system J0737$-$3039 and extract the sense of rotation of first born recycled pulsar PSR J0737$-$3039A. We find that this pulsar is rotating prograde in its orbit. This is the first direct measurement of the sense of rotation of a pulsar with respect to its orbit and a direct confirmation of the rotating lighthouse model for pulsars. This result confirms that the spin angular momentum vector is closely aligned with the orbital angular momentum, suggesting that kick of the supernova producing the second born pulsar J0737$-$3039B was small.
  • We present a search for optical bursts from the repeating fast radio burst FRB 121102 using simultaneous observations with the high-speed optical camera ULTRASPEC on the 2.4-m Thai National Telescope and radio observations with the 100-m Effelsberg Radio Telescope. A total of 13 radio bursts were detected, but we found no evidence for corresponding optical bursts in our 70.7-ms frames. The 5-sigma upper limit to the optical flux density during our observations is 0.33 mJy at 767nm. This gives an upper limit for the optical burst fluence of 0.046 Jy ms, which constrains the broadband spectral index of the burst emission to alpha < -0.2. Two of the radio pulses are separated by just 34 ms, which may represent an upper limit on a possible underlying periodicity (a rotation period typical of pulsars), or these pulses may have come from a single emission window that is a small fraction of a possible period.
  • We present a long-term spectral monitoring of the unique Double Pulsar binary PSR J0737-3039 corresponding to two "Large Programs" performed by XMM-Newton in 2006 and 2011. Spectral variability of pulsar emission in soft X-rays is not evident over 5 years, despite the significant relativistic spin precession in the considered time span (about 25deg). We provide, for the first time, evidence of hard X-ray emission from the system in the 5-8 keV energy band. The standard spectral analysis was coupled to the energy dependent spatial analysis to confirm this excess, most likely ascribed to iron line emission. The Fe K-alpha emission line at 6.4-6.97 keV was previously unheard-of in non-accreting binary systems and could testify to the presence of a relic disk that survived the supernova explosions that terminated the lives of the Double Pulsar's stellar progenitors. The existence of a relic disk in this system reinforces speculation about the presence of similar structures around other peculiar classes of isolated neutron stars.
  • Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, errors in the planetary ephemeris, the solar wind and instrumentation issues. All these effects can induce significant false detections of gravitational waves. We test mitigation routines to account for clock errors, ephemeris errors and the solar wind. We demonstrate that it is non-trivial to find an effective mitigation routine for the planetary ephemeris and emphasise that other spatially correlated signals may be present in the data.
  • 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 an analysis of X-ray, Ultraviolet and optical/near-IR photometric data of the transitional millisecond pulsar binary XSSJ12270-4859, obtained at different epochs after the transition to a rotation-powered radio pulsar state. The observations, while confirming the large-amplitude orbital modulation found in previous studies after the state change, also reveal an energy dependence of the amplitudes as well as variations on time scale of months. The amplitude variations are anti-correlated in the X-ray and the UV/optical bands. The average X-ray spectrum is described by a power law with \Gamma index of 1.07(8) without requiring an additional thermal component. The power law index \Gamma varies from 1.2 to 1.0 between superior and inferior conjunction of the neutron star. We interpret the observed X-ray behaviour in terms of synchrotron radiation emitted in an extended intrabinary shock, located between the pulsar and the donor star, which is eclipsed due to the companion orbital motion. The G5 type donor dominates the UV/optical and near-IR emission and is similarly found to be heated up to ? 6500K as in the disc state. The analysis of optical light curves gives a binary inclination 46 < i < 65deg and a mass ratio 0.11 < q <0.26. The donor mass is found to be 0.15 < M2 < 0.36Msun for a neutron star mass of 1.4Msun. The variations in the amplitude of the orbital modulation are interpreted in terms of small changes in the mass flow rate from the donor star. The spectral energy distribution from radio to gamma-rays is composed by multiple contributions that are different from those observed during the accretion-powered state.
  • 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.
  • The Square Kilometre Array (SKA) will use pulsars to enable precise measurements of strong gravity effects in pulsar systems, which yield tests of gravitational theories that cannot be carried out anywhere else. The Galactic census of pulsars will discover dozens of relativistic pulsar systems, possibly including pulsar -- black hole binaries which can be used to test the "cosmic censorship conjecture" and the "no-hair theorem". Also, the SKA's remarkable sensitivity will vastly improve the timing precision of millisecond pulsars, allowing probes of potential deviations from general relativity (GR). Aspects of gravitation to be explored include tests of strong equivalence principles, gravitational dipole radiation, extra field components of gravitation, gravitomagnetism, and spacetime symmetries.
  • Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged]
  • With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.
  • The pulsar PSR J1756$-$2251 resides in a relativistic double neutron star (DNS) binary system with a 7.67-hr orbit. We have conducted long-term precision timing on more than 9 years of data acquired from five telescopes, measuring five post-Keplerian parameters. This has led to several independent tests of general relativity (GR), the most constraining of which shows agreement with the prediction of GR at the 4% level. Our measurement of the orbital decay rate disagrees with that predicted by GR, likely due to systematic observational biases. We have derived the pulsar distance from parallax and orbital decay measurements to be 0.73$_{-0.24}^{+0.60}$ kpc (68%) and < 1.2 kpc (95% upper limit), respectively; these are significantly discrepant from the distance estimated using Galactic electron density models. We have found the pulsar mass to be 1.341$\pm$0.007 M$_\odot$, and a low neutron star (NS) companion mass of 1.230$\pm$0.007 M$_\odot$. We also determined an upper limit to the spin-orbit misalignment angle of 34{\deg} (95%) based on a system geometry fit to long-term profile width measurements. These and other observed properties have led us to hypothesize an evolution involving a low mass loss, symmetric supernova progenitor to the second-formed NS companion, as is thought to be the case for the double pulsar system PSR J0737$-$3039A/B. This would make PSR J1756$-$2251 the second compact binary system providing concrete evidence for this type of NS formation channel.
  • We present X-ray observations of the "redback" eclipsing radio millisecond pulsar and candidate radio pulsar/X-ray binary transition object PSR J1723-2837. The X-ray emission from the system is predominantly non-thermal and exhibits pronounced variability as a function of orbital phase, with a factor of ~2 reduction in brightness around superior conjunction. Such temporal behavior appears to be a defining characteristic of this variety of peculiar millisecond pulsar binaries and is likely caused by a partial geometric occultation by the main-sequence-like companion of a shock within the binary. There is no indication of diffuse X-ray emission from a bow shock or pulsar wind nebula associated with the pulsar. We also report on a search for point source emission and $\gamma$-ray pulsations in Fermi Large Area Telescope data using a likelihood analysis and photon probability weighting. Although PSR J1723-2837 is consistent with being a $\gamma$-ray point source, due to the strong Galactic diffuse emission at its position a definitive association cannot be established. No statistically significant pulsations or modulation at the orbital period are detected. For a presumed source detection, the implied $\gamma$-ray luminosity is $\lesssim$5% of its spin-down power. This indicates that PSR J1723-2837 is either one of the least efficient $\gamma$-ray producing millisecond pulsars or, if the detection is spurious, the $\gamma$-ray emission pattern is not directed towards us.
  • Intermediate mass black holes (IMBHs) are among the most elusive objects in contemporary astrophysics. Both theoretical and observational evidence of their existence is subject of debate. Conversely, both theory and observations confirm the presence of a large population of millisecond pulsars (MSPs) with low mass companions residing in globular cluster (GC) centers. If IMBHs are common in GC centers as well, then dynamical interactions will inevitably break up many of these binaries, causing the ejection of several fast MSPs in the Galactic halo. Such population of fast halo MSPs, hard to produce with 'standard' MSP generation mechanisms, would provide a strong, albeit indirect, evidence of the presence of a substantial population of IMBHs in GCs. In this paper we study in detail the dynamical formation and evolution of such fast MSPs population, highlighting the relevant observational properties and assessing detection prospects with forthcoming radio surveys.
  • We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.
  • We present the evolution of the radio emission from the 2.8-s pulsar of the double pulsar system PSR J0737-3039A/B. We provide an update on the Burgay et al. (2005) analysis by describing the changes in the pulse profile and flux density over five years of observations, culminating in the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases by 0.177 mJy/yr at the brightest orbital phases and the pulse profile evolves from a single to a double peak, with a separation rate of 2.6 deg/yr. The pulse profile changes are most likely caused by relativistic spin precession, but can not be easily explained with a circular hollow-cone beam as in the model of Clifton & Weisberg (2008). Relativistic spin precession, coupled with an elliptical beam, can model the pulse profile evolution well. This particular beam shape predicts geometrical parameters for the two bright orbital phases which are consistent and similar to those derived by Breton et al. (2008). However, the observed decrease in flux over time and B's eventual disappearance cannot be easily explained by the model and may be due to the changing influence of A on B.
  • As part of a survey for radio pulsars with the Parkes 64-m telescope we have discovered PSR J1622-4950, a pulsar with a 4.3-s rotation period. Follow-up observations show that the pulsar has the highest inferred surface magnetic field of the known radio pulsars (B ~ 3e14 G), exhibits significant timing noise and appears to have an inverted spectrum. Unlike the vast majority of the known pulsar population, PSR J1622-4950 appears to switch off for many hundreds of days and even in its on-state exhibits extreme variability in its flux density. Furthermore, the integrated pulse profile changes shape with epoch. All of these properties are remarkably similar to the only two magnetars previously known to emit radio pulsations. The position of PSR J1622-4950 is coincident with an X-ray source that, unlike the other radio pulsating magnetars, was found to be in quiescence. We conclude that our newly discovered pulsar is a magnetar - the first to be discovered via its radio emission.
  • In 2003 a previously unpulsed Einstein and ROSAT source cataloged as soft and dim (Lx of few 10^33 ergs) thermal emitting object, namely XTE J1810-197, was identified as the first unambiguous transient Anomalous X-ray Pulsar. Two years later this source was also found to be a bright highly polarized transient radio pulsar, a unique property among both AXPs and radio pulsars. In September 2006 the Swift Burst Alert Telescope (BAT) detected an intense burst from the candidate AXP CXOU J164710.2-455216, which entered in an outburst state reaching a peak emission of at least a factor of 300 higher than quiescence. Here, we briefly outline the recent results concerning the outburst phenomena observed in these two AXPs. In particular, XTE J1810-197 has probed to be a unique laboratory to monitor the timing and spectral properties of a cooling/fading AXP, while new important information have been inferred from X-ray and radio band simultaneous observations. CXOU J164710.2-455216 has been monitored in X-rays and radio bands since the very beginning of its outbursting state allowing us to cover the first phases of the outburst and to study the timing and spectral behavior during the first months.
  • We explore the origin of the Galactic high latitude black hole X-ray binary XTE J1118+480, and in particular its birth location and the magnitude of the kick received by the black hole upon formation in the supernova explosion. We constrain the age of the companion to the black hole using stellar evolution calculations between 2 Gyr and 5 Gyr, making an origin in a globular cluster unlikely. We therefore argue that the system was born in the Galactic disk and the supernova propelled it in its current high latitude orbit. Given the current estimates on its distance, proper motion and radial velocity, we back-trace the orbit of XTE J1118+480 in the Galactic potential to infer the peculiar velocity of the system at different disk crossings over the last 5 Gyr. Taking into account the uncertainties on the velocity components, we infer an average peculiar velocity of 183 \pm 31 km/s. The maximum velocity which the binary can acquire by symmetric supernova mass loss is about 100 km/s, which is 2.7 sigma away from the mean of the peculiar velocity distribution. We therefore argue that an additional asymmetric kick velocity is required. By considering the orientation of the system relative to the plane of the sky, we derive a 95% probability for a non null component of the kick perpendicular to the orbital plane of the binary. The distribution of perpendicular velocities is skewed to lower velocities with an average of 93^{+55}_{-60} km/s.
  • We have used deep high-resolution multiband images taken at the ESO Very Large Telescope to identify the optical binary companion to the millisecond pulsar (PSR J1911-5958A) located in the halo of the Galactic Globular Cluster NGC6752. The object turns out to be a blue star whose position in the Color Magnitude Diagram is consistent with the cooling sequence of a low mass (M~0.17-0.20Mo), low metallicity Helium white dwarf (He-WD) at the cluster distance. This is the second He-WD which has been found to orbit a millisecond pulsar in GGCs. Curiously both objects have been found to lie on the same mass He-WD cooling sequence. The anomalous position of PSR J1911-5958A with respect to the globular cluster center (~6') suggested that this system has recently (<1 Gyr) been ejected from the cluster core as the result of a strong dynamical interaction. The data presented here allows to constrain the cooling age of the companion within a fairly narrow range (~1.2-2.8 Gyr), therefore suggesting that such dynamical encounter must have acted on an already recycled millisecond pulsar.
  • The five millisecond pulsars that inhabit NGC 6752 display locations or accelerations remarkably different with respect to all other pulsars known in globular clusters. This may reflect the occurrence of an uncommon dynamics in the cluster core that could be attributed to the presence of a massive perturber. We here investigate whether a single intermediate-mass black hole, lying on the extrapolation of the mass versus sigma relation observed in galaxy spheroids, or, a less massive black hole binary could play the requested role.
  • NGC6752 hosts in its halo PSR J1911-5958A, a newly discovered binary millisecond pulsar which is the most distant pulsar ever known from the core of a globular cluster. Interestingly, its recycling history seems in conflict with a scenario of ejection resulting from ordinary stellar dynamical encounters. A scattering event off a binary system of two black holes with masses in the range of 3-50 solar masses that propelled PSR J1911-5958A into its current peripheral orbit seems more likely. It is still an observational challenge to unveil the imprint(s) left from such a dark massive binary on cluster's stars: PSR J1911-5958A may be the first case.
  • Was the black hole in XTE J1118+480 ejected from a globular cluster or kicked away from the galactic disk?
  • We have investigated different scenarios for the origin of the binary millisecond pulsar PSR J1911-5958A in NGC 6752, the most distant pulsar discovered from the core of a globular cluster to date. The hypothesis that it results from a truly primordial binary born in the halo calls for accretion-induced collapse and negligible recoil speed at the moment of neutron star formation. Scattering or exchange interactions off cluster stars are not consistent with both the observed orbital period and its offset position. We show that a binary system of two black holes with (unequal) masses in the range of 3-100 solar masses can live in NGC 6752 until present time and can have propelled PSR J1911-5958A into an eccentric peripheral orbit during the last ~1 Gyr.
  • The empirical relation between the X-ray luminosity (in the 2-10 keV band) and the rate of spin-down energy loss of a sample of 37 pulsars is re-examined considering recent data from ASCA, RXTE, BeppoSAX, Chandra, and XMM-Newton. The data show a significant scatter around an average correlation computed including statistical and systematic errors. All the X-ray luminosities are found to lie below a critical line corresponding to a maximum efficiency of conversion of rotational energy into X-rays varying between 0.1 and 80%. The large dispersion of the X-ray luminosity below the critical line indicates that other physical parameters uncorrelated with the period and period derivative need to be included to account for the observed emission at X-ray energies. We indicate a few possibilities that all conspire to lower the luminosity.