EDP Sciences
Free Access

This article has an erratum: [erratum]

Issue
A&A
Volume 547, November 2012
Article Number A27
Number of page(s) 5
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201219978
Published online 23 October 2012

© ESO, 2012

1. Introduction

The supergiant fast X-ray transient (SFXT) is a new class of massive X-ray binaries discovered by INTEGRAL observations (Sguera et al. 2005; Negueruela et al. 2006), which contains a magnetized neutron star that accretes material from an OB supergiant companion star. They show short and sporadic outbursts (Romano et al. 2007), span 2–5 orders of magnitude in dynamic range (e.g. in’t Zand 2005; Leyder et al. 2007), and represent the extreme case of X-ray variability in high mass X-ray binaries (HMXBs) with early-type supergiant companions (Masetti et al. 2006; Pellizza et al. 2006). The quiescent state is characterized by a low luminosity in the range of  ~1031−1033 ergs/s (e.g. in’t Zand 2005; Leyder et al. 2007).

Some SFXTs are believed to be X-ray pulsars with long spin periods, and their measured pulse periods are located in a wide range, from 4.7 s (AX J1841.0-0536/IGR J18410-0535, Bamba et al. 2001) to 1240 s (IGR J16418-4532, Walter et al. 2006; Sidoli 2011a). Generally speaking, a newly born neutron star is presumed to be a fast rotator with a spin period of a fraction of one second. If a neutron star is in a binary system, it may experience an ejector state (the conventional spin-powered energy-loss mechanism) and propeller state (the interaction between magnetosphere of the neutron star and the stellar wind of the companion) before accreting material from its optical companion (for a review see Bhattacharya & van den Heuvel 1991, and references therein). During the ejector and propeller regimes, the neutron star spins down and forms a slow rotator (van den Heuvel 1977; Davies et al. 1979; Davies & Pringle 1981) at the end of these two phases. Under certain assumptions (see Urpin et al. 1998a; Li & van den Heuvel 1999), the spin period can be limited to a few hundred seconds (~500 s). If the system undergoes a spherical accretion and three epochs, i.e. the ejector, supersonic propeller, and subsonic propeller epoch, the conditions under which a longer period can be achieved may be relaxed (Ikhsanov 2001, 2007).

Long-period pulsars are often identified as accretors and their periods gradually decrease during the accretion in the prograde scenario. However, observations indicate that dynamic torque reversals are present in some wind-fed persistent HMXBs (Bildsten et al. 1997; Nelson et al. 1997, and references therein). The spin-down torque due to the accretion reversals has the same order of magnitude as the spin-up torque, but is opposite in sign (Nelson et al. 1997). Accretion torque reversals can occur in an inhomogeneous wind (Ruffer 1997, 1999).

In this paper, we suggest that the long spin periods are associated with a retrograde accretion scenario. In Sect. 2, we review the theory of wind accretion in SGXBs/HMXBs and derive the spin evolution and the spin-down timescale in different accretion regimes. We apply our model to SFXTs in Sect. 3. Section 4 contains a short summary.

2. Retrograde wind accretion

OB supergiants (Reig 2011) can produce a powerful stellar wind that is driven by radiation, which is approximately isotropic, that is, the wind material leaves the surface uniformly in all directions, with a wind velocity of vw ~ 1000−2000 km s-1. The OB companions in SFXTs typically have a mass of Mc ~ 30   M, mass loss rate of w ~ 10-6   M   yr-1, and luminosity of log (Lc/L) ~ 5−6. Here, M is the solar mass, and L is the solar luminosity.

The magnetized neutron star in SGXBs accretes matter from an OB companion via the intense stellar wind (Verbunt & van den Heuvel 1995). We consider the case of a spherically symmetric stellar wind radiated from the OB companions that is disordered and inhomogeneously and is captured by relativistic neutron stars. The accreting neutron star can be either a prograde accretor or a retrograde accretor (Bildsten et al. 1997; Shakura et al. 2012). If the neutron star behaves like a retrograde accretor, it may spin down and rotate with a very long period because of the inverted accretion torque.

2.1. Characteristic radii

In the theory of wind accretion by a magnetized rotating neutron star, three characteristic radii (Davidson & Ostriker 1973; Lamb et al. 1973; Stella et al. 1986) can be defined as follows. Only part of the wind material from the OB companion will be captured by the gravitational field of the neutron star. The radiated wind flowing toward the neutron star within a certain radius will be accreted, whereas the flowing material outside that radius will escape. This radius, called the accretion radius racc (Bondi & Hoyle 1944), can be defined by noting that material will only be accreted if it has a kinetic energy lower than the potential energy in the gravitational field of the neutron star, that is, (1)where M is the neutron star mass, and vw,8 is a velocity in units of 1000 km s-1 (vw,8 = vw/(108 cm s-1)).

