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 X-ray emission of gamma-ray bursts (GRBs) is often characterized by an initial steep decay, followed by a nearly constant emission phase (so called "plateau") which can extend up to thousands of seconds. While the steep decay is usually interpreted as the tail of the prompt gamma-ray flash, the long-lasting plateau is commonly associated to the emission from the external shock sustained by energy injection from a long lasting central engine. A recent study proposed an alternative interpretation, ascribing both the steep decay and the plateau to high-latitude emission (HLE) from a "structured jet" whose energy and bulk Lorentz factor depend on the angular distance from the jet symmetry axis. In this work we expand over this idea and explore more realistic conditions: (a) the finite duration of the prompt emission, (b) the angular dependence of the optical depth and (c) the lightcurve dependence on the observer viewing angle. We find that, when viewed highly off-axis, the structured jet HLE lightcurve is smoothly decaying with no clear distinction between the steep and flat phase, as opposed to the on-axis case. For a realistic choice of physical parameters, the effects of a latitude-dependent Thomson opacity and finite duration of the emission have a marginal effect on the overall lightcurve evolution. We discuss the possible HLE of GW170817, showing that the emission would have faded away long before the first Swift-XRT observations. Finally, we discuss the prospects for the detection of HLE from off-axis GRBs by present and future wide-field X-ray telescopes and X-ray surveys, such as eROSITA and the mission concept THESEUS.

The first multi-messenger detection of a binary neutron star merger, GW170817, brought to the forefront the structured jet model as a way to explain multi-wavelength observations taken more than a year after the event. Here we show that the high-latitude emission from a structured jet can naturally produce an X-ray plateau in gamma-ray burst (GRB) light curves, independent of the radiation from an external shock. We calculate the radiation from a switched-off shell featuring an angular structure in both its relativistic bulk motion and intrinsic brightness. Our model is able to explain the shallow decay phase (plateau) often observed in GRB X-ray light curves. We discuss the possible contribution of the structured jet high-latitude emission to other distinctive features of GRB X-ray light curves, and its capability to explain the chromatic optical/X-ray light curve properties.

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars.

The X-ray afterglow plateau emission observed in many Gamma-ray Bursts (GRBs) has been interpreted as either being fueled by fallback onto a newly formed black hole, or by the spin-down luminosity of an ultra-magnetized millisecond neutron star. If the latter model is assumed, GRB X-ray afterglow light curves can be analytically reproduced. We fit a sample of GRB X-ray plateaus, interestingly yielding a distribution in the magnetic field versus spin period (B-P) diagram consistent with $B\propto P^{7/6}$. This is expected from the well-established physics of the spin-up line minimum period for Galactic millisecond pulsars. The normalisation of the relation we obtain perfectly matches spin-up line predictions for the expected masses ($\sim 1 M_{\odot}$) and radii ($\sim 10 {\rm ~km}$) of newly born magnetars, and mass accretion rates consistent with GRB expectations of $10^{-4} M_{\odot}/{\rm s} <\dot{M}< 10^{-1} M_{\odot}/{\rm s}$. Short GRBs with extended emission (SEE) appear towards the high period end of the distribution, while the long GRBs (LGRBs) towards the short period end. This result is consistent with spin-up limit expectations where the total accreted mass determines the position of the neutron star in the B-P diagram. The P-B distribution for LGRBs and SEE are statistically different, further supporting the idea that the fundamental plane relation \citep{dainotti16c,Dainotti2017} is a powerful discriminant among those populations. Our conclusions are robust against suppositions regarding the GRB collimation angle and magnetar breaking index, which shifts the resulting magnetar properties parallel to the spin-up line, and strongly support a magnetar origin for GRBs presenting X-ray plateaus.

The viscosity-driven "spin-flip" instability in newborn magnetars with interior toroidal magnetic fields is re-examined. We calculate the bulk viscosity coefficient ($\zeta$) of cold, $npe \mu$ matter in neutron stars (NS), for selected values of the nuclear symmetry energy and in the regime where $\beta$-equilibration is slower than characteristic oscillation periods. We show that: i) $\zeta$ is larger than previously assumed and the instability timescale correspondingly shorter; ii) for a magnetically-induced ellipticity $\epsilon_B \lesssim 4 \times 10^{-3}$, typically expected in newborn magnetars, spin-flip occurs for initial spin periods $\lesssim 2-3$ ms, with some dependence on the NS equation of state (EoS). We then calculate the detectability of GW signals emitted by newborn magnetars subject to "spin-flip", by accounting also for the reduction in range resulting from realistic signal searches. For an optimal range of $\epsilon_B \sim (1-5) \times 10^{-3}$, and birth spin period $\lesssim 2$ ms, we estimate an horizon of $\gtrsim 4$ Mpc, and $\gtrsim 30$ Mpc, for Advanced and third generation interferometers at design sensitivity, respectively. A supernova (or a kilonova) is expected as the electromagnetic counterpart of such GW events. Outside of the optimal range for GW emission, EM torques are more efficient in extracting the NS spin energy, which may power even brighter EM transients.

