The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was $L_{\rm burst} \sim2\times 10^{39}$ erg/s (at 4.8 kpc). Simultaneous observations with XMM-Newton and NuSTAR three days after the burst provided a source spectrum well fit by an absorbed blackbody ($N_{\rm H} = (1.13\pm0.03) \times 10^{23}$ cm$^{-2}$ and $kT = 1.16\pm0.03$ keV) plus a power-law ($\Gamma=0.0\pm1.3$) in the 1-20 keV band, with a luminosity of $\sim$$8\times10^{34}$ erg/s, dominated by the blackbody emission. From our timing analysis, we derive a dipolar magnetic field $B \sim 7\times10^{14}$ G, spin-down luminosity $\dot{E}_{\rm rot} \sim 1.4\times10^{36}$ erg/s and characteristic age of 240 yr, the shortest currently known. Archival observations led to an upper limit on the quiescent luminosity $<$$5.5\times10^{33}$ erg/s, lower than the value expected from magnetar cooling models at the source characteristic age. A 1 hr radio observation with the Sardinia Radio Telescope taken about 1 week after the X-ray burst detected a number of strong and short radio pulses at 1.5 GHz, in addition to regular pulsed emission; they were emitted at an average rate 0.9 min$^{-1}$ and accounted for $\sim$50% of the total pulsed radio fluence. We conclude that Swift ,J1818.0-1607 is a peculiar magnetar belonging to the small, diverse group of young neutron stars with properties straddling those of rotationally and magnetically powered pulsars. Future observations will make a better estimation of the age possible by measuring the spin-down rate in quiescence.

We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).

The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts.

