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

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

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

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

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

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

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

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

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

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

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

We propose a simple dynamical model that may account for the observed spectral and temporal properties of GRS 1915+105 and XTE J1550-5634. The model is based on the assumption that a fraction of the radiation emitted by a hot spot lying on the accreting disk is dynamically Comptonized by the relativistic jet that typically accompanies the microquasar phenomenon. We show that scattering by the jet produces a detectable modulation of the observed flux. In particular, we found that the phase lag between hard and soft photons depends on the radial position of the hot spot and, if the angle between the jet and the line of sight is sufficiently large, the lags of the fundamental and its harmonics may be either positive or negative.

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

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

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

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

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

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

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