We study how the viscosity of neutron star (NS) matter affects the distribution of tilt angles ($\chi$) between the spin and magnetic axes in young pulsars. Under the hypothesis that the NS shape is determined by the magnetically-induced deformation, and that the toroidal component of the internal magnetic field exceeds the poloidal one, we show that the dissipation of precessional motions by bulk viscosity can naturally produce a bi-modal distribution of tilt angles, as observed in radio/$\gamma$-ray pulsars, with a low probability of achieving $\chi \sim (20^\circ - 70^\circ)$ if the interior B-field is $\sim (10^{11} - 10^{15})$~G and the birth spin period is $\sim 10 - 300$~ms. As a corollary of the model, the idea that the NS shape is solely determined by the poloidal magnetic field, or by the centrifugal deformation of the crust, is found to be inconsistent with the tilt angle distribution in young pulsars. When applied to the Crab pulsar, with $\chi \sim 45^\circ - 70^\circ$ and birth spin $\sim$ 20 ms, our model implies that: (i) the magnetically-induced ellipticity is $\epsilon_B \gtrsim 3 \times 10^{-6}$; (ii) the measured positive $\dot{\chi} \sim 3.6 \times 10^{-12}$ rad s$^{-1}$ requires an additional viscous process, acting on a timescale $\lesssim 10^4$ yrs. We interpret the latter as crust-core coupling via mutual friction in the superfluid NS interior. One critical implication of our model is a GW signal at (twice) the spin frequency of the NS, due to $\epsilon_B \sim 10^{-6}$. This could be detectable by Advanced LIGO/Virgo operating at design sensitivity.

Flaring activity following gamma-ray bursts (GRBs), observed in both long and short GRBs, signals a long-term activity of the central engine. However, its production mechanism has remained elusive. Here we develop a quantitative model of the idea proposed by Perna et al. of a disc whose outer regions fragment due to the onset of gravitational instability. The self-gravitating clumps migrate through the disc and begin to evolve viscously when tidal and shearing torques break them apart. Our model consists of two ingredients: theoretical bolometric flare lightcurves whose shape (width, skewness) is largely insensitive to the model parameters, and a spectral correction to match the bandpass of the available observations, that is calibrated using the observed spectra of the flares. This simple model reproduces, with excellent agreement, the empirical statistical properties of the flares as measured by their width-to-arrival time ratio and skewness (ratio between decay and rise time). We present model fits to the observed lightcurves of two well-monitored flares, of GRB 060418 and of GRB 060904B. To the best of our knowledge, this is the first quantitative model able to reproduce the flare lightcurves and explain their global statistical properties.

We present the detection of persistent soft X-ray radiation with L_x ~ 10^41-10^42 erg/s at the location of the extremely luminous, double-humped transient ASASSN-15lh as revealed by Chandra and Swift. We interpret this finding in the context of observations from our multiwavelength campaign, which revealed the presence of weak narrow nebular emission features from the host-galaxy nucleus and clear differences with respect to superluminous supernova optical spectra. Significant UV flux variability on short time-scales detected at the time of the re-brightening disfavors the shock interaction scenario as the source of energy powering the long-lived UV emission, while deep radio limits exclude the presence of relativistic jets propagating into a low-density environment. We propose a model where the extreme luminosity and double-peaked temporal structure of ASASSN-15lh is powered by a central source of ionizing radiation that produces a sudden change of the ejecta opacity at later times. As a result, UV radiation can more easily escape, producing the second bump in the light-curve. We discuss different interpretations for the intrinsic nature of the ionizing source. We conclude that, IF the X-ray source is physically associated with the optical-UV transient, ASASSN-15lh most likely represents the tidal disruption of a main-sequence star by the most massive spinning black hole detected to date. In this case, ASASSN-15lh and similar events discovered in the future would constitute the most direct probes of very massive, dormant, spinning, supermassive black holes in galaxies. Future monitoring of the X-rays may allow us to distinguish between the supernova and TDE hypothesis.