Magnetars are believed to host the strongest magnetic fields in the present universe ($B\gtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We present the results of an XMM-Newton observation of the slowly rotating ($P = 3.4$ s), highly magnetized ($B \approx 3 \times 10^{13}$ G) radio pulsar PSR J0726-2612. A previous X-ray observation with the Chandra satellite showed that some of the properties of PSR J0726-2612 are similar to those of the X-ray Dim Isolated Neutron Stars (XDINSs), a small class of nearby slow pulsars characterized by purely thermal X-ray spectra and undetected in the radio band. We confirm the thermal nature of the X-ray emission of PSR J0726-2612, which can be fit by the sum of two blackbodies with temperatures $kT_1 = 0.074_{-0.011}^{+0.006}$ keV and $kT_2 = 0.14_{-0.02}^{+0.04}$ keV and emitting radii $R_1 = 10.4_{-2.8}^{+10.8}$ km and $R_2 = 0.5_{-0.3}^{+0.9}$ km, respectively (assuming a distance of 1 kpc). A broad absorption line modeled with a Gaussian profile centred at $0.39_{-0.03}^{+0.02}$ keV is required in the fit. The pulse profile of PSR J0726-2612 is characterized by two peaks with similar intensity separated by two unequal minima, a shape and pulsed fraction that cannot be reproduced without invoking magnetic beaming of the X-ray emission. The presence of a single radio pulse suggests that in PSR J0726-2612 the angles that the dipole axis and the line of sight make with the rotation axis, $\xi$ and $\chi$ respectively, are similar. This geometry differs from that of the two radio-silent XDINSs with a double peaked pulse profile similar to that of PSR J0726-2612, for which $\xi \sim 90^\circ$ and $\chi \sim 45^\circ$ have been recently estimated. These results strengthen the similarity between PSR J0726-2612 and the XDINSs and support the possibility that the lack of radio emission from the latter might simply be due to an unfavourable viewing geometry.

Magnetic field geometry is expected to play a fundamental role in magnetar activity. The discovery of a phase-variable absorption feature in the X-ray spectrum of SGR 0418+5729, interpreted as cyclotron resonant scattering, suggests the presence of very strong non-dipolar components in the magnetic fields of magnetars. We performed a deep XMM-Newton observation of pulsar 1E 2259+586, to search for spectral features due to intense local magnetic fields. In the phase-averaged X-ray spectrum, we found evidence for a broad absorption feature at very low energy (0.7 keV). If the feature is intrinsic to the source, it might be due to resonant scattering/absorption by protons close to star surface. The line energy implies a magnetic field of ~ 10^14 G, roughly similar to the spin-down measure, ~ 6x10^13 G. Examination of the X-ray phase-energy diagram shows evidence for a further absorption feature, the energy of which strongly depends on the rotational phase (E >~ 1 keV ). Unlike similar features detected in other magnetar sources, notably SGR 0418+5729, it is too shallow and limited to a small phase interval to be modeled with a narrow phase-variable cyclotron absorption line. A detailed phase-resolved spectral analysis reveals significant phase-dependent variability in the continuum, especially above 2 keV. We conclude that all the variability with phase in 1E 2259+586 can be attributed to changes in the continuum properties which appear consistent with the predictions of the Resonant Compton Scattering model.

We observed the slowly revolving pulsar 1E 161348-5055 (1E 1613, spin period of 6.67 h) in the supernova remnant RCW 103 twice with XMM-Newton and once with the Very Large Telescope (VLT). The VLT observation was performed on 2016 June 30, about a week after the detection of a large outburst from 1E 1613. At the position of 1E 1613, we found a near-infrared source with K_S = 20.68 +/- 0.12 mag that was not detected (K_S > 21.2 mag) in data collected with the same instruments in 2006, during X-ray quiescence. Its position and behavior are consistent with a counterpart in the literature that was discovered with the Hubble Space Telescope in the following weeks in adjacent near-IR bands. The XMM-Newton pointings were carried out on 2016 August 19 and on 2018 February 14. While the collected spectra are similar in shape between each other and to what is observed in quiescence (a blackbody with kT~0.5 keV plus a second, harder component, either another hotter blackbody with kT ~ 1.2 keV or a power law with photon index ~3), the two pointings caught 1E 1613 at different luminosity throughout its decay pattern: about 4.8E34 erg/s in 2016 and 1.2E34 erg/s in 2018 (0.5-10 keV, for the double-blackbody model and for 3.3 kpc), which is still almost about ten times brighter than the quiescent level. The pulse profile displayed dramatic changes, apparently evolving from the complex multi-peak morphology observed in high-luminosity states to the more sinusoidal form characteristic of latency. The inspection of the X-ray light curves revealed two flares with unusual properties in the 2016 observation: they are long (~1 ks to be compared with 0.1-1 s of typical magnetar bursts) and faint (~1E34 erg/s, with respect to 1E38 erg/s or more in magnetars). Their spectra are comparatively soft and resemble the hotter thermal component of the persistent emission.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the campaign, Swift BAT registered the occurrence of multiple bursts, accompanied by two other enhancements of the X-ray persistent flux. The long time span covered by our observations allowed us to study the spectral as well as the timing evolution of the source. After 11 months since the 2017 May outburst onset, the observed flux was about 15 times higher than its historical minimum level and a factor of 3 higher than the level reached after the 2006 outburst. This suggests that the crust has not fully relaxed to the quiescent level, or that the source quiescent level has changed following the multiple outburst activities in the past 10 years or so. This is another case of multiple outbursts from the same source on a yearly time scale, a somehow recently discovered behaviour in magnetars.

PSR B0943+10 is a mode-switching radio pulsar characterized by two emission modes with different radio and X-ray properties. Previous studies, based on simple combinations of blackbody and power law models, showed that its X-ray flux can be decomposed in a pulsed thermal plus an unpulsed non-thermal components. However, if PSR B0943+10 is a nearly aligned rotator seen pole-on, as suggested by the radio data, it is difficult to reproduce the high observed pulsed fraction unless magnetic beaming is included. In this work we reanalyze all the available X-ray observations of PSR B0943+10 with simultaneous radio coverage, modeling its thermal emission with polar caps covered by a magnetized hydrogen atmosphere or with a condensed iron surface. The condensed surface model provides good fits to the spectra of both pulsar modes, but, similarly to the blackbody, it can not reproduce the observed pulse profiles, unless an additional power law with an ad hoc modulation is added. Instead, the pulse profiles and phase-resolved spectra are well described using the hydrogen atmosphere model to describe the polar cap emission, plus an unpulsed power law. For the X-ray brighter state (Q-mode) we obtain a best fit with a temperature kT~0.09 keV, an emitting radius R~260 m, a magnetic field consistent with the value of the dipole field of 4x10^12 G inferred from the timing parameters, and a small angle between the magnetic and spin axis, $\xi$=5. The corresponding parameters for the X-ray fainter state (B-mode) are kT~0.08 keV and R~170 m.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

In 2003, the magnetar XTE J1810-197 started an outburst that lasted until early 2007. In the following 11 years, the source stayed in a quiescent/low activity phase. XTE J1810-197 is one of the closest magnetars, hence its X-ray properties can be studied in detail even in quiescence and an extended monitoring has been carried out to study its long term timing and spectral evolution. Here, we report the results of new X-ray observations, taken between September 2017 and April 2018, with XMM-Newton, Chandra and Nicer. We derived a phase-connected timing solution yielding a frequency derivative of -9.26(6)x10^-14 Hz s-1. This value is consistent with that measured between 2009 and 2011, indicating that the pulsar spin-down rate remained quite stable during the long quiescent period. A spectral analysis of all the X-ray observations taken between 2009 and 2018 does not reveal significant spectral and/or flux variability. The spectrum of XTE J1810-197 can be described by the sum of two thermal components with temperatures of 0.15 and 0.3 keV, plus a power law component with photon index 0.6. We also found evidence for an absorption line at ~1.2 keV and width of 0.1 keV. Thanks to the long exposure time of the summed XMM-Newton observations, we could also carry out a phase-resolved spectral analysis for this source in quiescence. This showed that the flux modulation can be mainly ascribed to the the warmer of the two thermal components, whose flux varies by ~45 per cent along the pulse phase.

The magnetosphere of strongly magnetized neutron stars, such as magnetars, can sustain large electric currents. The charged particles return to the surface with large Lorentz factors, producing a particle bombardment. We investigate the transport of radiation in the atmosphere of strongly magnetized neutron stars, in the presence of particle bombardment heating. We solve the radiative transfer equations for a gray atmosphere in the Eddington approximation, accounting for the polarization induced by a strong magnetic field. The solutions show the formation of a hot external layer and a low (uniform) temperature atmospheric interior. This suggests that the emergent spectrum may be described by a single blackbody with the likely formation of a optical/infrared excess (below ~ 1 eV). We also found that the polarization is strongly dependent on both the luminosity and penetration length of the particle bombardment. Therefore, the thermal emission from active sources, such as transient magnetars, in which the luminosity decreases by orders of magnitude, may show a substantial variation in the polarization pattern during the outburst decline. Our results may be relevant in view of future X-ray polarimetric missions such as IXPE and eXTP.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array.

Isolated Neutron Stars are some of the most exciting stellar objects known to astronomers: they have the most extreme magnetic fields, with values up to $10^{15}$ G, and, with the exception of stellar-mass black holes, they are the most dense stars, with densities of $\approx 10^{14}$ g cm$^{-3}$. As such, they are perfect laboratories to test theories of electromagnetism and nuclear physics under conditions of magnetic field and density unattainable on Earth. In particular, the interaction of radiation with strong magnetic fields is the cause of the {\em vacuum birefringence}, an effect predicted by quantum electrodynamics in 1936 but that lacked an observational evidence until now. Here, we show how the study of the polarisation of the optical radiation from the surface of an isolated neutron star yielded such an observational evidence, opening exciting perspectives for similar studies at other wavelengths.

New simultaneous X-ray and radio observations of the archetypal mode-switching pulsar PSR B0943+10 have been carried out with XMM-Newton and the LOFAR, LWA and Arecibo radio telescopes in November 2014. They allowed us to better constrain the X-ray spectral and variability properties of this pulsar and to detect, for the first time, the X-ray pulsations also during the X-ray-fainter mode. The combined timing and spectral analysis indicates that unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap are present during both radio modes and vary in a correlated way.

Recently, numerical calculations of the magnetic field evolution in neutron stars demonstrated the possible existence of a Hall attractor, a stage at which the evolution of the field driven by the Hall cascade ends. The existence of such a stage in neutron star magnetic evolution is very important, and can be potentially probed by observations. Here we discuss three types of objects which could have reached this stage. First, we briefly describe the evolution of normal radio pulsars with ages about a few hundred thousand years. Then we analyse in more detail observations of RX J1856.5-3754, one of the Magnificent Seven, focusing on the surface temperature distribution and comparing model predictions with the temperature map inferred from X-ray observations. Finally, we discuss the necessity of the Hall attractor stage to explain the hypothetical existence of accreting magnetars. We conclude that at the moment there is no direct confirmation of the Hall attractor stage in known sources. However, more detailed observations in the near future can demonstrate existence (or absence) of this stage of the crustal magnetic field evolution.

We report on 3.5 years of Chandra monitoring of the Galactic Centre magnetar SGR J1745-2900 since its outburst onset in April 2013. The magnetar spin-down has shown at least two episodes of period derivative increases so far, and it has slowed down regularly in the past year or so. We observed a slightly increasing trend in the time evolution of the pulsed fraction, up to about 55 per cent in the most recent observations. SGR J1745-2900 has not reached the quiescent level yet, and so far the overall outburst evolution can be interpreted in terms of a cooling hot region on the star surface. We discuss possible scenarios, showing in particular how the presence of a shrinking hot spot in this source is hardly reconcilable with internal crustal cooling and favors the untwisting bundle model for this outburst. Moreover, we also show how the emission from a single uniform hot spot is incompatible with the observed pulsed fraction evolution for any pair of viewing angles, suggesting an anisotropic emission pattern.

The recent detection of a large polarization degree in the optical emission of an isolated neutron star led to the suggestion that this has been the first evidence of vacuum polarization in a strong magnetic field, an effect predicted by quantum electrodynamics but never observed before. This claim was challanged in a paper by Capparelli, Maiani & Polosa (2017), according to whom a much higher polarization degree would be necessary to positively identify vacuum polarization. Here we show that their conclusions are biased by several inadequate assumptions and have no impact on the original claim.

e-ASTROGAM (`enhanced ASTROGAM') is a breakthrough Observatory mission dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. In the largely unexplored MeV-GeV domain, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on Galactic ecosystems. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and the promise of eLISA. Keywords: High-energy gamma-ray astronomy, High-energy astrophysics, Nuclear Astrophysics, Compton and Pair creation telescope, Gamma-ray bursts, Active Galactic Nuclei, Jets, Outflows, Multiwavelength observations of the Universe, Counterparts of gravitational waves, Fermi, Dark Matter, Nucleosynthesis, Early Universe, Supernovae, Cosmic Rays, Cosmic antimatter.

Bursts and flares are among the distinctive observational manifestations of magnetars, isolated neutron stars endowed with an ultra-strong magnetic field ($B \approx 10^{14}$--$10^{15}$ G). It is believed that these events arise in a hot electron-positron plasma which remains trapped within the closed magnetic field lines. We developed a simple radiative transfer model to simulate magnetar flare emission in the case of a steady trapped fireball. After dividing the fireball surface in a number of plane-parallel slabs, the local spectral and polarization properties are obtained integrating the radiative transfer equations for the two normal modes. We assume that magnetic Thomson scattering is the dominant source of opacity, and neglect contributions from second-order radiative processes, although double-Compton scattering is accounted for in establishing local thermal equilibrium in the fireball atmospheric layers. The observed spectral and polarization properties as measured by a distant observer are obtained summing the contributions from the patches which are visible for a given viewing geometry by means of a ray-tracing code. The spectra we obtained in the 1-100 keV energy range are thermal and can be described in terms of the superposition of two blackbodies. The blackbody temperature and the emitting area ratio are in broad agreement with the available observations. The predicted linear polarization degree is in general greater than 80% over the entire energy range and should be easily detectable by new-generation X-ray polarimeters, like IXPE, XIPE and eXTP.

We report on the results of a detailed phase-resolved spectroscopy of archival XMM--Newton observations of X-ray Dim Isolated Neutron Stars (XDINSs). Our analysis revealed a narrow and phase-variable absorption feature in the X-ray spectrum of RX J1308.6+2127. The feature has an energy of $\sim$740 eV and an equivalent width of $\sim$15 eV. It is detected only in $\sim$ 1/5 of the phase cycle, and appears to be present for the entire timespan covered by the observations (2001 December - 2007 June). The strong dependence on the pulsar rotation and the narrow width suggest that the feature is likely due to resonant cyclotron absorption/scattering in a confined high-B structure close to the stellar surface. Assuming a proton cyclotron line, the magnetic field strength in the loop is B$_{loop} \sim 1.7 \times 10^{14}$ G, about a factor of $\sim$5 higher than the surface dipolar magnetic field (B$_{surf} \sim 3.4 \times 10^{13}$ G). This feature is similar to that recently detected in another XDINS, RX J0720.4-3125, showing (as expected by theoretical simulations) that small scale magnetic loops close to the surface might be common to many highly magnetic neutron stars (although difficult to detect with current X-ray instruments). Furthermore, we investigated the available XMM--Newton, data of all XDINSs in search for similar narrow phase-dependent features, but could derive only upper limits for all the other sources.

Ultraluminous x-ray sources (ULXs) in nearby galaxies shine brighter than any X-ray source in our Galaxy. ULXs are usually modeled as stellar-mass black holes (BHs) accreting at very high rates or intermediate-mass BHs. We present observations showing that NGC5907 ULX is instead an x-ray accreting neutron star (NS) with a spin period evolving from 1.43~s in 2003 to 1.13~s in 2014. It has an isotropic peak luminosity of about 1000 times the Eddington limit for a NS at 17.1~Mpc. Standard accretion models fail to explain its luminosity, even assuming beamed emission, but a strong multipolar magnetic field can describe its properties. These findings suggest that other extreme ULXs (x-ray luminosity > 10^{41} erg/s) might harbor NSs.

The "Magnificent Seven" (M7) are a group of radio-quiet Isolated Neutron Stars (INSs) discovered in the soft X-rays through their purely thermal surface emission. Owing to the large inferred magnetic fields ($B\approx 10^{13}$ G), radiation from these sources is expected to be substantially polarised, independently on the mechanism actually responsible for the thermal emission. A large observed polarisation degree is, however, expected only if quantum-electrodynamics (QED) polarisation effects are present in the magnetised vacuum around the star. The detection of a strongly linearly polarised signal would therefore provide the first observational evidence of QED effects in the strong-field regime. While polarisation measurements in the soft X-rays are not feasible yet, optical polarisation measurements are within reach also for quite faint targets, like the M7 which have optical counterparts with magnitudes $\approx 26$--$28$. Here, we report on the measurement of optical linear polarisation for the prototype, and brightest member, of the class, RX\, J1856.5$-$3754 ($V\sim 25.5$), the first ever for one of the M7, obtained with the Very Large Telescope. We measured a polarisation degree $\mathrm{P.D.} =16.43\% \pm5.26\%$ and a polarisation position angle $\mathrm{P.A.}=145\fdg39\pm9\fdg44$, computed east of the North Celestial Meridian. The $\mathrm{P.D.}$ that we derive is large enough to support the presence of vacuum birefringence, as predicted by QED.

The evolution of magnetic field in isolated neutron stars is one of the most important ingredients in the attempt to build a unified description of these objects. A prediction of field evolution models is the existence of an equilibrium configuration, in which the Hall cascade vanishes. Recent calculations have explored the field structure in this stage, called the Hall attractor. We use X-ray data of near-by, cooling neutron stars to probe this prediction, as these sources are surmised to be close to or at Hall attractor phase. We show that the source RX J1856.5-3754 might be closer to the attractor than other sources of its class. Our modelling indicates that the properties of surface thermal emission, assuming that the star is in the Hall attractor, are in contradiction with the spectral data of RX J1856.5-3754.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

We report on simultaneous X-ray and radio observations of the mode-switching pulsar PSR B0943+10 obtained with the XMM-Newton satellite and the LOFAR, LWA and Arecibo radio telescopes in November 2014. We confirm the synchronous X-ray/radio switching between a radio-bright (B) and a radio-quiet (Q) mode, in which the X-ray flux is a factor ~2.4 higher than in the B-mode. We discovered X-ray pulsations, with pulsed fraction of 38+/-5% (0.5-2 keV), during the B-mode, and confirm their presence in Q-mode, where the pulsed fraction increases with energy from ~20% up to ~65% at 2 keV. We found marginal evidence for an increase in the X-ray pulsed fraction during B-mode on a timescale of hours. The Q-mode X-ray spectrum requires a fit with a two-component model (either a power-law plus blackbody or the sum of two blackbodies), while the B-mode spectrum is well fit by a single blackbody (a single power-law is rejected). With a maximum likelihood analysis, we found that in Q-mode the pulsed emission has a thermal blackbody spectrum with temperature ~3.4x10^6 K and the unpulsed emission is a power-law with photon index ~2.5, while during B-mode both the pulsed and unpulsed emission can be fit by either a blackbody or a power law with similar values of temperature and photon index. A Chandra image shows no evidence for diffuse X-ray emission. These results support a scenario in which both unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap (~1500 m^2) with a strong non-dipolar magnetic field (~10^{14} G), are present during both radio modes and vary in intensity in a correlated way. This is broadly consistent with the predictions of the partially screened gap model and does not necessarily imply global magnetospheric rearrangements to explain the mode switching.

We have analyzed XMM-Newton and Chandra observations of the transient magnetar XTE J1810-197, spanning more than 11 years, from the initial phases of the 2003 outburst to the current quiescent level. We investigated the evolution of the pulsar spin period and we found evidence for two distinct regimes: during the outburst decay, the spin derivative (nu_dot) was highly variable in the range -(2-4.5)x10^-13 Hz/s, while during quiescence the spin-down rate was more stable at an average value of -1x10^{-13} Hz/s. Only during ~3000 days (from MJD 54165 to MJD 56908) in the quiescent stage, it was possible to find a phase-connected timing solution, with nu_dot = -4.9x10^{-14} Hz/s, and a positive second frequency derivative, nu_dotdot = 1.8x10^{-22} Hz/s. These results are in agreement with the behavior expected if the outburst of XTE J1810-197 was due to a strong magnetospheric twist.

The observed polarization properties of thermal radiation from isolated, cooling neutron stars depend on both the emission processes at the surface and the effects of the magnetized vacuum which surrounds the star. Here we investigate the polarized thermal emission from X-ray Dim Isolated Neutron Stars, taking RX J1856.5-3754 as a representative case. The physical conditions of the star outermost layers in these sources is still debated, and so we consider emission from a magnetized atmosphere and a condensed surface, accounting for the effects of vacuum polarization as the radiation propagates in the star magnetosphere. We have found that, for a significant range of viewing geometries, measurement of the phase-averaged polarization fraction and phase-averaged polarization angle at both optical and X-ray wavelengths allow us to determine whether this neutron star has an atmosphere or a condensed surface. Our results may therefore be relevant in view of future developments of soft X-ray polarimeters.

We report on the discovery of a new member of the magnetar class, SGR J1935+2154, and on its timing and spectral properties measured by an extensive observational campaign carried out between July 2014 and March 2015 with Chandra and XMM-Newton (11 pointings). We discovered the spin period of SGR J1935+2154 through the detection of coherent pulsations at a period of about 3.24s. The magnetar is slowing-down at a rate of 1.43(1)x10^{-11} s/s and with a decreasing trend due to a negative second period derivative of -3.5(7)x10^{-19} s/s^2. This implies a surface dipolar magnetic field strength of about 2.2x10^{14} G, a characteristic age of about 3.6kyr and, a spin-down luminosity L_{sd} of about 1.7x10^{34} erg/s. The source spectrum is well modelled by a blackbody with temperature of about 500eV plus a power-law component with photon index of about 2. The source showed a moderate long-term variability, with a flux decay of about 25\% during the first four months since its discovery, and a re-brightening of the same amount during the second four months. The X-ray data were also used to study the source environment. In particular, we discovered a diffuse emission extending on spatial scales from about 1" up to at least 1' around SGR J1935+2154 both in Chandra and XMM-Newton data. This component is constant in flux (at least within uncertainties) and its spectrum is well modelled by a power-law spectrum steeper than that of the pulsar. Though a scattering halo origin seems to be more probable we cannot exclude that part, or all, of the diffuse emission is due to a pulsar wind nebula.

Giant $\gamma$-ray flares comprise the most extreme radiation events observed from magnetars. Developing on (sub)millisecond timescales and generating vast amounts of energy within a fraction of a second, the initial phase of these extraordinary bursts present a significant challenge for candidate trigger mechanisms. Here we assess and critically analyse the linear growth of the relativistic tearing instability in a globally twisted magnetosphere as the trigger mechanism for giant $\gamma$-ray flares. Our main constraints are given by the observed emission timescales, the energy output of the giant flare spike, and inferred dipolar magnetic field strengths. We find that the minimum growth time of the linear mode is comparable to the $e$-folding rise time, i.e. $\sim10^{-1}$ ms. With this result we constrain basic geometric parameters of the current sheet. We also discuss the validity of the presumption that the $e$-folding emission timescale may be equated with the growth time of an MHD instability.

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (10^14 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (~3 x 10^13 G) magnetic dipole field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

New-generation X-ray polarimeters currently under development promise to open a new window in the study of high-energy astrophysical sources. Among them, neutron stars appear particularly suited for polarization measurements. Radiation from the (cooling) surface of a neutron star is expected to exhibit a large intrinsic polarization degree due to the star strong magnetic field ($\approx 10^{12}-10^{15}$ G), which influences the plasma opacity in the outermost stellar layers. The polarization fraction and polarization angle as measured by an instrument, however, do not necessary coincide with the intrinsic ones derived from models of surface emission. This is due to the effects of quantum electrodynamics in the highly magnetized vacuum around the star (the vacuum polarization) coupled with the rotation of the Stokes parameters in the plane perpendicular to the line of sight induced by the non-uniform magnetic field. Here we revisit the problem and present an efficient method for computing the observed polarization fraction and polarization angle in the case of radiation coming from the entire surface of a neutron star, accounting for both vacuum polarization and geometrical effects due to the extended emitting region. Our approach is fairly general and is illustrated in the case of blackbody emission from a neutron star with either a dipolar or a (globally) twisted magnetic field.

The isolated neutron star (INS) 2XMM J104608.7-594306 is one of the only two to be discovered through their thermal emission since the ROSAT era. In a first dedicated XMM-Newton observation of the source, we found intriguing evidence of a very fast spin period. We re-observed 2XMM J104608.7-594306 with XMM-Newton to better characterise the spectral energy distribution of the source, confirm the candidate spin period, and possibly constrain the pulsar spin-down. Statistically acceptable spectral fits and meaningful physical parameters for the source are only obtained when the purely thermal spectrum is modified by at least one line in absorption. The implied distance is consistent with a location in (or in front of) the Carina nebula, and radiation radii are compatible with emission originating on most of the surface. Non-thermal X-ray emission is ruled out at levels above 0.5% of the source luminosity. Unfortunately, the second XMM-Newton observation proved inconclusive in terms of confirming (discarding) the fast candidate spin, providing an upper limit on the pulsed fraction of the source that is very close to the limiting sensitivity for detecting the modulation found previously. In the absence of an unambiguous period determination and an estimate of the magnetic field, the nature of the source remains open to interpretation. Its likely association with the Carina cluster and its overall spectral properties (only shared by a handful of other peculiar INSs) disfavour a standard evolutionary path, or one in which the source was previously recycled by accretion in a binary system. The INS 2XMM J104608.7-594306 may be similar to Calvera (1RXS J141256.0+792204), a neutron star for which the scenario of an evolved `anti-magnetar' has been discussed. A better age estimate and deeper radio and gamma-ray limits are required to further constrain the evolutionary state of the neutron star.

Magnetars are the strongest magnets in the present universe and the combination of extreme magnetic field, gravity and density makes them unique laboratories to probe current physical theories (from quantum electrodynamics to general relativity) in the strong field limit. Magnetars are observed as peculiar, burst--active X-ray pulsars, the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma Repeaters (SGRs); the latter emitted also three "giant flares," extremely powerful events during which luminosities can reach up to 10^47 erg/s for about one second. The last five years have witnessed an explosion in magnetar research which has led, among other things, to the discovery of transient, or "outbursting," and "low-field" magnetars. Substantial progress has been made also on the theoretical side. Quite detailed models for explaining the magnetars' persistent X-ray emission, the properties of the bursts, the flux evolution in transient sources have been developed and confronted with observations. New insight on neutron star asteroseismology has been gained through improved models of magnetar oscillations. The long-debated issue of magnetic field decay in neutron stars has been addressed, and its importance recognized in relation to the evolution of magnetars and to the links among magnetars and other families of isolated neutron stars. The aim of this paper is to present a comprehensive overview in which the observational results are discussed in the light of the most up-to-date theoretical models and their implications. This addresses not only the particular case of magnetar sources, but the more fundamental issue of how physics in strong magnetic fields can be constrained by the observations of these unique sources.

We present a long-term phase-coherent timing analysis and pulse-phase resolved spectroscopy for the two outbursts observed from the transient anomalous X-ray pulsar CXOU J164710.2-45521. For the first outburst we used 11 Chandra and XMM-Newton observations between September 2006 to August 2009, the longest baseline yet for this source. We obtain a coherent timing solution with $P=10.61065583(4)$ s, $\dot{P} = 9.72(1) \times 10^{-13}\;$s s$^{-1}$ and $\ddot{P} = -1.05(5)\times10^{-20}\; $s s$^{-2}$. Under the standard assumptions this implies a surface dipolar magnetic field of $\sim 10^{14}$ G, confirming this source as a standard-$B$ magnetar. We also study the evolution of the pulse profile (shape, intensity and pulsed fraction) as a function of time and energy. Using the phase-coherent timing solution we perform a phase-resolved spectroscopy analysis, following the spectral evolution of pulse-phase features, which hints at the physical processes taking place on the star. The results are discussed from the perspective of magneto-thermal evolution models and the untwisting magnetosphere model. Finally, we present similar analysis for the second, less intense, 2011 outburst. For the timing analysis we used Swift data together with 2 XMM-Newton and Chandra pointings. The results inferred for both outbursts are compared and briefly discussed in a more general framework.

We present the discovery of a strongly phase-variable absorption feature in the X-ray spectrum of the nearby, thermally-emitting, isolated neutron star RX J0720.4-3125. The absorption line was detected performing detailed phase-resolved spectroscopy in 20 XMM-Newton observations, covering the period May 2000 - September 2012. The feature has an energy of ~750eV, an equivalent width of ~30eV, and it is significantly detected for only ~20% of the pulsar rotation. The absorption feature appears to be stable over the timespan covered by the observations. Given its strong dependence on the pulsar rotational phase and its narrow width, a plausible interpretation is in terms of resonant proton cyclotron absorption/scattering in a confined magnetic structure very close to the neutron star surface. The inferred field in such a magnetic loop is B_loop ~ 2 x 10^{14} G, a factor of ~7 higher than the surface dipolar magnetic field.

In 2013 April a new magnetar, SGR 1745-2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the Centre of the Milky Way, Sagittarius A*. SGR 1745-2900 has a surface dipolar magnetic field of ~ 2x10^{14} G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity L_X ~ 5x10^{35} erg s^{-1}. Here we report on the long-term Chandra (25 observations) and XMM-Newton (8 observations) X-ray monitoring campaign of SGR 1745-2900, from the onset of the outburst in April 2013 until September 2014. This unprecedented dataset allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ~2 around June 2013, and reveals that a further increase occurred between 2013 Oct 30 and 2014 Feb 21. We find that the period derivative changed from 6.6x10^{-12} s s^{-1} to 3.3x10^{-11} s s^{-1} in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ~200 days of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of magnetospheres of isolated neutron stars. For a summary, we refer to the paper.

Diffuse radio emission was detected around the soft gamma-ray repeater SGR 1806-20, after its 2004 powerful giant flare. We study the possible extended X-ray emission at small scales around SGR 1806-20, in two observations by the High Resolution Camera Spectrometer (HRC-S) on board of the Chandra X-ray Observatory: in 2005, 115 days after the giant flare, and in 2013, during quiescence. We compare the radial profiles extracted from data images and PSF simulations, carefully considering various issues related with the uncertain calibration of the HRC PSF at sub-arcsecond scales. We do not see statistically significant excesses pointing to an extended emission on scales of arcseconds. As a consequence, SGR 1806-20 is compatible with being point-like in X-rays, months after the giant flare, as well as in quiescence.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We report on the quiescent state of the Soft Gamma Repeater SGR 0501+4516 observed by XMM-Newton on 2009 August 30. The source exhibits an absorbed flux ~75 times lower than that measured at the peak of the 2008 outburst, and a rather soft spectrum, with the same value of the blackbody temperature observed with ROSAT back in 1992. This new observation is put into the context of all existing X-ray data since its discovery in August 2008, allowing us to complete the study of the timing and spectral evolution of the source from outburst until its quiescent state. The set of deep XMM-Newton observations performed during the few-years timescale of its outburst allows us to monitor the spectral characteristics of this magnetar as a function of its rotational period, and their evolution along these years. After the first ~10 days, the initially hot and bright surface spot progressively cooled down during the decay. We discuss the behaviour of this magnetar in the context of its simulated secular evolution, inferring a plausible dipolar field at birth of 3x10^14 G, and a current (magneto-thermal) age of ~10 kyr.

The study of magnetars is of particular relevance since these objects are the only laboratories where the physics in ultra-strong magnetic fields can be directly tested. Until now, spectroscopic and timing measurements at X-ray energies in soft gamma-repeaters (SGRs) and anomalous X-ray pulsar (AXPs) have been the main source of information about the physical properties of a magnetar and of its magnetosphere. Spectral fitting in the ~ 0.5-10 keV range allowed to validate the "twisted magnetosphere" model, probing the structure of the external field and estimating the density and velocity of the magnetospheric currents. Spectroscopy alone, however, may fail in disambiguating the two key parameters governing magnetospheric scattering (the charge velocity and the twist angle) and is quite insensitive to the source geometry. X-ray polarimetry, on the other hand, can provide a quantum leap in the field by adding two extra observables, the linear polarization degree and the polarization angle. Using the bright AXP 1RXS J170849.0-400910 as a template, we show that phase-resolved polarimetric measurements can unambiguously determine the model parameters, even with a small X-ray polarimetry mission carrying modern photoelectric detectors and existing X-ray optics. We also show that polarimetric measurements can pinpoint vacuum polarization effects and thus provide an indirect evidence for ultra-strong magnetic fields.

It is now widely accepted that soft gamma repeaters and anomalous X-ray pulsars are the observational manifestations of magnetars, i.e. sources powered by their own magnetic energy. This view was supported by the fact that these `magnetar candidates' exhibited, without exception, a surface dipole magnetic field (as inferred from the spin-down rate) in excess of the electron critical field (~4.4E+13 G). The recent discovery of fully-qualified magnetars, SGR 0418+5729 and Swift J1822.3-1606, with dipole magnetic field well in the range of ordinary radio pulsars posed a challenge to the standard picture, showing that a very strong field is not necessary for the onset of magnetar activity (chiefly bursts and outbursts). Here we summarize the observational status of the low-magnetic-field magnetars and discuss their properties in the context of the mainstream magnetar model and its main alternatives.

The variety of the observational appearance of young isolated neutron stars must find an explanation in the framework of some unifying approach. Nowadays it is believed that such scenario must include magnetic field decay, the possibility of magnetic field emergence on a time scale $\lesssim 10^4$--$10^5$ yrs, significant contribution of non-dipolar fields, and appropriate initial parameter distributions. We present our results on the initial spin period distribution, and suggest that inconsistences between distributions derived by different methods for samples with different average ages can uncover field decay or/and emerging field. We describe a new method to probe the magnetic field decay in normal pulsars. The method is a modified pulsar current approach, where we study pulsar flow along the line of increasing characteristic age for constant field. Our calculations, performed with this method, can be fitted with an exponential decay for ages in the range $8\times 10^4$--$3.5 \times 10^5$ yrs with a time scale $\sim 5 \times 10^5$ yrs. We discuss several issues related to the unifying scenario. At first, we note that the dichotomy, among local thermally emitting neutron stars, between normal pulsars and the Magnificent Seven remains unexplained. Then we discuss the role of high-mass X-ray binaries in the unification of neutron star evolution. We note, that such systems allow to check evolutionary effects on a time scale longer than what can be probed with normal pulsars alone. We conclude with a brief discussion of importance of discovering old neutron stars accreting from the interstellar medium.

The center of our Galaxy hosts a supermassive black hole, Sagittarius (Sgr) A*. Young, massive stars within 0.5 pc of SgrA* are evidence of an episode of intense star formation near the black hole a few Myr ago, which might have left behind a young neutron star traveling deep into SgrA*'s gravitational potential. On 2013 April 25, a short X-ray burst was observed from the direction of the Galactic center. Thanks to a series of observations with the Chandra and the Swift satellites, we pinpoint the associated magnetar at an angular distance of 2.4+/-0.3 arcsec from SgrA*, and refine the source spin period and its derivative (P=3.7635537(2) s and \dot{P} = 6.61(4)x10^{-12} s/s), confirmed by quasi simultaneous radio observations performed with the Green Bank (GBT) and Parkes antennas, which also constrain a Dispersion Measure of DM=1750+/-50 pc cm^{-3}, the highest ever observed for a radio pulsar. We have found that this X-ray source is a young magnetar at ~0.07-2 pc from SgrA*. Simulations of its possible motion around SgrA* show that it is likely (~90% probability) in a bound orbit around the black hole. The radiation front produced by the past activity from the magnetar passing through the molecular clouds surrounding the Galactic center region, might be responsible for a large fraction of the light echoes observed in the Fe fluorescence features.

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 old pulsar PSR B0943+10 (P=1.1 s, characteristic age tau=5 Myr) is the best example of mode-switching radio pulsar. Its radio emission alternates between a highly organized state with regular drifting subpulses (B mode) and a chaotic emission pattern (Q mode). We present the results of XMM-Newton observations showing that the X-ray properties of PSR B0943+10 depend on its radio state (Hermsen et al. 2013). During the radio fainter state (Q mode) the X-ray flux is more than a factor two larger than during the B-mode and X-ray pulsations with about 50% pulsed fraction are detected. The X-ray emission of PSR B0943+10 in the B-mode is well described by thermal emission with blackbody temperature kT=0.26 keV coming from a small hot spot with luminosity of 7x10^28 erg/s, in good agreement with the prediction of the partially screened gap model, which also explains the properties of the radio emission in this mode. We derived an upper limit of 46% on the X-ray pulsed fraction in the B-mode, consistent with the geometry and viewing angle of PSR B0943+10 inferred from the radio data. The higher flux observed during the Q-mode is consistent with the appearance of an additional component with a power-law spectrum with photon index 2.2. We interpret it as pulsed non-thermal X-rays produced in the star magnetosphere. A small change in the beaming pattern or in the efficiency of acceleration of the particles responsible for the non-thermal emission can explain the reduced flux of this component during the radio B-mode.

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729, using all the available X-ray data obtained with RXTE, SWIFT, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: \dot{P}=4(1)x10^{-15} s/s, significant at ~3.5 sigma confidence level. This leads to a surface dipolar magnetic field of B_dip ~6x 10^{12} G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ~1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ~0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4x10^{-11} to 1.2x10^{-14} erg/s/cm^2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ~550 kyr, and a dipolar magnetic field at birth of ~10^{14} G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.

The problem of computing the pulse profiles from thermally emitting spots on the surface of a neutron star in general relativity is reconsidered. We show that it is possible to extend Beloborodov (2002) approach to include (multiple) spots of finite size in different positions on the star surface. Results for the pulse profiles are expressed by comparatively simple analytical formulas which involve only elementary functions.

X-ray observations unveiled various types of radio-silent Isolated Neutron Stars (INSs), phenomenologically very diverse, e.g. the Myr old X-ray Dim INS (XDINSs) and the kyr old magnetars. Although their phenomenology is much diverse, the similar periods (P=2--10 s) and magnetic fields (~10^{14} G) suggest that XDINSs are evolved magnetars, possibly born from similar populations of supermassive stars. One way to test this hypothesis is to identify their parental star clusters by extrapolating backward the neutron star velocity vector in the Galactic potential. By using the information on the age and space velocity of the XDINS RX J1856.5-3754, we computed backwards its orbit in the Galactic potential and searched for its parental stellar cluster by means of a closest approach criterion. We found a very likely association with the Upper Scorpius OB association, for a neutron star age of 0.42+/-0.08 Myr, a radial velocity V_r^NS =67+/- 13$ km s^{-1}, and a present-time parallactic distance d_\pi^NS = 123^{+11}_{-15} pc. Our result confirms that the "true" neutron star age is much lower than the spin-down age (tau_{sd}=3.8 Myrs), and is in good agreement with the cooling age, as computed within standard cooling scenarios. The mismatch between the spin-down and the dynamical/cooling age would require either an anomalously large breaking index (n~20) or a decaying magnetic field with initial value B_0 ~ 10^{14} G. Unfortunately, owing to the uncertainty on the age of the Upper Scorpius OB association and the masses of its members we cannot yet draw firm conclusions on the estimated mass of the RX J1856.5-3754 progenitor.

We present a long-term study of the 2011 outburst of the magnetar Swift J1834.9-0846 carried out using new Chandra observations, as well as all the available Swift, RXTE, and XMM-Newton data. The last observation was performed on 2011 November 12, about 100 days after the onset of the bursting activity that had led to the discovery of the source on 2011 August 07. This long time span enabled us to refine the rotational ephemeris and observe a downturn in the decay of the X-ray flux. Assuming a broken power law for the long-term light curve, the break was at ~46 d after the outburst onset, when the decay index changed from alpha ~ 0.4 to ~4.5. The flux decreased by a factor ~2 in the first ~50 d and then by a factor ~40 until November 2011 (overall, by a factor ~70 in ~100 d). At the same time, the spectrum, which was well described by an absorbed blackbody all along the outburst, softened, the temperature dropping from ~1 to ~0.6 keV. Diffuse X-ray emission extending up to 20" from the source was clearly detected in all Chandra observations. Its spatial and spectral properties, as well as its time evolution, are consistent with a dust-scattering halo due to a single cloud located at a distance of $\approx$200 pc from Swift J1834.9-0846, which should be in turn located at a distance of ~5 kpc. Considering the time delay of the scattered photons, the same dust cloud might also be responsible for the more extended emission detected in XMM-Newton data taken in September 2011. We searched for the radio signature of Swift J1834.9-0846 at radio frequencies using the Green Bank Radio Telescope and in archival data collected at Parkes from 1998 to 2003. No evidence for radio emission was found, down to a flux density of 0.05 mJy (at 2 GHz) during the outburst and ~0.2-0.3 mJy (at 1.4 GHz) in the older data.

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra, and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3-1606 (SGR 182-1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery (July 2011), through the first stages of its outburst decay (April 2012). Our X-ray timing analysis finds the source rotating with a period of P = 8.43772016(2) s and a period derivative Pdot = 8.3(2) x 10e-14 s s-1, which entails an inferred dipolar surface magnetic field of 2.7 x 10e13 G at the equator. This measurement makes Swift J1822.3-1606 the second lowest magnetic field magnetar (after SGR 0418+5729; Rea et al. 2010). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3+1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of Bp ~ 1.5 x 10e14 G and Btor ~ 7 x 10e14 G, respectively, and if its current age is ~550 kyr (more details in Rea et al. 2012).

We report on a detailed spectral analysis of all the available XMM-Newton data of RX J1856.5-3754, the brightest and most extensively observed nearby, thermally emitting neutron star. Very small variations (~1-2%) in the single-blackbody temperature are detected, but are probably due to an instrumental effect, since they correlate with the position of the source on the detector. Restricting the analysis to a homogeneous subset of observations, with the source at the same detector position, we place strong limits on possible spectral or flux variations from March 2005 to present-day. A slightly higher temperature (kT~61.5 eV, compared to the average value kT~61 eV) was instead measured in April 2002. If this difference is not of instrumental origin, it implies a rate of variation of about 0.15 eV/yr between April 2002 and March 2005. The high-statistics spectrum from the selected observations is well fit by the sum of two blackbody models, which extrapolate to an optical flux level in agreement with the observed value.

SXP 1062, a newly discovered Be/X-ray binary in the Small Magellanic Cloud, provides the first example of a robust association with a supernova remnant (SNR). The short age estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, ~ 10^4 yr. Here we discuss possible evolutionary scenarios for SXP 1062 in the attempt to reconcile its long spin period, P=1062 s, and short age. Our results indicate that SXP 1062 may host a neutron star born with a large initial magnetic field, typically in excess of ~ 10^14 G, which then decayed to ~ 10^13 G.

A subpopulation of neutron stars (NSs), known as central compact objects (CCOs) in supernova remnants, are suspected to be low-field objects basing on $P$-$\dot P$ measurements for three of them. The birth rate of low-field NSs is probably comparable with the birth rate of normal radio pulsars. However, among compact objects in High-Mass X-ray Binaries (HMXBs) we do not see robust candidates for low-field NSs. We propose that this contradiction can be solved if magnetic fields of CCOs was buried due to strong fall-back, and then the field emerges on the time scale $10^4$-$10^5$ yrs.

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) -- which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) -- represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306. Its evolutionary state is investigated by means of deep dedicated observations obtained with XMM-Newton, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi and AGILE missions. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula, has a soft spectrum with absorption features and no magnetospheric emission. The optical counterpart is fainter than V=27 and no gamma-ray emission is significantly detected. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P=18.6 ms. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3-1606 (SGR 1822-1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery, through the first stages of its outburst decay (covering the time-span from July 2011, until end of April 2012). We also report on archival ROSAT observations which witnessed the source during its likely quiescent state, and on upper limits on Swift J1822.3-1606's radio-pulsed and optical emission during outburst, with the Green Bank Telescope (GBT) and the Gran Telescopio Canarias (GTC), respectively. Our X-ray timing analysis finds the source rotating with a period of P=8.43772016(2) s and a period derivative \dot{P}=8.3(2)x10^{-14} s s^{-1} , which entails an inferred dipolar surface magnetic field of B~2.7x10^{13} G at the equator. This measurement makes Swift J1822.3-1606 the second lowest magnetic field magnetar (after SGR 0418+5729; Rea et al. 2010). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3-1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of B_{p}~1.5x10^{14} G and B_{tor}~7x10^{14} G, respectively, and if its current age is ~550 kyr.

We present estimates of initial spin periods, $P_0$, for radio pulsars associated with supernova remnants. By using the published data on 30 objects, we were able to derive a reliable estimate for the initial spin period, assuming standard magneto-dipole spin-down (braking index n=3), in many cases. Our set of estimates is still not sufficient to infer the exact shape of the initial period distribution. However, we show that a gaussian distribution with mean and deviation $\sim 0.1$ s is consistent with our results, while flat, wide distributions and very narrow ones are disfavored.

RX J0720.4-3125 is the most peculiar object among a group of seven isolated X-ray pulsars (the so-called "Magnificent Seven"), since it shows long-term variations of its spectral and temporal properties on time scales of years. This behaviour was explained by different authors either by free precession (with a seven or fourteen years period) or possibly a glitch that occurred around $\mathrm{MJD=52866\pm73 days}$. We analysed our most recent XMM-Newton and Chandra observations in order to further monitor the behaviour of this neutron star. With the new data sets, the timing behaviour of RX J0720.4-3125 suggests a single (sudden) event (e.g. a glitch) rather than a cyclic pattern as expected by free precession. The spectral parameters changed significantly around the proposed glitch time, but more gradual variations occurred already before the (putative) event. Since $\mathrm{MJD\approx53000 days}$ the spectra indicate a very slow cooling by $\sim$2 eV over 7 years.

High magnetic fields are a distinguishing feature of neutron stars and the existence of sources (the soft gamma repeaters and the anomalous X-ray pulsars) hosting an ultra-magnetized neutron star (or magnetar) has been recognized in the past few decades. Magnetars are believed to be powered by magnetic energy and not by rotation, as with normal radio pulsars. Until recently, the radio quietness and magnetic fields typically above the quantum critical value (Bq~4.4x10^{13} G), were among the characterizing properties of magnetars. The recent discovery of radio pulsed emission from a few of them, and of a low dipolar magnetic field soft gamma repeater, weakened further the idea of a clean separation between normal pulsars and magnetars. In this Letter we show that radio emission from magnetars might be powered by rotational energy, similarly to what occurs in normal radio pulsars. The peculiar characteristics of magnetars radio emission should be traced in the complex magnetic geometry of these sources. Furthermore, we propose that magnetar radio activity or inactivity can be predicted from the knowledge of the star's rotational period, its time derivative and the quiescent X-ray luminosity.

Using a large set of XMM-Newton observations we searched for long term spectral and flux variability of the isolated neutron star RX J1856.5-3754 in the time interval from April 2002 to October 2011. This is the brightest and most extensively observed source of a small group of nearby, thermally emitting isolated neutron stars, of which at least one member (RX J0720.4-3125, Hohle et al., 2010) has shown long term variability. A detailed analysis of the data obtained with the EPIC-pn camera in the 0.15-1.2 keV energy range reveals small variations in the temperature derived with a single blackbody fit (of the order of 1% around kT^inf \sim 61 eV). Such variations are correlated with the position of the source on the detector and can be ascribed to an instrumental effect, most likely a spatial dependence of the channel to energy relation. For the sampled instrumental coordinates, we quantify this effect as variations of \sim 4% and \sim 15 eV in the gain slope and offset, respectively. Selecting only a homogeneous subset of observations, with the source imaged at the same detector position, we find no evidence for spectral or flux variations of RX J1856.5-3754 from March 2005 to present-day, with limits of Delta kT^inf < 0.5% and Delta f_X < 3% (0.15-1.2 keV), with 3sigma confidence. A slightly higher temperature (kT^inf \sim 61.5 eV, compared to kT^\inf \sim 61 eV) was instead measured in April 2002. If this difference is not of instrumental origin, it implies a rate of variation \sim -0.15 eV yr^-1 between April 2002 and March 2005. The high-statistics spectrum from the selected observations is best fitted with the sum of two blackbody models, with temperatures kT_h^inf = 62.4_{-0.4}^{+0.6} eV and kT_s^\inf = 38.9_{-2.9}^{+4.9} eV, which account for the flux seen in the optical band. No significant spectral features are detected, with upper limits of 6 eV on their equivalent width.

The newly discovered Be/X-ray binary in the Small Magellanic Cloud, SXP 1062, provides the first example of a robust association with a supernova remnant (SNR). The short age estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, tau ~ 1e4 yr. As such, it allows to test current models of magneto-rotational evolution of neutron stars ina still unexplored regime. Here we discuss possible evolutionary scenarios for SXP 1062 in the attempt to reconcile its long spin period, P=1062 s, and short age. Although several options can be considered, like an anomalously long initial period or the presence of a fossil disc, our results indicate that SXP 1062 may host a neutron star born with a large initial magnetic field, typically in excess of ~ 1e14 G, which then decayed to ~ 1e13 G.

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.

Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions in a previously developed MC code to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angle-dependence of the spectra.

SGR 0418+5729 is a transient Soft Gamma-ray Repeater which underwent a major outburst in June 2009, during which the emission of short bursts was observed. Its properties appeared quite typical of other sources of the same class until long-term X-ray monitoring failed to detect any period derivative. The present upper limit on $\dot P$ implies that the surface dipole field is $B_p\lesssim 7.5\times 10^{12}\ {\rm G}$ (Rea et al 2010), well below those measured in other Soft Gamma-ray Repeaters (SGRs) and in the Anomalous X-ray Pulsars (AXPs), a group of similar sources. Both SGRs and AXPs are currently believed to be powered by ultra-magnetized neutron stars (magnetars, $B_p\approx 10^{14}$--$10^{15}\ {\rm G}$). SGR 0418+5729 hardly seems to fit in such a picture. We show that the magneto-rotational properties of SGR 0418+5729 can be reproduced if this is an aged magnetar, $\approx 1\ {\rm Myr}$ old, which experienced substantial field decay. The large initial toroidal component of the internal field required to match the observed properties of SGR 0418+5729 ensures that crustal fractures, and hence bursting activity, can still occur at present time. The thermal spectrum observed during the outburst decay is compatible with the predictions of a resonant compton scattering model (as in other SGRs/AXPs) if the field is low and the magnetospheric twist moderate.

The X-ray source 1E 161348-5055 lies at the centre of the 2-kyr-old supernova remnant RCW 103. Owing to its 24-ks modulation, orders-of-magnitude flux variability over a few months/years, and lack of an obvious optical counterpart, 1E 161348-5055 defies assignment to any known class of X-ray sources. Starting from April 2006, Swift observed 1E 161348-5055 with its X-ray telescope for 2 ks approximately once per month. During the five years covered, the source has remained in a quiescent state, with an average observed flux of 1.7e-12 erg/cm^2/s (1-10 keV), about 20 times lower than the historical maximum attained in its 1999-2000 outburst. The long time-span of the Swift data allows us to obtain an accurate measure of the period of 1E 161348-5055 [P = 24030.42(2) s] and to derive the first upper limit on its period derivative (|dP/dt| < 1.6e-9 s/s at 3 sigma).

We present high-speed optical photometry of the soft gamma repeater SGR 0501+4516, obtained with ULTRACAM on two consecutive nights approximately 4 months after the source was discovered via its gamma-ray bursts. We detect SGR 0501+4516 at a magnitude of i' = 24.4+/-0.1. We present the first measurement of optical pulsations from an SGR, deriving a period of 5.7622+/-0.0003 s, in excellent agreement with the X-ray spin period of the neutron star. We compare the morphologies of the optical pulse profile with the X-ray and infrared pulse profiles; we find that the optical, infrared and harder X-rays share similar double-peaked morphologies, but the softer X-rays exhibit only a single-peaked morphology, indicative of a different origin. The optical pulsations appear to be in phase with the X-ray pulsations and exhibit a root-mean-square pulsed fraction of 52+/-7%, approximately a factor of two greater than in the X-rays. Our results find a natural explanation within the context of the magnetar model for SGRs.

SGR 1833-0832 was discovered on 2010 March 19 thanks to the Swift detection of a short hard X-ray burst and follow-up X-ray observations. Since then, it was repeatedly observed with Swift, Rossi X-ray Timing Explorer, and XMM-Newton. Using these data, which span about 225 days, we studied the long-term spectral and timing characteristics of SGR 1833-0832. We found evidence for diffuse emission surrounding SGR 1833-0832, which is most likely a halo produced by the scattering of the point source X-ray radiation by dust along the line of sight, and we show that the source X-ray spectrum is well described by an absorbed blackbody, with temperature kT=1.2 keV and absorbing column nH=(10.4+/-0.2)E22 cm^-2, while different or more complex models are disfavoured. The source persistent X-ray emission remained fairly constant at about 3.7E-12 erg/cm^2/s for the first 20 days after the onset of the bursting episode, then it faded by a factor 40 in the subsequent 140 days, following a power-law trend with index alpha=-0.5. We obtained a phase-coherent timing solution with the longest baseline (225 days) to date for this source which, besides period P=7.5654084(4) s and period derivative dP/dt=3.5(3)E-12 s/s, includes higher order period derivatives. We also report on our search of the counterpart to the SGR at radio frequencies using the Australia Telescope Compact Array and the Parkes radio telescope. No evidence for radio emission was found, down to flux densities of 0.9 mJy (at 1.5 GHz) and 0.09 mJy (at 1.4 GHz) for the continuum and pulsed emissions, respectively, consistently with other observations at different epochs.

X-ray observations performed with ROSAT led to the discovery of a group (seven to date) of X-ray dim and radio-silent middle-aged isolated neutron stars (a.k.a. XDINSs), which are characterised by pure blackbody spectra (kT~40-100 eV), long X-ray pulsations (P=3-12 s), and appear to be endowed with relatively high magnetic fields, (B~10d13-14 G). RBS 1774 is one of the few XDINSs with a candidate optical counterpart, which we discovered with the VLT. We performed deep observations of RBS 1774 in the R band with the VLT to disentangle a non-thermal power-law spectrum from a Rayleigh-Jeans, whose contributions are expected to be very much different in the red part of the spectrum. We did not detect the RBS 1774 candidate counterpart down to a 3 sigma limiting magnitude of R~27. The constraint on its colour, (B-R)<0.6, rules out that it is a background object, positionally coincident with the X-ray source. Our R-band upper limit is consistent with the extrapolation of the B-band flux (assuming a 3 sigma uncertainty) for a set of power-laws F_nu ~nu^alpha with spectral indeces alpha<0.07. If the optical spectrum of RBS 1774 were non-thermal, its power-law slope would be very much unlike those of all isolated neutron stars with non-thermal optical emission, suggesting that it is most likely thermal. For instance, a Rayleigh-Jeans with temperature T_O = 11 eV, for an optically emitting radius r_O=15 km and a source distance d=150 pc, would be consistent with the optical measurements. The implied low distance is compatible with the 0.04 X-ray pulsed fraction if either the star spin axis is nearly aligned with the magnetic axis or with the line of sight, or it is slightly misaligned with respect to both the magnetic axis and the line of sight by 5-10 degrees

We present the first detailed joint modelling of both the timing and spectral properties during the outburst decay of transient anomalous X-ray pulsars. We consider the two sources XTE J1810-197 and CXOU J164710.2-455216, and describe the source decline in the framework of a twisted magnetosphere model, using Monte Carlo simulations of magnetospheric scattering and mimicking localized heat deposition at the NS surface following the activity. Our results support a picture in which a limited portion of the star surface close to one of the magnetic poles is heated at the outburst onset. The subsequent evolution is driven both by the cooling/varying size of the heated cap and by a progressive untwisting of the magnetosphere.

With two consecutive outbursts recorded in four months (October 2008 and January 2009), and a possible third outburst in 2007, 1E 1547.0-5408 is one of the most active transient anomalous X-ray pulsars known so far. Thanks to extensive X-ray observations, obtained both in the quiescent and active states, 1E 1547.0-5408 represents a very promising laboratory to get insights into the outburst properties and magnetar emission mechanisms. We performed a detailed timing and spectral analysis of four Chandra, three INTEGRAL, and one XMM observations collected over a two week interval after the outburst onset in January 2009. Several Swift pointings, covering a 1.5 year interval, were also analyzed in order to monitor the decay of the X-ray flux. We compare the characteristics of the two outbursts, as well as those of the active and quiescent states. We also discuss the long-term X-ray flux history of 1E 1547.0-5408 since its first detection in 1980, and show that the source displays three flux levels: low, intermediate and high.

Soft gamma-ray repeaters and anomalous X-ray pulsars are a small (but growing) group of X-ray sources characterized by the emission of short bursts and by a large variability in their persistent flux. They are believed to be magnetars, i.e. neutron stars powered by extreme magnetic fields 1E14-1E15 G). We found evidence for a magnetar with a low magnetic field, SGR 0418+5729, recently detected after it emitted bursts similar to those of soft gamma-ray repeaters. New X-ray observations show that its dipolar magnetic field cannot be greater than 8E12 G, well in the range of ordinary radio pulsars, implying that a high surface dipolar magnetic field is not necessarily required for magnetar-like activity. The magnetar population may thus include objects with a wider range of magnetic-field strengths, ages and evolutionary stages than observed so far.

Soft gamma repeaters and anomalous x-ray pulsars form a rapidly increasing group of x-ray sources exhibiting sporadic emission of short bursts. They are believed to be magnetars, i.e. neutron stars powered by extreme magnetic fields, B~10^{14}-10^{15} Gauss. We report on a soft gamma repeater with low magnetic field, SGR 0418+5729, recently detected after it emitted bursts similar to those of magnetars. X-ray observations show that its dipolar magnetic field cannot be greater than 7.5x10^{12} Gauss, well in the range of ordinary radio pulsars, implying that a high surface dipolar magnetic field is not necessarily required for magnetar-like activity. The magnetar population may thus include objects with a wider range of B-field strengths, ages and evolutionary stages than observed so far.

We report on a multi-wavelength study of the compact object candidate 1RXS J141256.0+792204 (Calvera). Calvera was observed in the X-rays with XMM/EPIC twice for a total exposure time of ~50 ks. The source spectrum is thermal and well reproduced by a two component model composed of either two hydrogen atmosphere models, or two blackbodies (kT_1~ 55/150 eV, kT_2~ 80/250 eV, respectively, as measured at infinity). Evidence was found for an absorption feature at ~0.65 keV; no power-law high-energy tail is statistically required. Using pn and MOS data we discovered pulsations in the X-ray emission at a period P=59.2 ms. The detection is highly significant (> 11 sigma), and unambiguously confirms the neutron star nature of Calvera. The pulse profile is nearly sinusoidal, with a pulsed fraction of ~18%. We looked for the timing signature of Calvera in the Fermi Large Area Telescope (LAT) database and found a significant (~5 sigma) pulsed signal at a period coincident with the X-ray value. The gamma-ray timing analysis yielded a tight upper limit on the period derivative, dP/dt < 5E-18 s/s (dE_rot/dt <1E33 erg/s, B<5E10 G for magneto- dipolar spin-down). Radio searches at 1.36 GHz with the 100-m Effelsberg radio telescope yielded negative results, with a deep upper limit on the pulsed flux of 0.05 mJy. Diffuse, soft (< 1 keV) X-ray emission about 13' west of the Calvera position is present both in our pointed observations and in archive ROSAT all-sky survey images, but is unlikely associated with the X-ray pulsar. Its spectrum is compatible with an old supernova remnant (SNR); no evidence for diffuse emission in the radio and optical bands was found. The most likely interpretations are that Calvera is either a central compact object escaped from a SNR or a mildly recycled pulsar; in both cases the source would be the first ever member of the class detected at gamma-ray energies.

Since the last phase coherent timing solution of the nearby radio-quiet isolated neutron star RX J0720.4-3125 six new XMM-Newton and three Chandra observations were carried out. The phase coherent timing solutions from previous authors were performed without restricting to a fixed energy band. However, we recently showed that the phase residuals are energy dependent, and thus phase coherent solutions must be computed referring always to the same energy band. We updated the phase coherent timing solution for RX J0720.4-3125 by including the recent XMM-Newton EPIC-pn, MOS1, MOS2 and Chandra ACIS data in the energy range 400-1000~eV. Altogether these observations cover a time span of almost 10~yrs. A further timing solution was obtained including the ROSAT pointed data. In this case, observations cover a time span of $\approx$16~yrs. To illustrate the timing differences between the soft band (120-400~eV) and the hard band (400-1000~eV) a timing solution for the soft band is also presented and the results are verified using a $\mathrm{Z_{n}^{2}}$ test. In contrast to previous work, we obtain almost identical solutions whether or not we include the ROSAT or Chandra data. Thanks to the restriction to the hard band, the data points from EPIC-pn are in better agreement with those from MOS1, MOS2 and Chandra than in previous works. In general the phase residuals are still large and vary with time. In particular, the latest XMM-Newton and Chandra data show that the phase residuals have attained relatively large and negative values.

There is a general consensus about the fact that the magnetar scenario provides a convincing explanation for several of the observed properties of the Anomalous X-ray Pulsars and the Soft Gamma Repeaters. However, the origin of the emission observed at low energies is still an open issue. We present a quantitative model for the emission in the optical/infrared band produced by curvature radiation from magnetospheric charges, and compare results with current magnetars observations.

In this paper, we discuss our first attempts to model the broadband persistent emission of magnetars within a self consistent, physical scenario. We present the predictions of a synthetic model that we calculated with a new Monte Carlo 3-D radiative code. The basic idea is that soft thermal photons (e.g. emitted by the star surface) can experience resonant cyclotron upscattering by a population of relativistic electrons threated in the twisted magnetosphere. Our code is specifically tailored to work in the ultra-magnetized regime; polarization and QED effects are consistently accounted for, as well different configurations for the magnetosphere. We discuss the predicted spectral properties in the 0.1-1000 keV range, the polarization properties, and we present the model application to a sample of magnetars soft X-ray spectra.

The high time resolution observations of the X-ray sky hold the key to a number of diagnostics of fundamental physics, some of which are unaccessible to other types of investigations, such as those based on imaging and spectroscopy. Revealing strong gravitational field effects, measuring the mass and spin of black holes and the equation of state of ultradense matter are among the goals of such observations. At present prospects for future, non-focused X-ray timing experiments following the exciting age of RXTE/PCA are uncertain. Technological limitations are unavoidably faced in the conception and development of experiments with effective area of several square meters, as needed in order to meet the scientific requirements. We are developing large-area monolithic Silicon Drift Detectors offering high time and energy resolution at room temperature, which require modest resources and operation complexity (e.g., read-out) per unit area. Based on the properties of the detector and read-out electronics that we measured in the lab, we developed a realistic concept for a very large effective area mission devoted to X-ray timing in the 2-30 keV energy range. We show that effective areas in the range of 10-15 square meters are within reach, by using a conventional spacecraft platform and launcher of the small-medium class.

Anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) are two small classes of X-ray sources strongly suspected to host a magnetar, i.e. an ultra-magnetized neutron star with $B\approx 10^14-10^15 G. Many SGRs/AXPs are known to be variable, and recently the existence of genuinely "transient" magnetars was discovered. Here we present a comprehensive study of the pulse profile and spectral evolution of the two transient AXPs (TAXPs) XTE J1810-197 and CXOU J164710.2-455216. Our analysis was carried out in the framework of the twisted magnetosphere model for magnetar emission. Starting from 3D Monte Carlo simulations of the emerging spectrum, we produced a large database of synthetic pulse profiles which was fitted to observed lightcurves in different spectral bands and at different epochs. This allowed us to derive the physical parameters of the model and their evolution with time, together with the geometry of the two sources, i.e. the inclination of the line-of-sight and of the magnetic axis with respect to the rotation axis. We then fitted the (phase-averaged) spectra of the two TAXPs at different epochs using a model similar to that used to calculate the pulse profiles ntzang in XSPEC) freezing all parameters to the values obtained from the timing analysis, and leaving only the normalization free to vary. This provided acceptable fits to XMM-Newton data in all the observations we analyzed. Our results support a picture in which a limited portion of the star surface close to one of the magnetic poles is heated at the outburst onset. The subsequent evolution is driven both by the cooling/varying size of the heated cap and by a progressive untwisting of the magnetosphere.

We report on the detailed study of the 2008 October outburst from the anomalous X-ray pulsar (AXP) 1E 1547.0-5408 discovered through the Swift/Burst Alert Telescope (BAT) detection of SGR-like short X-ray bursts on 2008 October 3. The Swift/X-ray Telescope (XRT) started observing the source after less than 100 s since the BAT trigger, when the flux (about 6E-11 erg/cm^2/s in the 2-10 keV range) was >50 times higher than its quiescent level. Swift monitored the outbursting activity of 1E 1547.0-5408 on a daily basis for approximately three weeks. This strategy allowed us to find a phase-coherent solution for the source pulsations after the burst, which, besides period and period derivative, requires a positive Period second derivative term (spin-down increase). The time evolution of the pulse shape is complex and variable, with the pulsed fraction increasing from 20% to 50% within the Swift observational window. The XRT spectra can be fitted well by means of a single component, either a power-law (PL) or a blackbody (BB). During the very initial phases of the outburst the spectrum is hard, with a PL photon index about 2 (or kT about 1.4 keV) which steepens to about 4 (or kT about 0.8 keV) within one day from the BAT trigger, though the two components are likely present simultaneously during the first day spectra. An INTEGRAL observation carried out five days after the trigger provided an upper limit of about 2E-11 erg/cm^2/s to the emission of 1E 1547.0-5408 in the 18-60 keV band.

Context: Isolated cooling neutron stars with thermal X-ray emission remain rarely detected objects despite many searches investigating the ROSAT data. Aims: We simulate the population of close-by young cooling neutron stars to explain the current observational results. Given the inhomogeneity of the neutron star distribution on the sky it is particularly interesting to identify promising sky regions with respect to on-going and future searches. Methods: Applying a population synthesis model the inhomogeneity of the progenitor distribution and the inhomogeneity of the X-ray absorbing interstellar medium are considered for the first time. The total number of observable neutron stars is derived with respect to ROSAT count rates. In addition, we present sky maps of neutron star locations and discuss age and distance distributions of the simulated neutron stars. Implications for future searches are discussed. Results: With our advanced model we can successfully explain the observed logN - logS distribution of close-by neutron stars. Cooling neutron stars will be most abundant in the directions of rich OB associations. New candidates are expected to be identified behind the Gould Belt, in particular in the Cygnus-Cepheus region. They are expected to be on average younger and then hotter than the known population of isolated cooling neutron stars. In addition, we propose to use data on runaway stars to search for more radio-quiet cooling neutron stars.

Emission of two short hard X-ray bursts on 2009 June 5 disclosed the existence of a new soft gamma-ray repeater, now catalogued as SGR 0418+5729. After a few days, X-ray pulsations at a period of 9.1 s were discovered in its persistent emission. SGR 0418+5729 was monitored almost since its discovery with the Rossi X-ray Timing Explorer (2-10 keV energy range) and observed many times with Swift (0.2-10 keV). The source persistent X-ray emission faded by a factor 10 in about 160 days, with a steepening in the decay about 19 days after the activation. The X-ray spectrum is well described by a simple absorbed blackbody, with a temperature decreasing in time. A phase-coherent timing solution over the 160 day time span yielded no evidence for any significant evolution of the spin period, implying a 3-sigma upper limit of 1.1E-13 s/s on the period derivative and of 3E+13 G on the surface dipole magnetic field. Phase-resolved spectroscopy provided evidence for a significant variation of the spectrum as a function of the stellar rotation, pointing to the presence of two emitting caps, one of which became hotter during the outburst. Finally, a deep observation of the field of SGR 0418+5729 with the new Gran Telescopio Canarias 10.4-m telescope allowed us to set an upper limit on the source optical flux of i'>25.1 mag, corresponding to an X-ray-to-optical flux ratio exceeding 10000, consistent with the characteristics of other magnetars.

On 2009 January 22 numerous strong bursts were detected from the anomalous X-ray pulsar 1E 1547.0-5408. Swift/XRT and XMM-Newton/EPIC observations carried out in the following two weeks led to the discovery of three X-ray rings centered on this source. The ring radii increased with time following the expansion law expected for a short impulse of X-rays scattered by three dust clouds. Assuming different models for the dust composition and grain size distribution, we fit the intensity decay of each ring as a function of time at different energies, obtaining tight constrains on the distance of the X-ray source. Although the distance strongly depends on the adopted dust model, we find that some models are incompatible with our X-ray data, restricting to 4-8 kpc the range of possible distances for 1E 1547.0-5408. The best-fitting dust model provides a source distance of 3.91 +/- 0.07 kpc, which is compatible with the proposed association with the supernova remnant G 327.24-0.13, and implies distances of 2.2 kpc, 2.6 kpc and 3.4 kpc for the dust clouds, in good agreement with the dust distribution inferred by CO line observations towards 1E 1547.0-5408. However, dust distances in agreement with CO data are also obtained for a set of similarly well-fitting models that imply a source distance of about 5 kpc. A distance of about 4-5 kpc is also favored by the fact that these dust models are already known to provide good fits to the dust-scattering halos of bright X-ray binaries.

Aims. Neutron stars (NSs) produced in the Milky Way are supposedly ten to the eighth - ten to the ninth, of which only $\sim 2 \times 10^{3}$ are observed. Constraining the phase space distribution of NSs may help to characterize the yet undetected population of stellar remnants. Methods. We perform Monte Carlo simulations of NS orbits, under different assumptions concerning the Galactic potential and the distribution of progenitors and birth velocities. We study the resulting phase space distributions, focusing on the statistical properties of the NS populations in the disk and in the solar neighbourhood. Results. It is shown that $\sim 80$ percent of NSs are in bound orbits. The fraction of NSs located in a disk of radius 20 kpc and width 0.4 kpc is $\lesssim 20$ percent. Therefore the majority of NSs populate the halo. Fits for the surface density of the disk, the distribution of heights on the Galactic plane and the velocity distribution of the disk, are given. We also provide sky maps of the projected number density in heliocentric Galactic coordinates (l, b). Our results are compared with previous ones reported in the literature. Conclusions. Obvious applications of our modelling are in the revisiting of accretion luminosities of old isolated NSs, the issue of the observability of the nearest NS and the NS optical depth for microlensing events. These will be the scope of further studies.

We report on recent results obtained thanks to Target of Opportunity observations of the two galactic sources SGR 1627-41 and 1E 1547-5408. These two transient sources present several similarities which support the interpretation of Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters as a single class of strongly magnetized neutron stars.

SGR 0526-66 was the first soft gamma-ray repeater (SGR) from which a giant flare was detected in March 1979, suggesting the existence of magnetars, i.e. neutron stars powered by the decay of their extremely strong magnetic field. Since then, very little information has been obtained on this object, mainly because it has been burst-inactive since 1983 and the study of its persistent X-ray emission has been hampered by its large distance and its location in a X-ray bright supernova remnant in the Large Magellanic Cloud. Here we report on a comprehensive analysis of all the available XMM-Newton observations of SGR 0526-66. In particular, thanks to a deep observation taken in 2007, we measured its pulsation period (P = 8.0544 +/- 0.0002 s) 6 years after its latest detection by Chandra. This allowed us to detect for the first time a significant reduction of its spin-down rate. From a comparison with two shorter XMM-Newton observations performed in 2000 and 2001, we found no significant changes in the spectrum, which is well modelled by an absorbed power-law with nH = 4.6E+21 cm^-2 and photon index = 3.27. The high luminosity (about 4E+35 erg/s, in the 1-10 keV energy band) still observed about 25 years after the latest detection of bursting activity places SGR 0526-66 in the group of bright and persistent magnetar candidates.

Using Chandra data taken on 2008 June, we detected pulsations at 2.59439(4) s in the soft gamma-ray repeater SGR 1627-41. This is the second measurement of the source spin period and allows us to derive for the first time a long-term spin-down rate of (1.9 +/- 0.4)E-11 s/s. From this value we infer for SGR 1627-41 a characteristic age of 2.2 kyr, a spin-down luminosity of 4E+34 erg/s (one of the highest among sources of the same class), and a surface dipole magnetic field strength of 2E+14 G. These properties confirm the magnetar nature of SGR 1627-41; however, they should be considered with caution since they were derived on the basis of a period derivative measurement made using two epochs only and magnetar spin-down rates are generally highly variable. The pulse profile, double-peaked and with a pulsed fraction of (13 +/- 2)% in the 2-10 keV range, closely resembles that observed by XMM-Newton in 2008 September. Having for the first time a timing model for this SGR, we also searched for a pulsed signal in archival radio data collected with the Parkes radio telescope nine months after the previous X-ray outburst. No evidence for radio pulsations was found, down to a luminosity level 10-20 times fainter (for a 10% duty cycle and a distance of 11 kpc) than the peak luminosity shown by the known radio magnetars.

We have carried out a search for radio emission at 820 MHz from six X-ray dim isolated neutron stars with the Robert C. Byrd Green Bank Radio Telescope. No transient or pulsed emission was found using fast folding, fast Fourier transform, and single-pulse searches. The corresponding flux limits are about 0.01 mJy for pulsed emission, depending on the integration time for the particular source and assuming a duty cycle of 2%, and 20 mJy for single dispersed pulses. These are the most sensitive limits to date on radio emission from X-ray dim isolated neutron stars. There is no evidence for isolated radio pulses, as seen in a class of neutron stars known as rotating radio transients. Our results imply that either the radio luminosities of these objects are lower than those of any known radio pulsars, or they could simply be long-period nearby radio pulsars with high magnetic fields beaming away from the Earth. To test the latter possibility, we would need around 40 similar sources to provide a 1 sigma probability of at least one of them beaming toward us. We also give a detailed description of our implementation of the Fast Folding Algorithm.

X-ray observations performed with ROSAT led to the discovery of 7 radio-silent isolated neutron stars (INSs) which are detected only through the relatively dim and purely thermal X-ray emission from the cooling star surface. A few of these INSs (a.k.a. X-ray Dim INSs, or XDINSs) have been also detected at optical wavelengths where they seem to feature thermal spectra. Optical studies of XDINSs thus play a crucial role in mapping the temperature distribution on the neutron star surface and in investigating the existence of an atmosphere around the neutron star. The aim of this work is to investigate the optical identification of the XDINS RX J0420.0-5022, tentatively proposed by Haberl et al. (2004) based on VLT observations. We re-analysed the original observations of Haberl et al. to assess the detection significance of the proposed counterpart and we performed deeper VLT observations aiming at a higher confidence detection. With a ~2 sigma detection significance and a re-computed flux of B=27.52+/-0.61, we can not rule out that the proposed counterpart was spurious and produced by the halo of a very bright nearby star. While we could not detect the proposed counterpart in our deeper VLT observations, we found evidence for a marginally significant (~3.9 sigma) detection of a similarly faint object (B= 27.5+/-0.3), ~ 0.5" north of it and coincident with the updated Chandra position of RX J0420.0-5022. Interestingly, the angular separation is consistent with the upper limit on the RX J0420.0-5022 proper motion (Motch et al. 2009), which suggests that we might have actually detected the Haberl et al. proposed counterpart. From the flux of the putative RX J0420.0-5022 counterpart we can rule out a > 7 optical excess with respect to the extrapolation of the XMM-Newton spectrum.

Within the magnetar scenario, the "twisted magnetosphere" model appears very promising in explaining the persistent X-ray emission from the Soft Gamma Repeaters and the Anomalous X-ray Pulsars (SGRs and AXPs). In the first two papers of the series, we have presented a 3D Monte Carlo code for solving radiation transport as soft, thermal photons emitted by the star surface are resonantly upscattered by the magnetospheric particles. A spectral model archive has been generated and implemented in XSPEC. Here we report on the systematic application of our spectral model to different XMM-Newton and Integral observations of SGRs and AXPs. We find that the synthetic spectra provide a very good fit to the data for the nearly all the source (and source states) we have analyzed.

In the past, the isolated, radio-quiet neutron star RX J0720.4-3125 showed variations in the spectral parameters (apparent radius, temperature of the emitting area and equivalent width of the absorption feature) seen in the X-ray spectra, not only during the spin period of 8.39s, but also over time scales of years. New X-ray observations of RX J0720.4-3125 with XMM Newton extend the coverage to about 7.5 years with the latest pointing performed in November 2007. Out of a total of fourteen available EPIC-pn datasets, eleven have been obtained with an identical instrumental setup (full frame read-out mode with thin filter), and are best suited for a comparative investigations of the spectral and timing properties of this enigmatic X-ray pulsar. We analysed the new XMM Newton observations together with archival data in order to follow the spectral and temporal evolution of RX J0720.4-3125 All XMM-Newton data were reduced with the standard XMM-SAS software package. A systematic and consistent data reduction of all these observations was emphasised in order to reduce systematic errors as far as possible. We investigate the phase residuals derived from data from different energy bands using different timing solutions for the spin period evolution and confirm the phase lag between hard and soft photons. The phase shift in the X-ray pulses between hard and soft photons varies with time and changes sign around MJD=52800 days, regardless of the chosen timing solution. The phase residuals[abridge]

We report here on the outburst onset and evolution of the new Soft Gamma Repeater SGR 0501+4516. We monitored the new SGR with XMM-Newton starting on 2008 August 23, one day after the source became burst-active, and continuing with 4 more observations, with the last one on 2008 September 30. Combining the data with the Swift-XRT and Suzaku data, we modelled the outburst decay over 160 days, and we found that the source flux decreased exponentially with a timescale of t_c=23.8 days. In the first XMM-Newton observation a large number of short X-ray bursts were observed, the rate of which decayed drastically in the following observations. We found large changes in the spectral and timing behavior of the source during the outburst, with softening emission as the flux decayed, and the non-thermal soft X-ray spectral component fading faster than the thermal one. Almost simultaneously to our XMM-Newton observations (on 2008 August 29 and September 2), we observed the source in the hard X-ray range with INTEGRAL, which clearly detected the source up to ~100keV in the first pointing, while giving only upper limits during the second pointing, discovering a variable hard X-ray component fading in less than 10 days after the bursting activation. We performed a phase-coherent X-ray timing analysis over about 160 days starting with the burst activation and found evidence of a strong second derivative period component (\ddot{P} = -1.6(4)x10^{-19} s/s^{-2}). Thanks to the phase-connection, we were able to study the the phase-resolved spectral evolution of SGR 0501+4516 in great detail. We also report on the ROSAT quiescent source data, taken back in 1992 when the source exhibits a flux ~80 times lower than that measured during the outburst, and a rather soft, thermal spectrum.

In January 2009, multiple short bursts of soft gamma-rays were detected from the direction of the anomalous X-ray pulsar 1E 1547.0-5408 by different satellites. Here we report on the observations obtained with the INTEGRAL SPI-ACS detector during the period with the strongest bursting activity. More than 200 bursts were detected at energies above 80 keV in a few hours on January 22. Among these, two remarkably bright events showed pulsating tails lasting several seconds and modulated at the 2.1 s spin period of 1E 1547.0-5408. The energy released in the brightest of these bursts was of a few 10^43 erg, for an assumed distance of 10 kpc. This is smaller than that of the three giant flares seen from soft gamma-ray repeaters, but higher than that of typical bursts from soft gamma-ray repeaters and anomalous X-ray pulsars.

We derive a strong upper bound on the amount of Primordial Black Holes (PBHs) that can still be present in the Universe. Gravitational capture of PBHs by the Milky Way stars during their formation and subsequent accretion would produce a dramatic depletion of disk stars and especially of white dwarfs, unless the average cosmic density and mass of PBHs are severely constrained. Our finding also helps to discriminate among the various production mechanisms of PBHs. Moreover, we show that a star becomes overluminous before its disappearance into a PBH for a time span independent of its mass, thereby providing a characteristic observational signature of the considered scenario. We stress that our result allows for the existence of stellar-mass black holes in a mass range that is forbidden by standard stellar evolution.

We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X and soft gamma-ray data of four magnetars, including anomalous X-ray pulsars, and soft gamma repeaters. In this scenario, non-thermal magnetar spectra in the soft X-rays result from resonant cyclotron scattering of the thermal surface emission by hot magnetospheric plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while an additional component is still required to model the hard X-ray persistent magnetar emission recently discovered by INTEGRAL. The latter is an important component in terms of magnetars' luminosity, and cannot be neglected when modelling the soft X-ray part of the spectrum.

Over the last decade, X-ray observations unveiled the existence of several classes of isolated neutron stars (INSs) which are radio-quiet or exhibit radio emission with properties much at variance with those of ordinary radio pulsars. The identification of new sources is crucial in order to understand the relations among the different classes and to compare observational constraints with theoretical expectations. A recent analysis of the 2XMMp catalogue provided less than 30 new thermally emitting INS candidates. Among these, the source 2XMM J104608.7-594306 appears particularly interesting because of the softness of its X-ray spectrum and of the present upper limits in the optical, which imply a logarithmic X-ray-to-optical flux ratio greater than 3.1, corrected for absorption. We present the X-ray and optical properties of 2XMM J104608.7-594306 and discuss its nature in the light of two possible scenarios invoked to explain the X-ray thermal emission from INSs: the release of residual heat in a cooling neutron star, as in the seven radio-quiet ROSAT-discovered INSs, and accretion from the interstellar medium. We find that the present observational picture of 2XMM J104608.7-594306 is consistent with a distant cooling INS with properties in agreement with the most up-to-date expectations of population synthesis models: it is fainter, hotter and more absorbed than the seven ROSAT sources and possibly located in the Carina Nebula, a region likely to harbour unidentified cooling neutron stars. The accretion scenario, although not entirely ruled out by observations, would require a very slow (~10 km/s) INS accreting at the Bondi-Hoyle rate.

There is an increasing theoretical and observational evidence that the external magnetic field of magnetars may contain a toroidal component, likely of the same order of the poloidal one. Such "twisted magnetospheres" are threaded by currents flowing along the closed field lines which can efficiently interact with soft thermal photons via resonant cyclotron scatterings (RCS). Actually, RCS spectral models proved quite successful in explaining the persistent ~1-10 keV emission from the magnetar candidates, the soft gamma-ray repeaters (SGRs) and the anomalous X-ray pulsars (AXPs). Moreover, it has been proposed that, in presence of highly relativistic electrons, the same process can give rise to the observed hard X-ray spectral tails extending up to ~200 keV. Spectral calculations have been restricted up to now to the case of a globally twisted dipolar magnetosphere, although there are indications that the twist may be confined only to a portion of the magnetosphere, and/or that the large scale field is more complex than a simple dipole. In this paper we investigate multipolar, force-free magnetospheres of ultra-magnetized neutron stars. We first discuss a general method to generate multipolar solutions of the Grad- Schluter-Shafranov equation, and analyze in detail dipolar, quadrupolar and octupolar fields. The spectra and lightcurves for these multipolar, globally twisted fields are then computed using a Monte Carlo code and compared with those of a purely dipolar configuration. Finally the phase-resolved spectra and energy-dependent lightcurves obtained with a simple model of a locally sheared field are confronted with the INTEGRAL observations of the AXPs 1RXS J1708-4009 and 4U 0142+61. Results support a picture in which the field in these two sources is not globally twisted.

We present the first Suzaku observation of the new Soft Gamma Repeater SGR 0501+4516, performed on 2008 August 26, four days after the onset of bursting activity of this new member of the magnetar family. The soft X-ray persistent emission was detected with the X-ray Imaging Spectrometer (XIS) at a 0.5-10 keV flux of 3.8E-11 erg/s/cm2, with a spectrum well fitted by an absorbed blackbody plus power-law model. The source pulsation was confirmed at a period of 5.762072+/-0.000002 s, and 32 X-ray bursts were detected by the XIS, four of which were also detected at higher energies by the Hard X-ray Detector (HXD). The strongest burst, which occurred at 03:16:16.9 (UTC), was so bright that it caused instrumental saturation, but its precursor phase, lasting for about 200 ms, was detected successfully over the 0.5-200 keV range, with a fluence of ~2.1E-7 erg/cm2 and a peak intensity of about 89 Crab. The entire burst fluence is estimated to be ~50 times higher. The precursor spectrum was very hard, and well modeled by a combination of two blackbodies. We discuss the bursting activity and X/gamma-ray properties of this newly discovered Soft Gamma Repeater in comparison with other members of the class.

Despite about a decade of observations, very little is known about the optical and infrared (IR) emission properties of the Soft Gamma-ray Repeaters (SGRs) and of the Anomalous X-ray Pulsars (AXPs), the magnetar candidates, and about the physical processes which drive their emission at these wavelengths. This is mainly due to the limited number of identifications achieved so far, five in total, and to the sparse spectral coverage obtained from multi-band optical/IR photometry. Aim of this work is to search for a likely candidate counterpart to the recently discovered transient radio AXP 1E 1547.0-5408. We performed the first deep near-IR (NIR) observations (Ks band) of 1E 1547.0-5408 with the VLT on three nights (July 8th, 12th, and August 17th), after the X-ray source rebrightening and during the subsequent decay reported around June 2007. We detected four objects within, or close to, the 3sigma radio position of 1E 1547.0-5408. The faintest of them (object 1) has a brightness K = 20.27 +/- 0.05, which would yield an unabsorbed X-ray to NIR flux ratio FX/FKs ~ 800 for 1E 1547.0-5408, i.e. on average lower than those derived for other magnetars. The non-detection of object 1 on the nights of July 8th and August 17th only allowed us to set an upper limit of Delta Ks ~ 0.2 on its NIR variability, which prevented us to search for correlations with the radio or X-ray flux. We detected no other object at the radio position down to a limit of Ks ~ 21.7 (at 5sigma), computed in our deepest VLT image (July 12th). From our observations we can not confidently propose a NIR counterpart to 1E 1547.0-5408. More NIR observations of object 1, e.g. to determine its colors and to monitor variability, would be conclusive to determine whether or not it can be considered a plausible candidate.

XTE J1810-197 is the first transient Anomalous X-ray Pulsar ever discovered. Its highly variable X-ray flux allowed us to study the timing and spectral emission properties of a magnetar candidate over a flux range of about two orders of magnitude. We analyzed nine XMM-Newton observations of XTE J1810-197 collected over a four years baseline (September 2003 - September 2007). EPIC PN and MOS data were reduced and used for detailed timing and spectral analysis. Pulse phase spectroscopic studies were also carried out for observations with sufficiently high signal to noise. We find that: (i) a three blackbodies model reproduces the spectral properties of the source over the entire outburst statistically better than the two blackbodies model previously used in the literature, (ii) the coldest blackbody is consistent with the thermal emission from the whole surface, and has temperature and radius similar to those inferred from ROSAT observations before the outburst onset, (iii) there is a spectral feature around 1.1 keV during six consecutive observations (since March 2005); if due to proton resonant cyclotron scattering, it would imply a magnetic field of around 2E14 G. This is in a very good agreement with the value from the spin period measurements.

After nearly a decade of quiescence, the soft gamma-ray repeater SGR 1627-41 reactivated on 2008 May 28 with a bursting episode followed by a slowly decaying enhancement of its persistent emission. To search for the still unknown spin period of this SGR taking advantage of its high flux state, we performed on 2008 September 27-28 a 120 ks long X-ray observation with the XMM-Newton satellite. Pulsations with P = 2.594578(6) s were detected at a higher than 6-sigma confidence level, with a double-peaked pulse profile. The pulsed fraction in the 2-12 keV range is 19% +/- 3% and 24% +/- 3% for the fundamental and the second harmonic, respectively. The observed 2-10 keV flux is 3.4E-13 erg/cm^2/s, still a factor of ~ 5 above the quiescent pre-burst-activation level, and the spectrum is well fitted by an absorbed power law plus blackbody model (photon index Gamma = 0.6, blackbody temperature kT = 0.5 keV, and absorption nH = 1.2E+23 cm^-2). We also detected a shell of diffuse soft X-ray emission which is likely associated with the young supernova remnant G337.0-0.1.

In May 2008 the soft gamma-ray repeater SGR 1627-41 resumed its bursting activity after nearly a decade of quiescence. After detection of a bright burst, Swift pointed its X-ray telescope in the direction of the source in less than five hours and followed it for over five weeks. In this paper we present an analysis of the data from these Swift observations and an XMM-Newton one performed when SGR 1627-41 was still in a quiescent state. The analysis of the bursts detected with Swift/BAT shows that their temporal and spectral properties are similar to those found in previous observations of SGR 1627-41 and other soft gamma-ray repeaters. The maximum peak luminosity of the bursts was about 2E+41 erg/s. Our data show that the outburst was accompanied by a fast flux enhancement and by a hardening of the spectrum with respect to the persistent emission.

High-energy observations have unveiled peculiar classes of isolated neutron stars which, at variance with radio pulsars, are mostly radio silent and not powered by the star rotation. Among these objects are the magnetars, hyper-magnetized neutron stars characterized by transient X-ray/gamma-ray emission, and neutron stars with purely thermal, and in most cases stationary, X-ray emission (a.k.a., X-ray dim isolated neutron stars or XDINSs). While apparently dissimilar in their high-energy behavior and age, both magnetars and XDINSs have similar periods and unusually high magnetic fields. This suggests a tantalizing scenario where the former evolve into the latter.Discovering so far uninvestigated similarities between the multi-wavelength properties of these two classes would be a further step forward to establish an evolutionary scenario. A most promising channels is the near infrared (NIR) one, where magnetars are characterized by a distinctive spectral flattening with respect to the extrapolation of the soft X-ray spectrum.We observed the two XDINSs RX J0420.0-5022 and RX J1856.5-3754 with the Multi-Conjugate Adaptive Optics Demonstrator (MAD) at the Very Large Telescope (VLT) as part of the instrument guaranteed time observations program, to search for their NIR counterparts. Both RX J1856.5-3754 and RX J0420.0-5022 were not detected down to K_s ~20 and Ks ~21.5, respectively. In order to constrain the relation between XDINSs and magnetars it would be of importance to perform deeper NIR observations. A good candidate is 1RXS J214303.7+065419 which is the XDINS with the highest inferred magnetic field.

Recent models of spectral formation in magnetars called renewed attention on electron-photon scattering in the presence of ultra-strong magnetic fields. Investigations presented so far mainly focussed on mildly relativistic particles and magnetic scattering was treated in the non-relativistic (Thomson) limit. This allows for consistent spectral calculations up to a few tens of keVs, but becomes inadequate in modelling the hard tails (<200 keV) detected by INTEGRAL from magnetar sources. In this paper, the second in a series devoted to model the X-/soft gamma-ray persistent spectrum of magnetar candidates, we present explicit, relatively simple expressions for the magnetic Compton cross-section at resonance which account for Landau-Raman scattering up to the second Landau level. No assumption is made on the magnetic field strength. We find that sensible departures from the Thomson regime can bealready present at B ~5E12 G. The form of the magnetic cross section we derived can be easily implemented in Monte Carlo transfer codes and a direct application to magnetar spectral calculations will be presented in a forthcoming study.

(Abridged) We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X-ray data of ten magnetars, including canonical and transient anomalous X-ray pulsars (AXPs), and soft gamma repeaters (SGRs). In this scenario, non-thermal magnetar spectra in the soft X-rays (i.e. below ~10 keV) result from resonant cyclotron scattering of the thermal surface emission by hot magnetospheric plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while using the same number of free parameters than the commonly used empirical blackbody plus power-law model. However, while the RCS model can alone reproduce the soft X-ray spectra of AXPs, the much harder spectra of SGRs below ~10 keV, requires the addition of a power-law component (the latter being the same component responsible for their hard X-ray emission). Although this model in its present form does not explain the hard X-ray emission of a few of these sources, we took this further component into account in our modeling not to overlook their contribution in the ~4-10 keV band. We find that the entire class of sources is characterized by magnetospheric plasma with a density which, at resonant radius, is about 3 orders of magnitudes higher than n_{GJ}, the Goldreich-Julian electron density. The inferred values of the intervening hydrogen column densities, are also in better agreement with more recent estimates inferred from the fit of single X-ray edges. For the entire sample of observations, we find indications for a correlation between the scattering depth and the electron thermal velocity, and the field strength. Moreover, in most transient anomalous X-ray pulsars the outburst state is characterized by a relatively high surface temperature which cools down during the decay.

Multiwavelength studies of the seven identified X-ray dim isolated neutron stars (XDINSs) offer a unique opportunity to investigate their surface thermal and magnetic structure and the matter-radiation interaction in presence of strong gravitational and magnetic fields. As a part of an ongoing campaign aimed at a complete identification and spectral characterization of XDINSs in the optical band, we performed deep imaging with the ESO Very Large Telescope (VLT) of the field of the XDINS RBS1774 (1RXS J214303.7 +065419). The recently upgraded FORS1 instrument mounted on the VLT provided the very first detection of a candidate optical counterpart in the B band. The identification is based on a very good positional coincidence with the X-ray source (chance probability ~2E-3). The source has B=27.4 +/- 0.2 (1 sigma confidence level), and the optical flux exceeds the extrapolation of the X-ray blackbody at optical wavelengths by a factor ~35 (+/- 20 at 3sigma confidence level). This is barely compatible with thermal emission from the neutron star surface, unless the source distance is d~200-300 pc, and the star is an almost aligned rotator or its spin axis is nearly aligned with the line of sight. At the same time, such a large optical excess appears difficult to reconcile with rotation-powered magnetospheric emission, unless the source has an extremely large optical emission efficiency. The implications and possible similarities with the optical spectra of other isolated NSs are discussed.

The anomalous X-ray pulsars and soft gamma-repeaters are peculiar high-energy sources believed to host a magnetar, i.e. an ultra-magnetized neutron star. Their persistent, soft X-ray emission (~1-10 keV)is usually modeled by the superposition of a blackbody and a power-law tail. It has been suggested that this spectrum forms as the thermal photons emitted by the star surface traverse the magnetosphere. Magnetar magnetospheres are likely different from those of ordinary radio-pulsars, since the external magnetic field may acquire a toroidal component as a consequence of the deformation of the star crust induced by the super-strong interior field. In turn, the magnetosphere will be permeated by currents that can boost primary photons through repeated scatterings. Here we present 3D Monte Carlo simulations of photon propagation in a twisted magnetosphere. Our model is based on a simplified treatment of the charge carriers velocity distribution which, however, accounts for the particle collective motion, in addition to the thermal one. Present treatment is restricted to conservative (Thomson) scattering in the electron rest frame. The code, nonetheless, is completely general and inclusion of the relativistic QED resonant cross section, which is required in the modeling of the hard (~20-200 keV) spectral tails observed in the magnetar candidates, is under way. The properties of emerging spectra have been assessed under different conditions, by exploring the model parameter space, including effects arising from the viewing geometry. Monte Carlo runs have been collected into a spectral archive. Two tabulated XSPEC spectral models, with and without viewing angles, have been produced and applied to the 0.1-10 keV XMM-Newton EPIC-pn spectrum of the AXP CXOU J1647-4552.

<Context>. We report on near-infrared (IR) observations of the three anomalous X-ray pulsars XTE J1810-197, 1RXS J1708-4009, 1E 1841-045 and the soft gamma-ray repeater SGR 1900+14, taken with the ESO-VLT, the Gemini, and the CFHT telescopes. <Aims>. This work is aimed at identifying and/or confirming the IR counterparts of these magnetars, as well as at measuring their possible IR variability. <Methods>. In order to perform photometry of objects as faint as Ks~20, we have used data taken with the largest telescopes, equipped with the most advanced IR detectors and in most of the cases with Adaptive Optics devices. The latter are critical to achieve the sharp spatial accuracy required to pinpoint faint objects in crowded fields. <Results>. We confirm with high confidence the identification of the IR counterpart to XTE J1810-197, and its IR variability. For 1E 1841-045 and SGR 1900+14 we propose two candidate IR counterparts based on the detection of IR variability. For 1RXS J1708-4009 we show that none of the potential counterparts within the source X-ray error circle can be yet convincingly associated with this AXP. <Conclusions>. The IR variability of the AXP XTE J1810-197 does not follow the same monotonic decrease of its post-outburst X-ray emission. Instead, the IR variability appears more similar to the one observed in radio band, although simultaneous IR and radio observations are crucial to draw any conclusion in this respect. For 1E 1841-045 and SGR 1900+14, follow-up observations are needed to confirm our proposed candidates with higher confidence.

The soft gamma-ray repeater SGR 1806-20 has been attracting a lot of attention owing to the fact that in December 2004 it emitted the most powerful giant flare ever observed. Here we present the results of the first Suzaku observation of SGR 1806-20, that seems to have reached a state characterized by a flux close to the pre-flare level and by a relatively soft spectrum. Despite this, the source remained quite active, as testified by several short bursts observed by Suzaku. We discuss the broadband spectral properties of SGR 1806-20 in the context of the magnetar model, considering its recent theoretical developments.

Magnetar's persistent emission above 10 keV was recently discovered thanks to the imaging capabilities of the IBIS coded mask telescope on board the INTEGRAL satellite. The only two sources that show some degree of long term variability are SGR 1806-20 and 1RXS J170849.0-400910. We find some indications that variability of these hard tails could be the driver of the spectral variability measured in these sources below 10 keV. In addition we report for the first time the detection at 2.8 sigma level of pulsations in the hard X-ray tail of SGR 1806-20.

In December 2004, the soft gamma-ray repeater SGR 1806-20 emitted the most powerful giant flare ever observed. This probably involved a large-scale rearrangement of the magnetosphere leading to observable variations in the properties of its X-ray emission. Here we present the results of the first Suzaku observation of SGR 1806-20, together with almost simultaneous observations with XMM-Newton and INTEGRAL. The source seems to have reached a state characterized by a flux close to the pre-flare level and by a relatively soft spectrum. Despite this, SGR 1806-20 remained quite active also after the giant flare, allowing us to study several short bursts observed by Suzaku in the 1-100 keV range. We discuss the broad-band spectral properties of SGR 1806-20, covering both persistent and bursting emission, in the context of the magnetar model, and consider its recent theoretical developments.

We have carried out a search for radio emission from six X-ray dim isolated neutron stars (XDINSs) observed with the Robert C. Byrd Green Bank Radio Telescope (GBT) at 820 MHz. No bursty or pulsed radio emission was found down to a 4sigma significance level. The corresponding flux limit is 0.01-0.04 mJy depending on the integration time for the particular source and pulse duty cycle of 2%. These are the most sensitive limits yet on radio emission from these objects.

We present our new advanced model for population synthesis of close-by cooling NSs. Detailed treatment of the initial spatial distribution of NS progenitors and a detailed ISM structure up to 3 kpc give us an opportunity to discuss the strategy to look for new isolated cooling NSs. Our main results in this respect are the following: new candidates are expected to be identified behind the Gould Belt, in directions to rich OB associations, in particular in the Cygnus-Cepheus region; new candidates, on average, are expected to be hotter than the known population of cooling NS. Besides the usual approach (looking for soft X-ray sources), the search in 'empty' $\gamma$-ray error boxes or among run-away OB stars may yield new X-ray thermally emitting NS candidates.

We report on a multi-band high-energy observing campaign aimed at studying the long term spectral variability of the Anomalous X-ray Pulsar (AXP) 1RXS J170849.0-400910, one of the magnetar candidates. We observed 1RXS J170849.0-400910 in Fall 2006 and Spring 2007 simultaneously with Swift/XRT, in the 0.1-10 keV energy range, and with INTEGRAL/IBIS, in the 20-200 keV energy range. Furthermore, we also reanalyzed, using the latest calibration and software, all the publicly available INTEGRAL data since 2002, and the soft X-ray data starting from 1999 taken using BeppoSAX, Chandra, XMM, and Swift/XRT, in order to study the soft and hard X-ray spectral variability of 1RXS J170849.0-400910. We find a long-term variability of the hard X-ray flux, extending the hardness-intensity correlation proposed for this source over 2 orders of magnitude in energy.

We report on two new XMM-Newton observations of the Anomalous X-ray Pulsar (AXP) 4U 0142+614 performed in March and July 2004, collecting the most accurate spectrum for this source to date. Furthermore, we analyse two short archival observations performed in February 2002 and January 2003 (the latter already reported by Gohler et al. 2005) in order to study the long term behaviour of this AXP. 4U 0142+614 appears to be relatively steady in flux between 2002 and 2004, and the phase-averaged spectrum does not show any significant variability between the four epochs. We derive the deepest upper limits to date on the presence of lines in the 4U 0142+614 spectrum as a function of energy: equivalent width in the 1-3 keV energy range < 4 eV and < 8 eV for narrow and broad lines, respectively. A remarkable energy dependence in both the pulse profile and the pulsed fraction is detected, and consequently pulse-phase spectroscopy shows spectral variability as a function of phase. By making use of XMM-Newton and INTEGRAL data, we successfully model the 1-250 keV spectrum of 4U 0142+614 with three models presented in Rea et al. (2007a), namely the canonical absorbed blackbody plus two power-laws, a resonant cyclotron scattering model plus one power-law and two log-parabolic functions.

We report on X-ray, optical, infrared and radio observations of the X-ray dim isolated neutron star (XDINS) 1RXS J214303.7+065419 (also known as RBS 1774). The X-ray observation was performed with the High Resolution Camera on board of the Chandra X-ray Observatory, allowing us to derive the most accurate position for this source (alpha = 21h43m3.38s, delta= +6deg54'17".53; 90% uncertainty of 0."6). Furthermore, we confirmed with a higher spatial accuracy the point-like nature of this X-ray source. Optical and infrared observations were taken in B, V, r', i', J, H and Ks filters using the Keck, VLT, Blanco and Magellan telescopes, while radio observations were obtained from the ATNF Parkes single dish at 2.9GHz and 708MHz. No plausible optical and/or infrared counterpart for RBS 1774 was detected within the refined sub-arsecond Chandra X-ray error circle. Present upper limits to the optical and infrared magnitudes are r'>25.7 and J>22.6 (5 sigma confidence level). Radio observations did not show evidence for radio pulsations down to a luminosity at 1.4 GHz of L < 0.02 mJy kpc^2, the deepest limit up to date for any XDINS, and lower than what expected for the majority of radio pulsars. We can hence conclude that, if RBS 1774 is active as radio pulsar, its non detection is more probably due to a geometrical bias rather than to a luminosity bias. Furthermore, no convincing evidence for RRAT-like radio bursts have been found. Our results on RBS 1774 are discussed and compared with the known properties of other thermally emitting neutron stars and of the radio pulsar population.

We present an empirical spectral modelling of the high energy emission of the anomalous X-ray pulsar 4U 0142+614, based on simultaneous Swift and INTEGRAL observations from X to gamma-ray energies. We adopted models contained in the XSPEC analysis package, as well as models based on recent theoretical studies, and restricted ourselves to those combinations of up to three components which produce a good fit while requiring the lowest number of free parameters. Only three models were found to fit satisfactorily the 0.5-250keV spectrum of 4U 0142+614: i) a ~0.4keV blackbody and two power-laws, ii) a resonant cyclotron scattering model plus a power-law and iii) two log-parabolic functions. We found that only the latter two models do not over-predict the infrared/optical emission observed simultaneously from this AXP, and only the log-parabolic functions can naturally account for the upper limits set by COMPTEL in the gamma-ray range. A possible interpretation of the two log-parabolae in terms of inverse Compton scattering of soft X-ray photons by very energetic particles is discussed.

We present an analysis of recent XMM-Newton and HST archive data of the ultraluminous X-ray source NGC 1313 X-2. Quasi-simultaneous observations taken with XMM-Newton, HST and VLT allow us to study both the X-ray light curve and its correlation with the optical emission of the two proposed ULX counterparts. At the end of December 2003 the source experienced a short, but intense flare, reaching a maximum luminosity of ~10^40 erg/s. At the same time, the optical flux of both the suggested counterparts did not show pronounced variations (<30%). Assuming that the ULX emission is isotropic and taking X-ray reprocessing into account, the optical data for one of the proposed counterparts are consistent with it being an early type, main sequence star of about 10-18 Msun losing matter through Roche-lobe overflow onto a ~120 Msun black hole at an orbital separation corresponding to a period of about 2 days.

Convincing evidence for long-term variations in the emission properties of the anomalous X-ray pulsar 1RXS J170849.0-400910 has been gathered in the last few years. In particular, and following the pulsar glitches of 1999 and 2001, XMM-Newton witnessed in 2003 a decline of the X-ray flux accompanied by a definite spectral softening. This suggested the existence of a correlation between the luminosity and the spectral hardness in this source, similar to that seen in the soft gamma-repeater SGR 1806-20. Here we report on new Chandra and Swift observations of 1RXS J170849.0-400910 performed in 2004 and 2005, respectively. These observations confirm and strengthen the proposed correlation. The trend appears to have now reversed: the flux increased and the spectrum is now harder. The consequences of these observations for the twisted magnetosphere scenario for anomalous X-ray pulsars are briefly discussed.

We report here on X-ray and IR observations of the Anomalous X-ray Pulsar 1RXS J170849-400910. First, we report on new XMM-Newton, Swift-XRT and Chandra observations of this AXP, which confirm the intensity-hardness correlation observed in the long term X-ray monitoring of this source by Rea et al. (2005). These new X-ray observations show that the AXP flux is rising again, and the spectrum hardening. If the increase in the source intensity is indeed connected with the glitches and a possible bursting activity, we expect this source to enter in a bursting active phase around 2006-2007. Second, we report on deep IR observations of 1RXS J170849-400910, taken with the VLT-NACO adaptive optics, showing that there are many weak sources consistent with the AXP position. Neither star A or B, previously proposed by different authors, might yet be conclusively recognised as the IR counterpart of 1RXS J170849-400910. Third, using Monte Carlo simulations, we re-address the calculation of the significance of the absorption line found in a phase-resolved spectrum of this source by Rea et al. (2003), and interpreted as a resonant scattering cyclotron feature.

We present a systematic analysis of all the BeppoSAX data of SGR1900+14. The observations spanning five years show that the source was brighter than usual on two occasions: ~20 days after the August 1998 giant flare and during the 10^5 s long X-ray afterglow following the April 2001 intermediate flare. In the latter case, we explore the possibility of describing the observed short term spectral evolution only with a change of the temperature of the blackbody component. In the only BeppoSAX observation performed before the giant flare, the spectrum of the SGR1900+14 persistent emission was significantly harder and detected also above 10 keV with the PDS instrument. In the last BeppoSAX observation (April 2002) the flux was at least a factor 1.2 below the historical level, suggesting that the source was entering a quiescent period.

The Magnificent Seven have all been discovered by their exceptional soft X-ray spectra and high ratios of X-ray to optical flux. They all are considered to be nearby sources. Searching for similar neutron stars with larger distances, one expects larger interstellar absorption resulting in harder X-ray counterparts. Current interstellar absorption treatment depends on chosen abundances and scattering cross-sections of the elements as well as on the 3D distribution of the interstellar medium. After a discussion of these factors we use the comprehensive 3D measurements of the Local Bubble by Lallement et al. 2003 to construct two simple models of the 3D distribution of the hydrogen column density. We test these models by using a set of soft X-ray sources with known distances. Finally, we discuss possible applications for distance estimations and population synthesis studies.

The X-ray spectra of the magnetar candidates are customarily fitted with an empirical, two component model: an absorbed blackbody and a power-law. However, the physical interpretation of these two spectral components is rarely discussed. It has been recently proposed that the presence of a hot plasma in the magnetosphere of highly magnetized neutron stars might distort, through efficient resonant cyclotron scattering, the thermal emission from the neutron star surface, resulting in the production of non-thermal spectra. Here we discuss the Resonant Cyclotron Scattering (RCS) model, and present its XSPEC implementation, as well as preliminary results of its application to Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters.

A ~50 ks XMM-Newton observation of SGR 1900+14 has been carried out in September 2005, after almost three years during which no bursts were detected from this soft gamma-ray repeater. The 0.8-10 keV spectrum was well fit by a power law plus blackbody model with photon index 1.9+/-0.1, temperature kT=0.47+/-0.02 keV and N_H = (2.12+/-0.08)x10^22 cm^-2, similar to previous observations of this source. The flux was ~5x10^{-12} erg/cm2/s, a factor 2 dimmer than the typical value and the smallest ever seen from SGR 1900+14. The long term fading of the persistent emission has been interrupted by the recent burst reactivation of the source. A target of opportunity XMM-Newton observation performed in April 2006 showed a flux ~15% higher. This variation was not accompanied by significant changes in the spectrum, pulsed fraction and light curve profile. We searched for emission and absorption lines in the spectra of the two observations, with negative results and setting tight upper limits of 50-200 eV (3 sigma), depending on the assumed line energy and width, on the equivalent width of lines in the 1-9 keV range.

In the last years, optical studies of Isolated Neutron Stars (INSs) have expanded from the more classical rotation-powered ones to other categories, like the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma-ray Repeaters (SGRs), which make up the class of the magnetars, the radio-quiet INSs with X-ray thermal emission and, more recently, the enigmatic Compact Central Objects (CCOs) in supernova remnants. Apart from 10 rotation-powered pulsars, so far optical/IR counterparts have been found for 5 magnetars and for 4 INSs. In this work we present some of the latest observational results obtained from optical/IR observations of different types of INSs.

The XMM-Newton spectra of the isolated neutron star RX J0720.4-3125 obtained over 4.5 years can be described by sinusoidal variations in the inferred blackbody temperature, the size of the emitting area and the depth of the absorption line with a period of 7.1 +/- 0.5 years, which we suggest to be the precession period of the neutron star. Precession of a neutron star with two hot spots of different temperature and size, probably not located exactly in antipodal positions, may account for the variations in the X-ray spectra, changes in the pulsed fraction, shape of the light curve and the phase-lag between soft and hard energy bands observed from RX J0720.4-3125. An independent sinusoidal fit to published and new pulse timing residuals from a coherent analysis covering ~12 years yields a consistent period of 7.7 +/- 0.6 years supporting the precession model.

We highlight similarities between recently discovered Rotating Radio Transients and X-ray Dim Isolated Neutron Stars. In particular, it is shown that X-ray Dim Isolated Neutron Stars have a birthrate comparable to that of Rotating Radio Transients. On the contrary, magnetars have too low a formation rate to account for the bulk of the radio transient population. The consequences of the recent detection of a thermal X-ray source associated with one of the Rotating Radio Transients on the proposed scenarios for these sources are also discussed.

The sky region containing the soft gamma-ray repeater SGR 1627-41 has been observed three times with XMM-Newton in February and September 2004. SGR 1627-41 has been detected with an absorbed flux of ~9x10^{-14} erg cm^{-2} s^{-1} (2-10 keV). For a distance of 11 kpc, this corresponds to a luminosity of ~3x10^{33} erg s^{-1}, the smallest ever observed for a Soft Gamma Repeater and possibly related to the long period of inactivity of this source. The observed flux is smaller than that seen with Chandra in 2001-2003, suggesting that the source was still fading and had not yet reached a steady quiescent level. The spectrum is equally well fit by a steep power law (photon index \~3.2) or by a blackbody with temperature kT~0.8 keV. We also report on the INTEGRAL transient IGR J16358-4726 that lies at ~10' from SGR 1627-41. It was detected only in September 2004 with a luminosity of ~4x10^{33} erg s^{-1} (for d=7 kpc), while in February 2004 it was at least a factor 10 fainter.

The study of thermal emission from isolated, cooling neutron stars plays a key role in probing the physical conditions of both the star crust and the core. The comparison of theoretical models for the star thermal evolution with the surface temperature derived from X-ray observations of sources of different age is one of the main tools to investigate the properties of the interior and constrain the equation of state. Here we propose to use population synthesis studies as an independent approach to test the physics governing the star cooling. Theoretical Log N - Log S distributions depend on the assumed neutron star thermal evolution. We have computed distributions for several different cooling scenarios and found that comparison with the observed Log N - Log S of isolated neutron stars is effective in discriminating among cooling models. Among the eleven cooling models considered in this paper, all of which may reproduce the observed temperature vs. age diagram, only at most three can explain the Log N - Log S distribution of close-by cooling neutron stars. The Log N - Log S test, being a ``global'' one and despite some limitations, appears capable of complementing the standard temperature vs. age test used up to now.

Recent Chandra and XMM-Newton observations of a number of X-ray ``dim'' pulsating neutron stars revealed quite unexpected features in the emission from these sources. Their soft thermal spectrum, believed to originate directly from the star surface, shows evidence for a phase-varying absorption line at some hundred eVs. The pulse modulation is relatively large (pulsed fractions in the range ~12%-35%), the pulse shape is often non-sinusoidal, and the hard X-ray color appears to be anti-correlated in phase with the total emission. Moreover, the prototype of this class, RX J0720.4-3125, has been found to undergo rather sensible changes both in its spectral and timing properties over a timescale of a few years. All these new findings seem difficult to reconcile with the standard picture of a cooling neutron star endowed with a purely dipolar magnetic field, at least if surface emission is produced in an atmosphere on top of the crust. In this paper we explore how a dipolar+quadrupolar star-centered field influence the properties of the observed lightcurves. The phase-resolved spectrum has been evaluated accounting for both radiative transfer in a magnetized atmosphere and general relativistic ray-bending. We computed over 78000 lightcurves varying the quadrupolar components and the viewing geometry. A comparison of the data with our model indicate that higher order multipoles are required to reproduce the observations.

We present archive ESO VLT photometric and spectroscopic data of the Ultraluminous X-ray source NGC 1313 X-2. The superb quality of the VLT images reveals that two distinct objects, with R magnitudes 23.7 and 23.6, are visible inside the Chandra error box. The two objects, separated by 0.75 arcsec, were unresolved in our previous ESO 3.6 m+EFOSC image. We show that both are stars in NGC 1313, the first a B0-O9 main sequence star of ~20 Msun, while the second a G supergiant of ~10 Msun. Irrespectively of which of the two objects the actual counterpart is, this implies that NGC 1313 X-2 is a high mass X-ray binary with a very massive donor.

XMM-Newton observed the soft gamma repeater SGR 1806-20 about two months after its 2004 December 27 giant flare. A comparison with the previous observations taken with the same instrument in 2003-2004 shows that the pulsed fraction and the spin-down rate have significantly decreased and that the spectrum slightly softened. These changes may indicate a global reconfiguration of the neutron star magnetosphere. The spectral analysis confirms that the presence of a blackbody component in addition to the power-law is required. Since this additional component is consistent with being constant with respect to the earlier observations, we explore the possibility of describing the long-term spectral evolution as only due to the power-law variations. In this case, the slope of the power-law does not significantly change and the spectral softening following the giant flare is caused by the increase of the relative contribution of the blackbody over the power-law component.

GCRT J1745-3009 is a peculiar transient radio-source in the direction of the Galactic Center. It was observed to emit a series of ~ 1 Jy bursts at 0.33 GHz, with typical duration ~ 10 min and at apparently regular intervals of ~ 77 min. If the source is indeed at the distance of the Galactic Center as it seems likely, we show that its observational properties are compatible with those expected from a double neutron star binary, similar to the double pulsar system J0737-3039. In the picture we propose the (coherent) radio emission comes from the shock originating in the interaction of the wind of the more energetic pulsar with the magnetosphere of the companion. The observed modulation of the radio signal is the consequence of an eccentric orbit, along which the separation between the two stars varies. This cyclically drives the shock inside the light cylinder radius of the less energetic pulsar.

We present a detailed analysis of XMM archival data of two Ultraluminous X-ray Sources (ULXs) in the nearby spiral galaxy NGC 1313: NGC 1313 X-1 and X-2. The spectral continuum of these sources was modeled with a soft thermal component plus a power-law. If the soft component originates from an accretion disk, the inferred mass of the compact remnant is > 100 M_sun, making it an Intermediate Mass Black Hole (IMBH). A detailed analysis of the residuals of the XMM EPIC-pn spectrum shows some evidence for the presence of an Oxygen emission line in NGC 1313 X-1. The simultaneous presence of an excess in emission, although at a much reduced significance level, at different energies in the X-ray spectra of NGC 1313 X-1 and X-2 is suggestive of typical emission lines from young supernova remnants. An optical counterpart for NGC 1313 X-2 was also identified. On an ESO 3.6 m image, the Chandra error box embraces a R \~ 23 mag stellar-like object and excludes a previously proposed optical counterpart.

We report on the results of a ~30 ks Chandra pointing of the soft gamma-ray repeater SGR 1806-20, the first X-ray observation with high spectral resolution performed after the 2004 December 27 Giant Flare. The source was found in a bursting active phase and with a significantly softer spectrum than that of the latest observations before the Giant Flare. The observed flux in the 2-10 keV range was ~2.2 x10^-11 erg/s/cm^2, about 20% lower than that measured three months before the event. This indicates that, although its giant flare was ~100 times more intense than those previously observed in two other soft gamma-ray repeaters, the post flare X-ray flux decay of SGR 1806-20 has been much faster. The pulsed fraction was about 3%, a smaller value than that observed before the flare. We discuss the different properties of the post-flare evolution of SGR 1806-20 in comparison to those of SGR 1900+14 and interpret the results as a strong evidence that a magnetospheric untwisting occurred (or is occurring) after the Giant Flare.

We report on the first XMM-Newton observation of the Anomalous X-ray Pulsar 1RXS J170849.0-400910. The source was observed in 2003 August and was found at a flux level a factor of about two lower than previous observations. Moreover, a significant spectral evolution appears to be present, the source exhibiting a much softer spectrum than in the past. Comparison of the present properties of 1RXS J170849.0-400910 with those from archival data shows a clear correlation between the X-ray flux and the spectral hardness. In particular, the flux and the spectral hardness reached a maximum level close to the two glitches the source experienced in 1999 and in 2001, and successively decreased. Although the excellent XMM-Newton spectral resolution should in principle allow us to detect the absorption line reported in a phase-resolved spectrum with BeppoSAX and interpreted as a cyclotron feature, we found no absorption features, neither in the phase averaged spectrum nor in the phase resolved spectra. We discuss in detail both the possibilities that the feature in the BeppoSAX data may have resulted from a spurious detection or that it is real and intrinsically variable. We then discuss a possible explanation for the glitches and for the softening of the source emission which followed the flux decrease, in the framework of the magnetar model.

The low energy (<10 keV) X-ray emission of the Soft Gamma-ray Repeater SGR1806-20 has been studied by means of four XMM-Newton observations carried out in the last two years, the latter performed in response to a strong sequence of hard X-ray bursts observed on 2004 October 5. The source was caught in different states of activity: over the 2003-2004 period the 2-10 keV flux doubled with respect to the historical level observed previously. The long term raise in luminosity was accompanied by a gradual hardening of the spectrum, with the power law photon index decreasing from 2.2 to 1.5, and by a growth of the bursting activity. The pulse period measurements obtained in the four observations are consistent with an average spin-down rate of 5.5x10e-10 s/s, higher than the values observed in the previous years. The long-term behavior of SGR1806-20 exhibits the correlation between spectral hardness and spin-down rate previously found only by comparing the properties of different sources (both SGRs and Anomalous X-ray Pulsars). The best quality spectrum (obtained on 6 September 2004) cannot be fitted by a single power law, but it requires an additional blackbody component (kT=0.79 keV, R_BB = 1.9 (d/15 kpc)^2 km), similar to the spectra observed in other SGRs and in Anomalous X-ray Pulsars. No spectral lines were found in the persistent emission, with equivalent width upper limits in the range 30-110 eV. Marginal evidence for an absorption feature at 4.2 keV is present in the cumulative spectrum of 69 bursts detected in September-October 2004.

We report the results of a recent (July 2004) XMM-Newton Target of Opportunity observation of the Anomalous X-ray pulsar 1E 1048.1-5937, together with a detailed re-analysis of previous observations carried out in 2000 and 2003. In July 2004 the source had a 2-10 keV flux of 6.2$\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ and a pulsed fraction P$_F$=0.68. The comparison of the three data sets shows the presence of an anti-correlation between flux and pulsed fraction, implying that previous estimates of the source energetics based on the assumption of a large and constant pulsed fraction might be significantly underestimated. The source spectrum is well described by a power law plus blackbody model (kT~0.63 keV, photon index $\Gamma$~2.7-3.5) or, alternatively, by the sum of two blackbodies of which the hotter is Comptonized by relativistic electrons. In this case the temperatures are kT${_1}$~0.2-0.3 keV and kT${_2}$~0.4-0.5 keV and the emitting area of the cooler component is consistent with the whole neutron star surface. The long term luminosity variation of a factor >~2 is accompanied by relatively small variations in the spectral shape. Phase resolved spectroscopy indicates a harder spectrum in correspondence of the pulse maximum. No spectral features have been detected with 4$\sigma$ limits on the equivalent width in the range ~10-220 eV, depending on line energy and width.

We report on the results of a deep XMM-Newton observation of RBS 1774, the most recent dim isolated neutron star candidate found in the ROSAT archive data. Spectral and timing analysis of the high-quality PN and MOS data confirm the association of this source with an isolated neutron star. The spectrum is thermal and blackbody-like, and there is evidence at a significance level > 4sigma that the source is an X-ray pulsar, with spin period of 9.437 s. Spectral fitting reveils the presence of an absorption feature at ~0.7 keV, but at this level data do not have enough resolution to allow us to discriminate between an absorption line or an edge. We compare the newly measured properties of RBS 1774 with those of other known dim isolated neutron stars, and discuss possible interpretations for the absorption feature.

In the last few years considerable observational resources have been devoted to study the thermal emission from isolated neutron stars. Detailed XMM and Chandra observations revealed a number of features in the X-ray pulse profile, like asymmetry, energy dependence, and possible evolution of the pulse profile over a time scale of months or years. Here we show that these characteristics may be explained by a patchy surface temperature distribution, which is expected if the magnetic field has a complex structure in which higher order multipoles contribute together with the dipole. We reconsider these effects from a theoretical point of view, and discuss their implications to the observational properties of thermally emitting neutron stars.

We present new population synthesis calculations of close young neutron stars. In comparison with our previous investigation we use a different neutron star mass spectrum and different initial spatial and velocity distributions. The results confirm that most of ROSAT dim radioquiet isolated neutron stars had their origin in the Gould Belt. Several tens of young neutron stars can be identified in future in ROSAT data at low galactic latitudes and some of them also can be EGRET unidentified sources.

We report on the analysis of new X-ray data obtained with XMM-Newton and Chandra from two ROSAT-discovered X-ray dim isolated neutron stars. RX J0806.4-4123 was observed with XMM-Newton in April 2003, 2.5 years after the first observation. The EPIC-pn data confirm that this object is an X-ray pulsar with 11.371 s neutron star spin period. The X-ray spectrum is consistent with absorbed black-body emission with a temperature kT = 96 eV and NH = 4x10^19 cm^-2 without significant changes between the two observations. Four XMM-Newton observations of RX J0420.0-5022 between December 2002 and July 2003 did not confirm the 22.7 s pulsations originally indicated in ROSAT data, but clearly reveal a 3.453 s period. A fit to the X-ray spectrum using an absorbed black-body model yields kT = 45 eV, the lowest value found from the small group of XDINs and NH = 1.0x10^20 cm^-2. Including a broad absorption line improves the quality of the spectral fits considerably for both objects and may indicate the presence of absorption features similar to those reported from RBS1223, RX J1605.3+3249 and RX J0720.4-3125. For both targets we derive accurate X-ray positions from the Chandra data and present an optical counterpart candidate for RX J0420.0-5022 with B = 26.6+/-0.3 mag from VLT imaging.

We briefly review selected results in astrophysics of neutron stars (NSs) obtained during the last two years, focusing on isolated radioquiet objects. We discuss in some details population synthesis of close-by isolated NSs, spectra of INSs, detection of spectral features in these sources (including cyclotron features), recent results on velocity distribution of NSs and accretion onto INSs.

In a cool neutron star (T < 1.e6 K) endowed with a rather high magnetic field (B > 1.e13 G), a phase transition may occur in the outermost layers. As a consequence the neutron star becomes `bare', i.e. no gaseous atmosphere sits on the top of the crust. The surface of a cooling, bare neutron star does not necessarily emit a blackbody spectrum because the emissivity is strongly suppressed at energies below the electron plasma frequency, omega_p. Since omega_p ~ 1 keV under the conditions typical of the denseelectron gas in the condensate, the emission from a T ~ 100 eV bare neutron star will be substantially depressed with respect to that of a perfect Planckian radiator at most energies. Here we present a detailed analysis of the emission properties of a bare neutron star. In particular, we derive the surface emissivity for a Fe composition in a range of magnetic fields and temperatures representative of cooling isolated neutron stars, like RX J1856.5-3754. We find that the emitted spectrum is strongly dependent on the electron conductivity in the solid surface layers. In the cold electron gas approximation (no electron-lattice interactions), the spectrum turns out to be a featureless depressed blackbody in the 0.1-2 keV band with a steeper low-energy distribution. When damping effects due to collisions between electrons and the ion lattice (mainly due to electron-phonon interactions) are accounted for, the spectrum is more depressed depending on the magnetic field strength. Details of the emitted spectrum are found, however, to be strongly dependent of the assumed treatment of the transition from the external vacuum to the metallic surface. The implications of out results to RX J1856.5-3754 and other isolated neutron stars are discussed.

We briefly discuss the evolutionary path and observational appearance of isolated neutron stars (INSs) focusing on radioquiet objects. There are many reasons to believe that these sources are extremely elusive once the star surface has cooled down: their high spatial velocities, the long propeller stage and/or the very low accretion efficiency. We describe recent population synthesis models of close-by young INSs, highlighting the major difficulties encountered in the past by these simulations in reproducing the observed properties of known sources. As we show, a likely possibility is that most of the INSs in the Solar proximity are young (less than few Myrs) neutron stars born in the Gould Belt. To stay hot enough and sustain X-ray emission for a time $\approx 1$ Myr, they probably need to be low- to medium-massive, with $M$ less than$\sim 1.35 M_\odot$.

We present new optical and {\it Chandra} observations of the field containing the ultraluminous X-ray source NGC 1313 X-2. On an ESO 3.6 m image, the {\it Chandra} error box embraces a $R=21.6$ point-like object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of NGC 1313 X-2 is $\sim 500$. The value of $f_X/f_{opt}$, the X-ray variability history and the spectral distribution derived from a re-analysis of the {\it ROSAT}, {\it ASCA} and {\it XMM} data indicate a luminous X-ray binary in NGC 1313 as a likely explanation for NGC 1313 X-2. If the X-ray soft component observed in the {\it XMM} EPIC spectrum originates from an accretion disk, the inferred mass of the compact remnant is $\approx 100 M_\odot$, making it an intermediate mass black hole. The derived optical luminosity ($L\approx 10^5 L_\odot$) is consistent with that of a $\approx 15-20 M_\odot$ companion. The properties of the environment of NGC 1313 X-2 are briefly discussed.

We present new optical and {\it Chandra} observations of the field containing the ultraluminous X-ray source NGC1313 X-2. On an ESO 3.6 m image, the {\it Chandra} error box embraces a $R=21.6$ mag stellar-like object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of NGC 1313 X-2 is $\sim 500$. The value of $f_X/f_{opt}$, the X-ray variability history and spectral distribution indicate a luminous X-ray binary in NGC 1313 as a likely explanation for NGC 1313 X-2. The inferred optical luminosity ($L\approx 10^5 L_\odot$) is consistent with that of a $\approx 10 M_\odot$ companion.

We model populations of isolated neutron stars in the Galaxy following their orbital and magneto-rotational evolution. It is shown that accretors become more abundant than coolers at fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1},$ and one can predict that about one accreting neutron star per square degree should be observed at the {\it Chandra} and {\it Newton} flux limits of $\sim 10^{-16}$ erg cm$^{-2}$ s$^{-1}.$ The soft ROSAT sources associated with isolated neutron stars can be relatively young cooling objects only if the neutron star birth rate in the Solar vicinity during the last $\sim 10^6$ yr is higher than that inferred from radiopulsar observations.

The radiative acceleration of particles and the electrostatic potential fields that arise in low density plasmas hit by radiation produced by a transient, compact source are investigated. We calculate the dynamical evolution and asymptotic energy of the charged particles accelerated by the photons and the radiation-induced electric double layer in the full relativistic, Klein-Nishina regime. For fluxes in excess of $10^{27}$ ${\rm erg} {\rm cm}^{-2} {\rm s}^{-1}$, the radiative force on a diluted plasma ($n\la 10^{11}$ cm$^{-3}$) is so strong that electrons are accelerated rapidly to relativistic speeds while ions lag behind owing to their larger inertia. The ions are later effectively accelerated by the strong radiation-induced double layer electric field up to Lorentz factors $\approx 100$, attainable in the case of negligible Compton drag. The asymptotic energies achieved by both ions and electrons are larger by a factor 2--4 with respect to what one could naively expect assuming that the electron-ion assembly is a rigidly coupled system. The regime we investigate may be relevant within the framework of giant flares from soft gamma-repeaters.

The origin of the local population of young, cooling neutron stars is investigated with a population synthesis model taking into account the contribution of neutron stars born in the Gould Belt, in addition to those originating in the Galactic disk. We estimate their emission in the soft X-ray band as a function of distance and age and construct the Log N -- Log S distribution. It is shown that the inclusion of neutron stars from the Gould Belt provides a good fit to the observed Log N -- Log S distribution. As the Sun is situated inside the Gould Belt, one can naturally explain the comparative local overabundance of massive progenitors and can remove the difficulty of the deficit of relatively bright ($\ga 0.1$ ROSAT PSPC cts s$^{-1}$) cooling neutron stars previously reported from models where only neutron stars from the Galactic disk were accounted for.

Recent Chandra observations have shown that the soft X-ray spectrum from the isolated neutron star RX J1856.5-3754 is featureless and best fitted by a blackbody, in contrast with the predictions of current neutron star atmospheric models. The star distance (~120-140 pc) implies a radiation radius of at most \~5-6 km, too small for any neutron star equation of state. Proposed explanations for such a small radius include a reduced X-ray emitting region (as a polar cap), or the presence of a more compact object, as a bare quark/strange star. However, both interpretations rely on the presumption that the quark star or the cap radiates a pure blackbody spectrum, and no justification for this assumption has been presented yet. Here we discuss an alternative possibility. Cool neutron stars (T < 10^6 K) endowed with high magnetic field (B > 10^13 G) may suffer a phase transition in the outermost layers. The neutron star is then left bare of the gaseous atmosphere. We computed spectra from naked neutron stars with a surface Fe composition. Depending on B, we found that the 0.1-2 keV emission can be featureless and indistinguishable from a blackbody. Owing to the reduced surface emissivity, the star only radiates ~30-50 % of the blackbody power, thus that the size of the emitting region is larger than for a perfect planckian emitter. When applied to RX J1856.5-3754 our model accounts for the observed X-ray properties and solves the paradox of the small radius. The apparent star radius is ~10-12 km, consistent with equations of state of a neutron star. The optical emission of RX J1856.5-3754 may be explained by the presence a thin gaseous shell on the top of the Fe condensate.

We present Log N -- Log S distribution for close-by young isolated neutron stars. On the basis of this distribution it is shown that the seven ROSAT isolated neutron stars (if they are young cooling objects) are genetically related to the Gould Belt. We predict, that there are about few tens unidentified close-by young isolated neutron stars in the ROSAT All-Sky Survey. The possibility that these seven peculiar sources contain a neutron star less massive and more magnetized then in ordinary radiopulsars is also discussed. In the aftermath of relatively close recent supernova explosions (1 kpc around the Sun, a few Myrs ago), a few black holes might have been formed, according to the local initial mass function. We thus discuss the possibility of determining approximate positions of close-by isolated black holes using data on runaway stars and simple calculations of binary evolution and disruption.

We report evidence of cyclotron resonance features from the Soft Gamma Repeater SGR 1806-20 in outburst, detected with the Rossi X-ray Timing Explorer in the spectrum of a long, complex precursor that preceded a strong burst. The features consist of a narrow 5.0 keV absorption line with modulation near its second and third harmonics (at 11.2 keV and 17.5 keV respectively). The line features are transient and are detected in the harder part of the precursor. The 5.0 keV feature is strong, with an equivalent width of ~ 500 eV and a narrow width of less than 0.4 keV. Interpreting the features as electron cyclotron lines in the context of accretion models leads to a large mass-radius ratio (M/R > 0.3 M_sun/km) that is inconsistent with neutron stars or that requires a low (5-7)x10^{11} G magnetic field that is unlikely for SGRs. The line widths are also narrow compared with those of electron cyclotron resonances observed so far in X-ray pulsars. In the magnetar picture, the features are plausibly explained as ion cyclotron resonances in an ultra-strong magnetic field that have recently been predicted from magnetar candidates. In this view, the 5.0 keV feature is consistent with a proton cyclotron fundamental whose energy and width are close to model predictions. The line energy would correspond to a surface magnetic field of 1.0x10^{15} G for SGR 1806-20, in good agreement with that inferred from the spin-down measure in the source.

We report the discovery of a new ultraluminous X-ray source (ULX) in the nearby galaxy NGC 4244 from Chandra archival data. The source, 1WGAJ1216.9+3743, is one of the least luminous and softest ULXs discovered so far. Its X-ray spectrum is the best available for a representative of the soft ULXs, a class of sources recently discovered by Chandra and XMM-Newton. 1WGAJ1216.9+3742 appears point-like in the Chandra image and has a complex spectral shape: a multicolor disk model, suitable for brighter ULXs, is inadequate for this source. 1WGAJ1216.9+3743 spectrum is heavily absorbed (NH ~ (1-4)E21 cm^-2),and very soft. The best-fit power-law model gives a photon index ~ 5 and implies a luminosity L(0.5-10 keV) ~ 1E39 erg/s. A comparison with previous detections shows that, despite the variability displayed by the source during the Chandra observation, 1WGAJ1216.9+3743 count rate, spectral shape and absorption are practically unchange over a 9-year period. We performed also deep optical imaging of the field containing the X-ray source, and found a possible R ~ 23.7 counterpart.

The effects of bulk motion comptonization on the spectral formation in a converging flow onto a black hole are investigated. The problem is tackled by means of both a fully relativistic, angle-dependent transfer code and a semi-analytical, diffusion-approximation method. We find that a power-law high-energy tail is a ubiquitous feature in converging flows and that the two approaches produce consistent results at large enough accretion rates, when photon diffusion holds. Our semi-analytical approach is based on an expansion in eigenfunctions of the diffusion equation. Contrary to previous investigations based on the same method we find that, although the power-law tail at really large energies is always dominated by the flatter spectral mode, the slope of the hard X-ray portion of the spectrum is dictated by the second mode and it approaches Gamma=3 at large accretion rate, irrespective of the model parameters. The photon index in the tail is found to be largely independent on the spatial distribution of soft seed photons when the accretion rate is either quite low (< 5 in Eddington units) or sufficiently high (> 10). On the other hand, the spatial distribution of source photons controls the photon index at intermediate accretion rates, when Gamma switches from the first to the second mode. Our analysis confirms that a hard tail with photon index Gamma <3 is produced by the up-scattering of primary photons onto infalling electrons if the central object is a black hole.

We suggest that the seven radio-quiet isolated neutron stars observed with ROSAT are young cooling objects associated to recent near-by supernova explosions which formed runaway stars and the Local Bubble, affecting the topology of the interstellar medium in the vicinity of the Sun (within a few hundred parsecs). In the aftermath of these explosions, a few black holes might have been formed, according to the local initial mass function. We thus discuss the possibility of determining approximate positions of close-by isolated black holes using data on runaway stars and simple calculations of binary evolution and disruption.

We model Galactic populations of accreting and cooling isolated neutron stars in the attempt to explore their link with a new class of dim soft X-ray sources revealed by ROSAT. For accretors we follow the magneto-rotational and dynamical evolution in the Galactic potential and a realistic large scale distribution of the interstellar medium is used. Under standard assumptions old neutron stars enter the accretor stage only if their magnetic field exceeds $\approx 10^{11}$--$10^{12}$ G. We predict about 1 source per square degree for fluxes $\approx 10^{-15}$--$10^{-16}$ erg cm$^{-2}$s$^{-1}$ in the energy range 0.5-2 keV. Cooling neutron stars are explored within a simpler model of local sources, including however interstellar absorption. They are found to be significantly less abundant at low fluxes, $<0.1$ sources per square degree, but dominate over accretors at higher flux levels ($\approx 10^{-12}$--$10^{-11}$ erg cm$^{-2}$s$^{-1}$). We suggest that the faint sources observed by ROSAT may be young cooling neutron stars with typical age $\la 10^6$ yrs, if the total number of young neutron stars in the Solar proximity is $\sim 10$ times higher than inferred from radiopulsars statistics.

We report on the identification of a new possible Isolated Neutron Star candidate in archival ROSAT observations. The source 1RXS J214303.7+065419, listed in the ROSAT Bright Survey as RBS 1774, is very soft, exhibits a thermal spectrum well fitted by a blackbody at ${T}\sim 90$ eV and has a low column density, ${N_H}\sim 5\times 10^{20}$ ${cm}^{-2}$. Catalogue searches revealed no known sources in other energy bands close to the X-ray position of RBS 1774. Follow-up optical observations with NTT showed no peculiar object within the X-ray error circle. The absence of any plausible optical counterpart down to ${m_R}\sim 23$ results in an X-ray to optical flux ratio in excess of 1000.

A great deal of interest has been recently raised in connection with the possibility that soft $\gamma$-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) contain {\em magnetars}, young neutron stars endowed with magnetic fields $\gtrsim 10^{14}$ G. In this paper we calculate thermal spectra from ultra-magnetized neutron stars for values of the luminosity and magnetic field believed to be relevant to SGRs and AXPs. Emergent spectra are found to be very close to a blackbody at the star effective temperature and exhibit a distinctive absorption feature at the proton cyclotron energy $E_{c,p}\simeq 0.63 (B/10^{14} {\rm G})$ keV. The proton cyclotron features (PCFs) are conspicuous (equivalent width of up to many hundreds eV) and relatively broad ($\Delta E/E \sim 0.05-0.2$). The detection of the PCFs is well within the capabilities of present X-ray spectrometers, like the HETGS and METGS on board Chandra. Their observation might provide decisive evidence in favor of the existence of magnetars.

