A&A 447, 813-825 (2006)
DOI: 10.1051/0004-6361:20052689
J. Ferreira1 - P.-O. Petrucci1 - G. Henri1 - L. Saugé1,2 - G. Pelletier1
1 - Laboratoire d'Astrophysique, Observatoire de Grenoble, BP 53,
38041 Grenoble Cedex 9, France
2 -
Institut de Physique Nucléaire de Lyon, 43 Bd 11 novembre 1918, 69622 Villeurbanne Cedex,
France
Received 12 January 2005 / Accepted 19 October 2005
Abstract
We present a new picture for the central regions
of Black Hole X-ray Binaries. In our view, these central regions have a
multi-flow configuration which consists in (1) an outer standard
accretion disc down to a transition radius rJ; (2) an inner magnetized
accretion disc below rJ driving (3) a non relativistic self-collimated
electron-proton jet surrounding, when adequate conditions for pair
creation are met; (4) a ultra relativistic electron-positron beam.
This accretion-ejection paradigm provides a simple explanation to the
canonical spectral states, from radio to X/-rays, by varying the
transition radius rJ and disc accretion rate
independently. Large values of rJ correspond to the
Quiescent state for low
and the Hard state for larger
.
These states are characterized by the presence of a steady electron-proton MHD jet emitted by the disc below rJ. The
hard X-ray component is expected to form at the jet basis. When rIbecomes smaller than the marginally stable orbit ri, the whole disc
resembles a standard accretion disc with no jet, characteristic of the
Soft state. Intermediate states correspond to situations where
.
At large
,
an unsteady pair cascade process is triggered
within the jet axis, giving birth to flares and ejection of relativistic
pair blobs. This would correspond to the luminous intermediate state,
sometimes referred to as the Very High state, with its associated
superluminal motions.
The variation of rJ independently of
is a necessary ingredient in this picture. It arises from the presence of a large scale vertical magnetic field threading the disc. Features such as possible hysteresis and the presence of quasi-periodic oscillations would naturally fit within this new framework.
Key words: black hole physics - accretion, accretion disks - magnetohydrodynamics (MHD) - ISM: jets and outflows - X-rays: binaries
Galactic Black Hole X-ray Binaries (hereafter BH XrBs) are binary systems that were first detected as X-ray sources. They harbor a massive compact object as a primary star, being therefore a black hole candidate, accreting matter from the companion. The X-ray emission probes the inner regions around the compact object and is interpreted as a signature of accretion (Tanaka & Lewin 1995; McClintock & Remillard 2003, hereafter McCR03, and references therein). There is now growing evidence that these objects also display ejection signatures: radio emission is commonly interpreted as the presence of compact steady jets or the sporadic ejection events also seen in infrared and X-ray bands (Buxton & Bailyn 2004; Fender & Belloni 2004). In fact, the correlation between radio luminosity (ejection) and X-ray luminosity (accretion) found in Active Galactic Nuclei (AGN) seems to be also consistent with XrBs (Robertson & Leiter 2004; Falcke et al. 2004; Gallo et al. 2003; Corbel et al. 2003; Choudhury et al. 2003). The similarity of the first detected galactic jet in 1E1740 with extragalactic jets gave birth to the name "microquasar'' (Mirabel & Rodriguez 1998).
What is so dramatic about microquasars is their multiple manifestations
through very different spectral states. They spend most of their
time in the Quiescent state which is characterized by a very low
accretion rate (
as low as
10-9). The multiwavelength
spectral energy distributions are thus very scarse but generally show a hard X-ray (2-10 keV)
spectrum with a power law photon index
and an optical/UV continuum reminiscent of a
104 K disc blackbody with strong emission lines (McCR03).
Occasionally, microquasars enter in outburst, resulting from a dramatic
increase of their accretion rate. During these outbursts, microquasars
show different canonical states, like the well known Hard and Soft
states. The hard state is characterized by a spectrum dominated above 2 keV by a hard power-law component (
)
with a cut-off around 100 keV (e.g. Grove et al. 1998; Zdziarski & Gierlinski 2003, and
references therein), and with a soft X-ray excess below 2 keV interpreted
as the presence of a cool (
0.01-0.2 keV) accretion disc. Strong
radio emission is observed during this state and some VLBI images
directly revealed spatially resolved structures. These are interpreted as non or only mildly relativistic (bulk Lorentz factor
)
steady jets (e.g. Dhawan et al. 2000; Stirling et al. 2001). On the contrary, the soft
state is dominated by a thermal blackbody-like component, typical of a
standard accretion disc of temperature ranging from 0.7 to 1.5 keV
(consistent with an inner disc radius
). A faint power-law component may still be present but with a steep photon index
(McCR03). This state is
also devoid of any radio emission which is interpreted as the absence of
jet (e.g. Corbel et al. 2000; Fender et al. 1999). Hard and soft state span a relatively large range in luminosity, i.e. from 10-2 to 1
.
Microquasars can also be observed in Intermediate states, generally during transitions between the hard and soft states. Intermediate states present relatively complex spectral and timing behaviors. Detailed studies have been done in the recent literature (see e.g. Belloni et al. 2005; Fender et al. 2004) and reveal a clear evolution with time between hard, variable (in X-ray) and radio-loud systems to softer, less variable and more radio-quiet ones. The transitions between the hard and soft "flavors'', which seems to correspond also to a transition between jet-producing and jet-free states, can be relatively abrupt especially at high luminosity level where they are apparently coincident with strong radio outbursts (Fender & Belloni 2004). These non steady ejection events display apparent superluminal velocities, indicative of a highly relativistic plasma (Dhawan et al. 2000; Mirabel & Rodríguez 1999). These different observational characteristics lead to the definition of the so-called "jet line'' by Fender & Belloni (2004), that separates hard/jet dominated states to soft/jet-free ones in hardness-intensity diagrams.
These various spectral states obviously carry a huge amount of
information about the physical processes behind accretion and
ejection. Moreover, since it is believed that microquasars are a scaled
down version of AGN, understanding the various states (along with their
transitions) in XrBs will certainly provide insights on the observed
differences between AGN (radio loud/radio quiet, FRI/FRII, blazars...).
![]() |
Figure 1:
A Standard Accretion Disc (SAD) is established down to a radius
rJ which marks the transition towards a low radiative Jet Emitting
Disc (JED), settled down to the last stable orbit. The JED is driving a
mildly relativistic, self-collimated electron-proton jet which, when
suitable conditions are met, is confining an inner ultra-relativistic
electron-positron beam. The MHD power
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The puzzle of the existence of these different spectral states is
enhanced by the fact that each state must correspond to a dynamically
steady state. Indeed, each state lasts a much longer time than the
inferred dynamical time scale. Let us consider the case of GRS1915+105,
which is the BH XrB with the most rapid time scales. It shows states
lasting 103 s with rising and decay times of the order of one
second, whereas the Keplerian orbit time scale is several milliseconds at
the inner radius (Belloni et al. 1997). One must therefore explain why a system
switches from one stationary state to another one. But before that, one
must first identify the relevant underlying dynamics describing each
state. Up to now, these canonical spectral states are still not fully
understood, let alone the transitions between them.
The most common paradigm used to interpret observations is based on the low radiative
efficiency of the ADAF model (Esin et al. 1997, and references
therein). Within this framework, the highest energy component is due to
the inner thick, low radiative disc whereas an outer standard disc is
responsible for the UV-soft X ray emission. By varying the transition
radius rtr between these two discs (as a function of )
one
gets reasonable fits to the spectral energy distributions (hereafter
SEDs, e.g. Narayan et al. 1996; Esin et al. 1998; Hameury et al. 1997; Esin et al. 2001). In the ADAF paradigm, the accretion power below the transition radius is
essentially stored as thermal energy of protons and eventually advected
below the black hole event horizon. Remarkably, energetic ejection events
appear to be always associated with hard and intermediate states and are
quenched during the soft phase (McCR03; and references therein). This has
led several authors to propose that ADAFs could indeed drive outflows or
its ADIOS extension (Blandford & Begelman 2004,1999). However, it must be realized
that there is only one source of energy, namely the release of
gravitational energy through accretion. Hence, whenever a disc is capable
of driving jets, these will carry away a fraction of the released
accretion energy. As a consequence, the disc luminosity will be
quenched. Observations tell us something consistent with this very simple
and unavoidable argument: whenever a steady jet is formed the disc may
observationally disappear. This does not mean that the disc is
really disappearing (i.e., inner region being depleted of its mass), only
that we do not see it anymore. Discs that drive jets are indeed radiating only a small fraction of the
accretion power released as shown in another class of accretion solution, namely Magnetised Accretion-Ejection Structures (hereafter MAES, see e.g. Ferreira & Pelletier 1993,1995; and Ferreira 2002, for a review). In this flow, a large scale magnetic field is threading the disc, exerts a torque leading to accretion and allows the production of self-confined jets. A logical consequence is that whenever a jet is
formed, the ADAF hypothesis is no more useful to explain the low
radiative efficiency of the accretion flow.
The goal of this paper is to provide an alternative view explaining the canonical spectral states observed in BH XrBs based on MAES. We do not intend to address the crucial issue of the transitions between these states, neither time scales involved nor possible hysteresis effects. This is delayed for future work. In this paper, we only expose the global physical picture and show that this framework is rich enough to explain all known spectral components. In a forthcoming paper, we will provide calculations of SEDs and compare them to observations.
The paper is then organized as follows. Section 2 describes in some detail the four physical components present within our paradigm, their dynamical properties as well as radiative processes governing their emission. The canonical spectral states of BH XrBs are then interpreted within this framework in Sect. 3. Section 4 is devoted to a discussion of some time scales and timing properties of our model that could be related to some observed features. Section 5 highlights questions opened by our framework as well as future developments.
We assume that the central regions of BH XrB are composed of four distinct flows: two discs, one outer "standard'' accretion disc (hereafter SAD) and one inner jet emitting disc (hereafter JED), and two jets, a non-relativistic, self-confined electron-proton MHD jet and, when adequate conditions for pair creation are fulfilled, a ultra-relativistic electron-positron beam. A sketch of our model is shown in Fig. 1 while the four dynamical components are discussed separately below. This is an extended version of the "two-flow'' model early proposed for AGN and quasars (Sol et al. 1989; Pelletier et al. 1988; Pelletier & Roland 1989; Henri & Pelletier 1991; Pelletier & Sol 1992) to explain the highly relativistic phenomena such as superluminal motions observed in these sources. This model provides a promising framework to explain the canonical spectral states of BH XrBs mainly by varying the transition radius rJ between the SAD and the JED. This statement is not new and has already been proposed in the past by different authors (e.g. Livio et al. 2003; King et al. 2004; Belloni et al. 1997; Esin et al. 1997) but our model distinguish itself from the others by the consistency of its disc-jet structure and by the introduction of a new physical component, the ultra-relativistic electron-positron beam, that appears during strong outbursts.
We believe that jets from BH XrBs are self-collimated because they follow the same accretion-ejection correlation as in AGN (Merloni et al. 2003; Fender et al. 2003; Corbel et al. 2003). This therefore implies the presence of a large scale vertical field anchored somewhere in the accretion disc. We think it is unlikely that such a field has a patchy distribution on the disc. Indeed, once a jet is launched, it exerts a torque on the underlying disc with a corresponding vertical flux of angular momentum and energy. This torque writes
where
is the toroidal field at the disc surface. This shows that the product Jr Bz must remain positive over the whole region driving jets. Such a condition is unlikely to be met with a vertical field changing its polarity from one zone to another because the currents (both radial and azimuthal) induced inside the disc will tend to cancel each other. The most efficient way to launch jets from an accretion disc is probably from a radial extension (the JED) with a large scale Bz of the same polarity.
This vertical field is therefore an unavoidable ingredient to produce self-collimated jets. So, where does this field come from? A first possible origin is in-situ generation of magnetic fields by dynamo. The huge difficulty is to provide an ordered large scale vertical field out of turbulent small scale seed fields. Opening of magnetic loops by the disc differential rotation might provide a mechanism
(Romanova et al. 1998). The second origin is advection of magnetic field by the
accreting material. We therefore assume that the whole accretion disc is pervaded by a large scale magnetic field Bz.
The presence of a large scale vertical field threading the disc is however not sufficient to drive super-Alfvénic jets. This field must be close to equipartition as shown by Ferreira & Pelletier (1995) and Ferreira (1997). The reason is twofold. On one hand, the magnetic field is vertically pinching the accretion disc so that a (quasi) vertical equilibrium is obtained only thanks to the gas and radiation pressure support. As a consequence, the field cannot be too strong. But on the other hand, the field must be strong enough to accelerate efficiently the plasma right at the disc surface (so that the slow-magnetosonic point is crossed smoothly). These two constraints can only be met with fields close to equipartition.
An important local parameter is therefore the disc magnetization
where
includes the plasma and radiation pressures. In our picture, a SAD is established down to a radius rJ where
becomes of order unity. Inside this radius, a JED with
is settled.
At any given time, the exact value of rJ depends on highly non-linear processes such as the interplay between the amount of new large scale magnetic field carried in by the accreting plasma (e.g. coming from the secondary) and turbulent magnetic diffusivity redistributing the magnetic flux already present. These processes are far to be understood. For the sake of simplicity, we will treat in the following rJ as a free parameter that may vary with time (see Sect. 3).
Accretion requires the presence of a negative torque extracting angular
momentum. In a SAD this torque is assumed to be of turbulent origin and
provides an outward transport of angular momentum in the radial
direction. It has been modeled as an "anomalous'' viscous torque of
amplitude
,
where
is a small
parameter (Shakura & Sunyaev 1973). The origin of this turbulence is now commonly
believed to arise from the magneto-rotational instability or MRI
(Balbus & Hawley 1991). The MRI requires the presence of a weak magnetic field
(
)
and is quenched when the field is close to equipartition. We
make the conjecture that a SAD no longer exists once
reaches
unity. We show below that this is very likely to occur in the innermost
regions.
The radial distribution Bz(r) is provided by the induction equation which describes the interplay between advection and diffusion. If we assume that, in steady state, the poloidal field is mostly vertical (no significant bending within the SAD) then this equation writes
![]() |
(1) |
![]() |
(3) |
![]() |
(4) |
To summarize, the accretion flow at r> rJ is a SAD with
fueled by the companion's mass flux and driving no outflow (constant accretion rate
). The global energy budget is
where
This inner region with
is fueled by the SAD at a rate
.
Since it undergoes mass loss,
the JED accretion rate is written as
![]() |
(6) |
![]() |
(8) |
This dynamical property has a tremendous implication on the JED
emissivity. The JED luminosity comes from the accretion power dissipated
within the disc by turbulence and transported away by photons, so
.
This dissipated power is very difficult to
estimate with precision because it requires a thorough description of the
turbulence itself. Thus, one usually uses crude estimates based on
"anomalous'' turbulent magnetic resistivity
(Joule heating) and
viscosity
(viscous heating). This translates into
where integration is made over the whole volume
occupied by the JED. The importance of local "viscous'' dissipation with
respect to the MHD Poynting flux leaving the disc is approximately given
by
![]() |
(9) |
![]() |
(10) |
To summarize, under quite general conditions on the turbulence within magnetized discs, most of the available accretion energy is powering the outflowing plasma (Ferreira & Pelletier 1993,1995). This is in strong contrast with ADAFs where the accretion power is stored as heat advected by the accreting plasma. In this case, low luminosity discs can be obtained as long as the central object possesses an event horizon. However, the power to magnetically drive jets is also missing. In the case of MAES, the JED is weakly dissipative while powerful jets are being produced regardless of the nature of the central object.
Complete calculations of MAES showed that isothermal or adiabatic
super-Alfvénic jets from Keplerian accretion discs were possible only
when a tiny fraction of the accreted mass is locally ejected. This
translates into a small ejection efficiency, typically
(Ferreira 1997; Casse & Ferreira 2000a). When some heat deposition occurs at the JED upper
layers a typical value of
becomes possible, even up to 0.5 but never reaching unity, in agreement with Eq. (7)
(Casse & Ferreira 2000b). This is a much lower mass loss than that assumed in ADIOS
models (Blandford & Begelman 1999).
The fact that the jet torque largely dominates the turbulent torque provides another striking difference between the internal structures of SADs and JEDs. Indeed, the angular momentum conservation provides a sonic Mach number measured at the disc midplane
Mass conservation in the JED writes
n | = | ![]() |
|
![]() |
![]() |
(12) |
![]() |
(14) |
![]() |
= | ![]() |
(15) |
![]() |
![]() |
![]() |
(16) |
![]() |
![]() |
![]() |
(17) |
![]() |
= | ![]() |
(18) |
The ejection to accretion rate ratio in a JED writes
.
In principle, the ejection efficiency
can be observationally deduced from the terminal jet speed. Indeed,
the maximum velocity reachable along a magnetic surface anchored on a
radius
(between ri and rJ) is
in the non-relativistic limit (see Ferreira 1997, for
relativistic estimates). Although a large power is provided to the
ejected mass (mainly electrons and protons), the mass loss (
)
is
never low enough to allow for speeds significantly relativistic required
by superluminal motions: MHD jets from accretion discs are basically non
or only mildly relativistic with
(Ferreira 1997). This is basically the reason why they can be efficiently
self-confined by the magnetic hoop stress. Indeed, in relativistic flows
the electric field grows so much that it counteracts the confining effect
due to the toroidal field. This dramatically reduces the self-collimation
property of jets (Pelletier 2004; Bogovalov 2001; Bogovalov & Tsinganos 2001).
Calculations of jets crossing the MHD critical points have been undergone under the self-similar ansazt (Ferreira & Casse 2004; Casse & Ferreira 2000b). In these calculations, the emission of the MHD jet has been neglected and all the available power is converted into ordered jet kinetic energy. However, a fraction of this power is always converted into heat and particle acceleration, leading to emission. In our case, jets from MAES have two distinct spectral components and the resultant SED may therefore be quite intricate. Producing a global SED is out of the scope of the present paper. It requires to fix several parameters, which is legitimate only by object fitting. This is postponed for future work.
We expect a small fraction of the jet power
to be converted
into particles, through first and/or second order Fermi acceleration,
populating the MHD jet with supra-thermal particles. These particles are
responsible for the bulk emission of the MHD jet. This is similar to
models of jet emission already proposed in the literature
(Vadawale et al. 2001; Falcke et al. 2004; Markoff et al. 2003; Markoff 2004; Markoff et al. 2001; Falcke & Biermann 1995).
In these models, the jet is assumed to be radiating self-absorbed
synchrotron emission in the radio band becoming then optically thin in
the IR-Optical bands and providing a contribution up to the X/
-rays.
A flat or even inverted spectrum index in the radio band is quite easily achieved by self-collimated jets for reasonable values of the parameters (e.g. the exponent p of the power law particle distribution).
Note that the MHD jet due to the MAES yields
.
In the framework of Heinz & Sunyaev (2003), this implies that such a jet would provide correlated radio and X-ray emissions close to the observed law, namely
Gallo et al. (2003); Corbel et al. (2003). However, the spectrum index in the X-ray band would not be steep enough, even taking into account the cooling of the particles. Moreover, the fundamental plane of BH activity of Merloni et al. (2003), namely the correlation between mass, radio and X-ray fluxes, cannot be explained by such synchrotron jets (Heinz 2004). Following these authors, we conclude that there must be another significant contribution to the X-ray emission in the Low/hard state. The same conclusion was independently reached by
Rodriguez et al. (2004) who suggested that, to explain the energy dependence of the
quasi-periodic oscillation (QPO) amplitude in GRS 1915+105, the high
energy spectrum of the source must be the sum of different emission
processes. Another argument comes from the study of the overall spectral energy distributions.
At least in some objects, the extrapolation of the X-ray power-law spectrum towards the optical and infrared bands is above the observed fluxes, which shows that hard X-rays cannot be direct synchrotron radiation from the jet (Kalemci et al. 2005).
Moreover, BH XrBs in the hard state generally exhibit a spectrum with a high energy cut-off around 100 keV (e.g. Grove et al. 1998; Zdziarski & Gierlinski 2003, and references therein). While naturally obtained if the emission is thermal (i.e. a comptonized corona), non-thermal emission requires a fine tuning of the parameters that we found questionable. However, our framework naturally provides another contribution to the high energy emission as long as a JED is present.
Jet production relies on a large scale magnetic field anchored on the
disc as much as on MHD turbulence triggered (and sustained) within
it. This implies that small scale magnetic fields are sheared by the disc
differential rotation, leading to violent release of magnetic energy at
the disc surface and related turbulent heat fluxes (e.g. Merloni & Fabian 2002; Galeev et al. 1979; Heyvaerts & Priest 1989; Stone et al. 1996). The energy released is actually tapping the MHD
Poynting flux flowing from the disc surface. We can safely assume that a
fraction f of it would be deposited at the jet basis,
with a total power
.
The dominant cooling term in this
optically thin medium is probably comptonization of soft photons emitted
by the outer SAD (with a small contribution from the underlying JED). These are
circumstances allowing a thermal plasma to reach a temperature as high as
100 keV, (Pietrini & Krolik 1995; Mahadevan 1997; Esin et al. 1997).
The computation of the exact spectral shape produced by this "corona''
through thermal comptonization requires sophisticated computations
(e.g. Haardt 1993; Poutanen & Svensson 1996) which are out of the scope of this
paper. Instead, a cut-off power law shape is generally used as zero-order
approximation. In this case, the high energy cut-off is rougly equal to
twice the plasma temperature (see e.g. Petrucci et al. 2000, for more
discussion). The power law photon index can also be approximated by a
simple fonction of the Compton amplification factor A, which is equal
to the ratio of the total luminosity outgoing from the jet basis to the soft
luminosity
entering in it (see e.g. Beloborodov 1999; Malzac et al. 2001):
In the most general case,
should include the SAD and JED
emissions but also the reprocessed radiation from the discs (both JED and SAD)
that are partly intercepted by the corona. In consequence, the photon
index depends implicitly on parameters like f,
,
rJ but
also on the bulk motion of the corona in a complex manner. This will be
precisely discussed in a forthcoming paper where calculations of SEDs will
be provided. Simple estimates can be given, however, in the case of large
rJ and
,
since in these conditions the SAD and JED emissions
become negligible compared to the reprocessed one. This situation has
been precisely studied by Malzac et al. (2001). The parameters C and
of
Eq. (19) obtained by these authors are equal to 2.19 and 2/15
respectively. Hard X-rays photon indexes in the range 1.4-2 are easily obtained
for different corona velocities and corona aspect ratios
(i.e. height/width). These values are in good agreement with what is
generally observed in the hard states where we expect large rJ(cf. Sect. 3 for a more detailed discussion). We can note
also that a decrease of rJ will result in a larger
,
due to
the increase of the SAD emission, and thus in a softening of the X-ray
spectrum.
Since the large scale magnetic field driving the self-confined jet is anchored onto the accretion disc which has a non zero inner radius, there is a natural hole on the axis above the central object with no baryonic outflow (this also holds for neutron stars). This hole provides a place for pair production and acceleration with the outer MHD jet acting as a sheath that confines and heats the pair plasma. This is the microquasar version of the "two flow'' model that has been successfully applied to the high energy emission of relativistic jets in AGNs (Renaud & Henri 1998; Marcowith et al. 1995; Henri & Pelletier 1991; Marcowith et al. 1998).
The
plasma is produced by
interaction, the
-ray photons being initially produced by a few relativistic
particles by Inverse Compton process, either on synchrotron photons (Synchrotron Self Compton or SSC) or on disc photons (External Inverse Compton or EIC). Detailed models for AGNs have shown that all processes can contribute, depending on the physical parameters of the system (magnetic field, disc luminosity, distance). We do not intend to build an explicit Spectral Energy Distribution of the pair plasma here, but we will just discuss the general mechanism of pair beam formation and the relevance of each process in explaining the various non thermal components.
It is well known that above 0.5 MeV photons can annihilate with themselves to produce an electron-positron pair. Usually, pairs are assume to cool once they are formed, producing at turn non thermal radiation. Some of this radiation can be absorbed to produce new pairs, but the overall pair yield never exceeds 10%. A key point of the two-flow model however is that the MHD jet launched from the
disc can carry a fair amount of turbulent energy, most probably through
its MHD turbulent waves spectrum. A fraction of this power can
be transferred to the pairs (
). Thus the freshly created pairs can be continuously reheated, triggering an efficient pair runaway process leading to a dense pair plasma (Henri & Pelletier 1991).
As we said, reacceleration is balanced by cooling through the combination of synchrotron, SSC and EIC processes. Synchrotron and SSC emission are quasi isotropic in the pair frame, but the external photon field is strongly anisotropic. The pair plasma will then experience a strong bulk acceleration due to the recoil term of EIC, an effect also known as the "Compton Rocket'' effect (Renaud & Henri 1998; O'Dell 1981). As was shown in previous works, this "rocket'' effect is the key process to explain relativistic motions (Renaud & Henri 1998; Marcowith et al. 1995).
At a given distance of the disc,
the bulk acceleration saturates at a characteristic Lorentz factor, depending only on
the radiation angular distribution. It is defined by the condition that the net radiation flux in the comoving frame (after Lorentz transform) vanishes. It can be shown easily that this characteristic Lorentz factor is approximately
on the axis of a standard accretion disc (Renaud & Henri 1998), as long as the distance z verifies
:
here the relevant disc inner radius is in fact the transition radius rJ. Noticeably this value does not depend on the disc luminosity (or accretion rate): it is only dependent on the angular distribution of the intensity, i.e. the radial dependency of the temperature
.
The photon field is in fact dominated by photons emitted at
.
The modification introduced by the JED in the central region is likely to be immaterial for two reasons. First although the luminosity decreases, the radial dependency remains almost unchanged. Second, as we argue below, the final bulk Lorentz factor depends only on the distance where the plasma decouples from the radiation, and not on the motion close to the core.
A pure pair plasma will thus experience a continuous bulk acceleration, the bulk velocity increasing slowly with the distance. At some point, the radiation field becomes too weak to act efficiently: this happens when the relaxation time towards the equilibrium Lorentz factor becomes larger than the dynamical time z/c. At this distance, the acceleration process stops and the plasma decouples from the external radiation field, moving on at a ballistic constant Lorentz factor
.
The asymptotic Lorentz factor depends essentially on the location of this critical distance.
In the original O'Dell's version of this process, the pair plasma was not
supposed to be reheated and this effect has been shown to be rather
inefficient, because cooling is always much faster than acceleration
(Phinney 1987). In fact, for a cold plasma, the mechanism reduces to the ordinary radiation pressure. Under these conditions, the critical distance is approximately
,where
is the soft photon compactness. For a near Eddington accreting disc,
and
(Phinney 1987). Although this is indeed a relativistic motion (an apparent superluminal motion is possible), this may not be high enough to account for high values around 5, as observed in microquasars.
In the two flow model however, continuous reheating of the pairs makes
the bulk acceleration more efficient, acting thus over a much larger distance: the radiation force is multiplied by
,
where
is the random (or relativistic temperature) of the pair plasma, (not to be confused with the bulk Lorentz factor
). Although the equilibrium Lorentz factor at a given distance is unchanged, the critical decoupling distance is much larger. The asymptotic
bulk Lorentz factor becomes
and values of 5 to 10 can be easily reached in near-Eddington accretion regime around stellar black
holes (Renaud & Henri 1998).
Producing this pair plasma requires thus altogether a strong MHD jet, a
radiative non-thermal component extending above the MeV range and a minimal
optical depth, namely
.
The non thermal component can indeed be associated with the steep power law observed during the intermediate states, given the fact that it seems to extend to MeV range without any break (McCR03). It is most probably due to Inverse Compton process on the disc photons. Indeed, the optical depth
for absorbing photons with energy
is approximately for a
spherical source of radius R filled by soft photons with density
by unit reduced energy:
![]() |
(20) |
It is noteworthy that the pair beam is intrinsically highly variable and
subject to an intermittent behavior. Indeed, once the pair creation is
triggered, a regulation mechanism must occur to avoid infinite power of
the pair plasma and limit the pair run-away. This is probably accomplished by
the quenching of the turbulence (
vanishes) when most of its
energy is suddenly tapped by the catastrophic number of newly created
pairs. These pairs will therefore simply expand freely, confined by the
heavier MHD jet. One would then expect a flare in the compact region,
followed by the ejection of a superluminal radio component, analoguous to those observed in AGNs (Saugé & Henri 2005, A&A, submitted). Such a situation can repeat itself as long as the
required physical conditions are met. Alternatively, it may also that the formation of a dense pair beam destroys the surrounding MHD jet, explaining the disappearance of the compact jet after a strong ejection event.
From Sect. 2, it is clear that the spectral appearance of a BH XrB critically depends
on the size of the JED relative to the SAD, namely rJ. As stated before, rJ is the
radius where the disc magnetization
becomes of order unity.
Thus, rJ depends on two quantities
and Bz(r,t). The total pressure is directly proportional to
since
.
As a consequence, any variation of the
accretion rate in the outer SAD implies also a change in the amplitude of the
total pressure. But we have to assume something about the time evolution of the large scale magnetic field threading the disc. Within our framework,
and Bz are two quantities that may vary independently with time.
Let us assume that Bz does not depend on the disc accretion rate. If, for instance,
undergoes a sudden increase (triggered by, e.g. some disc instability at the outer radii Lasota et al. 1996),
then there is a corresponding increase of
which is propagating
inwards, eventually reaching the inner JED. Since no magnetic flux is being
simultaneously added, the region where the vertical magnetic field
is close to equipartition shrinks, namely rJ decreases. If, on the
contrary, the disc accretion rate decreases (without a corresponding
decrease of the magnetic flux threading the disc), then the decrease of
requires some diffusion of the vertical magnetic field in the
inner regions in order to maintain equipartition, hence the JED expands and rJ increases.
According to this simplistic argument, one would expect an anti-correlation between rJ
and
,
namely a larger
implies a smaller rJ (and vice-versa).
Alternatively, one could also argue that a larger
implies more plasma within the disc and that a larger Bz would then be locally generated by dynamo. If such a process provides
,
then rJ would always remain unchanged, whatever
.
Another alternative could be advection of the companion's magnetic field along with the flow. Now, because of the stellar dynamo, such a field could also change with a time scale very different from that related to changes in
.
In any case, the amount and polarity advected along strong accretion phases would be an unknown function
.
The processes governing the amplitude and time scales of these
adjustments of rJ to a change in
are far too complex to be
addressed here. They depend on the nature of the magnetic diffusivity
within the disc but also on the radial distribution of the vertical magnetic field. We will
simply assume in the following that rJ and
are two independent parameters. In that respect, our view is very different from that of Mahadevan (1997); Esin et al. (1997) who considered only the dependency of
to explain the different spectral states of BH XrBs.
![]() |
Figure 2:
The canonical spectral states of BH X-ray binaries. a)
Quiescent state obtained with a low ![]() ![]() ![]() |
Open with DEXTER |
This state is characterized by a very low accretion rate (
as low
as
10-9) with a hard X-ray component. The ADAF model has been successfully applied to
some systems with a large transition radius between the ADAF and the
outer standard disc, namely
(e.g. Narayan et al. 1996; Hameury et al. 1997). However, such a model does not account for jets
and their radio emission, even though XrBs in quiescence seem also to
follow the radio/X-ray correlation (e.g. Gallo et al. 2004; Fender et al. 2003; Gallo et al. 2005).
Within our framework, a BH XrB in quiescence has a large rJ, so that a large zone in the whole disc is driving jets (Fig. 2a). The low
provides a low synchrotron jet luminosity, while the JED is optically thin, producing a SED probably very similar to that of an ADAF. We thus expect
.
The weak MHD Poynting flux prevents the ignition of the pair cascade process and no pair beam is produced.
It must be noted that slightly more complicated situations can arise depending on the actual value of .
For instance, the innermost denser regions of the JED could become optically thick for larger values of
,
i.e closer to the Hard state level.
Within our framework, the JED is now more limited radially than in the
Quiescent state, namely
(Fig. 2b). This transition radius corresponds to the inner
disc radius
as obtained within the SAD framework
(Zdziarski et al. 2004). Due to the higher
,
the JED may become optically
thick, which is required to explain the broad iron lines observed in some
binary systems (Sidoli & Mereghetti 2002; Frontera et al. 2001; Nowak et al. 2002). The low velocity of the plasma
expected at the jet basis is in good agreement with recent studies of
XrBs in Hard state (Gallo et al. 2003; Maccarone 2003). It can also explain the apparent
weakness of the Compton reflection (Zdziarski et al. 1999; Gilfanov et al. 1999) as already
suggested by Markoff et al. (2003, see also Beloborodov 1999; Malzac et al. 2001# and tested by
Markoff & Nowak (2004). In any case, the JED intrinsic emission is weak with
respect to that of the outer standard disc: most of the accretion power
flows out of the JED as an MHD Poynting flux. Nevertheless, the threshold
for pair creation is still not reached and there is no pair beam, hence
no superluminal motion. The MHD power is therefore shared between the jet
basis, whose temperature increases (the thermal "corona'') producing
X-rays, and the large-scale jet seen as the persistent (synchrotron)
radio emission.
Our interpretation of the Soft state relies on the disappearance of the
JED, i.e. when rJ becomes smaller or equal to ri (Fig. 2c). Depending on the importance of the magnetic flux in the disc, this may occur at different accretion rates. Thus, the threshold in
where there is no region anymore in the disc with equipartition fields may vary. The
whole disc adopts therefore a radial structure akin to the standard disc
model. Since no MHD jet is produced, all associated spectral signatures
disappear. Even if pair production may take place (when
is large), the absence of the confining MHD jet forbids the pairs to get warm enough and be accelerated: no superluminal motion should be
detected.
Note however that the presence of magnetic fields may be the cause of particle acceleration responsible for the weak hard-energy tail (McCR03, Zdziarski & Gierlinski 2003, and references therein).
This state has been first identified at large luminosities (
)
and was initially called Very High state. However, high
luminosity appeared to not be a generic feature since it has be observed at
luminosities as low as
(McCR03; Zdziarski & Gierlinski 2003). Therefore, the most prominent
feature is that these states are generally observed
during transitions between Hard and Soft states. Within our framework,
they correspond to geometrical situations where rJ is small but remains larger than ri(Fig. 2d). The flux of the outer standard disc is thus
important while the JED is occupying a smaller volume. The consequences on the spectral
shape are not straightforward since the importance of the different
spectral components relative to each other depends on the precise values
of rJ and
.
Such study is out of the scope of the present
paper and will be detailed elsewhere.
The crucial point however is that, in our framework, luminous
intermediate states (the so-called Very High State or VHS) with high
provide the best
conditions for the formation of the ultra-relativistic pair beam, as
described in details in Sect. 2.5: (1) a high luminosity; (2) a high energy steep power law spectrum extended up to the
-ray
bands; and (3) the presence of the MHD jet. The two first characteristics
enable a
opacity larger than unity (cf. Eq. (21) of
Sect. 2.5), while the MHD jet allows to confine the pair beam
and maintain the pair warm, a necessary condition to trigger a pair
runaway process. The total emission would be then dominated by the
explosive behavior of the pairs, with the sudden release of blobs. Each
blob produced in the beam first radiates in X and
-ray,
explaining the hard tail present in this state, and then, after a rapid
expansion, produces the optically thin radio emission. This pair beam would also explain the
superluminal ejections observed during this state in different
objects (e.g. Hannikainen et al. 2001; Sobczak et al. 2000). We conjecture that the exact moment where this occurs corresponds to the crossing of the "jet line'' recently proposed by (Fender et al. 2004, see also
Corbel et al. 2004). This corresponds to a transition from the "hard''
intermediate state to the "soft'' one.
The rapid increase of the pair beam pressure in the inner region of the MHD jet, during a strong outburst, may dramatically perturb the MHD jet production. Indeed, a huge pair pressure at the axis may enforce the magnetic surfaces to open dramatically, thereby creating a magnetic compression on the JED so that no more ejection is feasible. Alternatively, it is also possible that the racing of the pair process completely wears out the MHD Poynting flux released by the JED, suppressing the jet emission or even the jet itself. Whatever occurs (i.e. jet destruction or jet fading), we expect a suppression of the steady jet emission when a large outburst sets in. Interestingly, the detailed spectral and timing study of the radio/X-ray emission of four different black hole binaries during a major radio outburst (Fender et al. 2004) shows a weakening and softening of the X-ray emission as well as a the quenching of the radio emission after the burst. This is in good agreement with our expectations since the cooling of the pair beam should indeed results in a flux decrease and a softening of its spectrum.
Since both SAD and JED are quasi Keplerian, the first obvious time scale is the Keplerian orbit time, namely
![]() |
(22) |
Let us assume an increase in
triggered at some outer radius
.
This information propagates towards the center with the accretion flow, as a front of increased total pressure. The time scale involved is therefore
.
If we assume that rJ decreases (because of a decrease in
), then the JED (and its associated MHD jets) will progressively disappear. This evolution from a Hard state to an Intermediate or Soft state will be controlled by the advance of this front, namely
![]() |
(23) |
![]() |
(24) |
![]() |
(25) |
In our description of the canonical spectral states of XrBs we did not
mention the important issue of quasi-periodic oscillations or
QPOs. Low-frequency (0.1-30 Hz) QPOs in the X-ray bands (2-30 keV) are
indeed observed in the Hard and luminous Intermediate (VHS) states. In the latter state, higher
frequencies (up to 300 Hz) are also present (e.g. Remillard et al. 2002; McCR03).
We do not offer yet any precise explanation to these phenomena. However, we note that QPOs are stronger in the hard X-ray bands (20-30 keV; e.g. Rodriguez et al. 2004) and are correlated with the radio flux (Migliari et al. 2005). Their emission must then be related to the dynamics involved in the ejection events (both steady and eruptive).
Interestingly, our framework provides a promising environment for the onset of instabilities leading to QPOs. The inner pair beam is an intermittent flow from the inner regions (several ri) and could therefore be responsible for some of the high frequency QPOs. On the other hand, the MHD jet may provide low frequency QPOs. For instance, if some disc material, continuously ejected just outside rJ, is failing to become super-Alfvénic, then one would expect waves going back and forth between the disc surface and the Alfvén surface (located at
and
in cylindrical coordinates) where a
shock is occurring. A crude estimate of the frequency gives
since
and
in magnetically driven jets (Ferreira 1997)
and
is the rotation rate of the magnetic surface. For a radius
,
this gives a
10 Hz QPO. This clearly deserves further investigation.
We present in this paper a new paradigm for the accretion-ejection properties of Galactic Black Hole X-ray binaries. We assume the existence of a large scale magnetic field of bipolar topology in the innermost disc regions. Such a field allows for several dynamical phenomena to occur whose relative importance determine the observed spectral state of the binary. The dynamical constituents are: (1) an outer standard accretion disc (SAD) for r> rJ; (2) an inner Jet Emitting Disc (JED) for r<rJ driving; (3) a self-collimated non-relativistic electron-proton surrounding, when adequate conditions are met; (4) a ultra-relativistic electron-positron beam. The dynamical properties of each constituent have been thoroughly analyzed in previous works (e.g. Shakura & Sunyaev 1973; Saugé & Henri 2004; Renaud & Henri 1998; Ferreira & Pelletier 1995; Marcowith et al. 1997; Saugé & Henri 2003; Henri & Pelletier 1991), but it is the first time where they are invoked altogether as necessary ingredients to reproduce the different spectral states of a same object.
We showed that the various canonical states can be qualitatively explained by varying
independently the transition radius rJ and the disc accretion rate .
In our view, the Quiescent and Hard states are dominated by non relativistic jet production from the JED, providing henceforth a persistent synchrotron jet emission. The Soft state is obtained when the transition radius rJ becomes smaller than the last marginally stable orbit ri, a SAD is established throughout the whole accretion disc. Intermediate states, between Hard and Soft, are expected to display quite intricate and variable spectral energy distributions. Luminous Intermediate states, obtained during the Hard-to-Soft transitions, are those providing the unique conditions for intermittent pair creation. These pairs give rise to a ultra relativistic beam propagating on the MHD jet axis, explaining both the observed superluminal motions and hard energy tail.
The dynamical structure presented here (JED, SAD, MHD jet and, occasionally, a pair beam) seems to be consistent with all available information about the canonical spectral states of BH XrBs. However, a more quantitative analysis is critical. In particular, we need to show that the base of the MHD jet can indeed provide a hot corona with the correct spectral signature. Then, a precise estimate of the radio/X-ray correlation predicted by our model and its comparison to observations will be a test of prime importance for its validity. This is postponed to a future work (Petrucci et al., in preparation).
In our view, the magnetic flux available at the inner disc regions is a
fundamental and unavoidable ingredient that most probably
varies from one system to another. Since changing the amount of magnetic
flux changes the transition radius rJ, the characteristic value of
(hence luminosity) associated with each spectral state is also
modified. Also, if accreting material is carrying magnetic flux of
opposite direction (Tagger et al. 2004), then this should lead to a major
readjustment of the whole magnetic structure. Clearly, taking into
account the advection of a large scale magnetic field within the disc
introduces a whole new set of variable phenomena.
Finally, we note that the typical values of the magnetic field required to steadily
launch jets from JEDs, given in Eq. (13), are consistent with observational estimates (Gliozzi et al. 1999; ; Gnedin & Natsvlishvili 1997). We expect that BH XrBs should radiate above the MeV range during luminous intermediate states, when pairs are produced.
A similar proposal has been developed recently by Bosch-Ramon & Paredes (2004b,a) but it is here a natural
outcome of our model. Very interestingly, such a high energy emission has been recently detected by the HESS instrument in the TeV range (Aharonian et al. 2005). Note also that two microquasars were already detected by EGRET, LS 5039 (Paredes et al. 2000) and LS I +61 303 (Massi et al. 2004, for a recent review), and there is possibly other unidentified galactic EGRET sources (Paredes 2004).
Besides, the -ray emission
of the pair beam can occur further away along the jet at the
-ray
photosphere, as proposed for AGNs (e.g. Marcowith et al. 1995). Noticeably,
-ray spectra of the possible EGRET counterparts seem to exhibit a
break in the MeV range, very similar to that observed in many AGNs. This
break could be explained by the transition from the optically thin X-ray
component to the optically thick, photosphere dominated,
-ray
component.
This prediction of
-ray emission of microquasars during very high
flaring states could be tested by future GLAST observations.
Acknowledgements
We thank S. Corbel for a careful reading of the manuscript and T. Belloni for having sent us a draft of its paper before it was completely accepted.