No Arabic abstract
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.
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.
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.
Magnetars are believed to host the strongest magnetic fields in the present universe ($Bgtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.
We present time-domain Monte Carlo simulations of radio emission from cosmic ray air showers in the scheme of coherent geosynchrotron radiation. Our model takes into account the important air shower characteristics such as the lateral and longitudinal particle distributions, the particle track length and energy distributions, a realistic magnetic field geometry and the shower evolution as a whole. The Monte Carlo approach allows us to retain the full polarisation information and to carry out the calculations without the need for any far-field approximations. We demonstrate the strategies developed to tackle the computational effort associated with the simulation of a huge number of particles for a great number of observer bins and illustrate the robustness and accuracy of these techniques. We predict the emission pattern, the radial and the spectral dependence of the radiation from a prototypical 10^17 eV vertical air shower and find good agreement with our analytical results (Huege & Falcke 2003) and the available historical data. Track-length effects in combination with magnetic field effects surprisingly wash out any significant asymmetry in the total field strength emission pattern in spite of the magnetic field geometry. While statistics of total field strengths alone can therefore not prove the geomagnetic origin, the predicted high degree of polarisation in the direction perpendicular to the shower and magnetic field axes allows a direct test of the geomagnetic emission mechanism with polarisation-sensitive experiments such as LOPES. Our code provides a robust, yet flexible basis for detailed studies of the dependence of the radio emission on specific shower parameters and for the inclusion of additional radiation mechanism in the future.
Blue shifted absorption lines from highly ionised iron are seen in some high inclination X-ray binary systems, indicating the presence of an equatorial disc wind. This launch mechanism is under debate, but thermal driving should be ubiquitous. X-ray irradiation from the central source heats disc surface, forming a wind from the outer disc where the local escape velocity is lower than the sound speed. The mass loss rate from each part of the disc is determined by the luminosity and spectral shape of the central source. We use these together with an assumed density and velocity structure of the wind to predict the column density and ionisation state, then combine this with a Monte Carlo radiation transfer to predict the detailed shape of the absorption (and emission) line profiles. We test this on the persistent wind seen in the bright neutron star binary GX 13+1, with luminosity L/LEdd ~ 0.5. We approximately include the effect of radiation pressure because of high luminosity, and compute line features. We compare these to the highest resolution data, the Chandra third order grating spectra, which we show here for the first time. This is the first physical model for the wind in this system, and it succeeds in reproducing many of the features seen in the data, showing that the wind in GX13+1 is most likely a thermal-radiation driven wind. This approach, combined with better streamline structures derived from full radiation hydrodynamic simulations, will allow future calorimeter data to explore the detail wind structure.