We analyze archival Chandra HRC observations of the ultra luminous accreting pulsar M82-X2 (NuSTAR J095551+6940.8), and determine an upper limit of $< 1.7\times 10^{38}$~erg/s to its luminosity at an epoch at which it was undetected. Combined with other recent measurements, this confirms that the source X-ray emission has been highly variable during the last 15 years, ranging from a maximum of $10^{40}$ erg/s through intermediate values $\sim$ a few $\times 10^{39}$ erg/s, and down to a minimum that must be below the current detection threshold $\sim (2-3) \times 10^{38}$ erg/s . We interpret these results by means of a magnetically-threaded disk model: when at peak luminosity, the neutron star (NS) is close to spin equilibrium, its inner disk edge r_m ~ 10^8 cm is approximately half the corotation radius r_{co}, and radiation pressure dominates the disk out to r_{tr} ~ 10^9 cm. In the radiation pressure-dominated regime, r_m grows very slowly as the mass inflow rate drops: as a result, r_m < r_{co} remains valid until the mass accretion rate becomes ~ the Eddington accretion rate, allowing a wide range of accretion luminosities to the NS. Once the mass accretion rate is below Eddington, accretion onto the NS is inhibited because r_m > r_{co}, and the source luminosity is expected to drop by a large factor. We conclude that a magnetically threaded, radiation pressure-dominated disk, around a highly magnetized NS (B~10^{13} G) offers the best intepretation for all the currently observed properties of NuSTAR J095551+6940.8. This source offers an unprecedented opportunity to study the disk-magnetosphere interaction in a new regime of supercritical accretion, and across the transition between-radiation pressure and gas-pressure dominance inside the disk.

The identification of the Ultraluminous X-ray source (ULX) X-2 in M82 as an accreting pulsar has shed new light on the nature of a subset of ULXs, while rising new questions on the nature of the super-Eddington accretion. Here, by numerically solving the torque equation of the accreting pulsar within the framework of the magnetically threaded-disk scenario, we show that three classes of solutions, corresponding to different values of the magnetic field, are mathematically allowed. We argue that the highest magnetic field one, corresponding to B $\sim 10^{13}$ G, is favoured based on physical considerations and the observed properties of the source. In particular, that is the only solution which can account for the observed variations in $\dot{P}$ (over four time intervals) without requiring major changes in $\dot{M}$, which would be at odds with the approximately constant X-ray emission of the source during the same time. For this solution, we find that the source can only accomodate a moderate amount of beaming, 0.5 $\lesssim b < 1$. Last, we show that the upper limit on the luminosity, L$_X < 2.5 \times 10^{38}$ erg s$^{-1}$ from archival observations, is consistent with a highly-magnetized neutron star being in the propeller phase at that time.

Binary neutron star (NS) mergers are among the most promising sources of gravitational waves (GWs), as well as candidate progenitors for short Gamma-Ray Bursts (SGRBs). Depending on the total initial mass of the system, and the NS equation of state (EOS), the post-merger phase can see a prompt collapse to a black hole, or the formation of a supramassive NS, or even a stable NS. In the case of post-merger NS (PMNS) formation, magnetic field amplification during the merger will produce a magnetar with a large induced mass quadrupole moment, and millisecond spin. If the timescale for orthogonalization of the magnetic symmetry axis with the spin axis is sufficiently short the NS will radiate its spin down energy primarily via GWs. Here we study this scenario for various outcomes of NS formation: we generalise the set of equilibrium states for a twisted torus magnetic configuration to include solutions that, at a fixed exterior dipole field, carry a larger magnetic energy reservoir; we hence compute their magnetic ellipticity and the strength of the expected GW signal as a function of the magnitude of the dipole and toroidal field. The relative number of GW detections from PMNSs and from binary NSs is a strong function of the NS equation of state (EOS), being higher (~ 1%) for the stiffest EOSs and negligibly small for the softest ones. For intermediate-stiffness EOSs, such as the n=4/7 polytrope recently used by Giacomazzo \& Perna or the GM1 used by Lasky et al., the relative fraction is ~0.3%; correspondingly we estimate a GW detection rate from stable PMNSs of ~ (0.1-1) yr$^{-1}$ with Advanced detectors, and of ~ (100-1000) yr$^{-1}$ with third generation detectors such as the Einstein Telescope. Measurement of such GW signal would provide strong constraints on the NS EOS and on the nature of the binary progenitors giving rise to SGRBs.

The system RX J0806.3+1527 (HM Cnc) is a pulsating X-ray source with 100 per cent modulation on a period of 321.5 s (5.4 min). This period reflects the orbital motion of a close binary consisting of two interacting white dwarfs. Here we present a series of simultaneous X-ray (0.2-10 keV) and near-ultraviolet (2600 angstrom and 1928 angstrom) observations that were carried out with the Swift satellite. In the near-ultraviolet, the counterpart of RX J0806.3+1527 was detected at flux densities consistent with a blackbody with temperature 27E+3 K. We found that the emission at 2600 angstrom is modulated at the 321.5-s period with the peak ahead of the X-ray one by 0.28 cycles and is coincident within 0.05 cycles with the optical. This phase-shift measurement confirms that the X-ray hot spot (located on the primary white dwarf) is at about 80-100 degrees from the direction that connects the two white dwarfs. Albeit at lower significance, the 321.5-s signature is present also in the 1928-angstrom data; at this wavelength, however, the pulse peak is better aligned with that observed at X-rays. We use the constraints on the source luminosity and the geometry of the emitting regions to discuss the merits and limits of the main models for RX J0806.3+1527.

