A few years after its discovery as a magnetar, SGR 1935+2154 started a new burst-active phase on 2020 April 27, accompanied by a large enhancement of its X-ray persistent emission. Radio single bursts were detected during this activation, strengthening the connection between magnetars and fast radio bursts. We report on the X-ray monitoring of SGR 1935+2154 from ~3 days prior to ~3 weeks after its reactivation, using Swift, NuSTAR, and NICER. We detected X-ray pulsations in the NICER and NuSTAR observations, and constrained the spin period derivative to |Pdot| < 6e-11 s/s (3 sigma c.l.). The pulse profile showed a variable shape switching between single and double-peaked as a function of time and energy. The pulsed fraction decreased from ~34% to ~11% (5-10 keV) over ~10 days. The X-ray spectrum was well fit by an absorbed blackbody model with temperature decreasing from kT ~ 1.6 to 0.6-0.7 keV, plus a non-thermal component (Gamma ~ 1.2) observed up to ~25 keV with NuSTAR. The 0.3-10 keV X-ray luminosity increased in less than four days from ~ 4e33 erg/s to about 2.5e35 erg/s and then decreased again to 1.4e34 erg/s over the following three weeks of the outburst. We also detected several X-ray bursts, with properties typical of short magnetar bursts.

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.

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.

We report on the first NuSTAR observation of the transitional millisecond pulsar binary XSS J12270-4859 during its current rotation-powered state, complemented with a 2.5yr-long radio monitoring at Parkes telescope and archival XMM-Newton and Swift X-ray and optical data. The radio pulsar is mainly detected at 1.4GHz displaying eclipses over about 40% of the 6.91h orbital cycle. We derive a new updated radio ephemeris to study the 3-79keV light curve that displays a significant orbital modulation with fractional amplitude of 28+/-3%, a structured maximum centred at the inferior conjunction of the pulsar and no cycle-to-cycle or low-high-flaring mode variabilities. The average X-ray spectrum, extending up to about 70keV without a spectral break, is well described by a simple power-law with photon index Gamma = 1.17+/-0.08 giving a 3-79keV luminosity of 7.6(-0.8;+3.8)x10**32 erg/s, for a distance of 1.37(-0.15;+0.69)kpc. Energy resolved orbital light curves reveal that the modulation is not energy dependent from 3keV to 25keV and is undetected with an upper limit of about 10% above 25keV. Comparison with previous X-ray XMM-Newton observations in common energy ranges confirms that the modulation amplitudes vary on timescales of a few months, indicative of a non-stationary contribution of the intrabinary shock formed by the colliding winds of the pulsar and the companion. A more detailed inspection of energy resolved modulations than previously reported gives hints of a mild softening at superior conjunction of the pulsar below 3keV, likely due to the contribution of the thermal emission from the neutron star. The intrabinary shock emission, if extending into the MeV range, would be energetically capable alone to irradiate the donor star.

We present the evolution of the X-ray emission properties of the magnetar 1E 1547.0-5408 since February 2004 over a time period covering three outbursts. We analyzed new and archival observations taken with the Swift, NuSTAR, Chandra and XMM-Newton X-ray satellites. The source has been observed at a relatively steady soft X-ray flux of $\approx$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (0.3-10 keV) over the last 9 years, which is about an order of magnitude fainter than the flux at the peak of the last outburst in 2009, but a factor of $\sim$ 30 larger than the level in 2006. The broad-band spectrum extracted from two recent NuSTAR observations in April 2016 and February 2019 showed a faint hard X-ray emission up to $\sim$ 70 keV. Its spectrum is adequately described by a flat power law component, and its flux is $\sim$ $7 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (10-70 keV), that is a factor of $\sim$ 20 smaller than at the peak of the 2009 outburst. The hard X-ray spectral shape has flattened significantly in time, which is at variance with the overall cooling trend of the soft X-ray component. The pulse profile extracted from these NuSTAR pointings displays variability in shape and amplitude with energy (up to $\approx$ 25 keV). Our analysis shows that the flux of 1E 1547.0-5408 is not yet decaying to the 2006 level and that the source has been lingering in a stable, high-intensity state for several years. This might suggest that magnetars can hop among distinct persistent states that are probably connected to outburst episodes and that their persistent thermal emission can be almost entirely powered by the dissipation of currents in the corona.