We report the results of X-ray (2-10 keV) observations of the giant molecular clouds SgrB, SgrC and SgrD in the Galactic Center region, together with the discovery of the point-like source SAXJ1748.2-2808. The data have been obtained with the MECS instrument on the BeppoSAX satellite. The core of SgrB2 has an X-ray luminosity of 6x10^34 erg/s and its spectrum is characterized by a strong Fe emission line at 6.5 keV with an equivalent width of 2 keV. Faint diffuse X-ray emission is detected from SgrC and from the SNR G1.05-0.15 (SgrD). A new, unresolved source with a strong Fe line has been discovered in the SgrD region. This source, SAXJ1748.2-2808, is probably associated with a SiO and OH maser source at the Galactic Center distance. If so, its luminosity is 10^34 erg/s. We propose that the X-ray emission from SAX J1748.2-2808 is produced either by protostars or by a giant molecular cloud core. Emission from sources similar to SAX J1748.2-2808 could have an impact on the expected contribution on the observed Fe line emission from the Galactic ridge.

Transonic discs with accretion rates relevant to intrinsically bright Galactic X-ray sources ($L\approx 10^{38}$-$10^{39} {\rm erg s}^{-1}$) exhibit a time dependent cyclic behaviour due to the onset of a thermal instability driven by radiation pressure. In this paper we calculate radiation spectra emitted from thermally-unstable discs to provide detailed theoretical predictions for observationally relevant quantities. The emergent spectrum has been obtained by solving self-consistently the vertical structure and radiative transfer in the disc atmosphere. We focus on four particular stages of the disc evolution, the maximal evacuation stage and three intermediate stages during the replenishment phase. The disc is found to undergo rather dramatic spectral changes during the evolution, emitting mainly in the 1-10 keV band during outburst and in the 0.1-1 keV band off-outburst. Local spectra, although different in shape from a blackbody at the disc effective temperature, may be characterized in terms of a hardening factor $f$. We have found that $f$ is rather constant both in radius and in time, with a typical value $\sim 1.65$.