Although the internal structure of white dwarfs is considered to be generally well understood, the source and entity of viscosity is still very uncertain. We propose here to study white dwarf viscous properties using short period (< 1 hr), detached white dwarf binaries, such as the newly discovered ~12.8 min system. These binaries are wide enough that mass transfer has not yet started but close enough that the least massive component is subject to a measurable tidal deformation. The associated tidal torque transfers orbital energy, which is partially converted into heat by the action of viscosity within the deformed star. As a consequence, its outer non-degenerate layers expand, and the star puffs up. We self-consistently calculate the fractional change in radius, and the degree of asynchronism (ratio of stellar to orbital spin) as a function of the viscous time. Specializing our calculations to J0651, we find that the discrepancy between the measured radius of the secondary star and He white dwarf model predictions can be interpreted as tidal inflation if the viscous timescale is either ~2 10^5 yr or ~10^4 yr. Such values point to a non-microscopic viscosity, possibly given by tidally induced turbulence, or by magnetic field stresses with a magnetic field strength of 10-100 Gauss. Fortunately, these two timescales produce very different degree of asynchronism, with the shortest one, bringing the system much closer to synchronisation. A measurement of the stellar spin can thus univocally determined the mean viscosity. Extrapolating the secondary's radial expansion, we predict that the star will fill is Roche lobe at a separation which is 1.2-1.3 smaller than the current one. Applying this method to a future sample of systems can allow us to learn whether viscosity changes with mass and/or nuclear composition.

The binary pulsar J0737-3039 is the only known system having two observable pulsars, thus offering a unique laboratory to test general relativity and explore pulsar physics. Based on the low eccentricity and the position within the galactic plane, Piran & Shaviv (2004, 2005) argued that pulsar B had a non-standard formation scenario with little or no mass ejection. They have also predicted that the system would have a very slow proper motion. Pulsar timing measurements (Kramer et al. 2006; Deller et al. 2009) confirmed this prediction. The recent observations of the alignment between the spin of pulsar A and the binary orbit is also in agreement with this scenario. Detailed simulations of the formation process of pulsar B enable us to show that its progenitor, just before the collapse, was a massive O-Ne-Mg white dwarf surrounded by a tenuous, 0.1-0.16 M_sun, envelope. This envelope was ejected when the white dwarf collapsed to form a neutron star. Pulsar B was born as a slow rotator (spin period ~ 1 s) and a kick received when the pulsar formed changed its spin direction to the current one. This realization sheds light on the angular momentum evolution of the progenitor star, a process which is strongly affected by interaction with the binary companion. The slow proper motion of the system also implies that the system must have undergone a phase of mass transfer in which Star A shed a significant fraction of its mass onto B.

As we observe in the moon-earth system, tidal interactions in binary systems can lead to angular momentum exchange. The presence of viscosity is generally regarded as the condition for such transfer to happen. In this paper, we show how the orbital evolution can cause a persistent torque between the binary components, even for inviscid bodies. This preferentially occurs at the final stage of coalescence of compact binaries, when the orbit shrinks successively by gravitational waves and plunging on a timescale shorter than the viscous timescale. The total orbital energy transferred to the secondary by this torque is ~0.01 of its binding energy. We further show that this persistent torque induces a differentially rotating quadrupole perturbation. Specializing to the case of a secondary neutron star, we find that this non equilibrium state has an associated free energy of 10^47-10^48 erg, just prior to coalescence. This energy is likely stored in internal fluid motions, with a sizable amount of differential rotation. By tapping this free energy reservoir, a preexisting weak magnetic field could be amplified up to a strength of ~10^15 Gauss. Such a dynamically driven tidal torque can thus recycle an old neutron star into a highly magnetized neutron star, with possible observational consequences at merger.

Differential rotation induced by the r-mode instability can generate very strong toroidal fields in the core of accreting, millisecond spinning neutron stars. We introduce explicitly the magnetic damping term in the evolution equations of the r-modes and solve them numerically in the Newtonian limit, to follow the development and growth of the internal magnetic field. We show that the strength of the latter can reach large values, $B \sim 10^{14} $ G, in the core of the fastest accreting neutron stars. This is strong enough to induce a significant quadrupole moment of the neutron star mass distribution, corresponding to an ellipticity $|\epsilon_B}| \sim 10^{-8}$. If the symmetry axis of the induced magnetic field is not aligned with the spin axis, the neutron star radiates gravitational waves. We suggest that this mechanism may explain the upper limit of the spin frequencies observed in accreting neutron stars in Low Mass X-Ray Binaries. We discuss the relevance of our results for the search of gravitational waves.

