A&A 421, 659-666 (2004)
DOI: 10.1051/0004-6361:20035955
B. F. Liu - F. Meyer - E. Meyer-Hofmeister
Max-Planck-Institut für Astrophysik, Karl- Schwarzschildstr. 1, 85740 Garching, Germany
Received 30 December 2003 / Accepted 25 March 2004
Abstract
We study the possibility of a cool disk existing in the Galactic
Center in the framework of the disk-corona evaporation/condensation model.
Assuming an inactive disk near the gravitational capture
distance left over from an earlier evolutionary stage,
a hot corona should form above the disk since there is a continuous
supply of hot gas from stellar winds of the close-by massive stars.
We study the interaction between the disk and the corona. Whether the
cool disk can survive depends on the mass exchange between disk and
corona which is determined by the energy and pressure balance.
If evaporation is the dominant process and the rate is larger than the
Bondi accretion rate in the Galactic Center, the disk will be depleted
within a certain time and no persistent disk will exist. On the other
hand, if the interaction results in hot gas steadily condensing into
the disk, an inactive cool disk with little gas accreting towards the
central black hole might survive in the Galactic Center. For this case
we further investigate the Bremsstrahlung radiation from the hot corona and
compare it with the observed X-ray luminosity. Our model
shows that, for standard viscosity in the corona (
), the mass
evaporation rate is much higher than the Bondi accretion rate and the
coronal density is much larger than that inferred from Chandra observations. An
inactive disk can not survive such strong evaporation. For
small viscosity (
)
we find condensation solutions.
But detailed coronal structure computations show that in this case
there is too much X-ray radiation from the corona to be in
agreement with the observations. From this modeling we conclude that
there should be no thin/inactive disk presently in the Galactic
Center. However we do not exclude that the alternative non-radiative
model of Nayakshin (2004) might instead be realized in
nature and shortly discuss this question.
Key words: accretion, accretion disks - black hole physics - Galaxy: center
The low luminosity was explained by an advection-dominated accretion flow (ADAF), with a spectral fit first presented by Narayan et al. (1995). In the following years important observational results on the emission of Sgr A* and theoretical work lead to an improved model for radiatively inefficient accretion flows (RIAFs). RIAF models for Sgr A* are discussed in the recent reviews of Yuan et al. (2003) and Quataert (2003). According to these models, most of the thermal energy released by viscosity and increased by compression is retained in the gas and advected to the central black hole. The RIAF model in addition assumes that very little mass from large radii actually accretes onto the black hole while a large part is lost through outflows during the accretion. These models naturally yield the observed spectra of Sgr A*. A key constraint on these models is that the fraction of gravitational energy heating the electrons must be very small and hence a two-temperature treatment of the plasma is required.
Another possibility to explain the low luminosity of Sgr A* might come from the existence of cold molecular gas in the parsec region of the Galactic Center in the form of an inactive/dead thin disk without accretion. Falcke & Melia (1997) had suggested such a "fossil'' disk (for a review see Melia & Falcke 2001). The model assumes that gas captured at Bondi accretion rate condenses onto such an inactive disk without mass accretion onto the black hole. This runs into difficulties because the inflowing gas produces a fair amount of luminosity in the infrared (Narayan 2002) as it clashes onto the disk and loses its thermal and kinetic energy. This radiation is not seen. Nayakshin (2004) revisits the concept of an inactive disk. He suggests that, in the case of an extraordinarily long mean free path (larger than the pressure scale of the corona) and extremely low viscosity of the hot gas, the energy can be conducted to a very thin transition layer by free streaming electrons and is then radiated in infrared to UV wavelengths. As this radiation from the thin layer would be observed edge-on a discrepancy between predicted and observed luminosities could be avoided. This interesting suggestion deserves further analysis of the transition between the hot gas and the cool layer and the coupling between the hot ions and the energy transferring electrons.
As a further contribution to the issue of a cool disk around the Galactic Center we here study the vertical structure of a hot corona in interaction with a disk below. We assume a cool disk in the outer region around the circularization radius where the free fall of a Bondi type accreting hot gas goes over into a circular motion around the gravitational center due to its specific angular momentum. Such a disk might not be unreasonable since the system must have been quite bright during an earlier evolutionary stage and angular momentum that was released by the high accretion rate should have moved disk mass into outward regions. With winds from young stars being captured by the gravitational field of the black hole, a corona unavoidably forms above the cool disk. The question is, how do disk and corona evolve? Does the hot gas condense to the disk with little mass actually accreting to the black hole, or does mass rather evaporate from the disk to the corona overwhelming the incoming hot gas and even finally depleting the cool disk underneath?
The answer depends mainly on the rate at which the hot gas is supplied from the capture radius. If gas is supplied to the corona at a sufficiently high rate coronal gas condenses to the cool disk. If no gas is supplied from the outside or if the outside gas supply is too small mass instead evaporates from the disk into the corona. Both processes are the consequence of pressure and energy equilibrium between the disk and the corona (Meyer et al. 2000; Liu et al. 2002). Here we study in detail the structure of the corona for the case of the Galactic Center in order to see what the dominant process between disk and corona is, condensation or evaporation. Can the cool disk survive if evaporation is dominant? Furthermore, the detailed computation allows to calculate the Bremsstrahlung luminosity of the corona and compare it with the observed X-ray luminosity. We show that these results exclude the existence of any cool disk in our Galactic Center.
In Sect. 2 we describe the physics of the interaction between the disk and the corona. In particular we discuss how a radial inflow of mass from the outside affects the mass and energy balance in such a corona. In Sect. 3 we present numerical results and show how the value of the viscosity affects evaporation or condensation. We discuss several aspects of our results in Sect. 4, including a comparison with a non-radiative condensation model of Nayakshin (2004) and the question of a past disk being evaporated now in the Galactic Center. A conclusion follows in Sect. 5.
For a hot corona lying above a cool disk, interaction between the disk and the corona occurs via energy and mass exchange. The hot corona conducts heat downward by electrons. At the bottom the heat is radiated away. If the density in the corona is too low, Bremsstrahlung is inefficient and the thermal conductive flux heats up some of the disk gas leading to mass evaporation from the disk into the corona. The resulting density increase in the corona raises the radiation loss and thereby counteracts further evaporation. If the coronal density is too high, radiative cooling is too strong and gas condenses into the disk. At the final equilibrium density, cold gas steadily evaporates from the disk into the corona if mass is drained continuously from the corona inward by diffuse flow, or hot gas steadily condenses to the disk if the corona continuously gains mass by mass flow. For example, when there is no hot gas coming in through the outer boundary, (case 1), mass is continuously lost from the corona by accretion towards the central object. This is resupplied by evaporation from the surface of the cool disk as the corona tries to restore the density to the equilibrium level. If there is more hot gas being fed in at the outer boundary than what flows inward towards the center, (case 2), hot gas continuously condenses to the cool disk.
We now discuss the upper boundary condition for the corona at such radii. The earlier investigations showed that in general wind loss from the corona is an integral part of the solution, where a sonic transition occurs at some height and the wind flow cross section flares out with the wind expansion. In an advanced multi-zone modeling (Meyer-Hofmeister & Meyer 2003) we found that the wind pressure is highest at the distance where the evaporation efficiency is highest, at a few hundred Schwarzschild radii. The pressure in the expanding wind from this region even dominates over the pressure at a sonic point of winds from farther out regions and prevents sonic transition and wind loss there altogether. Thus for evaporation solutions in these outer regions we apply the condition of zero vertical velocity at a height of z=R. The actual height at this point is not important as long as it includes the lower down region where most of the coronal action occurs.
In the case of condensation solutions the incoming mass flow from the outside accretion anyhow dominates the coronal pressure and prevents any free wind expansion. We thus can apply the same boundary condition. In our equations we can also leave out the so-called flaring terms which are only important in the wind expansion geometry.
In the following we list the four ordinary differential equations describing the coronal flows above a disk in the Galactic Center.
Continuity of mass flow
Here
are density, pressure and temperature, vR and vzthe radial and vertical velocity, M is the black hole mass, G the gravitational constant and
the rotational frequency,
and
are electron and ion particle densities,
is the
Bremsstrahlung cooling rate, and
the ratio of specific heats.
is the viscosity parameter (ratio of viscous stress to
pressure).
The terms
and
account for radial mass and energy
flows with
,
explained in the next section.
The difference between
and
results from the fact that
the specific energy which the mass flow carries scales radially as
.
We consider stationary azimuthally symmetric flows.
Compton cooling is negligible for the coronal
structure at distances of 104 to 105 Schwarzschild radii.
Three possible contributions to Compton cooling have to be considered.
(1) The cooling by the radiation from the disk surface caused by mass flow
in the disk was investigated by Liu et al. (2002) and shown to be negligible
even for mass flow rates in the disk of
.
(2) We have estimated that the energy loss by Compton cooling from
radiation caused by reprocessing of coronal X-rays is always less than
of that by bremsstrahlung. (3) The Compton effect from the X-rays
from the central source and the surrounding gas is negligible
in our context because of the very low observed radiation.
At the lower boundary z0 we start our calculations at the
temperature
T=106.5 K.
This value is in the steep temperature profile in a thin transition
zone of nearly constant pressure. Its physics can be described by the
balance between gain of heat by thermal conduction and radiation loss.
Other effects like frictional heating and energy transport by the vertical
mass flow are negligible here (Meyer et al. 2000). This establishes a
relation between temperature and heat flux which can be scaled according
to the pressure (Smeleva & Syrovatskii 1973). If one temperature in the
profile is selected then, due to the scaling, one obtains a unique
relation between thermal heat flux and pressure (Liu et al. 1995).
We use this relation
![]() |
(5) |
![]() |
(6) |
Integration of Eq. (1) along the z-direction gives,
In our Galactic Center, hot gas in the central parsec produced by winds from massive stars is gravitationally captured by the black hole at a distance around 0.04 parsec which corresponds to 1 arcsec at the sky. If there is a cool (inactive) disk in this region, a corona forms above the disk, fed by the captured hot gas. The situation is a bit different from the one discussed in our previous evaporation model, where we assumed that cold gas flows inward from the outer region (e.g. a secondary star) via the thin disk and then feeds the corona through evaporation. There it was assumed that essentially no significant hot gas comes from a further outward located hot corona.
Now the situation is turned around, hot gas accretes through the corona without any cold gas being fed into the cool disk from the outer disk boundary. As explained in the last section, if much mass is fed into the corona, the pressure in the corona is high. The pressure and energy balance between disk and corona leads to mass condensing from the corona to the disk. On the other hand, if the feeding rate to the corona is low, mass evaporation from the disk to the corona is the dominant process. Strictly speaking, how large the mass exchange rate is and whether mass evaporates or condenses depend on the rate of net mass inflow into the coronal region, i.e., the difference between the mass feeding rate at the outer boundary and the mass loss rate to the central object at the inner boundary.
If the equilibrium between disk and corona leads to strong mass evaporation, the cool disk will eventually be depleted and mass then accretes via a hot corona/ADAF or RIAF. If condensation dominates, the hot gas captured by the gravitational field of the black hole finally flows to the disk and is deposited there with little mass actually accreting onto the black hole. In the following we investigate these coronal features for the case of our Galactic Center.
![]() |
Figure 1:
Radial mass flow rate in the corona for
given parameters
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In our terminology any value of
means there is more gas
flowing out of the coronal region than coming in. The net mass flow
out of the corona is contributed by mass evaporation from the disk.
The case of "standard evaporation'', considered in the
one-zone-model corresponds to
,
which means that the
typical mass flow in the corona is the same as the mass
evaporation into the corona at that distance.
Figure 1 shows the radial mass flow rate in the corona
(on both sides of the disk) for distances
to
.
The mass flow rate is
,
with
the typical radial drift velocity
(
isothermal sound
speed,
Kepler angular frequency, z0 and z1lower and upper boundary of the corona respectively).
The lower solid line represents the mass accretion rate
for the standard one-zone evaporation model.
We see that the mass flow rates in the coronal region from
to
are around
,
about 10 times higher than the
incoming mass flow estimated from Bondi accretion and Chandra
observations as
.
This indicates that the dominant process in the disk-corona system is that
gas evaporates from the disk to the corona and then flows towards the
central black hole. The hot gas captured at the Bondi radius is only
a minor contribution and hardly affects the coronal structure.
A strong corona with high density and high temperature
is built up above the disk by the mass evaporation.
![]() |
Figure 2:
Vertical structure of the corona with mass
evaporation from the disk to the corona (
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Figure 2
shows the coronal structure at
for both cases, evaporation and condensation. From the figure
we see that in the case of evaporation the typical temperature of the
coronal gas is
K (virial temperature
,
gas constant,
molecular
weight, taken as 0.62). The particle number density which follows
from pressure and temperature is larger than
.
Temperature and density both are much larger than
the values observed by Chandra at the Galactic Center.
A disk might have existed earlier but it could have been
eventually depleted by the mass evaporation since there is no mass
supply for the cool disk. Therefore, the
standard disk-corona evaporation model appears to exclude the continued
existence of a cool disk in the Galactic Center. Note that this
argument is independent on whether the cool disk is completely inert or
whether it allows some mass flow itself.
The condensation rate is far too high compared to the Bondi accretion rate in the Galactic Center. The density in such a corona is also much higher than the value inferred from the observations. Therefore, we cannot expect that, if there is a cool disk in the Galactic Center, hot gas captured by the black hole at the Bondi radius mainly condenses into the disk and is deposited there. Instead, the disk gas will evaporate into the corona, increasing the coronal flow inward significantly. As noted above, this process can finally deplete the cool disk. From then on no cool disk exists anymore in the Galactic Center.
![]() |
Figure 3:
Radial
mass flow rate in the corona ( upper curve) together with
condensation
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
More generally, we study the condensation and
evaporation solutions for a sequence of values of
.
Figure 3 shows how the mass flow rate and the
evaporation/condensation rate at
change with the
value of the parameter
.
The results show that condensation
solutions are connected with a high mass flow rate through the
coronal region. At
all mass that flows in the corona has
come from the outside and condenses into the cool disk. With increasing
a smaller and smaller part of the coronal mass flow settles
into the disk until at
no mass condenses and all mass that
comes from the outside continues inward. For increasing positive
mass evaporation contributes a growing fraction to the mass flow in
the corona until at
all mass flowing inward in the corona
comes from evaporation of the cool disk. (The difference between upper
and lower curve at
and 1 is
an artifact of the one-zone model approximation and results from the
different way in which averages over the one-zone area are defined.)
One may note that as the characteristic mass flow rate is proportional to the
pressure in the corona,
,
the results in Fig. 3
reflect the intrinsic relation: more pressure means more radiative
cooling which supports condensation, less pressure means less
radiative cooling and supports evaporation.
Thus, if we know the pressure (or mass content) in a coronal region,
our computations for the equilibrium between disk and corona allow
to derive how much mass condenses or evaporates. This is not directly
determined by the incoming mass flow alone, but depends on the net mass flow
into the region, i.e.,
the difference of incoming mass flow at the outer boundary and the
outgoing mass flow at the inner boundary. Therefore, when we
know both the inner and outer boundary conditions as in multi-zone
modeling (Meyer-Hofmeister & Meyer 2003), we
are able to determine the vertical and the radial structure of the
corona. Otherwise,
will be an open parameter in the range
of -1 and 1.
Combining Figs. 1 and 3, we find that, for
,
the mass flow rate in the corona is between the
two solid curves of Fig. 1. The same is true for all
distances of
.
In other words,
no matter how large a fraction of the coronal mass flow
condenses to the disk (
between -1 and 0),
or how large a fraction of the coronal mass flow originates from
disk evaporation (
between 0 and 1),
the mass flow in the corona would be much larger
than the Bondi accretion rate of
inferred from
observations. The detailed calculations also show that the hot gas
density is much larger than that obtained from the observations.
As a consequence of such a large mass flow rate, the luminosity of
coronal radiation in X-rays (note that it is not in the infrared),
estimated from
,
is very high. For instance, at
,
,
the
radiation from the corona would then be
,
much larger than the observed quiescent X-ray luminosity
.
The Bremsstrahlung radiation is
confirmed by detailed computations of the coronal structure.
Therefore, for standard ,
these consistent disk-corona model
calculations exclude a disk at distances
from the Galactic Center by
comparison of mass flow rate, density, and luminosity of the hot
gas predicted by theory with those actually observed.
In the last section we showed that for standard
no
condensation solution exists that is compatible with the observed
mass accretion rate in the Galactic Center.
We investigate whether agreement between
observation and a disk-corona analysis can be achieved for smaller
values of
.
In the models the viscosity parameter enters directly into the radial
drift velocity and the viscous release of heat, both of which become
smaller with smaller .
Since radiative cooling is not affected
by this change in
the balance between heating and cooling
requires less pressure in the corona. As a consequence one expects
that pressure, temperature, and radial mass flow in the corona
decrease when
becomes smaller (Meyer-Hofmeister & Meyer 2001).
This effect is shown in Fig. 1 where the dashed lines give
the mass flow rates in the corona for a smaller viscosity
.
The upper and lower dashed lines show the condensation
(
)
and the evaporation solutions (
)
for
.
For other values of
,
,
the mass flow rates lie between these two lines. Indeed, a small
viscosity in the corona results in a significantly decreased mass flow
rate. Density and pressure in the corona also are smaller.
For very small
,
the mass flow rate in the corona can then, in
principle, become comparable with the capture rate
of
.
If
is the Bondi accretion rate in the
Galactic Center,
,
consistency with a condensation solution (
)
can
only be obtained by assuming a small viscosity parameter
(for example,
for
at
and an even smaller value at smaller distances).
However, in either of the extreme cases,
or
,
or any case between
those,
,
our detailed computations of the coronal
structure show that the X-ray luminosity produced by Bremsstrahlung
radiation in the corona by far exceeds the luminosity observed. The
existence of a cool disk in the Galactic Center seems thus to be
excluded. For a condensation picture this was already pointed out
by Narayan (2002).
In addition to this result for the situation in the Galactic Center with the presently observed mass accretion rate the results of our investigation allow to discuss whether a formerly existing accretion disk in the context of gravitational instability and star formation in this disk has disappeared due to evaporation until now.
Should a presumptive cool disk however be cut off already at about 104
Schwarzschild radii the coronal solution has to allow for the escape
of a wind (if not suppressed by the accreting gas ram pressure).
Such a solution has already been obtained earlier (Meyer et al. 2000, for
the case
,
"the weakest" coronal radiation case). In this case
the corona in the upper layers is more tenuous and and the predicted X-ray
radiation is somewhat less but with
erg/s
still far above the value of 1033.6 erg/s allowed by Chandra
observations (Baganoff et al. 2003).
We find that for our solutions, the thermal flux
given by the Spitzer formula never exceeds the saturated value for free
streaming electrons and the mean free path remains always less than
the thermal scale height (
)). Since
is required
at the upper boundary,
as the density becomes smaller and
smaller in the upper layers, the temperature gradient approaches zero.
This results in a very small heat flux (from the Spitzer formula), which
is smaller than the "free-streaming'' heat flux. Though the mean free
path is quite large at these upper layers it also remains smaller than
the temperature scale height which approaches infinity.
We also compared the time needed for equipartition temperature between electrons and ions with the thermal timescale and confirmed that the former is always shorter than the latter. This means that in such a corona, far from the central black hole, electrons and ions are well coupled. Therefore, these models for the corona above a cool disk around the Galactic Center are self-consistent.
Recently, Nayakshin (2004) discussed an alternative model which suggests that all the hot gas captured by the black hole is deposited into an inactive disk (see also Nayakshin 2003). In that model, the corona is extremely tenuous so that the electron mean free path becomes large compared to the pressure scale height of the corona. Thus, a "free-streaming electron'' thermal conduction law is applied instead of the classic Spitzer formula. The corona is approximated as a homogeneous column on top of a very thin transition region to the cool disk, in contrast to the vertically layered corona discussed in this paper. The free-streaming particles of long mean free path are thought to transport the thermal energy into the cool disk at the saturation speed of the order of the sound velocity in the corona, thus the coronal gas would be able to sink down into the cool disk at a fraction of that speed and condense as tenuous gas without significant radiation in the X-ray band, low enough not to get into conflict with the observed low X-ray luminosity.
This model lies in a quite different parameter space as the the one discussed in our paper. It differs in the assumption of a very thin transition between the tenuous hot corona and the dense cool disk from the extended transition layer obtained in our model.
Assuming that both these different models are self-consistent the question arises which one might be realized in nature. We have seen that the addition of the hot tenuous accreting gas to the already existing layered corona discussed in this paper only constitutes a minor contribution to the dominant coronal evaporation process with its significant radiation in the X-ray band. Such coronal evaporation models were successfully used to explain the formation of inner holes in quiescent accretion disks and the soft/hard transition in spectra of soft transients and high-mass X-ray binaries. On the other hand if the tenuous non-radiative accretion would be set up before the standard evaporation could have established itself it might have prevailed until now. This question requires further investigation.
The evaporation rate depends on the value of the viscosity assumed.
For a standard value of
and dominant evaporation,
,
one has
corresponding to
/yr. Observations of
stars at the Galactic Center and especially spectroscopy of
one such star, S0-2, suggests that these are main sequence O/B stars
(Eisenhauer et al. 2003; Ghez et al. 2003).
The O/B stars close to the Galactic Center could not have formed
longer before their main sequence lifetime, of order of
106.5-107 yr (Maeder & Meynet 1989). Thus a disk that remained after the stars had
been formed but has evaporated by now should have contained gas not more
than 300 to 1000
.
This value is in the range of mass of the presently
observed bright O/B stars close to the Galactic Center. Interestingly,
but, perhaps by accident, it is also close to the stability limit of a
disk against self-gravitation. For an evaporating disk
with effective temperature of 50 K
at distance 1016.9 cm, the mass of the disk of the stability
limit (see e.g., Gammie 2001) is about
.
Thus a disk that had become unstable by self-gravitation and formed the presently observed young massive stars around the Galactic Center until the gravitational instability had ceased could perhaps have now completely disappeared by the process of coronal evaporation. But this remains a rather speculative question until we know more about the stellar population and its origin so close to the Galactic Center.
We have developed a model for a cool disk around the Galactic Center that has a corona above it, allowing for a coronal accretion of gas captured at the Bondi radius from stellar winds of massive stars. Our models incorporate hydrostatic layering, thermal heat conduction, friction, and radiative energy loss and allow for mass exchange between disk and corona by condensation or evaporation of gas, as well as mass gain or loss by radial flow in the corona (the latter a generalization of the former one-zone model).
We find that the solutions we obtain depend on the value of the
viscosity parameter .
For standard
values,
,
so much mass evaporates from the disk
into the hot corona that the additional mass flow from the outside is
a negligible contribution. The evaporation then with time could lead to
the complete disappearance of the original disk. Only for values
there might be solutions in which the incoming accreted gas
condenses into the disk. However, the vertical structure computations
show that in all these solutions the calculated
Bremsstrahlung of the corona in the X-ray band by far exceeds the
observed luminosity.
From this we conclude that, if our modeling is correct and applicable, at present no cool disk around the Galactic Center exists. We shortly discussed what limit this puts on an inert disk that might have originally remained from a phase of star formation in which the young bright stars presently seen close to the Galactic Center were formed. We also discuss the difference to the alternative non-radiative condensation model of Nayakshin (2004).
Acknowledgements
B. F. Liu would like to thank the Alexander von Humboldt-Foundation for support.