Thermal instability driven by radiation pressure might be relevant for intrinsically bright accreting sources. The most promising candidate where this instability seems to be at work is one of the two known galactic superluminal sources, GRS 1915+105 (Belloni et al. 1997). In spite of being of relevance, this scenario has not yet been confirmed by proper time-dependent modelling. Non-linear time-dependent calculations performed by Szuszkiewicz and Miller (1998) show that thermally unstable discs undergo limit-cycle behaviour with successive evacuation and refilling of the central parts of the disc. This evolution is very similar to the one proposed by Belloni et al. (1997) in their phenomenological model. Further investigations are needed to confirm the thermal instability being operational in this source. First of all the spectra emitted from the disc during its evolution should be calculated and compared with observations. Here such spectra are computed assuming local blackbody emission from the best studied transonic disc model.

In this article we briefly review our recent results on evolution and properties of isolated neutron stars (INSs) in the Galaxy. As the first step we calculate a {\it census} of INSs in our Galaxy. We infer a lower bound for the mean kick velocity of NSs, $ <V> \sim $(200-300) ${\rm km s^{-1}}$. The same conclusion is reached for both a constant magnetic field ($B\sim 10^{12}$ G) and for a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. These results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to $\sim$few percents. Then we show that the range of minimum value of magnetic moment, $\mu_b$: $\sim 10^{29.5}\ge \mu_b \ge 10^{28} {\rm G} {\rm cm}^3$, and the characteristic decay time, $t_d$: $\sim 10^8\ge t_d \ge 10^7 {\rm yrs}$, can be excluded assuming the standard initial magnetic momentum, $\mu_0=10^{30} {\rm G} {\rm cm}^3$, if accreting INSs are observed. For these parameters an INS would never reach the stage of accretion from the interstellar medium even for a low space velocity of the star and high density of the ambient plasma. The range of excluded parameters increases for lower values of $\mu_0$. It is shown that old accreting INSs become more abundant than young cooling INSs at X-ray fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. We can predict that about one accreting INS per square degree should be observed at the {\it Chandra} and {\it Newton} flux limits of $\sim 10^{-16}$ erg cm$^{-2}$ s$^{-1}.$ The weak {\it ROSAT} sources, associated with INSs, can be young cooling objects, if the NSs birth rate in the solar vicinity during the last $\sim 10^6$ yr was much higher than inferred from radiopulsar observations.