We report on a detailed spectral characterization of the non-thermal X-ray emission for a large sample of gamma-ray pulsars in the second Fermi-LAT catalogue. We outline the criteria adopted for the selection of our sample, its completeness, and critically describe different approaches to estimate the spectral shape and flux of pulsars. We perform a systematic modelling of the pulsars' X-ray spectra using archival observations with XMM-Newton, Chandra, and NuSTAR and extract the corresponding non-thermal X-ray spectral distributions. This set of data is made available online and is useful to confront with predictions of theoretical models.

We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $\Gamma = 2.06\pm 0.14$ and the second has $\Gamma = 2.03\pm 0.27$, compared to $\Gamma = 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $\nu \approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region.

We report on a NIR, optical and X-ray campaign performed in 2017 with the XMM-Newton and Swift satellites and the VLT/HAWK-I instrument on the transitional MSP PSR J1023+0038. NIR observations were performed in fast-photometric mode in order to detect any fast variation of the flux and correlate them with the optical and X-ray light curves. The optical curve shows the typical sinusoidal modulation at the orbital period (4.75hr). No flaring or flickering is found in the optical, neither signs of transitions between active and passive states. On the contrary, the NIR curve displays a bimodal behaviour, showing strong flares in the first part of the curve, and an almost flat trend in the rest. The X-ray curves show a few low/high mode transitions, but no flaring activity. One of the low/high mode transition is found to happen at the same time as the emission of an infrared flare. This can be interpreted as the emission of a jet: the NIR flare could be due to the evolving spectrum of the jet, which possesses a break frequency that moves from higher (NIR) to lower (radio) frequencies after the launching, that has to happen at the low/high mode transition. We also present the cross correlation function between the optical and near infrared curves. Due to the bimodality of the NIR curve, we divided it in two parts (flaring and quiet). While the CCF of the quiet part is found to be flat, the one referring to the flaring part shows a narrow peak at ~10s, which indicates a delay of the NIR emission with respect to the optical. This lag can be interpreted as reprocessing of the optical emission at the light cylinder radius with a stream of matter spiraling around the system due to a phase of radio-ejection. This strongly supports a different origin of the NIR flares observed for PSR J1023+0038 with respect to the optical and X-ray flaring activity reported in other works on the same source.

We apply a synchro-curvature spectral emission model based on characterizing the dynamics of magnetospheric particles to fit the phase-average spectra of the most extended database for the non-thermal spectra of pulsars. We consider 36 pulsars with well-determined non-thermal spectra from X-rays to gamma-rays. The sample includes Crab and Crab twin, for which the spectra extends even to the optical/ultraviolet and infrared energies. We find that the model --with just three physical parameters and a global scaling-- can fit the observations well across eight orders of magnitude for 18 of the 36 pulsars studied. Additionally, we find a set of 8 pulsars for which the model still provides arguably good fits and another set of 10 pulsars for which the model fails in reproducing the spectra. We discuss why, propose and provide physical interpretations for a simple model extension (related to the geometry of the accelerating system with regards to the observer) that allows dealing with all such cases, ultimately providing very good fits for all pulsars. The extended model is still austere, adding only two additional parameters to the former set, of the same kind of the ones previously used. We use these fits to discuss issues going from the observed spectral origin, to the extent of the dominance of synchrotron or curvature regimes, the use of a model as predictor for searching new non-thermal pulsars starting from gamma-ray surveys, and how the model offers a setting where phase shifts between X-ray and gamma-ray light curves would naturally arise.

