We consider the largest observed sample including all intermediate-frequency peaked (IBL) and high-frequency peaked (HBL) flaring blazars above 100 GeV up to redshift $z = 0.6$. We show that the best-fit regression line of the emitted spectral indices $\Gamma_{\rm em} (z)$ is a concave parabola decreasing as $z$ increases, thereby implying a statistical correlation between the $\{\Gamma_{\rm em} (z) \}$ distribution and $z$. This result contradicts our expectation that such a distribution should be $z$-independent. We argue that the above correlation does not arise from any selection bias. We show that our expectation naturally emerges provided that axion-like particles (ALPs) are put into the game. Moreover, ALPs can also explain why flat spectrum radio quasars emit up to 400 GeV, in sharp contradiction with conventional physics. So, the combination of the two very different but consistent results -- taken at face value -- leads to a hint at an ALP with mass $m = {\cal O} (10^{-10} \, {\rm eV})$ and two-photon coupling in the range $2.94 \times 10^{- 12} \, {\rm GeV}^{- 1} < g_{a \gamma \gamma} < 0.66 \times 10^{- 10} \, {\rm GeV}^{- 1}$. As a bonus, the Universe would become considerably more transparent above energies $E \gtrsim 1 \, {\rm TeV}$ than dictated by conventional physics. Our prediction can be checked not only by the new generation of observatories like CTA, HAWC, GAMMA-400, LHAASO, TAIGA-HiSCORE and HERD, but also thanks to the planned laboratory experiments ALPS II (upgraded), STAX, IAXO and with other techniques now being developed by Avignone and collaborators.

Soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are slowly rotating, isolated neutron stars that sporadically undergo episodes of long-term flux enhancement (outbursts) generally accompanied by the emission of short bursts of hard X-rays. This behaviour can be understood in the magnetar model, according to which these sources are mainly powered by their own magnetic energy. This is supported by the fact that the magnetic fields inferred from several observed properties of AXPs and SGRs are greater than - or at the high end of the range of - those of radio pulsars. In the peculiar case of SGR 0418+5729, a weak dipole magnetic moment is derived from its timing parameters, whereas a strong field has been proposed to reside in the stellar interior and in multipole components on the surface. Here we show that the X-ray spectrum of SGR 0418+5729 has an absorption line, the properties of which depend strongly on the star's rotational phase. This line is interpreted as a proton cyclotron feature and its energy implies a magnetic field ranging from 2E14 gauss to more than 1E15 gauss.

The middle-aged PSR J0357+3205 is a nearby, radio-quiet, bright gamma-ray pulsar discovered by the Fermi mission. Our previous Chandra observation revealed a huge, very peculiar structure of diffuse X-ray emission, originating at the pulsar position and extending for > 9' on the plane of the sky. To better understand the nature of such a nebula, we have studied the proper motion of the parent pulsar. We performed relative astrometry on Chandra images of the field spanning a time baseline of 2.2 yr, unveiling a significant angular displacement of the pulsar counterpart, corresponding to a proper motion of 0.165"+/-0.030" yr^(-1). At a distance of ~500 pc, the space velocity of the pulsar would be of ~390 km s^(-1) assuming no inclination with respect to the plane of the sky. The direction of the pulsar proper motion is perfectly aligned with the main axis of the X-ray nebula, pointing to a physical, yet elusive link between the nebula and the pulsar space velocity. No optical emission in the H_alpha line is seen in a deep image collected at the Gemini telescope, which implies that the interstellar medium into which the pulsar is moving is fully ionized.