In this Letter we propose a simple thermal comptonization model to account for the observed properties of the phase lags associated to the ``plateau'' intervals of GRS 1915+105. By invoking a temperature stratification in a corona and assuming that the optical depth of the comptonizing region increases as the disk inner radius moves inward, we are able to reproduce both the observed colors and time lags in the continuum.

In some models of $\gamma$-ray bursts super-strong electric fields ($E\sim 10^{14} {\rm statvolt cm}^{-1}$) have been surmized. Such large fields may provide copious pair production through vacuum polarization. Here we examine various astrophysical scenarios where huge electric field can be present. It is shown that when the conditions for quantistic vacuum breakdown are met, pair creation through conventional processes is likely to be always more important. Although of great physical interest, quantum pair production seems therefore of little astrophysical relevance.

We present an algorithm for two-dimensional radiative transfer in axisymmetric, circumstellar media. The formal integration of the transfer equation is performed by a generalization of the short characteristics (SC) method to spherical coordinates. Accelerated Lambda Iteration (ALI) and Ng's algorithm are used to converge towards a solution. By taking a logarithmically spaced radial coordinate grid, the method has the natural capability of treating problems that span several decades in radius, in the most extreme case from the stellar radius up to parsec scale. Flux conservation is guaranteed in spherical coordinates by a particular choice of discrete photon directions and a special treatment of nearly-radially outward propagating radiation. The algorithm works well from zero up to very high optical depth, and can be used for a wide variety of transfer problems, including non-LTE line formation, dust continuum transfer and high temperature processes such as compton scattering. In this paper we focus on multiple scattering off dust grains and on non-LTE transfer in molecular and atomic lines. Line transfer is treated according to an ALI scheme for multi-level atoms/molecules, and includes both random and systematic velocity fields. The algorithms are implemented in a multi-purpose user-friendly radiative transfer program named RADICAL. We present two example computations: one of dust scattering in the Egg Nebula, and one of non-LTE line formation in rotational transitions of HCO$^{+}$ in a flattened protostellar collapsing cloud.