Transitional pulsars provide us with a unique laboratory to study the physics of accretion onto a magnetic neutron star. PSR J1023+0038 (J1023) is the best studied of this class. We investigate the X-ray spectral properties of J1023 in the framework of a working radio pulsar during the active state. We modelled the X-ray spectra in three modes (low, high, and flare) as well as in quiescence, to constrain the emission mechanism and source parameters. The emission model, formed by an assumed pulsar emission (thermal and magnetospheric) plus a shock component, can account for the data only adding a hot dense absorber covering ~30% of the emitting source in high mode. The covering fraction is similar in flaring mode, thus excluding total enshrouding, and decreases in the low mode despite large uncertainties. This provides support to the recently advanced idea of a mini-pulsar wind nebula (PWN), where X-ray and optical pulsations arise via synchrotron shock emission in a very close (~100 km, comparable to a light cylinder), PWN-like region that is associated with this hot absorber. In low mode, this region may expand, pulsations become undetectable, and the covering fraction decreases.

PSR J1023+0038 is the first millisecond pulsar discovered to pulsate in the visible band; such a detection took place when the pulsar was surrounded by an accretion disk and also showed X-ray pulsations. We report on the first high time resolution observational campaign of this transitional pulsar in the disk state, using simultaneous observations in the optical (TNG, NOT, TJO), X-ray (XMM-Newton, NuSTAR, NICER), infrared (GTC) and UV (Swift) bands. Optical and X-ray pulsations were detected simultaneously in the X-ray high intensity mode in which the source spends $\sim$ 70% of the time, and both disappeared in the low mode, indicating a common underlying physical mechanism. In addition, optical and X-ray pulses were emitted within a few km, had similar pulse shape and distribution of the pulsed flux density compatible with a power-law relation $F_{\nu} \propto \nu^{-0.7}$ connecting the optical and the 0.3-45 keV X-ray band. Optical pulses were detected also during flares with a pulsed flux reduced by one third with respect to the high mode; the lack of a simultaneous detection of X-ray pulses is compatible with the lower photon statistics. We show that magnetically channeled accretion of plasma onto the surface of the neutron star cannot account for the optical pulsed luminosity ($\sim 10^{31}$ erg/s). On the other hand, magnetospheric rotation-powered pulsar emission would require an extremely efficient conversion of spin-down power into pulsed optical and X-ray emission. We then propose that optical and X-ray pulses are instead produced by synchrotron emission from the intrabinary shock that forms where a striped pulsar wind meets the accretion disk, within a few light cylinder radii away, $\sim$ 100 km, from the pulsar.

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.

In a recent work, we numerically studied the radiative properties of the reverberation phase of pulsar wind nebulae (PWNe), i.e., when the reverse shock created by the supernova explosion travels back towards the pulsar, compressing the wind bubble. We focused on several well-characterized PWNe and used them as examples for introducing the concept of superefficiency. The latter is a period of the PWN evolution, happening within reverberation, where the luminosity in a given band exceeds the spin-down power at the time. Here, we explore a broad range of PWN models to study their reverberation and superefficiency phases in a systematic way. Armed with these models we consider two aspects: On the one hand, we analyze via Monte Carlo simulations how many Galactic PWNe are expected to be reverberating or in a superefficiency stage at any given time, providing the first such estimations. On the other hand, we focus on searching for observational signatures of such periods. We analyze archival observations and check for the existence of possible candidates for superefficient PWNe. We also provide predictions for the future evolution of the magnetar nebula J1834.9-0846 (which we consider to be starting its reverberation period) along the next 50 years. Using our simulations as input we study how sensitive current and future X-ray satellites (like eXTP or Athena) will be to observe such evolution, concluding that they will be able to track it in detail.

