A&A 387, 918-930 (2002)
DOI: 10.1051/0004-6361:20020294
H. C. Spruit - B. Deufel
Max-Planck-Institut für Astrophysik, Postfach 1317, 85741 Garching, Germany
Received 31 August 2001 / Accepted 26 February 2002
Abstract
We show that the inner regions of a cool accretion disk in an X-ray binary can
transform into an advective, ion supported accretion flow (an optically thin ADAF,
here called ISAF),
through events involving only the known properties of the Coulomb interaction
in a two-temperature plasma, standard radiation processes, and viscous heating.
The optically thin inner edge of the disk is heated to a few 100 keV by the
strong flux of hot ions from the surrounding hot ISAF. We show that he
resident ions in this "warm'' disk are thermally unstable due to internal viscous
heating, and heat up to their virial temperature. The innermost disk
regions thus evaporate and feed the ISAF. These processes are
demonstrated with time dependent calculations of a two-temperature
plasma in vertical hydrostatic equilibrium, including heating by
external ions, internal proton-electron energy exchange, and viscous
heating. The process complements the "coronal'' evaporation mechanism
which operates at larger distances from the central object.
Key words: accretion, accretion disks - X-rays: binaries - black hole physics - radiations mechanisms: general, radiative transfer
Accreting galactic black holes (BHC) and active galactic nuclei (AGN) are often observed with two different spectral components: a soft component which is probably due to a multi-color blackbody from an optically thick, geometrically thin standard disk (SSD Shakura & Sunyaev 1973), and a hard component which is linked to an optically thin and geometrically thick flow. The hard component (an approximate power law with high energy cut off at keV) is most likely produced by inverse Compton scattering of soft photons on a hot thermal plasma (Sunyaev & Titarchuk 1980). Shapiro et al. 1976 showed that accretion can take place in the form of a two-temperature plasma, with the properties needed to produce such Comptonized radiation. Stable, optically thin, two-temperature flows were studied by Ichimaru (1977) and Rees et al. (1982). Extensive theoretical work on these accretion flows, stressing the role of the advection of internal energy was done by Narayan & Yi (1994, 1995a,b), for a review see Narayan et al. (1998), see also Esin et al. (1997) and references therein.
The effects of advection are the same for the optically thin two-temperature flows and optically thick, radiation supported flows, and both types of flow are now customarily called ADAF's. The distinction between these types, which was explicit in the older labels "ion supported'' and "radiation supported'', is still needed in many applications, however. For this reason, we propose here the name "ion supported accretion flow'' or ISAF as a means of distinguishing the optically thin type of ADAF from the radiation supported ("RSAF'') type.
Figure 1: Energy channels in a cool accretion disk embedded in a hot region. The energy dissipation in the cool disk is assumed small compared the atmosphere (a corona or an ion supported ADAF, here called ISAF). Squared boxes show physical processes, round boxes the particles involved. Heavy arrows and boxes show the main energy channel: viscous dissipation in the ISAF heats the protons there, which illuminate the cool disk below. By Coulomb interactions the ISAF protons loose their energy mainly to electrons, producing a warm (80 keV) layer which radiates this energy by Compton-upscattering of soft photons from the cool disk below. The protons in the warm surface layer are largely outside this main energy channel. Near the transition radius the cool part disappears and the warm part heats up to several hundred keV. Then the energy exchange of disk protons with the electrons is slow (shown by the saw tooth line) and viscous heating of protons becomes important. | |
Open with DEXTER |
ISAFs have attracted attention because of their potential to explain the spectra of X-ray transients (Esin et al. 1997 and references therein). The observations are consistent with an accretion flow that consists of two zones: an interior ISAF that extends from the black hole horizon to a transition radius , followed by an optically thick, geometrically thin and cool standard disk outside . A partial overlap between the two regions is probable since observations show evidence for the close vicinity of hot and cold matter in the central regions of BHCs and AGNs. This is indicated by a K iron fluorescence line at 6.4 keV and a Compton reflection component between 10-30 keV.
A critical element of such an accretion geometry is the change from the geometrically thin SSD to the hot ISAF flow at the transition radius. An alternative interpretation of the spectra is given by the "magnetic flare'' or disk corona model (Haardt & Maraschi 1991; Maraschi & Haardt 1997 and references therein; di Matteo et al. 1999; Merloni et al. 2000 and references therein), in which this transition is absent and the X-ray emitting plasma is a layer on top of a disk, heated by mechanical energy transfer from the disk.
Signatures of a disk transition radius have recently been found in the power density spectra of Cyg X-1 (Churazov et al. 2001; Gilfanov et al. 2000). From an analysis of the reflection component Di Salvo et al. (2001) showed that in Cyg X-1 the transition radius is located between and , if the observed reflected spectrum is due to a smeared component, or if the reflection is unsmeared, e.g. from the companion star or the outer disk.
We take these observations as reasonable indications that a transition from a cool optically thick to a hot optically thin accretion flow does in fact take place. But how the SSD-ISAF transition works is still under debate (cf. Manmoto et al. 2000 and references therein). Meyer et al. (2000) propose that the transition from the cold disk to the optically thin flow is due to a heat flow by electron conduction from a hot, friction-heated corona to the cold disk below (see also Meyer & Meyer-Hofmeister 1994). This model has its maximum evaporation efficiency at a large distance from the hole (a few 100 Schwarzschild radii). If, due to a high accretion rate in the cool disk, not all material is evaporated until that distance, the cool disk will survive until the last stable orbit. A transition radius further in than 100 Schwarzschild radii is inconsistent with this picture. Rózanska & Czerny (2000) investigate conductive and radiative coupling of an accretion powered corona with an underlying cool disk. For low accretion rates they find that the disk completely evaporates whereas high accretion rates prevent the SSD-ISAF transition as in Meyer et al. (2000). From a mathematical point of view Abramowicz et al. (1998) show that, if the transition region is not too wide, the region must rotate with super Keplerian orbital speed. Based on this property Kato & Manmoto (2000) demonstrate that trapped low-frequency oscillations are possible in the transition region.
In this paper we show how the inner disk regions, where the coronal evaporation process does not work, can evaporate into an ISAF. Only a few, well-known ingredients need to be invoked: the coupling of protons and electrons by the Coulomb interaction in a fully ionized plasma, standard radiation processes, and viscous heating.
The starting point is the view, supported by the observations mentioned, that an ISAF and a cool disk can coexist. That is, there is a partial overlap between the cold disk and the ISAF. In the overlap region, there is a very strong interaction between the two, since the ISAF consists of energetic ions (10-100 MeV) that penetrate the cool disk to a significant depth. The goal in the next sections is to determine the nature of the energy and mass exchange in this interaction region, and to show that it will lead to evaporation at the inner edge of the cool disk. The argument is then closed by determining the conditions under which these processes can consistently lead to the coexistence of the disk and the ISAF that was assumed at the outset.
Due to their low temperature, the conductive energy flux carried from the ISAF to the disk by the electrons is negligible. This is one of the reasons why the coronal evaporation mechanism that functions well at larger distances from the central mass fails for a two-temperature plasma. The energy flux to the cool disk is carried almost entirely by the ions. This ion energy flux can not take over the role played by electron conduction in the coronal evaporation process, for two reasons. One is that the mean free path of the ions is not negligible as it is in the case of electrons. The stopping length of 10-100 MeV ions penetrating into a cool disk corresponds to a Thomson optical depth of order unity. This prevents the development of the thin energy deposition layer that is needed to heat the plasma to high temperatures. More important is the fact that the energy carried by these ions primarily heats the electrons of the disk. These in turn radiate it very effectively by inverse Compton scattering and limit the plasma temperature in the interaction region to values of 100 keV, well below the virial temperature. Instead of evaporation, the loss of the ISAF ions to the cool disk is a very effective condensation process. This condensation is an important sink to an ISAF flow generated by evaporation at larger distances, and has to be overcome by a sufficiently powerful evaporation process at some location in the disk. This location will turn out to be a region near the inner edge of the disk.
The surface layer radiates its energy by Compton-upscattering of soft photons from the cool disk below. The energy of the ISAF is thus carried to the disk by the hot ISAF protons, and radiated away by the disk electrons. The protons in the warm surface layer are largely outside this main energy channel.
In the present analysis we concentrate on ion heating, which is necessarily very strong in the overlap region between an ion supported flow and a cool disk. A radiative flux from the ISAF can exist additionally, but in order to separate these contributions we assume here that the ISAF itself is radiatively inefficient, so that its direct contribution to the radiation from the accretion flow negligible. Finally, the viscous energy release in the warm layer turns out to be small compared with incident energy flux from the ISAF. Viscous heating is thus energetically unimportant in the surface layer, but it will turn out to be crucial to the evaporation process.
Finally there is the cool disk (1 keV) below the warm layer. This region is outside of the reach of the impinging hot protons. The cool disk serves as a thermalizer for the downward directed radiation of the warm layer. It is the source of the soft photons which are Comptonized in the warm region and keep it at moderate temperatures. The viscous energy release in the cool disk is assumed to be negligible as a source of soft photons.
In Paper II we have shown that close to the inner edge of an accretion disk, where the surface density (and the optical depth) of the disk gets small, the penetration of virialized protons heats the entire vertical disk structure to temperatures of several 100 keV (which is equivalent to a disappearance of the cold part in Fig. 1). This corresponds to the first step in our model and is at the same time the starting point of the present investigation. At high temperatures the time scale for establishing thermal equilibrium between the disk protons and electrons is not short compared to the dynamical time scale any more. The viscous energy channel (due to internal heating of protons by friction) gets important now because coupling to the electrons is weak.
It will turn out that the physics described depends almost only on the dimensionless radius from the hole , the dimensionless accretion rate , and the dimensionless viscosity parameter . The results are thus scaleable between AGN and BHC cases. Where explicit values of the physical parameters are needed, we take those of a typical BHC case.
In the next three subsections, we briefly review the physics associated with the penetration of the protons into the disk. This has been discussed before in detail in the references given in the previous subsection. Here, we address a few conceptual issues, such as the validity of the approximation that the energy of the incident protons is transferred mainly to the electrons in the disk, the accuracy of a nonrelativistic treatment, and the charge balance between ISAF and disk.
At the typical energies of the protons incident on the cool disk, the energy
loss is mostly by long-range Coulomb interactions with the electrons in the
disk (small-angle scattering on the large number of electrons in a Debye
sphere). This is opposite to the case of protons with a temperature near
that of the plasma in which they move. In the latter case, the equilibration
among the protons is faster than between electrons and protons, by a
factor of order
.
To see how this apparent
contradiction is resolved, consider the basic result for the energy loss of a
charged particle moving in a fully ionized, charge-neutral plasma. This was
derived by Spitzer (1962) (making use of Chandrasekhar's (1942)
earlier results on dynamical friction). Introduce as a measure of distance in the
plasma the Thomson optical depth ,
i.e.
,
where
is the
Thomson cross section, n the electron density and l the distance. The
rate of change of energy
of a proton moving with velocity
v in a field of particles with charge e and mass
(the "field
particles'' in Spitzer's nomenclature) is then given by Spitzer's Eqs. (5)-(15). In
our notation, this can be written in terms of the energy loss length
for interaction with the field particles f,
(2) |
(3) |
(4) |
For incoming protons of high energy, however, the result is different because
is not sufficiently small any more. In the high-v limit,
,
and one has
(5) |
(6) |
(7) |
Spitzer's treatment is non-relativistic, while virialized ISAF protons near the hole can reach sub-relativistic temperatures. A fully relativistic treatment of the Coulomb interactions in a plasma has been given by Stepney & Guilbert (1983). We have compared the classical treatment according to Spitzer's theory with this relativistic result in Deufel et al. (2001, 2002), and found it to be accurate to better than 5% for proton temperatures <100 MeV. The classical approximation in Spitzer's analysis therefore does not introduce a significant error for the problem considered here.
For high energies, the Coulomb energy loss becomes so small that loss by direct nuclear collisions becomes competitive. This happens (cf. Stepney & Guilbert 1983) at MeV, an energy that is not reached by virialized protons except in the tail of their distribution. We ignore these direct nuclear collisions. Note, however, that a gradual nuclear processing by such collisions can be important (Aharonian & Sunyaev 1984), in particular for the production of the Lithium. The Lithium overabundances seen in the companions of LMXB (Martín et al. 1994a), may in fact be a characteristic signature of the interaction of an ADAF and a disk described here (Martín et al. 1994b; Spruit 1997).
As the protons slow down, they eventually equilibrate with the field protons. This last part of the process is not accurately described by the energy loss formula (1). In addition to the simple energy loss of a particle moving on a straight path through the plasma, one has to take into account the random drift in direction and energy resulting from the interaction with the fluctuating electric field in the plasma. This drift can be ignored to first order (end of Sect. 5.2 in Spitzer 1962), but takes over in the final process of equilibration with the plasma. This last phase involves negligible energy transfer compared with the initial energy of the protons in the present calculations, and can be ignored here.
The protons penetrating into the disk imply a current that has to be balanced by a "return current''. As in all such situations, this return current results from the electric field that builds up due to the proton current. As this electric field develops, it drives a flow of electrons from the ADAF to the disk which maintains the charge balance. Since the electron density in the disk is high, the return current does not involve a high field strength.
Any mechanism heating the protons in the disk, even if only small, can potentially give rise to a thermal instability since the disk protons do not loose energy efficiently. Their energy loss by radiation is negligible on account of their high mass, and the transfer of energy to the disk electrons by Coulomb interaction is inefficient in a sufficiently hot plasma. In addition, by the nature of the Coulomb interaction, the time scale for the disk protons to equilibrate with the disk electrons increases with temperature (cf. Sect. 1.3). Depending on the strength of the mechanism heating the protons and its temperature dependence, a runaway may occur in which the disk protons continue to heat up as their cooling rate continues to decrease.
One heating process to be taken into account is internal viscous energy release (due to friction). At the low surface densities in the region of interest, the amount of energy released by friction is small compared with the heating flux by the external ISAF protons. If we may assume that a reasonable fraction of the viscous energy release goes into the protons, this small amount can still be quite important for the energy balance of the disk protons because of their low energy loss rate, and must therefore be taken into account. A second heating process is energy transfer from the incident ISAF protons to protons in the disk by Coulomb interaction. Though we have seen (Sect. 1.3) that the incident protons loose their energy mainly to the electrons, the small amount transferred to the disk protons may conceivably be relevant for their energy balance. In the following we investigate these heating mechanisms in turn, an conclude that viscous heating can give rise to a runaway, while heating of the disk protons by Coulomb interaction with the incident protons does not.
The energy balance of the protons in the warm disk is conveniently described in terms of heating and cooling time scales due to the processes involved. First we consider the cooling time scale by transfer of energy from protons to electrons in the disk.
The energy exchange timescale by Coulomb interactions between nonrelativistic
protons and electrons with Maxwellian distributions at temperatures
and ,
respectively, is given by Spitzer's (1962)
classical result:
The relativistic correction to this result is small for the conditions of interest here. In Deufel et al. (2002) we have shown that Spitzer's formalism is accurate to better than 5% compared to the relativistic treatment by Stepney & Guilbert (1983), for proton energies below 100 MeV and electron temperatures keV.
Next, the viscous heating timescale in an accretion disk is given by
Comparing Eqs. (8), (9) yields a
critical electron density
The amount of mass involved in the instability is the mass in the
atmosphere above the level where the density has dropped to
.
Let this mass be
(per unit surface area of the disk).
The time scale on which the protons heat up in the unstable region is
just the viscous heating time scale
.
The rate at which
the atmosphere evaporates is thus approximately
(14) |
Using Eq. (12) to replace
in Eq. (11)
this yields, with Eq. (10):
An evaporating part exists in every disk atmosphere, but at low temperatures only the highest layers evaporate. The mass loss is unimportant in such a case and can not change the properties of the disk. Equation (16) shows that the evaporation rate is a relatively strong function of R, and . Whether situations exist where this evaporation rate becomes important is investigated in Sect. 3.
In this subsection we show that Coulomb interaction of the incident protons with the resident disk protons is relatively unimportant for the conditions encountered. It is included in the numerical results reported in the next section, however.
The main energy channel is from the hot protons to the electrons in the warm surface layer (see Sect. 1.3 and Paper II). The properties of the Coulomb interactions at the temperature of the incoming protons are such that only a small fraction of their energy is transferred to the protons in the layer (which we will call "field protons'' here following the terminology in Spitzer 1962). The energy budget of the field protons is also small, however (cf. Fig. 1), so the effect of this channel on the proton temperature in the warm layer has to be considered.
The interaction of a fast proton with the protons in a much cooler plasma is given by Spitzer's result discussed in Sect. 1. In terms of the rate of change of the energy of the incoming proton, it can be expressed as
The heating rate
per unit volume due to the interaction
of hot protons with the field protons is then
We conclude that heating of the ions in the warm disk by the incident ions alone is not sufficient to drive the ion temperature away from equilibrium with the electrons. For this to happen, viscous heating has to be included.
We test the inferred heating instability numerically, starting with a warm disk of moderate optical depth and investigating its temporal evolution due to internal viscous heating, and p'-p and p-e energy exchange as described in the previous subsection.
We use a plane-parallel, one-dimensional model geometry. The vertical density distribution through the atmosphere is found from the equation of hydrostatic equilibrium as in Paper II, except that here we do not account for the pressure exerted by the decelerating protons.
The temporal evolution of the temperature profile is computed with a fourth order Runge-Kutta method. The timestep of each integration is set to the shortest timescale of the heating and cooling processes involved. The change of temperature per timestep within a certain volume is due to the rate of change of the enthalpy there, .
The temperature increase per timestep as a function of optical depth
due to viscous heating of the protons is then given by
The temperature change of the field protons due to the impinging hot
protons is computed similar to the proton illumination method
described in detail in Papers I and II. Thus we follow protons from a
Maxwellian distribution with virial temperature through the warm atmosphere
of the disk and record their energy losses. The energy flux in the illuminating
protons is related to the mass condensation rate
(see Eq. (26) and Sect. 3 for details).
The kinetic energy decrease of the fast protons (p') due to their
interactions with the field protons (p) as a function of
optical depth can be derived from Eq. (17). This gives the
heating rate
and the contribution of p'p exchange
to the temperature change per time step is
Figure 2: Evolution of the temperature profile due to the viscous instability of the protons. The model layer has optical depth and initial temperature keV (dashed line). Dashed-dotted lines show the profile after 10 thermal time scales, solid lines after t=20. Values of the viscosity parameter are 0.1, 0.3 and 0.5 (numbers at the lines). For comparison the temperature profile after due to interactions alone is also shown (dotted line). | |
Open with DEXTER |
When the disk protons in our simulation have reached their local virial temperature, we do not allow a further temperature increase. In a more realistic calculation, this limit on the temperature would come about through the advection of internal energy. Such a calculation requires a more detailed model of the ISAF and its sources and sinks of mass and energy, which is beyond the scope of the present calculation.
As an example we consider a warm disk with an optical
depth (from the surface to the midplane)
,
an
initially isothermal temperature profile
keV,
illuminated by an ISAF with accretion rate
,
where
is the Eddington accretion rate.
Figure 3: Local evaporation and condensation rates from a warm (100-500 keV) layer of optical depth , formed by proton illumination of an initially cool disk. Solid line shows the condensation rate onto the layer from an ISAF accreting at , as a function of distance from the hole. Dashed lines show the evaporation rates from the layer, for three values of its temperature. In the innermost region, the net effect is evaporation, feeding the ISAF, further out the ISAF condenses onto the layer. At lower ISAF accretion rates ( , dotted) the evaporating region becomes wider (see text for details). | |
Open with DEXTER |
Figure 2 shows the temperature profiles for various values of the viscosity parameter after (dashed-dotted lines) and (solid lines), where is the thermal time scale of the disk. For all disk protons heat up to the virial temperature, whereas for and only a part of the disk is subject to the instability. The extent of the unstable part is in excellent agreement with our estimate in Sect. 2.1. In test calculations in which only heating and energy coupling are included, the temperature increase of the disk protons is considerably smaller than with the viscous heating included. This confirms our conclusion from Sect. 2.2 that interactions alone can not heat the bulk of the internal protons considerably above the ambient electron temperature (cf. Eq. (20)).
Summarizing our investigation of the energy channels in a warm disk we conclude that viscous energy release in a hot plasma causes a runaway temperature increase in the unstable upper part of the disk atmosphere, until the protons have reached their local virial temperature. The process affects a significant part of the stratification if the electron temperature is above 100 keV.
In Sect. 2 we have shown that the internal viscous heating of the disk protons leads to a mass evaporation rate according to Eq. (16).
At this point we do not yet know whether the mass loss from a warm disk region is high enough to completely evaporate the disk. At the same time as the upper atmosphere of the warm disk evaporates, the hot protons from ISAF condense into it and increase the surface density. For an effective evaporation of the disk the mass loss rate must be higher than the condensation rate. To compute the condensation rate we need an estimate of the density in the ISAF.
In our previous numerical simulations of warm disks (Paper II) we parameterized proton mass flux from the ISAF by scaling the energy flux of the incident protons with the local energy dissipation rate in the ISAF. Here we adopt a slightly more realistic mass flux rate in an ISAF.
In a thin disk approximation the surface density
of the ISAF,
with accretion rate ,
is
Figure 4: Decrease of the proton mass flux with depth in a thin, hydrostatic, isothermal disk with total optical depth , at disk temperatures of 1, 50 and 200 keV, for and . At temperatures above 100 keV, only a fraction of the incident protons is captured by the disk. | |
Open with DEXTER |
Figure 3 shows this comparison for different values for the viscosity parameter , different temperatures and accretion rates. For values of the viscosity parameter and evaporation dominates when the temperature of the warm state exceeds 300 keV. As the ISAF density (i.e. the accretion rate) decreases, the condensation rate of the protons into the the disk decreases and evaporation dominates over the condensation rates over a wider range of radii.
The condensation rate given by Eq. (26) is actually an overestimate, since it assumes that all incident protons are stopped in the disk. While this is correct for cool disks, for a high temperature plasma the rate of the electron-proton energy exchange is small and a disk with low surface density gets optically thin for the penetrating hot protons. This is demonstrated by Fig. (5) where we show how the incident proton flux changes with depth into a warm disk. The Thomson optical depth in this example is ( ). At the temperatures of a cool standard disk (1 keV) almost all protons are absorbed in this layer. But at the high temperatures of the warm state the disk is optically thin and practically all ISAF protons fly through the disk without being absorbed. The penetrating protons do not add to the surface density in this case, and evaporation should therefore be possible even at lower values of .
At temperatures below 100 keV, the evaporation rate from the warm state into the ISAF is quite small and can not balance the loss from it by condensation. This is roughly the temperature of the warm surface layer on a cool optically thick disk heated by proton illumination (the situation sketched in Fig. 1). Thus the relative importance of evaporation and condensation reverses just at the point where the cool component disappears. As long as a cool disk is present, the thermal instability in its warm surface layer is relatively weak, while it effectively absorbs all incoming protons. Once the cool component is gone, the temperature and evaporation rate goes up, while at the same time mass condensation by stopping of protons in the disk becomes ineffective.
Figures 3 and 4 show that net evaporation takes place close to the hole preferentially at low accretion rates in the ISAF, which may sound counterintuitive. It is a consequence of the fact that the evaporation rate does not depend on the flux of incident protons (cf. Eq. (16)). This again is a result of the fact that the field protons in the warm disk lie "outside the main energy channel'': the incident proton energy flux sets the electron temperature (through the Comptonization balance) but does not affect the proton temperature directly. The temperature of the field protons is determined by the secondary balance between viscous heating of the field protons and their energy loss to the electrons.
At low accretion rates, however, the hot proton flux eventually becomes insufficient to keep the layer "warm'': it will cool down to low (1 keV) temperatures by bremsstrahlung losses. We have studied this transition in Paper II, where we found that a layer of optical depth can just be kept in the warm state for an accretion rate in the ISAF of about .
At lower accretion rates, a warm disk is possible only when the bremsstrahlung losses are lower, at lower electron densities. Since the electron density is proportional to the optical depth of the disk, and the brems losses proportional to , the minimum ISAF accretion rate needed to maintain a warm disk state by proton illumination scales as . The optical depth of the disk vanishes towards its inner edge, so we expect that there is always a region close to the inner edge of the disk where evaporation takes place, even at very low accretion rates.
Summarizing this argument, the optical depth
of a warm disk
at the point where it matches onto the cool disk, depends on the accretion rate
in the ISAF as
In the above we have started with the assumption that an ISAF and a cool disk coexist, and found the conditions under which the disk can feed mass into the ISAF. We have done this by considering the conditions at each distance from the hole separately. To turn the ingredients into a consistent picture, we have to consider the mass flux through the system, so that conditions as a function of distances from the hole are connected, in the way sketched in Fig. 5. At large distance, we have a standard cool disk. Closer in, an ISAF surrounds it and condenses onto it, producing what we have called here the warm, proton-illuminated layer. Even further in, the vertical optical depth of the disk becomes too low to sustain cooling by brems losses, and the whole disk transforms to a "warm disk'' state (Paper II) of nearly uniform temperature. The upper layers of this warm disk evaporate to feed the ISAF assumed at the outset, and at the inner edge of the disk, at some radius , the entire mass flux through the disk has evaporated into the ISAF. We now investigate what the conditions are for such a radial structure to be possible in a steady state.
The disk, whether in a cool (1 keV) or warm (300 keV) state, is still very cold compared to the local virial temperature, so that the standard thin disk approximation is valid. The difference with a standard steady disk is that the mass flux is now a function of distance, due to the condensation and evaporation from/to the ISAF. We first consider the modifications to the -disk diffusion equation that this causes.
Figure 5: Evaporation from the warm disk: a cool disk (black, 1 keV) partly extends into an ISAF and is exposed to the hot virialized protons from the hot torus (open arrows, light grey). Above and below the cool disk a heated surface layer is produced (80 keV, dark grey) due to proton illumination. At the inner edge at R0 the surface density of the cool disk gets small and a warm disk develops (300-500 keV, grey). Net mass loss into the ISAF evaporates the warm region between (filled arrows). Due to the high temperatures in the warm disk this region is a source for hard X-rays. | |
Open with DEXTER |
The surface mass density of the disk material is
The change of the mass accretion rate with radius through the cool
disk due to evaporation can be expressed by
As in the standard derivation, the thin disk diffusion equation follows from the angular momentum equation, which now includes a term for the angular momentum carried with the evaporating/condensing material. For the present purpose, it is sufficient to assume that condensing and evaporating material just has the same specific angular momentum, as the disk.
The equation for the angular momentum balance in an evaporating disk
then is
We use the usual
prescription for the viscosity
(Shakura & Sunyaev 1973),
Multiplying the continuity equation (Eq. (31)) with
yields after subtraction from Eq. (32) an
expression for the mass flux
in the cool disk in the
stationary case (
):
Thus we set
at the inner edge of the disk, as in
standard accretion theory. With
Eqs. (30), (37), (38) and
we get
an expression for the second integration constant C2,
This expression for the surface density, though strictly derived for
steady conditions, is still approximately valid if the position of the
inner edge changes slowly. We are interested in a true stationary
case, however. In this case the accretion rate in the disk will be
exactly equal to the total mass loss due to evaporation in the disk,
,
or in other words all matter
drifting inward through the cool disk has evaporated when
is reached. Eq. (36) then simplifies to
We can now estimate the distance over which the process of evaporation
into the ISAF takes place. Let R0 be the innermost radius where the
cool disk component exists. Inside this, there is only a warm disk
(see Fig. 5). Evaporation takes place both from the
warm layer on top of the cool disk at R>R0 and from the warm disk
region
,
but the warm disk region is expected to
contribute most, since its temperature is significantly higher. Thus
we equate, for the present approximate purpose, the evaporating region
with the warm disk region. Assume that the relative extent
of the warm disk is small. The evaporation rate
Eq. (16), with
,
depends only
on temperature. The temperature of the warm disk is relatively fixed
(Paper II), so we can set
constant as well.
Equation (40) for the surface density as a function of
distance from the inner edge is then
We can now make an estimate of the relative extent
of the warm
disk. R0 is the maximum radius where the warm disk can
exist (at larger surface density, it develops a cool disk component).
The optical depth at R0 is thus given by Eq. (27). Computing
the optical depth from Eq. (43) and equating this to
we have
Our model describes a situation which naturally develops if a cold standard Shakura-Sunyaev disk is truncated within a hot, optically thin flow (ISAF) in the inner regions of the accreting system. We have proposed a SSD-ISAF transition based on few, well-known physical processes: Spitzer's theory of the energy exchange in a fully ionized plasma, and standard viscous heating due to friction in the accretion disk. We have shown that the transition of the cool disk material into an ISAF is the logical and inevitable consequence of these basic interactions. The process involves two steps:
(i) At the inner edge of the disk the surface density of the cool disk gets low. Virial protons penetrating from the ISAF heat the cool disk electrons. The electron temperature is limited ( keV) since they can radiate their energy efficiently via bremsstrahlung and Comptonization. Once the disk is too thin, proton heating overcomes the radiative losses everywhere in the disk. The disk heats up, expands and radiative losses become even more inefficient. With increasing temperature pair production sets in, but also the proton heating gets less efficient at higher temperature. Finally a new equilibrium for the thin disk is found at several 100 keV. We label this the "warm state'', since its temperature is intermediate between that of the cool disk and the virial temperature. An important aspect in this process is that protons in the disk are outside of the main energy channel. The main energy is transferred from the external ISAF protons to the disk electrons, which loose this energy via radiation. For the formation of the warm disk state the internal viscous heating of the protons is completely unimportant, but not so in the second step of the process:
(ii) At the temperatures of the warm state the disk protons and electrons are not coupled very tightly any more, as the timescale for establishing thermal equilibrium is not short compared to the thermal timescale. Now the minor energy channel due to viscous heating in the warm disk becomes important for the energy budget of the protons, because the viscously released energy can not be exchanged very efficiently with the ambient electrons. The upper part of the warm disk, where the densities are lowest and the Coulomb exchange time scales longest, is subject to a thermal instability. The size of the unstable region depends mainly on the viscosity parameter and the temperature of the warm disk. The higher the temperatures and , the deeper is the unstable region. The protons there are heated to their local virial temperature, and become part of the ISAF: the warm state evaporates. In the warm state no effective mass condensation takes place, since it at the same time becomes optically thin for penetrating hot protons. The mass evaporation due to the viscous instability of the warm state is therefore the main process which determines the mass budget in this region.
On the basis of this picture, we have also looked at the radial structure of an accretion disk in which mass exchange with an ISAF takes place. Using the evaporation and condensation rates derived, we have investigated the conditions under which accretion is possible in such a way that the entire accretion flow eventually evaporates. In this case a steady state is possible with an inner edge to the disk at some finite radius outside the last stable orbit. We find that such a steady state is indeed possible for plausible values for the accretion rate and viscosity parameter. The steady state condition determines a relation between the accretion rate and the value of . Also, it determines the width of the warm, evaporating disk region; we find for of order 0.3 and accretion at a tenth of Eddington.
A potentially important factor which we have not been able to include in our picture of the warm disk region is the cooling effect of soft photons from the cool disk extending just outside the warm disk. If such photons can enter the warm region, they will cause a cooling by inverse Compton scattering on the warm electrons. Since both the cool and warm disk are quite thin (), the radial optical depth of the warm disk is large, and the angle subtended by the cool disk as seen from the interior of the warm disk is small (to visualize this, see Fig. 5). It thus seems likely that the effect of such cooling on the energy balance in the warm disk can be small, but this point requires closer scrutiny.
The extent of the region (in distance from the hole) where the physical conditions assumed here apply is limited. At large radii the proton temperature in the optically thin, hot region decreases (R-1). At lower proton temperatures the proton penetration depth into the cool disk also gets smaller. This limits the region where a warm state can be produced. Without the warm state our mechanism will probably not work efficiently enough to transfer all material from a cool SSD into material from an ISAF. Therefore we do not expect our mechanism to work at large radii from the black hole. But a combination of the coronal evaporation flow model as suggested by Meyer et al. (2000) and the two stage model proposed here could in principle cover a large range in the radial direction for the SSD-ISAF transition to occur.
Acknowledgements
This work was done in the research network "Accretion onto black holes, compact stars and proto stars'' funded by the European Commission under contract number ERBFMRX-CT98-0195.