We present a detailed investigation of atmospheres around accreting neutron stars with high magnetic field ($B\gtrsim 10^{12}$ G) and low luminosity ($L\lesssim 10^{33}$ erg/s). We compute the atmospheric structure, intensity and emergent spectrum for a plane-parallel, pure hydrogen medium by solving the transfer equations for the normal modes coupled to the hydrostatic and energy balance equations. The hard tail found in previous investigations for accreting, non-magnetic neutron stars with comparable luminosity is suppressed and the X-ray spectrum, although still harder than a blackbody at the star effective temperature, is nearly planckian in shape. Spectra from accreting atmospheres, both with high and low fields, are found to exhibit a significant excess at optical wavelengths above the Rayleigh-Jeans tail of the X-ray continuum.

In this Letter we investigate the properties of advection-dominated accretion flows (ADAFs) fed by the evaporation of a Shakura-Sunyaev accretion disk (SSD). In our picture the ADAF fills the central cavity evacuated by the SSD and extends beyond the transition radius into a coronal region. We find that, because of global angular momentum conservation, a significant fraction of the hot gas flows away from the black hole forming a transsonic wind, unless the injection rate depends only weakly on radius (if $r^2\dot\sigma\propto r^{-\xi}$, $\xi< 1/2$). The Bernoulli number of the inflowing gas is negative if the transition radius is $\lesssim 100$ Schwarzschild radii, so matter falling into the hole is gravitationally bound. The ratio of inflowing to outflowing mass is $\approx 1/2$, so in these solutions the accretion rate is of the same order as in standard ADAFs and much larger than in advection-dominated inflow/outflow models (ADIOS). The possible relevance of evaporation-fed solutions to accretion flows in black hole X-ray binaries is briefly discussed.