We report on a multi-wavelength study of the unclassified X-ray source CXOU J110926.4-650224 (J1109). We identified the optical counterpart as a blue star with a magnitude of $\sim$20.1 (3300-10500 $\require{mediawiki-texvc} \AA$). The optical emission was variable on timescales from hundreds to thousands of seconds. The spectrum showed prominent emission lines with variable profiles at different epochs. Simultaneous XMM-Newton and NuSTAR observations revealed a bimodal distribution of the X-ray count rates on timescales as short as tens of seconds, as well as sporadic flaring activity. The average broad-band (0.3-79 keV) spectrum was adequately described by an absorbed power law model with photon index of $\Gamma$=1.63$\pm$0.01 (at 1$\sigma$ c.l.), and the X-ray luminosity was (2.16$\pm$0.04)$\times$10$^{34}$ erg s$^{-1}$ for a distance of 4 kpc. Based on observations with different instruments, the X-ray luminosity has remained relatively steady over the past $\sim$15 years. J1109 is spatially associated with the gamma-ray source FL8Y J1109.8-6500, which was detected with Fermi at an average luminosity of (1.5$\pm$0.2)$\times$10$^{34}$ erg s$^{-1}$ (assuming the distance of J1109) over the 0.1-300 GeV energy band between 2008 and 2016. The source was undetected during ATCA radio observations that were simultaneous with NuSTAR, down to a 3$\sigma$ flux upper limit of 18 $\mu$Jy/beam (at 7.25 GHz). We show that the phenomenological properties of J1109 point to a binary transitional pulsar candidate currently in a sub-luminous accretion disk state, and that the upper limits derived for the radio emission are consistent with the expected radio luminosity for accreting neutron stars at similar X-ray luminosities.

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.

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.

Based on a theoretical selection of pulsars as candidates for detection at X-ray energies, we present an analysis of archival X-ray observations done with Chandra and XMM-Newton of PSR J1747-2958 (the pulsar in the "Mouse" nebula), PSR J2021+3651 (the pulsar in the "Dragonfly" nebula), and PSR J1826-1256. X-ray pulsations from PSR J1747-2958 and PSR J1826-1256 are detected for the first time, and a previously reported hint of an X-ray pulsation from PSR J2021+3651 is confirmed with a higher significance. We analyze these pulsars' spectra in regards to the theoretically predicted energy distribution, finding a remarkable agreement, and provide here a refined calculation of the model parameters taking into account the newly derived X-ray spectral data.

We report on the results of optical, near-infrared (NIR) and mid-infrared observations of the black hole X-ray binary candidate (BHB) MAXI J1535-571 during its 2017/2018 outburst. During the first part of the outburst (MJD 58004-58012), the source shows an optical-NIR spectrum that is consistent with an optically thin synchrotron power-law from a jet. After MJD 58015, however, the source faded considerably, the drop in flux being much more evident at lower frequencies. Before the fading, we measure a de-reddened flux density of $\gtrsim$100 mJy in the mid-infrared, making MAXI J1535-571 one of the brightest mid-infrared BHBs known so far. A significant softening of the X-ray spectrum is evident contemporaneous with the infrared fade. We interpret it as due to the suppression of the jet emission, similar to the accretion-ejection coupling seen in other BHBs. However, MAXI J1535-571 did not transition smoothly to the soft state, instead showing X-ray hardness deviations, associated with infrared flaring. We also present the first mid-IR variability study of a BHB on minute timescales, with a fractional rms variability of the light curves of $\sim 15-22 \%$, which is similar to that expected from the internal shock jet model, and much higher than the optical fractional rms ($\lesssim 7 \%$). These results represent an excellent case of multi-wavelength jet spectral-timing and demonstrate how rich, multi-wavelength time-resolved data of X-ray binaries over accretion state transitions can help refining models of the disk-jet connection and jet launching in these systems.

