A variable absorption feature in the X-ray spectrum of a magnetar
Soft--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[mereghetti08, rea11]. This behaviour can be understood in the magnetar model[thompson95, thompson96, tlk02], 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[kouveliotou98, thompson01, vietri07] 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[rea13], whereas a strong field has been proposed to reside in the stellar interior[rea10, turolla11] and in multipole components on the surface[guver11]. 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 210 gauss to more than 10 gauss.
Istituto Universitario di Studi Superiori, piazza della Vittoria 15, I-27100 Pavia, Italy.
Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, INAF, via Bassini 15, I-20133 Milano, Italy.
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, via Bassi 6, I-27100 Pavia, Italy.
Dipartimento di Fisica e Astronomia, Università di Padova, via Marzolo 8, I-35131 Padova, Italy.
Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK.
AIM CEA/Irfu/Service d’Astrophysique, Orme des Merisiers, F-91191 Gif-sur-Yvette, France.
Osservatorio Astronomico di Roma, INAF, via Frascati 33, I-00040 Monteporzio Catone, Italy.
Institut de Ciències de l’Espai (IEEC–CSIC), Campus UAB, Torre C5, 2a planta, E-08193 Barcelona, Spain.
On 2009 June 5 two short bursts of hard X-rays, detected by Fermi and other satellites, revealed the previously unknown source SGR 0418+5729[vanderhorst10]. Subsequent observations with the Rossi X-ray Timing Explorer (RXTE), Swift, Chandra and X-ray Multi-mirror Mission (XMM) Newton satellites found the new SGR to be an X-ray pulsar with a period of 9.1 s and a luminosity of erg s (in the 0.5–10 keV band and for a distance of 2 kpc)[vanderhorst10, esposito10]. During the three years after the onset of the outburst, the spectrum softened and the luminosity declined by three orders of magnitude, but remained still too high to be powered by rotational energy[esposito10, rea10, rea13]. The measured spin-down rate of s s translates (under the assumption of rotating magnetic dipole in vacuo) into a magnetic field G at the magnetic equator[rea13], a value well in the range of normal radio pulsars. However, the presence of high-order multipolar field components of G close to the surface has been invoked to interpret the spectrum of the source in the framework of atmosphere models[guver11]. In any case, a strong crustal magnetic field ( G) seems to be required to explain the overall properties of SGR 0418+5729 within the magnetar model[turolla11, rea13].
Hints of the presence of an absorption feature at 2 keV in the spectrum of SGR 0418+5729 were found in the phase-resolved analysis of data (with relatively low-count statistics) from the Swift X-ray Telescope (XRT) taken during 2009 July 12–16[esposito10]. Thanks to the large collecting area and good spectral resolution of the European Photon Imaging Camera (EPIC), we were able to perform a more detailed investigation using data collected by XMM-Newton during a 67-ks long observation performed on 2009 August 12, when the source flux was still high ( erg cm s in the 2–10 keV band).
To examine the spectral variations as a function of the star’s rotational phase without making assumptions about the X-ray spectral energy distribution of SGR 0418+5729, we produced a phase–energy image by binning the EPIC source counts into energy and rotational phase channels and then normalising to the phase-averaged energy spectrum and pulse profile. The normalised phase–energy image (Fig.1) shows a prominent V-shaped feature in the phase interval 0.1–0.3. This is produced by a lack of counts in a narrow energy range with respect to nearby energy channels, that is, an absorption feature at a phase-dependent energy. The regular shape of the feature in the phase–energy plane as well as its presence in the three independent EPIC detectors (see Supplementary Fig. 5) exclude the possibility that it results from statistical fluctuations in the number of counts or from an instrumental effect. Another absorption feature is visible at low energies at phase 0.5–0.6.
We extracted from the EPIC data the phase-averaged spectrum of SGR 0418+5729, as well as the spectra from 50 phase intervals of width 0.02 rotational cycles, as described in the Supplementary Information. The phase-averaged spectrum can be adequately fit by either a two-blackbody model ( for 196 degrees of freedom, d.f.) or a blackbody plus power-law model ( for 196 d.f.), corrected for interstellar absorption (see refs 11 and 12 for other models that can fit the spectrum).
The 15 spectra extracted from the phase intervals 0.1–0.3 and 0.5–0.6, unlike those of the remaining phases, cannot be fitted by a renormalisation of the phase-averaged best-fit model, which gives in most cases null hypothesis probabilities in the range – (see Supplementary Fig. 4). They are instead well fitted (null hypothesis probability 0.03) by the addition of a narrow absorption line component, which can be equally well modelled with a Gaussian profile or a cyclotron absorption line model[makishima90] (the improvement obtained by adding a cyclotron component in the phase intervals 0.1–0.3 and 0.5–0.6 can be seen in Supplementary Fig. 4). The best-fit line parameters as a function of phase are shown in Fig. 2 and an example of phase-resolved spectrum is displayed in Fig. 3.
We searched for the phase-dependent absorption feature in all the available X-ray observations of SGR 0418+5729 and found that it was present in the phase interval 0–0.3, and up to higher energies than in XMM-Newton, in RXTE data taken during the first two months of the outburst (see Supplementary Fig. 6).
Absorption features have been observed in the X-ray spectra of various classes of neutron stars[truemper78, heindl04, turolla09, haberl06, vankerkwijk07, bignami03, gotthelf13, kdm12] and interpreted as being due to either cyclotron absorption (by electrons or protons) or bound–bound atomic transitions. However, variations in the line energy as a function of the rotational phase as large as in SGR 0418+5729 (by a factor in one-tenth of a cycle) have not been seen in any source.
In a neutron star atmosphere, different atomic transitions might be responsible for a phase-variable absorption feature if temperature, elemental abundance or magnetic field vary strongly on the surface. The line energies observed in SGR 0418+5729 (1–5 keV) rule out transitions in magnetised H and He, which occur below 1 keV[potekhin98, medin08]. On the other hand, the absorption spectra of heavier elements are much more complex (see, for example, ref. 26 for C, O and Ne) and some lines could occur at high-enough energies. However, to explain the phase resolved spectra of SGR 0418+5729, the physical conditions of a heavy-element atmosphere are forced to vary in such a way that a single transition should dominate the opacity at each of the phases where the absorption line is detected.
A more straightforward explanation for the line variability can instead be given if the feature is due to cyclotron resonant scattering. The cyclotron energy (in keV) for a particle of charge and mass in magnetic field (in gauss) is given by
where (which is 0.8 at the star surface for typical neutron star mass and radius and km, respectively) accounts for the gravitational redshift at distance from the neutron star centre, is the mass of the electron, and is the velocity of light. In this case, the phase variability of the feature energy would simply be due to the different fields experienced by the charged particles interacting with the photons directed towards us as the neutron star rotates.
If the absorbers/scatterers are electrons hovering near the star surface, the expected line energy is 70 keV for the dipole field at the equator of SGR 0418+5729 ( G); this line energy is more than 10 times higher than that observed. A possible way to explain this large discrepancy might be to assume that the electrons producing the line are located higher up in the magnetosphere in a dipolar geometry, where the magnetic field is smaller ( to have ). Moreover, such an electron population should also be nearly mono-energetic, or subrelativistic, in order to prevent Compton scattering from washing out the feature, which would require a mechanism to maintain slowly moving electrons confined in a small volume high in the magnetosphere.
If the particles responsible for the cyclotron scattering are protons, the energy range of the SGR 0418+5729 spectral feature requires a magnetic field G (it would be even larger for heavier ions). In the framework of the magnetar model, the unprecedented phase-variability of the line energy can be explained by the complex topology of the magnetar magnetospheres, in which global and/or localised twists play an important part[tlk02]. This is particularly true for SGR 0418+5729, which has a weak dipolar component, as testified by the small spin-down value, whereas a much stronger internal magnetic field has been advocated to explain its X-ray luminosity and burst activity[rea10, turolla11]. Furthermore, the presence of small-scale, strong, multipolar components of the surface field has been inferred by fitting its phase-averaged X-ray spectrum with models of magnetized neutron star atmospheres[guver11].
In this context, the observed line variability might be due to the presence of either strong magnetic field gradients along the surface or vertical structures (with a spatially dependent field) emerging from the surface. To work out how the dynamic magnetosphere of a magnetar should look, an analogy with the solar corona in the proximity of sunspots has been proposed (see, for example, ref. 27). In particular, localised, baryon-rich, magnetic structures, in the form of rising flux tubes, or ‘prominences’, produced by magnetic reconnection or the emergence of the internal field near a crustal fault, have been proposed to explain some of the observed properties of the giant flare emitted in 2004 by SGR 1806–20[gelfand05, masada10]. If a similar scenario, albeit on a reduced scale, occurred during the outburst of SGR 0418+5729, a spectral feature might arise as thermal photons from the hot spot (a small hot region on the neutron star surface, responsible for most of the X-ray emission, which could be itself related to the prominence) cross the plasma threading the magnetic loop. A proton density cm is needed to produce a resonant scattering depth of order unity[tlk02]. Protons, being heavy, do not rise much above the surface and move subrelativistically[tlk02], so resonant scattering in the prominence is likely to produce a narrow feature instead of an extended tail. As the star rotates, photons emitted in different directions pass through portions of the prominence with different magnetic field, density and size, giving rise to the observed variations of the line centroid and width. A simple quantitative model based on this picture is presented in Supplementary Information. Results, obtained with a geometry consistent with the constraints derived from the X-ray pulsed fraction of SGR 0418+5729, are in good agreement with the observed variations of the feature with phase (Fig. 1).
We thank G. Goggi and C. Paizis for useful discussion. This research is based on data and software provided by the the ESA XMM-Newton Science Archive (XSA) and the NASA/GSFC High Energy Astrophysics Science Archive Research Center (HEASARC). The Authors acknowledge partial funding from INAF through a PRIN 2010 grant and ASI through contract I/032/10/0.
The authors declare that they have no competing financial interests. Correspondence and requests for materials should be addressed to A.T. (email: firstname.lastname@example.org).
1 XMM-Newton data analysis
The data have been processed with version 12 of the
Scientific Analysis Software (SAS) and we used the
most recent (2012 October) calibration files available for the
EPIC instrument. EPIC consists of two MOS and one
pn CCD cameras sensitive to photons with energy
between 0.2 and 10 keV. During the 2009 observation (see also
Refs. 32, 33), the two MOS and the pn cameras
were set in Small Window mode (time resolution of 0.3 s and 5.7
ms, respectively); all detectors were operated with the thin
optical blocking filter. Periods in which the particle background
was unusually high because of soft proton flares were excluded
using an intensity filter. This reduced the net exposure time to
30.7 ks, 50.0 ks and 50.7 ks for pn, MOS1 and MOS2, respectively.
Again according to standard procedure, photon event grades higher than 12
for the MOS cameras and 4 for the pn were filtered out. Photon
arrival times were converted to the Solar System barycentre
reference frame, by using the coordinates ,
(J2000) and the spin phases were computed with the timing
parameters of Ref. 32 (, Epoch 54993 MJD,
valid over the range MJD 54993–55463). The phase–energy images
extracted for each of the three EPIC cameras are shown in Fig. 5.
The V-shaped feature indicating the presence of a phase-dependent
absorption line is present in the three independent datasets.
We also analysed in the same way the three XMM-Newton observations
of SGR 0418+5729 performed from September 2010 to August 2012 (the
observation settings and the source flux are reported in
Ref. 35). However, due to the lower flux, they were
not sensitive enough to detect the absorption feature, even if it
were still present.
For the source spectra we used the pn counts extracted from a circular region with radius of 35; the background spectra were extracted from source-free regions on the same chip as the target. The ancillary response files and the spectral redistribution matrices were generated with the SAS tasks arfgen and rmfgen, respectively. For the spectral analysis, performed in the 0.3–10 keV energy range, we used the XSPEC fitting package version 12.4. The abundances adopted were those of Ref. 37 and for the photoelectric absorption we used the cross-sections from Ref. 38.
2 RXTE and Swift data analysis
The RXTE and Swift observations used in this research have been
already presented in Ref. 39, to which we refer
for more information. All data were reprocessed and analysed with
version 6.12 of the HEAsoft package and the CALDB
calibration data base available in 2013 February. Apart from this,
the Swift data were reduced exactly as described in
For each RXTE/PCA data-set, we ran the xenon2fits script to
combine the GoodXenon files into science event tables with 256
energy bins. We used only data from the Proportional Counter Unit
PCU-2, since it is the best-calibrated unit of the PCA instrument,
and selected photons from the layer 1 of the detector. We
corrected in the event tables the photon arrival times to the
Solar System barycentre using the script fxbary.
Fig. 6 shows the phase–energy images (rescaled to the
phase-averaged spectrum) obtained from PCA observations of SGR 0418+5729 in three different epochs.
phase-dependent absorption feature is apparent in the phase
interval 0–0.3, and up to higher energies (10 keV) than in
XMM-Newton, in the data taken in the first two months after the
onset of the outburst, and marginally visible in August 2009 due
to the low signal-to-noise ratio. The possible presence of the
line at even higher energies cannot be tested due to the very few
source photons detected at these energies.
The phase-resolved spectra, with 50 phase bins, were extracted
with fasebin (which also barycentres the data) and combined
with fbadd. Consistent phase-averaged spectra were used by
the script pcarsp to make the 256-bin response matrices,
which were combined with addrmf. Using the 2009
June–August spectra, we performed a similar analysis as described
for the XMM-Newton data. While the results were consistent, owing
to the limited spectral capabilities of the PCA instrument the
analysis added no new information on the
characteristics of the absorption line.
The phase–energy image obtained from the Swift/XRT data taken on 2009 July 12–16 is consistent with the EPIC one. Although the limited spectral capabilities and high instrumental background of the PCA instrument and the poor counting statistics of the XRT data are not adequate for a detailed characterisation of the feature properties, these data indicate that the absorption feature observed with EPIC had been present in the spectrum of SGR 0418+5729 for at least two months.
3 A simple proton cyclotron resonance model
We present here a simple model to illustrate how resonant proton
cyclotron scattering can produce a phase-variable absorption
feature in the X-ray spectrum of SGR 0418+5729. Let us assume that thermal
emission from the star surface comes from a
small hot spot
SGR 0418+5729 the actual angular size is assuming a
distance of 2 kpc and a star radius km; see e.g. Ref. 33) and that an
ultra-strong, small-scale magnetic field is present above the spot
(the large scale field is a dipole with surface strength
G). The field is taken to thread a
plasma-loaded magnetic loop of radius ,
containing (non-relativistic) protons of density . For the sake of simplicity,
we neglect the loop extent in the radial direction, while
retaining a finite transverse width, and assume that the B-field
lines are along the loop. Although many different geometries may
be envisaged, in the following the loop is assumed to have the
shape of a spherical lune (the round surface of a spherical wedge)
with the diameter on the star surface and the spot at its centre,
dihedral angle and inclination with respect to
the surface normal (see Fig. 7).
Photons of frequency emitted by the spot may undergo resonant cyclotron scattering as they traverse the loop and the optical depth is
where the expression for the resonant cross section (e.g. Ref. 40) was used, is the angle between the photon direction and , is the radial coordinate counted from the spot and is the particle cyclotron frequency. If is above G, photons in the –10 keV range can only resonantly scatter on protons. Since we are assuming that matter is confined in a very thin layer at and , equation (1) reduces to
so a proton density is required
to make the loop thick to resonant scattering. If this occurs, the
flux emitted by the spot will be reduced by a factor
in traversing the baryon-loaded loop. Because
only photons with do scatter, this will
produce a monochromatic absorption line at . This is
clearly an effect of our approximations. The finite radial extent
of both the loop and the emitting spot will result in a broadening
of the absorption feature.
Thermal photons are emitted from the spot with a given angular pattern, but only those propagating along the unit vector reach the observer; is characterized by the angle it makes with the surface normal and by the associated azimuth . Because of general relativistic ray bending, differs from , the angle between the line-of-sight (LOS) and the spot normal. We relate them using Beloborodov’s approximation (see also Ref. 42)
where is the star Schwarzschild radius.
When the star rotates, both and change with the rotational phase , where is the star spin period. By introducing the angles and that the LOS and the spot normal form with the rotation axis, respectively, it is
the latter follows from geometrical considerations (see again
Fig. 7; note that the azimuths associated to and
coincide) and the last equality holds for
oriented along the meridian which passes through the spot
and pointing north. The angle is counted
from and changes by in a cycle if it is
while it oscillates between a minimum and a maximum
value (which depend on and ) in the opposite case.
Because and change during a cycle, the ray
direction may or may not intersect the loop, so an
absorption feature may be present at certain phases only.
What actually happens depends on the loop geometry (here
its transverse angular width and inclination
) and on the angles and .
and are not known for SGR 0418+5729, the only
constraint coming from the measured pulsed fraction (),
which implies that their sum is about .
When the ray directed towards the observer crosses the loop, the intersection will occur at different positions according to the phase, giving rise to features at different energies, proportional to the local magnetic field intensity. To obtain a phase–line energy relation to compare with observations, we need to introduce how the magnetic field depends on position in the loop. For simplicity, we assume that the field varies only along the loop, linearly decreasing its intensity as its transverse width increases:
where is the magnetic field intensity at the base of the loop (where the transverse width vanishes and the field is strongest), is a multiplicative factor measuring how fast the field
varies along the loop and is the angle between the lune
diameter and the meridian passing through the spot.
Selecting convenient combinations of the parameters, our simple analytical model can well reproduce the variation of the line energy with phase observed in SGR 0418+5729, including its relative phase with respect to the maximum of the pulse profile, which occurs at (see Fig. 2) and in our model corresponds to the phase where is maximum. The red line in Fig. 1 shows the result obtained for a loop with (i.e., parallel to the equator), , (consistent with the 50% pulsed fraction), and (see Table 1 for a list of all the model parameters). The magnetic field in the portion of the loop swept by the line of sight varies from G to G. Although this simple model cannot explain with a single symmetric loop also the low energy line observed at phase 0.55, the latter might originate from a second loop with a smaller angular size located on the other side of the hot spot.
- 2. Turner, M. J. L. et al. The European Photon Imaging Camera on XMM-Newton: The MOS cameras. A&A 365, L27–L35 (2001).
- 3. Strüder, L. et al. The European Photon Imaging Camera on XMM-Newton: The pn-CCD camera. A&A 365, L18–L26 (2001).
- 4. Rea, N. et al. A Low-Magnetic-Field Soft Gamma Repeater. Science 330, 944 (2010).
- 5. Turolla, R., Zane, S., Pons, J. A., Esposito, P. & Rea, N. Is SGR 0418+5729 Indeed a Waning Magnetar? ApJ 740, 105 (2011).
- 6. van der Horst, A. J. et al. Discovery of a New Soft Gamma Repeater: SGR J0418+5729. ApJ 711, L1–L6 (2010).
- 7. Rea, N. et al. The Outburst Decay of the Low Magnetic Field Magnetar SGR 0418+5729. ApJ 770, 65 (2013).
- 8. Arnaud, K. A. XSPEC: The First Ten Years. In Jacoby, G. H. & Barnes, J. (eds.) Astronomical Data Analysis Software and Systems V, vol. 101 of Astronomical Society of the Pacific Conference Series, 17–20 (ASP, San Francisco, 1996).
- 9. Wilms, J., Allen, A. & McCray, R. On the Absorption of X-Rays in the Interstellar Medium. ApJ 542, 914–924 (2000).
- 10. Balucinska-Church, M. & McCammon, D. Photoelectric absorption cross sections with variable abundances. ApJ 400, 699 (1992).
- 11. Esposito, P. et al. Early X-ray and optical observations of the soft gamma-ray repeater SGR0418+5729. MNRAS 405, 1787–1795 (2010).
- 12. Thompson, C., Lyutikov, M. & Kulkarni, S. R. Electrodynamics of Magnetars: Implications for the Persistent X-Ray Emission and Spin-down of the Soft Gamma Repeaters and Anomalous X-Ray Pulsars. ApJ 574, 332–355 (2002).
- 13. Beloborodov, A. M. Gravitational Bending of Light Near Compact Objects. ApJ 566, L85–L88 (2002).
- 14. Turolla, R. & Nobili, L. Pulse Profiles from Thermally Emitting Neutron Stars. ApJ 768, 147 (2013).
- 15. Makishima, K. et al. Observations of the peculiar hard X-ray transient X 0331+53 (V 0332+53). PASJ 42, 295–315 (1990).
|Angle between the spot normal and the rotation axis|
|Angle between the LOS and the rotation axis|
|Angle between the lune diameter and the spot meridian|
|Semi-amplitude of the lune transverse angle|
|Lune inclination angle|
|Gravitational redshift at the neutron star surface|
|Magnetic field intensity at the base of the loop|
|Linear term in Eq. 6|