The recent discovery of the "weak field, old magnetar", the soft gamma repeater SGR 0418+5729, whose dipole magnetic field is less than 7.5 \times 10^{12} G, has raised perplexing questions: How can the neutron star produce SGR-like bursts with such a low magnetic field? What powers the observed X-ray emission when neither the rotational energy nor the magnetic dipole energy are sufficient? These observations, that suggest either a much larger energy reservoir or a much younger true age (or both), have renewed the interest in the evolutionary sequence of magnetars. We examine, here, a phenomenological model for the magnetic field decay: B_dip} \propto (B_dip)^{1+a} and compare its predictions with the observed period, P,the period derivative, \dot{P}, and the X-ray luminosity, L_X, of magnetar candidates. We find a strong evidence for a dipole field decay on a timescale of \sim 10^3 yr for the strongest (\sim 10^{15} G) field objects, with a decay index within the range 1 \leq a < 2 and more likely within 1.5\lesssim a \lesssim 1.8. The decaying field implies a younger age than what is implied by the spinown age. Surprisingly, even with the younger age, the energy released in the dipole field decay is insufficient to power the X-ray emission, suggesting the existence of a stronger internal field, B_int. Examining several models for the internal magnetic field decay we find that it must have a very large (> 10^{16} G) initial value. Our findings suggest two clear distinct evolutionary tracks -- the SGR/AXP branch and the transient branch, with a possible third branch involving high-field radio pulsars that age into low luminosity X-ray dim isolated neutron stars.

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.

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.

We investigate a model for the shallow decay phases of Gamma-ray Burst (GRB) afterglows discovered by Swift/XRT in the first hours following a GRB event. In the context of the fireball scenario, we consider the possibility that long-lived energy injection from a millisecond spinning, ultramagnetic neutron star (magnetar) powers afterglow emission during this phase. We consider the energy evolution in a relativistic shock subject to both radiative losses and energy injection from a spinning down magnetar in spherical symmetry. We model the energy injection term through magnetic dipole losses and discuss an approximate treatment for the dynamical evolution of the blastwave. We obtain an analytic solution for the energy evolution in the shock and associated lightcurves. To fully illustrate the potential of our solution we calculate lightcurves for a few selected X-ray afterglows observed by Swift and fit them using our theoretical lightcurves. Our solution naturally describes in a single picture the properties of the shallow decay phase and the transition to the so-called normal decay phase. In particular, we obtain remarkably good fits to X-ray afterglows for plausible parameters of the magnetar. Even though approximate, our treatment provides a step forward with respect to previously adopted approximations and provides additional support to the idea that a millisecond spinning (1-3 ms), ultramagnetic (B$\sim 10^{14}-10^{15}$ G) neutron star loosing spin energy through magnetic dipole radiation can explain the luminosity, durations and shapes of X-ray GRB afterglows.

This contribution summarizes, as of early 2008, the observational and theoretical understanding of the physics, and emission properties of two Double Degenerate Binaries with the shortest orbital period known to date. In particular, the Unipolar Inductor Model and its coupling to GW emission have been introduced to explain a number of puzzling features that these two sources have in common and that are thought difficult to reconcile with models of mass transfer in such systems. Emphasis was put on the relevant new physical features that characterize the model. The role of spin-orbit coupling through the Lorentz torque and the role of GW emission in keeping the electric interaction active at all times has been thoroughly discussed in all their implications. Application of the model to both HM Cnc and V407 Vul accounts in a natural way for their main observational properties. Constraints on physical parameters are derived in order for the model to work, and can be verified by future observations.

A 321.5 s modulation was discovered in 1999 in the X-ray light curve of HM Cnc. In 2001 and 2002, optical photometric and spectroscopic observations revealed that HM Cnc is a very blue object with no intrinsic absorptions but broad (FWHM 1500 km s^-1) low equivalent width emission lines (EW 1-6A), which were first identified with the HeII Pickering series. The combination of X-ray and optical observations pictures HM Cnc as a double degenerate binary hosting two white dwarfs, and possibly being the shortest orbital period binary discovered so far. The present work is aimed at studying the orbital motion of the two components by following the variations of the shape, centroid and intensity of the emission lines through the orbit. In February 2007, we carried out the first phase resolved optical spectroscopic study with the VLT/FORS2 in the High Time Resolution (HIT) mode, yielding five phase bins in the 321 s modulation. Despite the low SNR, the data show that the intensity of the three most prominent emission lines, already detected in 2001, varies with the phase. These lines are detected at phases 0.2-0.6 where the optical emission peaks, and marginally detected or not detected at all elsewhere. Moreover, the FWHM of the emission lines in the phase resolved spectra is smaller, by almost a factor 2, than that in the the phase-averaged 2001 spectrum. Our results are consistent with both the pulsed optical component and emission lines originating in the same region which we identify with the irradiated surface of the secondary. Moreover, regardless of the origin of the un-pulsed optical continuum, we note that the EWs of the emission lines might be up to -15 / -25A, larger than thought before; these values are more similar to those detected in cataclysmic variables. All the findings further confirm that the 321s modulation observed in HM Cnc is the orbital period of the system, the shortest known to date.

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.