We report on the first continuous, 80 day optical monitoring of the transitional millisecond pulsar PSR J1023+0038 carried out in mid-2017 with Kepler in the K2 configuration, when an X-ray subluminous accretion disk was present in the binary. Flares lasting from minutes to 14 hr were observed for 15.6% of the time, which is a larger fraction than previously reported on the basis of X-ray and past optical observations, and more frequently when the companion was at the superior conjunction of the orbit. A sinusoidal modulation at the binary orbital period was also present with an amplitude of ~16%, which varied by a few percent over timescales of days, and with a maximum that took place 890 +/- 85 s earlier than the superior conjunction of the donor. We interpret these phenomena in terms of reprocessing of the X-ray emission by an asymmetrically heated companion star surface and/or a non-axisymmetric outflow possibly launched close to the inner Lagrangian point. Furthermore, the non-flaring average emission varied by up to ~ 40% over a time scale of days in the absence of correspondingly large variations of the irradiating X-ray flux. The latter suggests that the observed changes in the average optical luminosity might be due to variations of the geometry, size, and/or mass accretion rate in the outer regions of the accretion disk.

We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $\Gamma$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected.

We report on the first simultaneous XMM-Newton, NuSTAR and Swift observations of the transitional millisecond pulsar PSR J1023+0038 in the X-ray active state. Our multi-wavelength campaign allowed us to investigate with unprecedented detail possible spectral variability over a broad energy range in the X-rays, as well as correlations and lags among emissions in different bands. The soft and hard X-ray emissions are significantly correlated, with no lags between the two bands. On the other hand, the X-ray emission does not correlate with the UV emission. We refine our model for the observed mode switching in terms of rapid transitions between a weak propeller regime and a rotation-powered radio pulsar state, and report on a detailed high-resolution X-ray spectroscopy using all XMM-Newton Reflection Grating Spectrometer data acquired since 2013. We discuss our results in the context of the recent discoveries on the system and of the state of the art simulations on transitional millisecond pulsars, and show how the properties of the narrow emission lines in the soft X-ray spectrum are consistent with an origin within the accretion disc.

We present the results of the systematic study of all magnetar outbursts observed to date, through a reanalysis of data acquired in about 1100 X-ray observations. We track the temporal evolution of the outbursts soft X-ray spectral properties and the luminosities of the single spectral components as well as of the total emission. We model empirically all outburst light curves, and estimate the characteristic decay time-scales as well as the energetics involved. We investigate the link between different parameters (e.g. the luminosity at the peak of the outburst and in quiescence, the maximum luminosity increase, the decay time-scale and energy of the outburst, the neutron star surface dipolar magnetic field and characteristic age, etc.), and unveil several correlations among these quantities. We discuss our results in the context of the internal crustal heating and twisted bundle models for magnetar outbursts. This study is complemented by the Magnetar Outburst Online Catalogue (http://magnetars.ice.csic.es), an interactive data base where the user can plot any combination of the parameters derived in this work, and download all data.

PSR J1023+0038 (J1023) is a binary system hosting a neutron star and a low mass companion. J1023 is the best studied transitional pulsar, alternating a faint eclipsing millisecond radio pulsar state to a brighter X-ray active state. At variance with other Low Mass X-ray binaries, this active state reaches luminosities of only ~$10^{34}$ erg s$^{-1}$, showing strong, fast variability. In the active state, J1023 displays: i) a high state ($L_X\sim7\times10^{33}$ erg s$^{-1}$, 0.3-80 keV) occurring ~80% of the time and during which X-ray pulsations at the neutron star spin period are detected (pulsed fraction ~8%); ii) a low state ($L_X~10^{33}$ erg s$^{-1}$) during which pulsations are not detected (~<3%); and iii) a flaring state during which sporadic flares occur in excess of ~$10^{34}$ erg s$^{-1}$, with no pulsation too. The transition between the high and the low states is very rapid, on a ~10 s timescale. Here we put forward a plausible physical interpretation of the high and low states based on the (fast) transition among the propeller state and the radio pulsar state. We modelled the XMM-Newton spectra of the high, low and radio pulsar states, finding a good agreement with this physical picture.

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.

