We present a statistical study of the glitch population and the behaviour of the glitch activity across the known population of neutron stars. An unbiased glitch database was put together based on systematic searches of radio timing data of 898 rotation-powered pulsars obtained with the Jodrell Bank and Parkes observatories. Glitches identified in similar searches of 5 magnetars were also included. The database contains 384 glitches found in the rotation of 141 of these neutron stars. We confirm that the glitch size distribution is at least bimodal, with one sharp peak at approximately $20\, \rm{\mu\,Hz}$, which we call large glitches, and a broader distribution of smaller glitches. We also explored how the glitch activity $\dot{\nu}_{\rm{g}}$, defined as the mean frequency increment per unit of time due to glitches, correlates with the spin frequency $\nu$, spin-down rate $|\dot{\nu}|$, and various combinations of these, such as energy loss rate, magnetic field, and spin-down age. It is found that the activity is insensitive to the magnetic field and that it correlates strongly with the energy loss rate, though magnetars deviate from the trend defined by the rotation-powered pulsars. However, we find that a constant ratio $\dot\nu_{\rm{g}}/|\dot\nu| = 0.010 \pm 0.001$ is consistent with the behaviour of all rotation-powered pulsars and magnetars. This relation is dominated by large glitches, which occur at a rate directly proportional to $|\dot{\nu}|$. The only exception are the rotation-powered pulsars with the highest values of $|\dot{\nu}|$, such as the Crab pulsar and PSR B0540$-$69, which exhibit a much smaller glitch activity, intrinsically different from each other and from the rest of the population. The activity due to small glitches also shows an increasing trend with $|\dot\nu|$, but this relation is biased by selection effects.

Braking index measurements offer the opportunity to explore the processes affecting the long-term spin evolution of pulsars and possible evolutionary connections between the various pulsar populations. For young pulsars the long-term trends are generally obscured by short term phenomena such as timing noise and the recoveries form large glitches. Here we present a new method to overcome the latter and report on braking index measurements for the Vela-like pulsars PSR B1800-21 and PSR B1823-13, an updated measurement for Vela and new estimates for four more glitching pulsars observed at Jodrell Bank Observatory. The values of braking indices describe the long-term evolution of the pulsars across the P-P' diagram. Despite some measurements being affected by considerable uncertainties, there is evidence for a common trend involving low braking indices (n<=2) among young glitching pulsars. Such values introduce a new variant in the evolution of young pulsars and their relationship with other populations in the P-P' diagram. Low braking indices also imply that these pulsars could be a few times older than their characteristic ages. We consider PSR B1757-24 and conclude that the pulsar could be old enough to be related to the supernova remnant G5.4-1.2. Between glitches, the short-term evolution of Vela-like pulsars is characterised by large inter-glitch braking indices n_ig > 10. We interpret both short and long term trends as signatures of the large glitch activity, and speculate that they are driven by short-term post-glitch re-coupling and a cumulative long-term decoupling of superfluid to the rotation of the star.

Context: GeV gamma-ray pulsations from over 140 pulsars have been characterized using the Fermi Large Area Telescope, enabling improved understanding of the emission regions within the neutron star magnetospheres, and the contributions of pulsars to high energy electrons and diffuse gamma rays in the Milky Way. The first gamma-ray pulsars to be detected were the most intense and/or those with narrow pulses. Aims: As the Fermi mission progresses, progressively fainter objects can be studied. In addition to more distant pulsars (thus probing a larger volume of the Galaxy), or ones in high background regions (thus improving the sampling uniformity across the Galactic plane), we detect pulsars with broader pulses or lower luminosity. Adding pulsars to our catalog with inclination angles that are rare in the observed sample, and/or with lower spindown power, will reduce the bias in the currently known gamma-ray pulsar population. Methods: We use rotation ephemerides derived from radio observations to phase-fold gamma rays recorded by the Fermi Large Area Telescope, to then determine the pulse profile properties. Spectral analysis provides the luminosities and, when the signal-to-noise ratio allows, the cutoff energies. We constrain the pulsar distances by different means in order to minimize the luminosity uncertainties. Results: We present six new gamma-ray pulsars with an eclectic mix of properties. Three are young, and three are recycled. They include the farthest, the lowest power, two of the highest duty-cycle pulsars seen, and only the fourth young gamma-ray pulsar with a radio interpulse. We discuss the biases existing in the current gamma-ray pulsar catalog, and steps to be taken to mitigate the bias.

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