We present a state-of-the-art scenario for newly born magnetars as strong sources of Gravitational Waves (GWs)in the early days after formation. We address several aspects of the astrophysics of rapidly rotating, ultramagnetized neutron stars (NSs), including early cooling before transition to superfluidity, the effects of the magnetic field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, oblique rotator and the strength of the electromagnetic torque on the newly born NS. We show that our scenario is consistent with recent studies of SNRs surrounding AXPs and SGRs in the Galaxy that constrain the electromagnetic energy input from the central NS to be <= 10^51 erg. We further show that if this condition is met, then the GW signal from such sources is potentially detectable with the forthcoming generation of GW detectors up to Virgo cluster distances where an event rate <= 1/yr can be estimated. Finally, we point out that the decay of an internal magnetic field in the 10^16 G range couples strongly to the NS cooling at very early stages, thus significantly slowing down both processes: the field can remain this strong for at least 10^3 yrs, during which the core temperature stays higher than several times 10^8 K.

We report on the intense burst ``forest'' recorded on 2006 March 29 which lasted for ~30s. More than 40 bursts were detected both by BAT and by XRT, seven of which are rare intermediate flares (IFs): several times 10^{42} ergs were released. The BAT data were used to carry out time-resolved spectroscopy in the 14-100keV range down to 8ms timescales. This unique dataset allowed us to test the magnetar model predictions such as the magnetically trapped fireball and the twisted magnetosphere over an unprecedented range of fluxes and with large statistics (in terms of both photons and IFs). We confirmed that a two blackbody component fits adequately the time-resolved and integrated spectra of IFs. However, Comptonization models give comparable good reduced chi^2. Moreover, we found: i) a change of behavior, around ~10^{41} erg/s, above which the softer blackbody shows a sort of saturation while the harder one still grows to a few times 10^{41} erg/s; ii) a rather sharp correlation between temperature and radii of the blackbodies (R^2 prop kT^{-3}), which holds for the most luminous parts of the flares (approximately for L_{tot} > 10^{41} erg/s). Within the magnetar model, the majority of these findings can be accounted for in terms of thermalised emission from the E-mode and O-mode photospheres. Interestingly, the maximum observed luminosity coming from a region of ~15km matches the magnetic Eddington luminosity at the same radius, for a surface dipole field of ~8 x 10^{14} G (virtually equal to the one deduced from the spindown of SGR 1900+14).

Previous studies of the X-ray flux and spectral properties of 1RXS J170849-400910 showed hints of a possible correlation with the spin glitches that occurred in 1999 and 2001. However, due to the sparseness of spectral measurements and the paucity of detected glitches no firm conclusion could be drawn. We retrieved and analysed archival XTE pointings of 1RXS J170849-400910 covering the time interval between January 2003 and June 2006 and carried out a detailed timing analysis by means of phase fitting techniques. We detected two relatively large glitches Delta nu / nu of 1.2 and 2.1 10^-6 occurred in January and June 2005. Interestingly, the occurrence times of these glitches are in agreement with the predictions made in our previous studies. This finding strongly suggests a connection between the flux, spectral and timing properties of 1RXS J170849-400910.

In 2003 a previously unpulsed Einstein and ROSAT source cataloged as soft and dim (Lx of few 10^33 ergs) thermal emitting object, namely XTE J1810-197, was identified as the first unambiguous transient Anomalous X-ray Pulsar. Two years later this source was also found to be a bright highly polarized transient radio pulsar, a unique property among both AXPs and radio pulsars. In September 2006 the Swift Burst Alert Telescope (BAT) detected an intense burst from the candidate AXP CXOU J164710.2-455216, which entered in an outburst state reaching a peak emission of at least a factor of 300 higher than quiescence. Here, we briefly outline the recent results concerning the outburst phenomena observed in these two AXPs. In particular, XTE J1810-197 has probed to be a unique laboratory to monitor the timing and spectral properties of a cooling/fading AXP, while new important information have been inferred from X-ray and radio band simultaneous observations. CXOU J164710.2-455216 has been monitored in X-rays and radio bands since the very beginning of its outbursting state allowing us to cover the first phases of the outburst and to study the timing and spectral behavior during the first months.