As many as $10^9$ neutron stars populate the Galaxy, but only $\approx 10^3$ are directly observed as pulsars or as accreting sources in X-ray binaries. In principle also the accretion of the interstellar medium may make isolated neutron stars shine, and their weak luminosity could be detected in soft X-rays. Recent ROSAT observations have convincingly shown that neutron stars accreting from the interstellar medium are extremely rare, if observed at all, in contrast with earlier theoretical predictions. Until now two possible explanations for their elusiveness have been proposed: their velocity distribution may peak at $\sim 200-400 {\rm km s}^{-1}$, as inferred from pulsar statistics, and this would severely choke accretion; the magnetic field may decay on timescales $\sim 10^8-10^9$ yr, preventing a large fraction of neutron stars from entering the accretor stage. The search for accreting neutron stars has produced up to now a handful of promising candidates. While little doubt is left that these objects are indeed isolated neutron stars, the nature of their emission is still controversial. In particular accreting objects can be confused with much younger, cooling neutron stars. However, a combination of observations and theoretical modeling may help in discriminating between the two classes.

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical sample of stars for different values of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean velocity $0\leq < V>\leq 550$ ${\rm km s^{-1}}$. The spin-down, induced by dipole losses and the interaction with the ambient medium, is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All Sky Survey an upper limit of $\sim 10$ accreting neutron stars within $\sim 140$ pc from the Sun, we infer a lower bound for the mean kick velocity, $ < V>\ga 200-300$ ${\rm km s^{-1}}$. The same conclusion is reached for both a constant ($B\sim 10^{12}$ G) and a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. Present results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to $\la 1%$.

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical sample of stars for different values of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean velocity $0\leq < V>\leq 550$ ${\rm km s^{-1}}$. The spin--down, induced by dipole losses and the interaction with the ambient medium, is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All Sky Survey an upper limit of $\sim 10$ accreting neutron stars within $\sim 140$ pc from the Sun, we infer a lower bound for the mean kick velocity, $ < V>\gtrsim 200-300$ ${\rm km s^{-1}},$ corresponding to a velocity dispersion $\sigma_V\gtrsim 125-190$ km s$^{-1}$. The same conclusion is reached for both a constant magnetic field ($B\sim 10^{12}$ G) and a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. Such high velocities are consistent with those derived from radio pulsar observations. Present results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to less than 1%.