Glitches are sudden spin-up events that punctuate the steady spin down of pulsars and are thought to be due to the presence of a superfluid component within neutron stars. The precise glitch mechanism and its trigger, however, remain unknown. The size of glitches is a key diagnostic for models of the underlying physics. While the largest glitches have long been taken into account by theoretical models, it has always been assumed that the minimum size lay below the detectability limit of the measurements. In this paper we define general glitch detectability limits and use them on 29 years of daily observations of the Crab pulsar, carried out at Jodrell Bank Observatory. We find that all glitches lie well above the detectability limits and by using an automated method to search for small events we are able to uncover the full glitch size distribution, with no biases. Contrary to the prediction of most models, the distribution presents a rapid decrease of the number of glitches below ~0.05 $\mu$Hz. This substantial minimum size indicates that a glitch must involve the motion of at least several billion superfluid vortices and provides an extra observable which can greatly help the identification of the trigger mechanism. Our study also shows that glitches are clearly separated from all the other rotation irregularities. This supports the idea that the origin of glitches is different to that of timing noise, which comprises the unmodelled random fluctuations in the rotation rates of pulsars.

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

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.

The study of high-magnetic-field pulsars is important for examining the relationships between radio pulsars, magnetars, and X-ray-isolated neutron stars (XINSs). Here we report on X-ray observations of three such high-magnetic-field radio pulsars. We first present the results of a deep XMM-Newton observation of PSR J1734-3333, taken to follow up on its initial detection in 2009. The pulsar's spectrum is well fit by a blackbody with a temperature of 300 +/- 60 eV, with bolometric luminosity L_bb = 2.0(+2.2 -0.7)e+32 erg/s = 0.0036E_dot for a distance of 6.1 kpc. We detect no X-ray pulsations from the source, setting a 1 sigma upper limit on the pulsed fraction of 60% in the 0.5-3 keV band. We compare PSR J1734-3333 to other rotation-powered pulsars of similar age and find that it is significantly hotter, supporting the hypothesis that the magnetic field affects the observed thermal properties of pulsars. We also report on XMM-Newton and Chandra observations of PSRs B1845-19 and J1001-5939. We do not detect either pulsar, setting 3 sigma upper limits on their blackbody temperatures of 48 and 56 eV, respectively. Despite the similarities in rotational properties, these sources are significantly cooler than all but one of the XINSs, which we attribute to the two groups having been born with different magnetic fields and hence evolving differently.

Pulsar glitches are traditionally viewed as a manifestation of vortex dynamics associated with a neutron superfluid reservoir confined to the inner crust of the star. In this Letter we show that the non-dissipative entrainment coupling between the neutron superfluid and the nuclear lattice leads to a less mobile crust superfluid, effectively reducing the moment of inertia associated with the angular momentum reservoir. Combining the latest observational data for prolific glitching pulsars with theoretical results for the crust entrainment we find that the required superfluid reservoir exceeds that available in the crust. This challenges our understanding of the glitch phenomenon, and we discuss possible resolutions to the problem.

Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered by young and energetic pulsars. HESS J1857+026 is a spatially extended gamma-ray source detected by H.E.S.S. and classified as a possible pulsar wind nebula candidate powered by PSR J1856+0245. We search for gamma-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1.5 GHz, we phase-fold 36 months of gamma-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete gamma-ray spectral and morphological analysis. No gamma-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident with the TeV source HESS J1857+026. The gamma-ray spectrum is well described by a simple power-law with a spectral index of 1.53 \pm 0.11_{\rm stat} \pm 0.55_{\rm syst}$ and an energy flux of $G(0.1$--100 GeV$)=(2.71 \pm 0.52_{\rm stat} \pm 1.51_{\rm syst}) \times 10^{-11}$ ergs cm$^{-2}$ s$^{-1}$. The $\gamma$-ray luminosity is $L_{PWN}^{\gamma} (0.1$--100 GeV$)=(2.5 \pm 0.5_{stat} \pm 1.5_{syst}) \times 10^{35} (\frac{d}{9 kpc})^2$ ergs s$^{-1}$, assuming a distance of 9 kpc. This implies a $\gamma-$ray efficiency of $\sim$ 5% for $\dot{E}=4.6 \times 10^{36}$ erg $s^{-1}$, in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new constraints on its pulsar wind nebula nature.

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind nebula HESS J1825-137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of sigma = 0.56{\deg} $\pm$ 0.07{\deg} for an assumed Gaussian model. The 1 - 100 GeV LAT spectrum of this source is well described by a power-law with a spectral index of 1.38 $\pm$ 0.12 $\pm$ 0.16 and an integral flux above 1 GeV of (6.50 $\pm$ 0.21 $\pm$ 3.90) x 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV-TeV emission in a cooling-limited nebula powered by the pulsar PSR J1826-1334.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.