RX J0822-4300 is the Central Compact Object associated with the Puppis A supernova remnant. Previous X-ray observations suggested RX J0822-4300 to be a young neutron star with a weak dipole field and a peculiar surface temperature distribution dominated by two antipodal spots with different temperatures and sizes. An emission line at 0.8 keV was also detected. We performed a very deep (130 ks) observation with XMM-Newton, which allowed us to study in detail the phase-resolved properties of RX J0822-4300. Our new data confirm the existence of a narrow spectral feature, best modelled as an emission line, only seen in the `Soft' phase interval - when the cooler region is best aligned to the line of sight. Surprisingly, comparison of our recent observations to the older ones yields evidence for a variation in the emission line component, which can be modelled as a decrease in the central energy from ~0.80 keV in 2001 to ~0.73 keV in 2009--2010. The line could be generated via cyclotron scattering of thermal photons in an optically thin layer of gas, or - alternatively - it could originate in low-rate accretion by a debris disk. In any case, a variation in energy, pointing to a variation of the magnetic field in the line emitting region, cannot be easily accounted for.

We present the AGILE gamma-ray observations in the energy range 50 MeV - 10 GeV of the supernova remnant (SNR) W44, one of the most interesting systems for studying cosmic-ray production. W44 is an intermediate-age SNR (20, 000 years) and its ejecta expand in a dense medium as shown by a prominent radio shell, nearby molecular clouds, and bright [SII] emitting regions. We extend our gamma-ray analysis to energies substantially lower than previous measurements which could not conclusively establish the nature of the radiation. We find that gamma-ray emission matches remarkably well both the position and shape of the inner SNR shocked plasma. Furthermore, the gamma-ray spectrum shows a prominent peak near 1 GeV with a clear decrement at energies below a few hundreds of MeV as expected from neutral pion decay. Here we demonstrate that: (1) hadron-dominated models are consistent with all W44 multiwavelength constraints derived from radio, optical, X-ray, and gamma-ray observations; (2) ad hoc lepton-dominated models fail to explain simultaneously the well-constrained gamma-ray and radio spectra, and require a circumstellar density much larger than the value derived from observations; (3) the hadron energy spectrum is well described by a power-law (with index s = 3.0 \pm 0.1) and a low-energy cut-off at Ec = 6 \pm 1 GeV. Direct evidence for pion emission is then established in an SNR for the first time.

The Large Area Telescope (LAT) onboard the Fermi satellite opened a new era for pulsar astronomy, detecting gamma-ray pulsations from more than 60 pulsars, ~40% of which are not seen at radio wavelengths. One of the most interesting sources discovered by LAT is PSR J0357+3205, a radio-quiet, middle-aged (tau_C ~0.5 Myr) pulsar standing out for its very low spin-down luminosity (Erot ~6x10^33 erg/s), indeed the lowest among non-recycled gamma-ray pulsars. A deep X-ray observation with Chandra (0.5-10 keV), coupled with sensitive optical/infrared ground-based images of the field, allowed us to identify PSR J0357+3205 as a faint source with a soft spectrum, consistent with a purely non-thermal emission (photon index Gamma=2.53+/-0.25). The absorbing column (NH=8+/-4x10^20 cm^-2) is consistent with a distance of a few hundred parsecs. Moreover, the Chandra data unveiled a huge (9 arcmin long) extended feature apparently protruding from the pulsar. Its non-thermal X-ray spectrum points to synchrotron emission from energetic particles from the pulsar wind, possibly similar to other elongated X-ray tails associated with rotation-powered pulsars and explained as bow-shock pulsar wind nebulae (PWNe). However, energetic arguments, as well as the peculiar morphology of the diffuse feature associated with PSR J0357+3205 make the bow-shock PWN interpretation rather challenging.

The well known Crab Nebula is at the center of the SN1054 supernova remnant. It consists of a rotationally-powered pulsar interacting with a surrounding nebula through a relativistic particle wind. The emissions originating from the pulsar and nebula have been considered to be essentially stable. Here we report the detection of strong gamma-ray (100 MeV-10 GeV) flares observed by the AGILE satellite in September, 2010 and October, 2007. In both cases, the unpulsed flux increased by a factor of 3 compared to the non-flaring flux. The flare luminosity and short timescale favor an origin near the pulsar, and we discuss Chandra Observatory X-ray and HST optical follow-up observations of the nebula. Our observations challenge standard models of nebular emission and require power-law acceleration by shock-driven plasma wave turbulence within a ~1-day timescale.

Prompted by the Fermi LAT discovery of a radio-quiet gamma-ray pulsar inside the CTA 1 supernova remnant, we obtained a 130 ks XMM-Newton observation to assess the timing behavior of this pulsar. Exploiting both the unprecedented photon harvest and the contemporary Fermi LAT timing measurements, a 4.7 sigma single peak pulsation is detected, making PSR J0007+7303 the second example, after Geminga, of a radio-quiet gamma-ray pulsar also seen to pulsate in X-rays. Phase-resolved spectroscopy shows that the off-pulse portion of the light curve is dominated by a power-law, non-thermal spectrum, while the X-ray peak emission appears to be mainly of thermal origin, probably from a polar cap heated by magnetospheric return currents, pointing to a hot spot varying throughout the pulsar rotation.

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

1E 1207.4-5209, the peculiar Central Compact object in the G296.5+10.0 supernova remnant, has been proposed to be an "anti-magnetar" - a young neutron star born with a weak dipole field. Accretion, possibly of supernova fallback material, has also been invoked to explain a large surface temperature anisotropy as well as the generation of peculiar cyclotron absorption features superimposed to its thermal spectrum. Interestingly enough, a faint optical/infrared source was proposed as a possible counterpart to 1E 1207.4-5209, but later questioned, based on coarse positional coincidence. Considering the large offset of 1E 1207.4-5209 with respect to the center of its host supernova remnant, the source should move at ~70 mas/yr. Thus, we tested the association by measuring the proper motion of the proposed optical counterpart. Using HST observations spanning 3.75 years, we computed a 3 sigma upper limit of 7 mas/yr. Absolute astrometry on the same HST data set also places the optical source significantly off the 99% confidence Chandra position. This allows us to safely rule out the association. Using the HST data set, coupled to ground-based observations collected at the ESO/VLT, we set the deepest limits ever obtained to the optical/infrared emission from 1E 1207.4-5209. By combining such limits to the constraints derived from X-ray timing, we rule out accretion as the source of the thermal anisotropy of the neutron star.

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

The population of stellar black holes (SBHs) in the Galaxy and galaxies generally is poorly known in both number and distribution. SBHs are the fossil record of the massive stars in galaxy evolution and may have produced some (if not all) of the intermediate mass (\gsim100\Msun) black holes (IMBHs) and, in turn, the central supermassive black holes (SMBHs) in galactic nuclei. For the first time, a Galaxy-wide census of accreting black holes, and their more readily recognizable tracer population, accreting neutron stars (NSs), could be measured with a wide-field hard X-ray imaging survey and soft X-ray and optical/IR prompt followup -- as proposed for the EXIST mission.

Using gamma-ray data collected by the Astrorivelatore Gamma ad Immagini LEggero (AGILE) satellite over a period of almost one year (from 2007 July to 2008 June), we searched for pulsed signals from 35 potentially interesting radio pulsars, ordered according to $F_{\gamma}\propto \sqrt{\dot{E}} d^{-2}$ and for which contemporary or recent radio data were available. AGILE detected three new top-ranking nearby and Vela-like pulsars with good confidence both through timing and spatial analysis. Among the newcomers we find pulsars with very high rotational energy losses, such as the remarkable PSR B1509-58 with a magnetic field in excess of 10^13 Gauss, and PSR J2229+6114 providing a reliable identification for the previously unidentified EGRET source 3EG 2227+6122. Moreover, the powerful millisecond pulsar B1821-24, in the globular cluster M28, is detected during a fraction of the observations. Four other promising gamma-ray pulsar candidates, among which is the notable J2043+2740 with an age in excess of 1 million years, show a possible detection in the timing analysis only and deserve confirmation.

Over the next decade, we can expect time domain astronomy to flourish at optical and radio wavelengths. In parallel with these efforts, a dedicated transient "machine" operating at higher energies (X-ray band through soft gamma-rays) is required to reveal the unique subset of events with variable emission predominantly visible above 100 eV. Here we focus on the transient phase space never yet sampled due to the lack of a sensitive, wide-field and triggering facility dedicated exclusively to catching high energy transients and enabling rapid coordinated multi-wavelength follow-up. We first describe the advancements in our understanding of known X-ray transients that can only be enabled through such a facility and then focus on the classes of transients theoretically predicted to be out of reach of current detection capabilities. Finally there is the exciting opportunity of revealing new classes of X-ray transients and unveiling their nature through coordinated follow-up observations at longer wavelengths.

In this letter we calculate the energy loss of highly magnetized neutron star due to friction with quantum vacuum, namely Quantum Vacuum Friction (QVF). Taking into account one-loop corrections in the effective Heisenberg-Euler Lagrangian of the light-light interaction, we derive an analytic expression for QVF allowing us to consider magnetic field at the surface of the star as high as $10^{11} $T. In the case of magnetars with high magnetic field above the QED critical field, we show that the energy loss by QVF dominates the energy loss process. This has important consequences, in particular on the inferred value of the magnetic field. This also indicates the need for independent measurements of magnetic field, energy loss rate, and of the braking index to fully characterize magnetars.

1E 161348-5055 (1E 1613) is a point-like, soft X-ray source originally identified as a radio-quiet, isolated neutron star, shining at the center of the 2000 yr old supernova remnant RCW103. 1E 1613 features a puzzling 6.67 hour periodicity as well as a dramatic variability over a time scale of few years. Such a temporal behavior, coupled to the young age and to the lack of an obvious optical counterpart, makes 1E 1613 a unique source among all compact objects associated to SNRs. It could either be the first low-mass X-ray binary system discovered inside a SNR, or a peculiar isolated magnetar with an extremely slow spin period. Analysis of archival IR observations, performed in 2001 with the VLT/ISAAC instrument, and in 2002 with the NICMOS camera onboard HST unveils a very crowded field. A few sources are positionally consistent with the refined X-ray error region that we derived from the analysis of 13 Chandra observations. To shed light on the nature of 1E 1613, we have performed deep IR observations of the field with the NACO instrument at the ESO/VLT, searching for variability. We find no compelling reasons to associate any of the candidates to 1E 1613. On one side, within the frame of the binary system model for the X-ray source, it is very unlikely that one of the candidates be a low-mass companion star to 1E 1613. On the other side, if the X-ray source is an isolated magnetar surrounded by a fallback disc, we cannot exclude that the IR counterpart be hidden among the candidates. If none of the potential counterparts is linked to the X-ray source, 1E 1613 would remain undetected in the IR down to Ks>22.1. Such an upper limit is consistent only with an extremely low-mass star (an M6-M8 dwarf) at the position of 1E 1613, and makes rather problematic the interpretation of 1E 1613 as an accreting binary system.

While waiting for new gamma-ray burst detections, the Burst Alert Telescope (BAT) on board Swift covers each day ~50% of the sky in the hard X-ray band (``Survey data''). The large field of view (FOV), high sensitivity and good angular resolution make BAT a potentially powerful all-sky hard X-ray monitor, provided that mask--related systematics can be properly accounted for. We have developed and tested a complete procedure entirely based on public Swift/BAT software tools to analyse BAT Survey data, aimed at assessing the flux and spectral variability of bright sources in the 15-150 keV energy range. Detailed tests of the capabilities of our procedure were performed focusing, in particular, on the reliability of spectral measurements over the entire BAT FOV. First, we analyzed a large set of Crab observations, spread over ~7 months. Next, we studied the case of GRO J1655-40, a strongly variable source, which experienced a 9-month long outburst, beginning on February 2005. Such an outburst was systematically monitored with the well-calibrated PCA and HEXTE instruments on board the RXTE mission. Thanks to the good BAT temporal coverage of the source, we have been able to cross-check BAT light-curves with simultaneous HEXTE ones. The Crab tests have shown that our procedure recovers both the flux and the source spectral shape over the whole FOV of the BAT instrument. Moreover, by cross-checking GRO J1655-40 light-curves obtained by BAT and HEXTE, we found the spectral and flux evolution of the outburst to be in very good agreement. Using our procedure, BAT reproduces HEXTE fluxes within a 10-15% uncertainty with a 3sigma sensitivity of ~20 mCrab for an on-axis source, thus establishing its capability to monitor the evolution of relatively bright hard X-rays sources.

PSR B0540-69 is the Crab twin in the Large Magellanic Cloud. Age, energetic and overall behaviour of the two pulsars are very similar. The same is true for the general appearance of their pulsar wind nebulae (PWNe). Analysis of Hubble Space Telescope images spanning 10 years unveiled significant variability in the PWN surrounding PSR B0540-69, with a hot spot moving at ~0.04c. Such behaviour, reminiscent of the variability observed in the Crab nebula along the counter-jet direction, may suggest an alternative scenario for the geometry of the system. The same data were used to assess the pulsar proper motion. The null displacement recorded over 10 y allowed us to set a 3sigma upper limit of 290 km/s to the pulsar velocity.

The limited angular resolution of gamma-ray telescopes prevents a direct identification of the majority of sources detected so far. This is particularly true for the low latitude, probably galactic, ones only 10 % of which has been identified. Most counterparts of the identified low-latitude gamma-ray sources are Isolated Neutron Stars (INS), both radio-loud and radio-quiet (Geminga-like) objects, which are characterised by an extremely high value of the X-ray-to-optical flux ratio f_X/f_opt. Therefore, the systematic X-ray and optical coverage of low-latitude unidentified gamma-ray sources aiming at high f_X/f_opt sources seems one of the most promising ways to spot INS candidate counterparts. Since low latitude sources are heavily affected by the interstellar absorption at both X-ray and optical wavelengths, we have focussed on two middle-latitude, probably galactic, GRO/EGRET sources: 3EG J0616-3310 and 3EG J1249-8330. These two sources, which could belong to a local galactic population, have been selected owing to their relatively good positional accuracy, spectral shape and lack of candidate extragalactic radio counterparts. Here we report on X-ray observations of the two gamma-ray error boxes performed with XMM-Newton and on their optical follow-up carried on with the Wide Field Imager at the ESO/MPG 2.2m telescope. Less than half of the ~300 sources detected by the X-ray coverage have no optical counterparts. Among those, we have selected few interesting sources with f_X/f_opt > 100, which we consider promising INS candidates.

Observations with the Newton X-ray Multimirror Mission (XMM) show a strong periodic modulation at 6.67+/-0.03 hours of the X-ray source at the centre of the 2,000-year-old supernova remnant RCW 103. No fast pulsations are visible. If genetically tied to the supernova remnant, the source could either be an X-ray binary, comprising a compact object and a low-mass star in an eccentric orbit, or an isolated neutron star. In the latter case, its age-period combination would point to a peculiar magnetar, dramatically slowed-down, possibly by a supernova debris disc. Both scenarios require non-standard assumptions on the formation and evolution of compact objects in supernova explosions.

The strange timing property of X-ray pulsar 1E 1207.4-5209 can be explained by the hypothesis that it is a member of an ultra-compact binary system. This paper confronts the ultra-compact assumption with the observed properties of this pulsar. The gravitational potential well of an ultra-compact binary can enlarge the corotation radius and thus make it possible for accreting material to reach the surface of the NS in the low accretion rate case. Thus the generation of the absorption features should be similar to the case of accreting pulsars. The close equality of the energy loss by fast cooling of the postsupernova neutron star and the energy dissipation needed for a wide binary evolving to an ultra-compact binary demonstrates that the ultra-compact binary may be formed in 10-100yr after the second supernova explosion. Moreover, the ultra-compact binary hypothesis can well explain the the absence of optical counterpart and the observed two black body emissions. We suggest a simple method which can test the binary nature directly with XMM-Newton and Chandra observations. We further predict that the temperature of the two black bodies should vary at different pulse periods.

A deep (100 ks) XMM-Newton observation of Geminga has shown two faint tails of diffuse X-ray emission, extending for ~2' behind the pulsar, well aligned with the proper motion (PM) direction. We report here on a recent ~20 ks Chandra observation, which unveils a new structure, ~25'' long and ~5'' thick, starting at the pulsar position and perfectly aligned with the PM direction, with a surface brightness ~40 times higher than that of the XMM Tails. The Chandra comet-like feature has a remarkably hard spectrum (photon index 0.9-1.4) and a luminosity of ~5.5x10^28 erg s^-1, comparable to the energetics of the larger XMM one. Geminga is thus the first neutron star to show a clear X-ray evidence of a large-scale, outer bow-shock as well as a short, inner cometary trail.

The radio quiet neutron star 1E1207.4-5209 has been the target of a 260 ks XMM-Newton observation, which yielded, as a by product, an harvest of about 200 serendipitous X-ray sources above a limiting flux of 2E-15 erg/cm2/s, in the 0.3-8 keV energy range. In view of the intermediate latitude of our field (b~10 deg), it comes as no surprise that the logN-logS distribution of our serendipitous sources is different from those measured either in the Galactic Plane or at high galactic latitudes. Here we shall concentrate on the analysis of the brightest sources in our sample, which unveiled a previously unknown Seyfert-2 galaxy.

We present numerical calculations of the photon-light-pseudoscalar-boson conversion in the recently discovered binary pulsar system J0737-3039. Light pseudoscalar bosons (LPBs) oscillate into photons in the presence of strong magnetic fields. In the context of this binary pulsar system, this phenomenon attenuates the light beam emitted by one of the pulsars, when the light ray goes through the magnetosphere of the companion pulsar. We show that such an effect is observable in the gamma-ray band since the binary pulsar is seen almost edge-on, depending on the value of the LPB mass and on the strenght of its two-photon coupling. Our results are surprising in that they show a very sharp and significant (up to 50%) transition probability in the gamma-ray ($>$ tens of MeV) domain. The observations can be performed by the upcoming NASA GLAST mission.

In this letter we study the propagation of light in the neighbourhood of magnetised neutron stars. Thanks to the optical properties of quantum vacuum in the presence of a magnetic field, light emitted by background astronomical objects is deviated giving rise to a phenomenon of the same kind as the gravitational one. We give a quantitative estimation of this effect and we discuss the possibility of its observation. We show that this effect could be detected monitoring the evolution of the recently discovered double neutron star system J0737-3039.

XMM-Newton observations of PSR B0656+14, PSR B1055-52 and Geminga have substantially increased the statistics available for these three isolated neutron stars, so apparently similar to deserve the nickname of "Three Musketeers" (Becker & Truemper, 1997). Here we shall take advantage of the EPIC statistics to perform phase resolved spectroscopy for all three objects. The phase-averaged spectrum of the three musketeers is best described by a three component model. This includes two blackbody components, a cooler one, possibly originating from the bulk of the star surface, and a hotter one, coming from a smaller portion of the star surface (a "hot spot"), plus a power law. The relative contributions of the three components are seen to vary as a function of phase, as the stars' rotation bring into view different emitting regions. The hot spots, which have very different apparent dimensions (in spite of the similarity of the three neutron stars polar cap radii) are responsible for the bulk of the phase variation. The amplitude of the observed phase modulation is also markedly different for the three sources. Another striking aspect of our phase-resolved phenomenology is the apparent lack of any common phase alignement between the observed modulation patterns for the two blackbody components. They are seen to vary in phase in the case of PSR B1055-52, but in anti-phase in the case of PSR B0656+14. These findings do not support standard and simplistic models of neutron star magnetic field configuration and surface temperature distribution.

Most counterparts of the identified low-latitude gamma-ray sources are Isolated Neutron Stars (INSs). They are characterized by an extremely high value of the X/optical flux ratio, therefore a systematic X-ray/optical coverage of the fields of unidentified low-latitude gamma-ray sources is the best way to unveil possible INS counterparts. Since the low-latitude sources are heavily affected by the interstellar absorption in both the X-ray and optical bands, we decided to apply the above strategy to two middle-latitude EGRET sources, which could belong to a local galactic population: 3EG J0616-3310 and 3EG J1249-8330. Here we report on the global X-ray characterisation of about 300 objects, on their candidate optical counterparts and on the preliminary results of their identification.

Isolated neutron stars are seen in x-rays through their nonthermal and/or surface thermal emissions. XMM-Newton observations of the Geminga pulsar show a 43 electron volt spectrum from the whole neutron star surface, as well as a power-law component above 2 kiloelectron volts. In addition, we have detected a hot (170 electron volts) thermal emission from an 60-meter-radius spot on the pulsar s surface. Such a thermal emission, only visible at selected phase intervals, may be coming from polar hot spot(s), long thought to exist as a result of heating from magnetospheric accelerated particles. It may provide the missing link between the x-ray and gamma-ray emission of the pulsar.

We present the results of a 50 ks long X-ray observation of the relativistic double pulsar system PSR J0737-3039 obtained with the XMM-Newton satellite in March 2004. The source has a soft spectrum (power law photon index = 3.5(+0.5,-0.3)) and a 0.2-10 keV luminosity of ~3x10^30 erg/s (assuming a distance of 500 pc), consistent with the values derived from a previous Chandra observation. No flux variations have been detected, implying the absence of large orbital modulations. The high time resolution of the EPIC instrument has allowed us to perform the first search for X-ray pulsations from this system. The result was negative, with an upper limit of 60% on the pulsed fraction of the 22 ms pulsar.

The limited angular resolution of gamma-ray telescopes prevents the straight identification of the majority of the sources detected so far. This is particularly true for the low latitude, probably galactic ones, only 10 % of which has been identified. Most of the counterparts are Isolated Neutron Stars, both radio-loud and radio-quiet, characterised by an extremely high value of X-ray-to-optical flux ratio. The best way to search for INSs in the error boxes of unidentified EGRET sources is to perform the X-ray coverage of the gamma-ray field, followed by the optical characterization of each X-ray source. We applied this procedure to two EGRET sources, which could belong to a local galactic population: 3EG J0616-3310 and 3EG J1249-8330. Here we report on the analysis of about 300 X-ray objects, as well as on their optical study.

We report on the deep (100 ksec) XMM-Newton/EPIC observation of the field of the Geminga pulsar. The unprecedented throughput of the instrument allowed to detect two elongated parallel X-ray tails trailing the neutron star. They are aligned with the object's supersonic motion, extend for ~2', and have a nonthermal spectrum produced by electron-synchrotron emission in the bow shock between the pulsar wind and the surrounding medium. Electron lifetime against synchrotron cooling matches the source transit time over the X-ray features' length. Such an X-ray detection of a pulsar bow shock allows us to gauge the pulsar electron injection energy and the shock magnetic field while constraining the angle of Geminga's motion and the local matter density. We give also preliminary results on the timing and spectral analysis of the \~63,000 photons collected from the neutron star.

The discovery of deep spectral features in the X-ray spectrum of 1E1207.4-5209 focussed the attention of the astronomical community on this radio-quiet NS, making it the most intensively observed INS ever. The harvest of X-ray photons, collected mainly by XMM-Newton, unveiled a number of unique characteristics, raising questions on this source very nature.

In August 2002, XMM-Newton devoted two full orbits to the observation of 1E 1207.4-5209, making this isolated neutron star the most deeply scrutinized galactic target of the mission. Thanks to the high throughput of the EPIC instrument, ~360,000 photons were collected from the source, allowing for a very sensitive study of the temporal and spectral behaviour of this object. The spectral data, both time-averaged and phase-resolved, yield one compelling interpretation of the observed features: cyclotron absorption from one fundamental (~0.7 keV) and three harmonics, at ~1.4, ~2.1 and ~2.8 keV. Possible physical consequences are discussed, also on the basis of the obvious phase variations of the features' shapes and depths. We also present deep VLT optical data which we have used to search for a counterpart, with negative results down to ~27.

The discovery of deep spectral features in the X-ray spectrum of 1E1207.4-5209 has pushed this Isolated Neutron Star (INS) out of the chorus line, since no other INS has shown significant features in its X-ray continuum. On August 2002, XMM-Newton devoted a two-orbit TOO observation to this target with the aim to better understand the nature of such spectral features, using much improved statistics. Indeed, the 260 ksec observation yielded 360,000 photons from 1E1207.4-5209, allowing for a very sensitive study of the temporal and spectral behaviour of this object.

SIMBOL-X is a hard X-ray mission, operating in the 0.5-70 keV range, which is proposed by a consortium of European laboratories for a launch around 2010. Relying on two spacecraft in a formation flying configuration, SIMBOL-X uses a 30 m focal length X-ray mirror to achieve an unprecedented angular resolution (30 arcsec HEW) and sensitivity (100 times better than INTEGRAL below 50 keV) in the hard X-ray range. SIMBOL-X will allow to elucidate fundamental questions in high energy astrophysics, such as the physics of accretion onto Black Holes, of acceleration in quasar jets and in supernovae remnants, or the nature of the hard X-ray diffuse emission. The scientific objectives and the baseline concepts of the mission and hardware design are presented.

Isolated neutron stars are higly magnetized, fast rotating end points of stellar evolution. They are now becoming directly observable through X-ray astronomy, owing to their high surface temperatures. In particular, features in their X-ray spectra could reveal the presence of atmospheres, or gauge their unknown magnetic fields through the cyclotron process, as in the classic case of X-ray binaries. All isolated neutron stars spectra observed so far, however, appear as featureless thermal continua. The unique exception is 1E 1207.4-5209. In its spectrum, previous observations had detected two deep absorption features, still too undefined for any unambiguous interpretation. Here we report on a much longer X-ray observation, in which the star's spectrum shows three distinct features, regularly spaced at 0.7, 1.4 and 2.1 keV, plus a fourth, lower significance one, at 2.8 keV. Such features vary in phase with the star rotation. Cyclotron resonant absorption is their logical interpretation, yielding a magnetic field strenght of 8x10^10 Gauss in the case of electrons. This is the first direct measure of an isolated neutron star magnetic field.

We have observed the radio-quiet X-ray pulsar 1E 1207.4-5209 with the high throughput EPIC cameras onboard XMM-Newton. The spectrum of this peculiar source is characterized by two broad absorption features which present significant substructures and show a clear phase-dependence. We believe that these features represent a strong evidence for the presence of a magnetized atmosphere containing heavy elements.

We present the results of an XMM-Newton observation of the radio-quiet X-ray pulsar 1E 1207-5209 located at the center of the shell-like supernova remnant G296.5+10.0. The X-ray spectrum is characterized by the presence of two phase-dependent absorption lines at energies ~0.7 keV and ~1.4 keV. Moreover, these broad spectral features have significant substructure, suggesting that they are due to the blending of several narrower lines. We interpret such features as evidence for an atmosphere containing metals and a magnetic field value of a few 10^12 G, consistent with the observed spin-down rate P_dot=(1.98+/-0.83)*10^(-14) s/s. Since 1E 1207-5209 is the only X-ray emitting pulsar showing evidence of such features, we tentatively link them to the unique combination of age and energetics that characterize this object. We suggest that a young age and a low level of magnetospheric activity are favorable conditions for the detection of atomic spectral features from metals in neutron star atmospheres, which would be either blanketed by a thin layer of accreted hydrogen in older objects or masked by non-thermal processes in young energetic pulsars.

The limited angular resolution of gamma-ray telescopes prevents the straight identification of the majority of the sources detected so far. While the only certified galactic gamma-ray emitters are associated with pulsars, more than 90% of the nearly 100 low latitude sources detected by EGRET on the GRO lack an identification. The best identification strategy devised over the years relies on the use of X-ray and optical data to single out possible counterparts, rating them on the basis of their extreme F_x/F_opt. Here we describe our multiwavelength programme based on the EPIC mapping of selected EGRET sources complemented by the optical coverage of the same fields obtained through the Wide Field Imager (WFI) of the ESO 2.2 m telescope. The field of view of the WFI, comparable to that of EPIC, will allow to directly cross-correlate X and optical data speeding up significantly the selection of interesting candidates worth follow-up optical (and X-ray) studies.

The distance to the Vela pulsar (PSR B0833-45) has been traditionally assumed to be 500 pc. Although affected by a significant uncertainty, this value stuck to both the pulsar and the SNR. In an effort to obtain a model free distance measurement, we have applied high resolution astrometry to the pulsar V~23.6 optical counterpart. Using a set of five HST/WFPC2 observations, we have obtained the first optical measurement of the annual parallax of the Vela pulsar. The parallax turns out to be 3.4 +/- 0.7 mas, implying a distance of 294(-50;+76) pc, i.e. a value significantly lower than previously believed. This affects the estimate of the pulsar absolute luminosity and of its emission efficiency at various wavelengths and confirms the exceptionally high value of the N_e towards the Vela pulsar. Finally, the complete parallax data base allows for a better measurement of the Vela pulsar proper motion (mu_alpha(cos(delta))=-37.2 +/- 1.2 mas/yr; mu_delta=28.2 +/- 1.3 mas/yr after correcting for the peculiar motion of the Sun) which, at the parallax distance, implies a transverse velocity of ~65 km/s. Moreover, the proper motion position angle appears specially well aligned with the axis of symmetry of the X-ray nebula as seen by Chandra. Such an alignment allows to assess the space velocity of the Vela pulsar to be ~81 km/s.

HST observations have contributed significantly to our knowledge on the behaviour of Isolated Neutron Stars (INSs) as optical emitters. First, HST has been instrumental both to discover new optical counterparts (PSR B1055-52, PSR B1929+10, PSR B0950+08) and to confirm proposed identifications (PSR B0656+14). Second, HST multicolor photometry provided useful information to characterize the optical emission mechanism(s) at work in middle-aged INSs like PSR B0656+14 and Geminga. Last, but not least, the superior angular resolution of the HST allowed both to perform a very accurate morphological study of the plerionic environments of young pulsars (e.g. the Crab and PSR B0540-69) and to perform very accurate astrometric measurements yielding proper motions (Crab, Vela, Geminga, PSR B0656+14) and parallaxes (Geminga).

High resolution images of the Crab-like SNR SNR0540-69 have been recently obtained by the HRC, on board the Chandra X-ray Observatory, and by the FORS1 instrument, attached to ESO/VLT1. A detailed comparison of these images with an archival Planetary Camera one will be presented.

Recent optical/UV data of Geminga (Mignani et al. 1998) have strengthen the significance of the emission feature present over the thermal continuum best fitting the EUV/soft X-ray data. Here we present a phenomenological model interpreting the data in terms of an ion cyclotron feature, arising from a thin outer layer of hot plasma covering Geminga's polar caps. The width of the feature favours for Geminga the strongly oblique rotator configuration, as already deduced from the soft X-ray data (Halpern & Ruderman 1993). Our interpretation implies a magnetic field of 3 - 5 10d11 G, consistent with the value deduced from the dynamical parameters of the pulsar. This would represent the first case of an independent measurement of the surface magnetic field of an isolated neutron star. The expected time modulation of Geminga's feature is also discussed.

We report on the optical observations of the gamma-ray pulsar PSR1706-44 performed with the Test Camera of the VLT-UT1 during the telescope Science Verification programme. With a limiting magnitude of V~ 27.5, these are the deepest optical observations ever performed for this pulsar, for which, however, no likely counterpart has been detected. The non-detection of PSR1706-44 sets an upper limit on its optical luminosity varying from 2 10^28 to 5 10^29 erg/ s, depending on the exact amount of interstellar absorption.

SNR0540-69 is a supernova remnant in the LMC, harbouring a young (~1600 yrs) radio/optical/X-ray pulsar (P=50 ms). Ground based Halpha imaging of the region has shown a unique spiral-like structure centered around the pulsar. In narrow band HST imaging, the feature seems resolved in a ring-like structure, probably ejected by the progenitor star in a pre-supernova phase (> 10^4 yrs ago).

New optical and near-UV HST observations of Geminga are presented. When compared with previous ground-based and HST imaging, the data confirm and better define the presence of an emission feature centered at 6,000 Angstrom and superimposed on the thermal continuum best fitting the extreme-UV/soft X-ray data. This feature may be interpreted in terms of cyclotron emission originated from a mixture of H/He ions in the neutron star's atmosphere. In the case of pure Hydrogen, the feature wavelength would imply a magnetic field of order 3-5 10d11 G, consistent with the value deduced from the dynamical parameters of the pulsar. If due to a cyclotron emission, the observation of this feature would represent the first case of an in situ measurement of the surface magnetic field of an isolated neutron star.

The EGRET catalogue of unidentified X-ray sources has more objects along the galactic disk than at high galactic latitude, where identifications are comparatively easier. On the other hand, the Egret/GRO mission has already identified several known radio pulsars as gamma-ray sources as well as discovering Geminga's nature as a pulsar. If Geminga is not a unique case, as it is very likely not to be, than other galactic sources could, in fact, be radio quiet isolated neutron stars. For these, the identification work is extremely difficult and should anyway start from high resolution X-ray/optical data.

In view of the data gathered in September 1997, we review the flux values collected so far for the "fuzziness" seen in the optical counterpart of GRB970228. Comparison between the ground based data collected in March and the data of September 1997 suggests a fading of the fuzz. Given the diversity of the data in hand, the magnitude of the effect and its significance are not easy to quantify. Only new images, both from the ground and with the Space Telescope, directly comparable to the old ones could settle this problem.

Accuracy in the absolute position in the sky is one of the limiting factors for pulsar timing, and timing parameters have a direct impact on the understanding of the physics of Isolated Neutron Stars (INS). We report here on a high-accuracy measurement of the optical position of Geminga (V=25.5), the only known radio-quiet INS. The procedure combines the Hipparcos and Tycho catalogues, ground-based astrometric data,and Hubble Space Telescope (HST) Wide Field Planetary Camera (WFPC2) images, to yield Geminga's absolute position to within ~40 mas (per coordinate). Such a positional accuracy, unprecedented for the optical position of a pulsar or an object this faint, is needed to combine in phase gamma-ray photons collected over more than 20 years, i.e. over 2.5 billions of star' revolutions. Although quite a difficult task, this is the only way to improve our knowledge of the timing parameters of this radio silent INS.

After a quarter of a century of gamma-ray burst (GRB) astronomy, the Italian-Dutch satellite BeppoSAX on Feb 28th, 1997 detected a soft X-ray afterglow from GRB 970228 and positioned it accurately. This made possible the successful detection of an optical transient. Two public Hubble Space Telescope (HST) images of the GRB/optical transient region were taken on March 26th and April 7th, 1997. They are analyzed here, with the purpose of understanding the nature of GRB 970228. We find that the position of the faint point-like object V ~ 26 seen at the transient location changed by 0.40 +/-0.10 pixels in 12 days, corresponding to a proper motion of ~ 550 mas/year. By comparison, four adjacent sources in the same field do not show any significant displacement, with astrometric residuals close to zero and average absolute displacements less than 0.09 pixels. If confirmed, this result would strongly support the galactic nature of GRB 970228.

After reviewing the multifrequency behaviour of the Isolated Neutron Stars detected so far, we analyze for each object the efficiency of conversion of the star's rotational energy loss into optical, X and gamma radiation. Although the number of pulsars detected at different wavelengths is rather limited, a pattern is seen to emerge from our analysis pointing towards the period derivative as the leading parameter to describe the multifrequency emission of INS. One object in particular, PSR 1509-58, stands out as archetypal for the Pdot-dependence of its luminosity at different wavelengths.