To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

Isolated Neutron Stars (INSs) were the first sources identified in the field of high-energy gamma-ray astronomy. At first, in the 70s, there were only two identified sources, the Crab and Vela pulsars. However, although few in number, these objects were crucial in establishing the very concept of a gamma-ray source. Moreover, they opened up significant discovery space both in the theoretical and phenomenological fronts. The need to explain the copious gamma-ray emission of these pulsars led to breakthrough developments in understanding the structure and physics of neutron star magnetospheres. In parallel, the 20-year-long chase to understand the nature of Geminga unveiled the existence of a radio-quiet, gamma-ray-emitting, INS, adding a new dimension to the INS family. Today we are living through an extraordinary time of discovery. The current generation of gamma-ray detectors has vastly increased the population of known of gamma-ray-emitting neutron stars. The 100 mark was crossed in 2011 and we are now approaching 150. The gamma-ray-emitting neutron star population exhibits roughly equal numbers of radio-loud and radio-quiet young INSs, plus an astonishing, and unexpected, group of isolated and binary millisecond pulsars (MSPs). The number of MSPs is growing so rapidly that they are on their way to becoming the most numerous members of the family of gamma-ray-emitting Neutron Stars (NSs) . Even as these findings have set the stage for a revolution in our understanding of gamma-ray-emitting neutron stars, long term monitoring of the gamma-ray sky has revealed evidence of flux variability in the Crab Nebula as well as in the pulsed emission from PSR J2021+4026, challenging a four-decade-old, constant-emission paradigm. Now we know that both pulsars and their nebulae can, indeed, display variable emission.

Using archival as well as freshly acquired data, we assess the X-ray behaviour of the Fermi/LAT gamma-ray pulsars listed in the First Fermi source catalog. After revisiting the relationships between the pulsars' rotational energy losses and their X and gamma-ray luminosities, we focus on the distance-indipendent gamma to X-ray flux ratios. When plotting our Fgamma/Fx values as a function of the pulsars' rotational energy losses, one immediately sees that pulsars with similar energetics have Fgamma/Fx spanning 3 decades. Such spread, most probably stemming from vastly different geometrical configurations of the X and gamma-ray emitting regions, defies any straightforward interpretation of the plot. Indeed, while energetic pulsars do have low Fgamma/Fx values, little can be said for the bulk of the Fermi neutron stars. Dividing our pulsar sample into radio-loud and radio-quiet subsamples, we find that, on average, radio-quiet pulsars do have higher values of Fgamma/Fx, implying an intrinsec faintness of their X-ray emission and/or a different geometrical configuration. Moreover, despite the large spread mentioned above, statistical tests show a lower scatter in the radio-quiet dataset with respect to the radio-loud one, pointing to a somewhat more constrained geometry for the radio-quiet objects with respect to the radio-loud ones.

2009 has been an extraordinary year for gamma-ray pulsar astronomy and 2010 promises to be equally good. Not only have we registered an extraordinary increase in the number of pulsars detected in gamma rays, but we have also witnessed the birth of new sub-families: first of all, the radio-quiet gamma pulsars and later an ever growing number of millisecond pulsars, a real surprise. We started with a sample of 7 gamma-ray emitting neutron stars (6 radio pulsars and Geminga) and now the Fermi-LAT harvest encompasses 24 "Geminga-like" new gamma-ray pulsars, a dozen millisecond pulsars and about thirty radio pulsars. Moreover, radio searches targeted to LAT unidentified sources yielded 18 new radio millisecond pulsars, several of which have been already detected also in gamma rays. Thus, currently the family of gamma-ray emitting neutron stars seems to be evenly divided between classical radio pulsars, millisecond pulsars and radio quiet neutron stars.

When it comes to identifying or to characterizing gamma-ray sources, X-ray observations are of paramount importance. Correlated X-and-gamma-ray flux variations are a powerful identification tool, if the gamma-ray source is unidentified, or an important diagnostic tool to understand the behaviour of an identified source. Moreover, X-ray observations of non-variable unidentified gamma-ray sources, both galactic and extragalactic, can unveil interesting candidate counterparts, narrowing down the search space and improving significantly the chances for a successful identification. Swift observations of Fermi gamma-ray-selected pulsar error boxes provide accurate positions of likely counterparts. This makes it possible both to confirm pulsations and to improve timing solution, opening up a new synergy between X and gamma ray observations.

Combining resolved spectroscopy with deep imaging, XMM-Newton is providing new insights on the particle acceleration processes long known to be at work in the magnetospheres of isolated neutron stars. According to a standard theoretical interpretation, in neutron stars' magnetospheres particles are accelerated along the B field lines and, depending on their charge, they can either move outward, to propagate in space, or be funnelled back, towards the star surface. While particles impinging on the neutron star surface should heat it at well defined spots, outgoing ones could radiate extended features in the neutron star surroundings. By detecting hot spots, seen to come in and out of sight as the star rotates, as well as extended features trailing neutron stars as they move in the interstellar medium, XMM-Newton provides the first end-to-end test to the particle acceleration process.