We consider a black hole accreting spherically from the surrounding medium. If accretion produces a luminosity close to the Eddington limit the hole acquires a net charge so that electrons and ions can fall with the same velocity. The condition for the electrostatic field to be large enough to break the vacuum near the hole horizon translates into an upper limit for the hole mass, $M\sim 6.6\times 10^{20} {\rm g}.$ The astrophysical conditions under which this phaenomenon can take place are rather extreme, but in principle they could be met by a mini black hole residing at the center of a star.

The existence of bimodal disks is investigated. Following a simple argument based on energetic considerations we show that stationary, bimodal accretion disk models in which a Shakura--Sunyaev disk (SSD) at large radii matches an advection dominated accretion flow (ADAF) at smaller radii are never possible using the standard slim disk approach, unless some extra energy flux is present. The same argument, however, predicts the possibility of a transition from an outer Shapiro--Lightman--Eardley (SLE) disk to an ADAF, and from a SLE disk to a SSD. Both types of solutions have been found.

Old neutron stars (ONSs) which have radiated away their internal and rotational energy may still shine if accreting the interstellar medium. Rather stringent limits from the analysis of ROSAT surveys indicate that most optimistic predictions on ONSs observability are in excess of a factor as large as $\sim 100$. Here we explore two possible evolutionary scenarios that may account for the paucity of ONSs. In the first it is assumed that the ONS population is not too fast ($V<100 km s^{-1}$) and that magnetic field decay guides the evolution. In the second, NSs move with high speed ($V>100$ km s$^{-1}$) and preserve their magnetic field at birth. We find that according to the former scenario most ONSs are now in the propeller phase, while in the latter nearly all ONSs are silent, dead pulsars.

Old neutron stars (ONSs) which have radiated away their internal and rotational energy may still shine if accreting the interstellar medium. Despite their large number, only two promising candidates have been detected so far and rather stringent limits on their observability follow from the analysis of ROSAT surveys. This contrasts with optimistic theoretical estimates that predicted a large number of sources in ROSAT fields. We have reconsidered the issue of ONSs observability, accounting for the spin and magnetic field evolution over the neutron star lifetime. In the framework of a spin-induced field decay model, we show that the total number of ONSs which are, at present, in the accretion stage is reduced by a factor ~5 over previous figures if the characteristic timescale for crustal current dissipation is ~ 10^8 - 10^9 yr. This brings theoretical predictions much closer to observational limits. Most ONSs should be at present in the propeller phase and, if subject to episodic flaring, they could be observable.

The structure of static atmospheres around unmagnetized neutron stars undergoing steady, spherical accretion is discussed. We focus on the ``hot'' configurations presented by Turolla et al. (1994) and calculate the radiation spectrum using a characteristics method. In particular, it is found that e$^+$--e$^-$ pair production may affect significantly the external atmospheric layers, where positron and proton number densities become of the same order. The consequent increase of the scattering opacity lowers the Eddington limit and this, in turn, may drive a dynamical instability if the accretion luminosity is large enough, ultimately producing a rapid expulsion of the envelope. If ``hot'' states are indeed accessible, this mechanism could give rise to transient phenomena in hard X--rays of potential great astrophysical interest.

As first shown by Shvartsman (1970), a neutron star accreting close to the Eddington limit must acquire a positive charge in order for electrons and protons to move at the same speed. The resulting electrostatic field may contribute to accelerating positrons produced near the star surface in conjunction with the radiative force. We reconsider the balance between energy gains and losses, including inverse Compton (IC), bremsstrahlung and non--radiative scatterings. It is found that, even accounting for IC losses only, the maximum positron energy never exceeds $\approx 400$ keV. The electrostatic field alone may produce energies $\approx 50$ keV at most. We also show that Coulomb collisions and annihilation with accreting electrons severely limit the number of positrons that escape to infinity.

The transport of photons in steady, spherical, scattering flows is investigated. The moment equations are solved analytically for accretion onto a Schwarzschild black hole, taking into full account relativistic effects. We show that the emergent radiation spectrum is a power law at high frequencies with a spectral index smaller (harder spectrum) than in the non--relativistic case. Radiative transfer in an expanding envelope is also analyzed. We find that adiabatic expansion produces a drift of injected monochromatic photons towards lower frequencies and the formation of a power--law, low--energy tail with spectral index $-3$.

Accretion onto old, isolated neutron stars is reviewed. The detection of their X--ray emission with ROSAT is discussed and a summary of the current status of observations is presented.

Results are presented from a time--dependent, numerical investigation of spherical accretion onto black holes, within the framework of relativistic radiation hydrodynamics. We have studied the stability of self--consistent, stationary solutions of black hole accretion with respect to thermal and radiative perturbations and also the non--linear evolution of unstable, high temperature models, heated by the hard radiation produced by the accretion flow itself in the inner region near to the horizon. In some cases, a hydrodynamic shock forms at around $10^3$--$10^4$ Schwarzschild radii, where Compton heating exceeds radiative cooling. The calculations were made using a suitably designed radiation hydrodynamics code, in which radiative transfer is handled by means of the PSTF moment formalism and which contains an original treatment of the radiation temperature equation.

The contribution of weakly--magnetized ($B\sim 10^9$ G) neutron stars accreting the interstellar medium to the diffuse X--ray emission observed in the Galactic Center is investigated. It is shown that, under rather conservative assumptions about the neutron stars and gas distributions, the accretion luminosity can account for a sizable fraction, possibly most, of the detected X--ray flux in the 2.5--7 keV band. In particular, model results are compared with {\it Granat\/} data and show a general agreement in both the flux energy and radial distributions.

In this paper we present a characteristic method for solving the transfer equation in differentially moving media in a curved spacetime. The method is completely general, but its capabilities are exploited at best in presence of symmetries, when the existence of conserved quantities allows to derive analytical expressions for the photon trajectories in phase space. In spherically--symmetric, stationary configurations the solution of the transfer problem is reduced to the integration of a single ordinary differential equation along the bi--parametric family of characteristic rays. Accurate expressions for the radiative processes relevant to continuum transfer in a hot astrophysical plasma have been used in evaluating the source term, including relativistic e--p, e--e bremsstrahlung and Compton scattering. A numerical code for the solution of the transfer problem in moving media in a Schwarzschild spacetime has been developed and tested. Some applications, concerning ``hot'' and ``cold'' accretion onto non--rotating black holes as well as static atmospheres around neutron stars, are presented and discussed.

Old isolated accreting neutron stars may show up among unidentified soft X--ray sources detected by the ROSAT All Sky Survey. We argue that the chances of identification are greater for ONSs located in the closest overdense regions of the solar neighbourhood. In particular, we consider the neutral hydrogen wall in the second Galactic quadrant ($15^\circ < l < 150^\circ$) which shrinks the estimated contour of the Local Bubble to $\approx 16-30$ pc. Due to their vicinity, we expect $\sim 10$ ONSs to be detectable at a relatively high flux limit ($\sim 0.1$ counts/s) in the 0.2--2.4 keV band. This implies that about 5 \% of unidentified sources above this threshold could be ONSs. No optical counterpart is present, but EUV emission from these objects could be detected using EUVE Lex filter at the highest sensitivity limits.

The issue of the observability of Old Isolated Neutron Stars (ONSs) accreting from the interstellar gas, is reconsidered using the spectra presented in Zampieri {\it et al.\/} (1995). In particular, we focus our attention on the overall soft X--ray emission of ONSs, which may provide a substantial contribution to the X--ray background in the 0.5--2 keV band. It is found that if accretion is channeled into polar caps by a magnetic field of $\sim 10^9$ G, up to $\sim$ $12-25 \%$ of the unresolved soft excess observed at high latitudes in the X--ray background can be explained in terms of emission from accreting isolated ONSs. However, this would imply the detection with ROSAT of $\sim 10$ sources deg$^{-2}$ above a threshold of $\sim 10^{-14} \, {\rm erg\, cm^{-2} s^{-1}}$. A brief reconsideration of the observability of isolated ONSs as individual X--ray sources is also presented.

General--relativistic, frequency--dependent radiative transfer in spherical, differentially--moving media is considered. In particular we investigate the character of the differential operator defined by the first two moment equations in the stationary case. We prove that the moment equations form a hyperbolic system when the logarithmic velocity gradient is positive, provided that a reasonable condition on the Eddington factors is met. The operator, however, may become elliptic in accretion flows and, in general, when gravity is taken into account. Finally we show that, in an optically thick medium, one of the characteristics becomes infinite when the flow velocity equals $\pm c/\sqrt 3$. Both high--speed, stationary inflows and outflows may therefore contain regions which are ``causally'' disconnected.

The spectral properties of X--ray radiation produced in a static atmosphere around a neutron star accreting at very low rates are investigated. Previous results by Alme \& Wilson (1973) are extended to the range $10^{-7}\leq L/L_{Edd}\leq 10^{-3}$ to include the typical luminosities, $L\sim 10^{31}-10^{32} \ {\rm ergs\, s^{-1}}$, expected from isolated neutron stars accreting the interstellar medium. The emergent spectra show an overall hardening with respect to the blackbody at the neutron star effective temperature in addition to a significant excess over the Wien tail. The relevance of present results in connection with the observability of low--luminosity X--ray sources is briefly discussed.

A possible sound origin for the photospheric oscillations in the X--ray bursting sources 1608-522 and 2127+119 is suggested. It is shown that standing sound waves in an expanding spherical envelope can have periods very close to the observed ones. The quite large ratio, $\sim$ 10, of the periods in the two sources is explained in terms of different wave regimes. The relevance of sound oscillations to the observed QPO in type II bursts of the Rapid Burster is also discussed.

We present estimates for the accretion rates in 13 X--ray bursting sources which exhibit photospheric expansion, basing on theoretical models of stationary, radiatively driven winds from neutron stars. The relatively high values obtained, $\Mdot_{acc}\magcir 10^{-9} \MS$, are in accordance with theoretical limits for unstable helium burning, and, at the same time, almost never exceed the ``dynamical'' limit for stationary accretion, $\sim 10 \Mdot_{Edd}$. The only exceptions are 1820-30, already known to be a very peculiar object, and 1608-522; there are indications, however, that in both sources, accretion could be non--stationary.

A model for stationary, radiatively driven winds from X--ray bursting neutron stars is presented. General relativistic hydrodynamical and radiative transfer equations are integrated from the neutron star surface outwards, taking into account for helium nuclear burning in the inner, dense, nearly hydrostatic shells. Radiative processes include both bremsstrahlung emission--absorption and Compton scattering; only the frequency--integrated transport is considered here. It is shown that each solution is characterized by just one parameter: the mass loss rate $\Mdot$, or, equivalently, the envelope mass $\Menv$. We found that, owing to the effects of Comptonization, steady, supersonic winds can exist only for $\Mdot$ larger than a limiting value $\Mdot_{min} \approx\Mdot_{E}$. Several models, covering about two decades in mass loss rate, have been computed for given neutron star parameters. We discuss how the sequence of our solutions with decreasing $\Menv$ can be used to follow the time evolution of a strong X--ray burst during the expansion/contraction phase near to the luminosity maximum. The comparison between our numerical results and the observational data of Haberl {\it et al.\/} (1987) for the bursts in 4U/MXB 1820-30 gives an estimate for both the spectral hardening factor and the accretion rate in this source.

Stationary, spherical accretion onto an unmagnetized neutron star is here reconsidered on the wake of the seminal paper by Zel'dovich \& Shakura (1969). It is found that new ``hot'' solutions may exist for a wide range of luminosities. These solutions are characterized by a high temperature, $10^{9}\div 10^{11}$ K, and arise from a stationary equilibrium model where the dominant radiative mechanisms are multiple Compton scattering and bremsstrahlung emission. For low luminosities, $\mincir 10^{-2} \ L_{E}$, only the ``cold'' (\`a la Zel'dovich and Shakura) solution is present.

We investigate the effects of dynamical Comptonization on the emergent radiation spectrum produced by near critical accretion onto a neutron star. The flow dynamics and the transfer of radiation are self--consistently solved in the case of a spherically symmetric, ``cold'', pure scattering flow, including general relativity. A sequence of models, each characterized by the value of the total observed luminosity, was obtained assuming that the spectrum at the star surface is black body in shape. It is found that when the luminosity approaches the Eddington limit dynamical effects become important shifting the spectrum to the blue and producing a power--law, high--energy tail. The relevance of these results in connection with the observed spectral properties of LMXBs, and of Cyg X--2 in particular, are discussed.