When the magnetic pressure of the neutron star can balance the ram pressure of the accreting material, a magnetosphere forms and the inflowing matter may stop at this preferred radius – the magnetospheric radius rmag (Elsner & Lamb 1977; Davies & Pringle 1981). It can be defined by (2)where is the magnetic field strength at rmag, and the ρ(rmag) and v(rmag) are wind density and wind velocity at rmag, respectively. μ is the magnetic moment of the neutron star. For rmag > racc, the inflows cannot experience a significant gravitational field. Therefore, the wind speed at rmag is v(rmag) = vw. Assuming a nonmagnetized spherically symmetric wind (Elsner & Lamb 1977), the wind density ρ(rmag) reads ρ(rmag) ≈ w/4πa2vw (Davies & Pringle 1981; Bozzo et al. 2008), where a is the orbital separation of the two components in the system and a ≫ rmag is assumed. Therefore, the magnetospheric radius is (3)where w, − 6 = w/(10-6   M/yr), μ33 = μ/(1033 G cm3), and a10 is related to the orbital separation as a ~ 1012a10 cm  cm. Pb,10 is the orbital period Pb in units of 10 days, and M30 is the companion mass in units of 30 solar masses. When rmag < racc, the gravitational field of the neutron star dominates the falling of wind matter. The wind density and wind velocity at rmag can therefore be taken as and , respectively. The accretion rate  depends on racc and a according to (Frank et al. 2002) (4)Then the magnetospheric radius is given by (5)A corotation radius rcor can be defined because the spin angular velocity (Ωs = 2π/Ps) of the neutron star is equal to the Keplerian angular velocity () of matter being accreted, (6)where Ps,3 is the spin period of the neutron star Ps in units of 103 s. The competition among the three characteristic radii leads to different accretion regimes, so do the accretion modes and the behavior of the neutron star.

2.2. Pure propeller

In a system in which the magnetospheric radius is larger than the accretion radius (rmag > racc), the stellar wind, flowing around the magnetosphere, may not experience a significant gravitational field and interacts directly with the magnetosphere. The neutron star behaves like a pure propeller and spins down via rotational energy dissipation that results from the interaction between the magnetosphere and stellar winds, similar as in radio pulsars with magnetic inhibition (Stella et al. 1986; Bozzo et al. 2008). The spin-down torque can be expressed as (Lynden-Bell & Pringle 1974; Wang 1981; Lipunov 1992) (7)where κt ≤ 1 is a dimensionless parameter of the order of unity (Davies & Pringle 1981; Lipunov 1992). The equation governing the spin evolution can be written as (Davidson & Ostriker 1973; Ghosh & Lamb 1978) (8)where I is the moment of inertia of the neutron star, and T is the total torque imposed on the neutron star. Therefore, the spin-down rate reads (9)where I45 is the moment of inertia in units of 1045 g cm2. Accordingly, if the neutron star spins down to a period of 1000 s, the spin-down timescale is given by (10)The luminosity in this regime results from the shock forming near rmag, (11)

2.3. Retrograde propeller – quasi-spherical capture

If the magnetospheric radius is smaller than the accretion radius, the stellar wind radiated from the companion can penetrate through the accretion radius and halt at the magnetopause. We consider the regime of rmag > rcor in this part. The inflow in the region between racc and rmag falls approximately in a spherical configuration. The neutron star behaves like a supersonic rotator, i.e., the linear velocity of rotation of magnetosphere is much higher than the velocity of sound in the wind clumps, and ejects some of that wind material out of the accretion radius, leading to the dissipation of some rotational energy and thus spinning down. This scenario corresponds to the propeller mechanism and imposes a spin-down torque on the neutron star. The torque resulting from the propeller mechanism reads (Lipunov 1992) (12)In addition, due to the high speed of strong stellar wind, the accumulation of accreted matter is faster than the ejection rate by the propeller mechanism. More and more wind clumps may be deposited outside the magnetosphere (Stella et al. 1986). The deposit is disordered by the supersonic rotation of the neutron star, forming turbulent motion around the magnetosphere. The velocity of turbulent motion of accumulated stellar wind at the magnetoaphere boundary is close to the velocity of sound (Lipunov 1992). Because of the disorder of the stellar wind and the turbulent motion (Davies & Pringle 1981), the flowing accumulated material can be either prograde or retrograde around the magnetosphere. If the wind flows retrogradely, it will impose an inverted accretion torque on the neutron star and spin it down. This accretion torque can be written as (Lipunov 1992) (13)Accordingly, the total torque imposed on the neutron star in this regime is (14)Substituting Eq. (14) into (8), we can obtain the spin evolutionary law in the regime of racc > rmag > rcor, (15)In this scenario, the neutron star can spin down to 1000 s during a timescale of (16)However, if the stellar wind is not strong enough, the accreting matter is approximately in radial free fall as it approaches the magnetopshere boundary. The captured material cannot accumulate near the magnetosphere due to the supersonic rotation of the neutron star. Consequently, the accretion torque contains only the spin-down torque imposed by the propeller mechanism. This alternative scenario is discussed in Urpin et al. (1998b,c). It is less relevant for SFXTs, since the OB giant companions emit strong stellar wind.