On September 21, 2006, an intense (~10^39 erg/s) and short (20 ms) burst was detected by Swift BAT at a position consistent with that of the candidate Anomalous X-ray Pulsar, CXOU J164710.2-455216, discovered by Chandra in 2005. Swift follow-up observations began about 13 hours after the event and found the source at a 1-10keV flux level of about 4.5 x 10^-11 erg/s/cm^2, i.e. ~300 times brighter than measured 5 days earlier by XMM. We report the results obtained from Swift BAT observations of the burst and subsequent Swift XRT observations carried out during the first four months after the burst. These data are complemented with those from two XMM observations (carried out just before and after the BAT event) and four archival Chandra observations carried out between 2005 and 2007. We find a phase coherent solution for the source pulsations after the burst. The evolution of the pulse phase comprises an exponential component decaying with timescale of 1.4d which we interpret as the recovery stage following a large glitch (Delta nu / nu about 6 x 10^-5). We also detect a quadratic component corresponding to a spin-down rate of Pdot ~ 9 x 10^-13, implying a magnetic field strength of 10^14 Gauss. During the first Swift XRT observation taken 0.6 days after the burst, the spectrum showed a kT = 0.65keV blackbody (R_BB = 1.5km) plus a Gamma=2.3 power-law accounting for about 60% of the 1-10 keV observed flux. Analysis of Chandra archival data, taken during 2005 when the source was in quiescence, reveal that the modulation in quiescence is 100% pulsed at energies above 4 keV and consistent with the (unusually small-sized) blackbody component being occulted by the neutron star as it rotates.

Soft Gamma Repeaters and the Anomalous X-ray Pulsars are believed to contain slowly spinning "magnetars". The enormous energy liberated in the 2004 Dece 27 giant flare from SGR 1806-20, together with the likely recurrence time of such events, points to an internal magnetic field strength ~ 10^{16} G. Such strong fields are expected to be generated by a coherent alpha-Omega dynamo in the early seconds after the Neutron Star formation, if its spin period is of a few milliseconds at most. A substantial deformation of the NS is caused by such fields and a newborn millisecond-spinning magnetar would thus radiate for a few days a strong gravitational wave signal. Such a signal may be detected with Advanced LIGO-class detectors up to the distance of the Virgo cluster, where ~ 1 magnetar per year are expected to form. Recent X-ray observations reveal that SNRs around magnetar candidates do not show evidence for a larger energy content than standard SNRs (Vink & Kuiper 2006). This is at variance with what would be expected if the spin energy of the young, millisecond NS were radiated away as electromagnetic radiation andd/or relativistic particle winds and, thus, transferred quickly to the expanding gas shell. We show here that these recent findings can be reconciled with the idea of magnetars being formed with fast spins, if most of their initial spin energy is radiated thorugh GWs. In particular, we find that this occurs for essentially the same parameter range that would make such objects detectable by Advanced LIGO-class detectors up to the Virgo Cluster.