In search for the counterpart to the Fermi-LAT source 3FGL J0838.8-2829, we performed a multi-wavelength campaign, in the X-ray band with Swift and XMM-Newton, performed infrared, optical (with OAGH, ESO-NTT and IAC80) and radio (ATCA) observations, as well as analysed archival hard X-ray data taken by INTEGRAL. We report on three X-ray sources consistent with the position of the Fermi-LAT source. We confirm the identification of the brightest object, RX J0838-2827, as a magnetic cataclysmic variable, that we recognize as an asynchronous system (not associated with the Fermi-LAT source). RX J0838-2827, is extremely variable in the X-ray and optical bands, and timing analysis reveals the presence of several periodicities modulating its X-ray and optical emission. The most evident modulations are interpreted as due to the binary system orbital period of ~1.64h and the white dwarf spin period of ~1.47h. A strong flux modulation at ~15h is observed at all energy bands, consistent with the beat frequency between spin and orbital periods. Optical spectra show prominent Hbeta, HeI and HeII emission lines Doppler-modulated at the orbital period and at the beat period. Therefore, RX J0838-2827, accretes through a diskless configuration and could be either a strongly asynchronous polar or a rare example of a pre-polar system in its way to reach synchronism. Among the other two X-ray sources, XMM J083850.4-282759 showed a variable X-ray emission, with a powerful flare lasting ~600s, similar to what is observed in transitional millisecond pulsars during the sub-luminous disc state: that would possibly associate this source with the Fermi-LAT source.

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.

The 6.67 hr periodicity and the variable X-ray flux of the central compact object (CCO) at the center of the SNR RCW 103, named 1E 161348-5055, have been always difficult to interpret within the standard scenarios of an isolated neutron star or a binary system. On 2016 June 22, the Burst Alert Telescope (BAT) onboard Swift detected a magnetar-like short X-ray burst from the direction of 1E 161348-5055, also coincident with a large long-term X-ray outburst. Here we report on Chandra, NuSTAR, and Swift (BAT and XRT) observations of this peculiar source during its 2016 outburst peak. In particular, we study the properties of this magnetar-like burst, we discover a hard X-ray tail in the CCO spectrum during outburst, and we study its long-term outburst history (from 1999 to July 2016). We find the emission properties of 1E 161348-5055 consistent with it being a magnetar. However in this scenario, the 6.67 hr periodicity can only be interpreted as the rotation period of this strongly magnetized neutron star, which therefore represents the slowest pulsar ever detected, by orders of magnitude. We briefly discuss the viable slow-down scenarios, favoring a picture involving a period of fall-back accretion after the supernova explosion, similarly to what is invoked (although in a different regime) to explain the "anti-magnetar" scenario for other CCOs.

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 the identification of the new Galactic Center (GC) transient Swift J174540.7-290015 as a likely low mass X-ray binary (LMXB) located at only 16 arcsec from Sgr A*. This transient was detected on 2016 February 6th during the Swift GC monitoring, and it showed long-term spectral variations compatible with a hard to soft state transition. We observed the field with XMM-Newton on February 26th for 35 ks, detecting the source in the soft state, characterised by a low level of variability and a soft X-ray thermal spectrum with a high energy tail (detected by INTEGRAL up to ~50 keV), typical of either accreting neutron stars or black holes. We observed: i) a high column density of neutral absorbing material, suggesting that Swift J174540.7-290015 is located near or beyond the GC and; ii) a sub-Solar Iron abundance, therefore we argue that Iron is depleted into dust grains. The lack of detection of FeK absorption lines, eclipses or dipping suggests that the accretion disc is observed at a low inclination angle. Radio (VLA) observations did not detect any radio counterpart to Swift J174540.7-290015. No evidence for X-ray or radio periodicity is found. The location of the transient was observed also in the near-IR with GROND at MPG/ESO La Silla 2.2m telescope and VLT/NaCo pre- and post-outburst. Within the Chandra error region we find multiple objects that display no significant variations.