In the retrograde propeller phase, both the interaction with magnetic fields by the propeller mechanism and inverted accretion contribute to the luminosity, which reads (17)In addition, the wind material halts at the magnetopause and forms a shock near rmag. Energy is also released through the shock, and the corresponding luminosity is given by Eq. (17) in Bozzo et al. (2008).

2.4. Retrograde accretor – quasi-disk accretion

When racc > rmag and rcor > rmag, the disordered stellar wind penetrates the accretion radius and falls directly toward the magnetosphere. In this case, the velocity of sound at the magnetosphere boundary is much higher than the linear velocity of the neutron star rotation, and the neutron star is a subsonic rotator. Owing to the relatively slow rotation of the neutron star, the rotating magnetosphere cannot eject the wind material. The inflows experience a significant gravitational field and accumulate near the magnetosphere. The accumulated material, with high initial velocity, can form an extended adiabatic, tenuous, and retrogradely rotating atmosphere around the magnetosphere, which impose an inverted accretion torque and spin the neutron star down. In addition, the reversal rotation between the adiabatic atmosphere and the magnetosphere may cause a flip-flop behavior of some clumps in the atmosphere, ejecting a part of material and leading to the loss of thermal energy.

The spin-down torque resulting from the flip-flop ejection has the same form as Eq. (7) (Davies & Pringle 1981; Lipunov 1992). The accretion spin-down torque imposed by the inverted atmosphere has the same expression as Eq. (13). Accordingly, the total spin-down torque imposed on the neutron star in the regime of racc > rmag and rcor > rmag is (18)Substituting the above expression into Eq. (8), we obtain the spin evolution equation, (19)The spin-down (to a period of 1000 s) timescale reads (20)The corresponding luminosity is (21)

3. Application to SFXTs