We further discuss the Unipolar Inductor Model (UIM) coupled to GW emission (Dall'Osso et al. 2005) and compare it to observed properties of the two candidate ultrashort period binaries RX J0806+15 and RX J1914+24 . We consider the measured orbital periods, period derivatives and inferred X-ray luminosities of these two sources and find constraints on system parameters in order for the model to account for them. We find that these properties point to the two sources being in different regimes of the UIM, with the requirement of low magnetic moment primaries ($\sim 10^{30} G cm^3$) for both. Given this weak magnetization, RX J0806+15 has a sufficiently low luminosity that it can be interpreted as having a primary spin almost synchronous to and just slightly slower than the orbital motion. Its measured orbital spin-up is only slightly affected by spin-orbit coupling and is mostly due to GW emission. RX J1914+24, on the other hand, is too bright in X-rays and has too slow an orbital spin-up for the same regime to apply. We suggest that this binary system may be emitting GWs at a significantly higher rate than implied by its measured $\dot{\omega}_o = 6 \times 10^{-17} rad s^{-2}$. The latter is explained, in this framework, by the primary spin being slightly faster than the orbital motion ($\alpha <or= 1.1$). In this case, the associated spin-orbit coupling transfers to the orbit a significant amount of angular momentum, thus partially balancing that lost to GW emission. All expectations can be tested in the near future to confirm the viability of the model.

There is growing evidence that two classes of high-energy sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars contain slowly spinning ``magnetars'', i.e. neutron stars whose emission is powered by the release of energy from their extremely strong magnetic fields (>10^15 G. We show here that the enormous energy liberated in the 2004 December 27 giant flare from SGR1806-20 (~5 10^46 erg), together with the likely recurrence time of such events, requires an internal field strength of > 10^16 G. Toroidal magnetic fields of this strength are within an order of magnitude of the maximum fields that can be generated in the core of differentially-rotating neutron stars immediately after their formation, if their initial spin period is of a few milliseconds. A substantial deformation of the neutron star is induced by these magnetic fields and, provided the deformation axis is offset from the spin axis, a newborn fast-spinning magnetar would radiate for a few weeks a strong gravitational wave signal the frequency of which (0.5-2 kHz range) decreases in time. The signal from a newborn magnetar with internal field > 10^16.5 G could be detected with Advanced LIGO-class detectors up to the distance of the Virgo cluster (characteristic amplitude h_c about 10^-21). Magnetars are expected to form in Virgo at a rate approx. 1/yr. If a fraction of these have sufficiently high internal magnetic field, then newborn magnetars constitute a promising new class of gravitational wave emitters.

We consider the Unipolar Inductor Model (Goldreich & Lynden-Bell 1969) applied to Double Degenerate Binaries (DDBs) with ultrashort periods (Wu et al. 2002). In this model a magnetized primary white dwarf has a slight asynchronism between its spin and orbital motion, so that the (non-magnetic) secondary experiences a motional electric field when moving through the primary field lines. This induces a current flow between the two stars and provides an electric spin-orbit coupling mechanism for the primary. We study the combined effect of Gravitational Wave emission and electric spin-orbit coupling on the evolution of the primary degree of asynchronism and the associated rate of electric current dissipation in such systems, assuming that the primary's spin is not affected by any other mechanisms. In particular, we show that in ultrashort period binaries the emission of GW pumps energy in the electric circuit as to keep it steadily active. This happens despite the fact that spin-orbit coupling can rapidly synchronize the primary, because GW represent a slow desynchronizing mechanism steadily substracting orbital angular momentum to the system. A slightly asynchronous steady-state is thus achieved, determined by the balance between these two competing effects. This can be shown to correspond to a condition where the total available electric energy is conserved, because of GW emission, while dissipation, synchronization and orbital shrinking continue.

We have discovered rapid Quasi Periodic Oscillations (QPOs) in RXTE/PCA measurements of the pulsating tail of the 27th December 2004 giant flare of SGR 1806-20. QPOs at about 92.5Hz are detected in a 50s interval starting 170s after the onset of the giant flare. These QPOs appear to be associated with increased emission by a relatively hard unpulsed component and are seen only over phases of the 7.56s spin period pulsations away from the main peak. QPOs at about 18 and 30Hz are also detected, 200-300s after the onset of the giant flare. This is the first time that QPOs are unambiguously detected in the flux of a Soft Gamma-ray Repeater, or any other magnetar candidate. We interpret the highest QPOs in terms of the coupling of toroidal seismic modes with Alfven waves propagating along magnetospheric field lines. The lowest frequency QPO might instead provide indirect evidence on the strength of the internal magnetic field of the neutron star.

RX J0806.3+1527 is thought to be a 321s orbital period (the shortest known) double white dwarf binary system. According to the double degenerate binary (DDB) scenario this source is expected to be one of the strongest gravitational wave (GW) emitter candidates. In the last years RX J0806.3+1527 has been studied in great details, through multiwavelength observational campaigns and from the point of view of data analysis result interpretations. We present here the timing results obtained thanks to a 3.5-year long optical monitoring campaign carried out by the Very Large Telescope (VLT) and the Telescopio Nazionale Galileo (TNG) which allowed us to detect and study the orbital period derivative (spin-up at a rate of about 10^-3 s/yr) of the 321s modulation, to detect the linear polarisation (at a level of about 2%), and to study the broad band energy spectrum. The VLT/TNG observational strategy we used allowed us to rely upon a P-Pdot coherent solution which we finally extended backward to the 1994 ROSAT observation of RX J0806.3+1527.

We report on a timing analysis of archival observations of the Anomalous X-ray Pulsar 1RXSJ170849.0-400910 made with the RXTE Proportional Counter Array. We detect a new large glitch (Delta nu / nu ~ 3 x 10^-6) which occurred between 2001 March 27 and 2001 May 6, with an associated large increase in the spin-down rate (Delta nudot/nudot ~ 0.3). The short time (1.5 yrs) elapsed from the previously detected glitch and the large amplitude of the new spin-up place this source among the most frequent glitchers, with large average glitch amplitudes, similar to those of the Vela pulsar. The source shows different recoveries after the glitches: in the first one it is well described by a long term linear trend similar to those seen in Vela-like glitches; in the second case the recovery is considerably faster and is better described by an exponential plus a fractional change in the long-term spin-down rate of the order of 1%. No recovery of the latter is detected but additional observations are necessary to confirm this result. Observed glitch properties are compared to those of radio pulsar glitches; current models are discussed in light of our results. It appears that glitches may represent yet another peculiarity of AXPs. Starquake-based models appear to be prefered on qualitative grounds. Alternative models can be applied to individual glitches but fail in explaining both. Thus the two events may as well arise from two different mechanisms.