We present the first optical polarimetric study of the two transitional pulsars PSR J1023+0038 and XSS J12270-4859. This work is focused on the search for intrinsical linear polarisation (LP) in the optical emission from the two systems. We carried out multiband optical and NIR photo-polarimetry of the two systems using the ESO NTT at La Silla (Chile), equipped with the EFOSC2 and the SOFI instruments. XSS J12270-4859 was observed during its radio-pulsar state; we did not detect LP in all bands, with 3 sigma upper limits of, e.g., 1.4% in the R-band. We built the NIR-optical averaged spectral energy distribution (SED) of the system, that could be well described by an irradiated black body with radius $R_{*} = 0.33\pm0.03\,R_{\odot}$ and albedo $\eta=0.32\pm0.05$, without the need of further components (thus excluding the visible presence of an extended accretion disc and/or of relativistic jets). The case was different for PSR J1023+0038, that was in its accretion phase during our campaign. We measured a LP of $1.09\pm0.27\%$ and $0.90\pm 0.17\%$ in the V and R bands, respectively. The phase-resolved polarimetric curve of the source in the R-band reveals a hint of a sinusoidal modulation at the source 4.75 hr orbital period, peaked at the same orbital phase as the light curve. The measured LP of PSR J1023+0038 could in principle be interpreted as scattering with free electrons (that can be found in the accretion disc of the system or even in the hot corona that surrounds the disc itself) or to synchrotron emission from a relativistic particles jet or outflow. However, the NIR-optical SED of the system built starting from our dataset did not suggest the presence of a jet. We conclude that the optical LP observed for PSR J1023+0038 is possibly due to Thomson scattering with electrons in the disc, as also suggested from the possible modulation of the R-band LP at the system orbital period.

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.

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.

The X-ray source RX J2015.6+3711 was discovered by ROSAT in 1996 and recently proposed to be a cataclysmic variable (CV). Here we report on an XMM-Newton observation of RX J2015.6+3711 performed in 2014, where we detected a coherent X-ray modulation at a period of 7196+/-11 s, and discovered other significant (>6sigma) small-amplitude periodicities which we interpret as the CV spin period and the sidebands of a possible ~12 hr periodicity, respectively. The 0.3-10 keV spectrum can be described by a power law (Gamma = 1.15+/-0.04) with a complex absorption pattern, a broad emission feature at 6.60+/-0.01 keV, and an unabsorbed flux of (3.16+/-0.05)x10^{-12} erg/s/cm^2. We observed a significant spectral variability along the spin phase, which can be ascribed mainly to changes in the density of a partial absorber and the power law normalization. Archival X-ray observations carried out by the Chandra satellite, and two simultaneous X-ray and UV/optical pointings with Swift, revealed a gradual fading of the source in the soft X-rays over the last 13 years, and a rather stable X-ray-to-optical flux ratio (F_X/F_V ~1.4-1.7). Based on all these properties, we identify this source with a magnetic CV, most probably of the intermediate polar type. The 2 hr spin period makes RX J2015.6+3711 the second slowest rotator of the class, after RX J0524+4244 ("Paloma", P_spin~2.3 hr). Although we cannot unambiguously establish the true orbital period with these observations, RX J2015.6+3711 appears to be a key system in the evolution of magnetic CVs.

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.

We report on multi-frequency, wideband radio observations of the Galactic Center magnetar (SGR 1745$-$2900) with the Green Bank Telescope for $\sim$100 days immediately following its initial X-ray outburst in April 2013. We made multiple simultaneous observations at 1.5, 2.0, and 8.9 GHz, allowing us to examine the magnetar's flux evolution, radio spectrum, and interstellar medium parameters (such as the dispersion measure (DM), the scattering timescale and its index). During two epochs, we have simultaneous observations from the Chandra X-ray Observatory, which permitted the absolute alignment of the radio and X-ray profiles. As with the two other radio magnetars with published alignments, the radio profile lies within the broad peak of the X-ray profile, preceding the X-ray profile maximum by $\sim$0.2 rotations. We also find that the radio spectral index $\gamma$ is significantly negative between $\sim$2 and 9 GHz; during the final $\sim$30 days of our observations $\gamma \sim -1.4$, which is typical of canonical pulsars. The radio flux has not decreased during this outburst, whereas the long-term trends in the other radio magnetars show concomitant fading of the radio and X-ray fluxes. Finally, our wideband measurements of the DMs taken in adjacent frequency bands in tandem are stochastically inconsistent with one another. Based on recent theoretical predictions, we consider the possibility that the dispersion measure is frequency-dependent. Despite having several properties in common with the other radio magnetars, such as $L_{\textrm{X,qui}}/L_{\textrm{rot}} \lesssim 1$, an increase in the radio flux during the X-ray flux decay has not been observed thus far in other systems.

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.