Some SFXTs that contain sporadically accreting neutron stars are believed to be X-ray pulsars and display very different pulse periods, from a few seconds (e.g. IGR J16479-4514, Ps = 2.14 s, http://hera.ph1.uni-koeln.de/~heintzma/Integral/SFXT.htm; IGR J18483-0311, Ps = 21 s, Sguera et al. 2007; IGR J17544-2619, Ps = 71.49 s, Drave et al. 2012; AX J1749.1-2733, Ps = 132 s, Karasev et al. 2007, 2008; IGR J111215-5952, Ps = 187 s, Swank et al. 2007; IGR J16465-4507, Ps = 228 s, Lutovinov et al. 2005) to a few thousand seconds (e.g. IGR J16418-4532, Ps = 1240 s, Walter et al. 2006; Sidoli et al. 2012; XTE J1739-302, reported for periodicity in the 1500–2000 s range, Sguera et al. 2006). Their orbital periods also span a very wide range, from 1.4 days to a few hundred days (e.g. Sidoli 2011b, and references therein).

Table 1

Different accretion regimes for eight SFXTs.

According to the definition of the three aforementioned characteristic radii (see Eqs. (1), (5) and (6)), they depend on some physical properties of the system, including the neutron star mass and magnetic field strength, stellar wind velocity, orbital period of the binary system, the mass loss rate of stellar wind, the mass of the companion, and the neutron star spin period. In this work we take the neutron star mass as 1.4   M and magnetic field as 1012 G, the mass of OB supergiant as 30   M, the mass loss rate and velocity of the stellar wind as 10-6   M yr-1 and 1000 km s-1, respectively. Then the accretion radius is about 1010 cm, and the corotation radius depends only on the spin period as . The magnetospheric radius depends also only on the orbital period according to in the region of rmag > racc and when rmag < racc. Therefore, we can determine the accretion regime in which a certain SFXT is located, according to the spin period and orbital period (see Table 1).

IGR J18483-0311 and IGR J17544-2619 display not-too-slow rotation and relatively short orbital periods (for IGR J18483-0311, Ps = 21 s, http://hera.ph1.uni-koeln.de/~heintzma/Integral/SFXT.htm, Porb = 18 days, Jain et al. 2009; for IGR J17544-2619, Ps = 71.49 s, Drave et al. 2012, Porb = 4.9 days, Clark et al. 2010). From the expressions of rmag and rcor (see Eqs. (5) and (6)), one can see that the relation among these three radii is racc > rmag > rcor. This means that IGR J18483-0311-like sources can behave like a propeller in the retrograde scenario.

Very long spin period SFXTs with narrow orbits, e.g. IGR J16418-4532 (Ps = 1240 s, Walter et al. 2006; Sidoli et al. 2011a; Porb = 3.74 days, Levine et al. 2011), typically have a large corotation radius. Moreover, the magnetospheric radius is one order of magnitude smaller than the accretion radius. Therefore, this class of systems can be an accretor with retrograde accretion. In addition, IGR J16479-4514, with Ps = 2.14 s (http://hera.ph1.uni-koeln.de/~heintzma/Integral/SFXT.htm) and Porb = 3.3 days (Jain et al. 2009), also present a somewhat larger corotation radius. Therefore, it can also behave like an retrograde accretor.

There are several slowly rotating systems with spin periods of a few hundred seconds and relatively wide orbital separation, such as IGR J111215-5952 (Ps = 187 s, Swank et al. 2007; Porb = 165 days, Sidoli et al. 2007), AX J1749.1-2733 (Ps = 132 s, Karasev et al. 2007, 2008; Porb = 185.5 days, Zurita & Chaty 2008), and IGR J16465-4507 (Ps = 228 s, Lutovinov et al. 2005; Porb = 30.32 days, Clark et al. 2010). As listed in Table 1, the relations among the three aforementioned characteristic radii of these sources are rmag > racc > rcor. Therefore, these three sources can be in the pure propeller regime.

The outbursts of XTE J1739-302 observed by several observatories show a characteristic timescale for variability of about 1500–2000 s (Sguera et al. 2006). This source has an orbital period of 51.47 days (Drave et al. 2010). Both the corotation radius and magnetospheric radius are slightly larger than the accretion radius (see Table 1 for details). However, due to the uncertainty in the spin period, it is difficult to determine its exact status. More accurate knowledge about the properties of the system is needed to determine whether the system is a pure propeller or a retrograde propeller.

It can be seen that the three characteristic radii in several sources have the same order of magnitude (see Table 1 for details). IGR J16418-4532 shows a dynamical range of two orders of magnitude and quasi-periodic flares (Sidoli 2011b). These features were taken to suggest that its X-ray emission is driven by a transitional accretion regime in a prograde accretion scenario (Sidoli 2011b). In this sense, we claim that this class of SFXTs frequently experience a transitional state between the propeller regime and the retrograde accretor phase.

The above distinction depends strongly on the assumption that these systems have the same properties except for spin and orbital periods. The particularities for each system should be investigated for a more detailed study.

4. Summary

Neutron stars in a wind-fed accreting system with disordered and inhomogeneous stellar winds can be either a prograde or a retrograde accretor. An accretor may exhibit different behaviors in different accretion regimes, depending on the competition among the three characteristic radii.

When the magnetospheric radius is located outside the accretion radius, the wind cannot experience significant gravitational focusing and forms a bow shock near the magnetosphere. The interaction between the stellar wind and the magnetosphere contributes to the loss of rotational energy and the spin-down of neutron star. This is the pure propeller regime. If the accretion radius is larger than the magnetospheric radius, the wind flowing from OB companions with a high velocity directly interacts with the magnetosphere and forms a shock at rmag (Harding & Leventhal 1992). In this case, when rmag is larger than the rcor, the neutron star can behave like a retrograde propeller. The shock, formed near the magnetosphere, flip-flops from side to side and makes the uniform wind disordered and inhomogeneous. Therefore, the neutron star accretes matter with a relatively high angular velocity, but with opposite signs (Livio 1992), consistent with the retrograde scenario. If both the accretion radius and corotation radius are larger than the magnetospheric radius, the inflowing matter can fall directly and halt near the magnetopause. Because of the high falling velocity, there is no effective cooling mechanism. Moreover, the dissipation of the rotational energy also produces additional heat. The spherically symmetric accreting material forms a quasi-disk, adiabatic and tenuous atmosphere around the neutron star, which can retrogradely rotate with the magnetosphere and spin the neutron star down. In the case of retrograde scenario, the neutron star can spin down to a period of 1000 s during a timescale of about 103–104 yr.

Acknowledgments

It is a pleasure to appreciate enlightening suggestions from S. N. Zhang. This work was supported by the National Science Council of the Republic of China under grant NSC 99-2112-M-007-017-MY3.

References

All Tables

Table 1

Different accretion regimes for eight SFXTs.

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.