The binary millisecond radio pulsar PSR J1023+0038 has been recently the subject of multiwavelength monitoring campaigns which revealed that an accretion disc has formed around the neutron star (since 2013 June). We present here the results of X-ray and UV observations carried out by the Swift satellite between 2013 October and 2014 May, and of optical and NIR observations performed with the REM telescope, the Liverpool Telescope, the 2.1-m telescope at the San Pedro M\'artir Observatory and the 1.52-m telescope at the Loiano observing station. The X-ray spectrum is well described by an absorbed power law, which is softer than the previous quiescent epoch (up to 2013 June). The strong correlation between the X-ray and the UV emissions indicates that the same mechanism should be responsible for part of the emission in these bands. Optical and infrared photometric observations show that the companion star is strongly irradiated. Double-peaked emission lines in the optical spectra provide compelling evidence for the presence of an outer accretion disc too. The spectral energy distribution from IR to X-rays is well modelled if the contributions from the companion, the disc and the intra-binary shock emission are all considered. Our extensive data set can be interpreted in terms of an engulfed radio pulsar: the radio pulsar is still active, but undetectable in the radio band due to a large amount of ionized material surrounding the compact object. X-rays and gamma-rays are produced in an intra-binary shock front between the relativistic pulsar wind and matter from the companion and an outer accretion disc. The intense spin-down power irradiates the disc and the companion star, accounting for the UV and optical emissions.

The Swift satellite monitored the quiescence of the low-mass X-ray binary transient Aql X-1 on a weekly basis during the March-November 2012 interval. A total of 42 observations were carried out in the soft X-ray (0.3-10 keV) band with the X-ray Telescope on board Swift. We investigated the X-ray variability properties of Aql X-1 during quiescence by tracking luminosity variations and characterising them with a detailed spectral analysis. The source is highly variable in this phase and two bright flares were detected, with peak luminosities of ~4x10^{34} erg s^{-1} (0.3-10 keV). Quiescent X-ray spectra require both a soft thermal component below ~2 keV and a hard component (a power law tail) above ~2 keV. Changes in the power law normalisation alone can account for the overall observed variability. Therefore, based on our data set, the quiescent X-ray emission of Aql X-1 is consistent with the cooling of the neutron star core and with mechanisms involving the accretion of matter onto the neutron star surface or magnetosphere.

Aql X-1 is the prototypical low mass X-ray binary transient. The Rossi X-ray Timing Explorer All Sky Monitor provided a ~16 yr coverage revealing 20 outbursts. This is by far the most extensive legacy of outbursts from the same source. We investigated the outbursts characteristics in terms of energetics, peak luminosities, durations, decays and recurrence times. We found that bright outbursts (peak luminosity >10^{37} erg s^-1) equal in number dimmer outbursts (<10^{36.6} erg s^-1). The peak luminosity does not correlate with outburst energetics, durations or quiescent times. We analysed the latest stages of the outbursts searching for exponential and/or linear decays. Light curve modeling led to constraints on the outer disk radius and enabled us to estimate the viscosity and the irradiation parameters. The former is larger than what has been obtained for other, shorter orbital period, transients, while the latter is somewhat smaller. This might be related to the longer orbital period of Aql X-1 with respect to other transient X-ray binaries.