Issue |
A&A
Volume 499, Number 1, May III 2009
|
|
---|---|---|
Page(s) | 349 - 355 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/200810956 | |
Published online | 27 March 2009 |
Thermal instability in X-ray photoionized media in active galactic nuclei
II. Role of the thermal conduction in warm absorber
B. Czerny1 - L. Chevallier1,2 - A. C. Gonçalves2,3,4 - A. Rózanska1 - A.-M. Dumont2
1 - Copernicus Astronomical Center, Bartycka 18,
00-716 Warsaw, Poland
2 -
LUTH, Observatoire de Paris, CNRS,
Université Paris Diderot, 5 Place Jules Janssen,
92190 Meudon, France
3 -
Observatoire Astronomique de Strabsourg, Université
Louis Pasteur, CNRS, 11 rue de l'Université,
67000 Strasbourg, France
4 -
CAAUL, Observatório Astronómico de Lisboa,
Tapada da Ajuda, 1349-018 Lisboa, Portugal
Received 12 September 2008 / Accepted 19 February 2009
Abstract
Context. A photoionized gas under constant pressure can display thermal instability with three or more solutions for possible thermal equilibrium. A unique solution of the structure of the irradiated medium is obtained only if electron conduction is considered.
Aims. The subject of our study is to estimate how the efect of thermal conduction affects the structure and transmitted spectrum of the warm absorber, computed by solving radiative transfer with the code TITAN.
Methods. We developed a new computational mode for the code TITAN to obtain several solutions for a given external conditions and we test a posteriori which solution is the closest one to the required integral condition based on conduction.
Results. We demonstrate that the automatic mode of the code TITAN provides the solution to the radiative transfer, that is generally consistent with the estimated exact solution within a few per cent accuracy, with larger errors for some line intensities (up to 20 per cent) for iron lines in an intermediate ionization state.
Key words: instabilities - radiation mechanisms: thermal - radiative transfer - galaxies: active - X-rays: general - X-rays: galaxies
1 Introduction
Recent modeling of the warm absorbing medium in the central regions of active galactic nuclei (AGN) has shown that the cloud models in constant total (gas plus radiation) pressure (Gonçalves et al. 2006; Rózanska et al. 2006; Chevallier et al. 2007) are more sophisticated than those in constant density (Kaastra et al. 2002; Kaspi et al. 2001; Netzer 1993), or those invoking a continuous wind (Steenbrugge et al. 2005; Chelouche & Netzer 2005). Constant pressure models, however, are not simple to deal with, as they require addressing the subject of thermal instability.
The existence of thermally unstable bubbles in gas affected by
external illumination was
studied by Field (1965) with applications to the interstellar
and intergalactic media.
More generally, thermal instabilities
develop in photoionized gas in thermal equilibrium, i.e. in gas
where the radiative heating is balanced by radiative cooling
(e.g. Krolik et al. 1981).
These authors have shown that for a specific range of the ionization
parameter, a medium illuminated by X-rays can achieve local
thermal equilibrium for three different values of the
temperature and the density.
Such multiple solutions provide a natural explanation for
the two-phase gas scenario where colder and denser clouds at
K are
embedded in a hotter medium at
K and
are at the origin of the description of the interstellar medium,
the warm absorber, and illuminated accretion disk atmospheres in terms of a
clumpy medium (Rózanska 1999; Rózanska & Czerny 1996; Torricelli-Ciamponi & Courvoisier 1998; Raymond 1993; Krolik et al. 1981; Smith & Raine 1988; Krolik & Kriss 2001).
The uniqueness of the equilibrium state can be restored when the electron conduction is taken into account (e.g. Begelman & McKee 1990; McKee & Begelman 1990). Mathematically, it happens because the energy balance equation with conduction becomes a differential instead of an algebraic equation. Physically, the equilibrium distribution is achieved either by stationary heat transport or by evaporation or condensation of the material.
Analytical estimations of the role of thermal conducion in planetary nebulae (Sage & Seaton 1973) have shown that thermal conduction is not an important process for heating filamentary regions in which the ionization is maintained by diffuse radiation.
The influence of thermal conduction in AGN and galactic black holes (GBH) has been studied in several papers, in the context of disk/corona boundary (Rózanska 1999; Dullemond 1999; Maciolek-Niedzwiecki et al. 1997), disk evaporation (Rózanska & Czerny 2000b; Meyer & Meyer-Hofmeister 1994; Liu et al. 1999), clumpy accretion inflow (Torricelli-Ciamponi & Courvoisier 1998), and clumpy outflow (Krolik 1998). But the efect has never been studied in the case of warm absorbers. This is important, since high-resolution spectra by X-ray statellites require proper models.
Presently, no complex radiative transfer code for regions illuminated by X-rays addresses the conduction term; therefore, codes tread thermal instabilities by two main approaches: i) by fixing the discontinuity position at the optical depth found during the first integration (Nayakshin et al. 2000), ii) by finding a continuous intermediate but approximate solution, as discussed first by Rózanska et al. (2002) and later in more detail by Gonçalves et al. (2007, hereafter Paper I).
In Paper I we estimated the maximum error that can be made due to arbitrariness in the discontinuity location when the conduction is not included. In this paper (hereafter Paper II), we develop a method determining the exact location of the discontinuity, by iteratively determining the warm absorber cloud structure and radiative transfer complete solution with a self-consistent location of the discontinuity and by comparing this spectrum to the spectrum produced in a standard way by the TITAN code. This position is crucial since, in the models of warm absorber under constant pressure, the amount of energy absorbed in lines strongly depends on the structure of the cloud.
In the next section, we comment on the model and the TITAN transfer-photoionization code, and we describe our approach to incorporating the effect of thermal conduction. Results are presented and compared to the analytical model of Sage & Seaton (1973) in Sect. 3. Section 4 is dedicated to the physical discussion of these models.
2 Methodology
It is well known that a numerical difficulty appears when trying to
estimate the strength of thermal conduction. The fundamental
quantity, i.e. the heat conductive flux, is proportional to the
temperature gradient. In the numerical approach to the radiative
transfer, the computations are performed using a certain grid,
defined by a given number of layers with a given thickness, that
are used to reproduce the total gas slab.
This grid imposes a limit to the temperature gradient. Even
if the temperature distribution clearly shows a discontinuity, i.e.
the temperature difference between two consecutive layers, ,
is a large fraction of the highest value, T, the estimated
temperature gradient,
may not be high since the
geometrical width of the layer
,
where the sharp
temperature drop occurs, is set by the coarse spatial grid.
To estimate accurately the heat conductive flux, the grid would
have to be refined too much for practical use, since the actual
spatial extension of the transition zone, is
close to the Field length (see: Rózanska 1999),
and its optical depth for
electron scattering,
,
can be estimated as
![]() |
(1) |
where



In this study, we use an integral formulation of the solution to the problem, which is not grid-dependent (Rózanska & Czerny 2000a). The method is based on the integral criterion for the location of the discontinuity. Therefore, we do not include the conduction flux in the energy balance but aim at determining the unique position of the discontinuity. To estimate this unique position, we solve the radiative transfer within an irradiated cloud under constant pressure using the TITAN code.
TITAN is a transfer-photoionization code developed by
Dumont et al. (2000) and Collin et al. (2004) to correctly
model optically thick (Thomson optical depth up to several tens),
ionized media; it can be applied equally to thinner media
(Thomson depth 0.001 to 0.1).
The code solves full radiative transfer
through the photoionized gas, together with the ionization equilibrium, the
thermal equilibrium, and the statistical equilibrium of all the levels of
each ion. Our atomic data include
lines from ions
and atoms of H, He, C, N, O, Ne, Mg, Si, S, and Fe.
The code relies on
the accelerated lambda iteration (ALI) method, which allows
for exact treatment of the transfer of both the continuum
and the lines. It includes all relevant physical processes
(e.g., photoionization, radiative and dielectronic
recombination, fluorescence and Auger processes, collisional
ionization, radiative and collisional excitation/de-excitation)
and all induced processes. As an output, it
gives the ionization, density, and temperature structures,
as well as the reflected, outward, and transmitted spectra.
The energy balance is ensured locally with a precision of
0.01%, and globally with a precision of 1%. For more details
on the code and its evolution, see Paper I and references therein.
The code presently offers three operational modes.
- (i)
- The first mode, used in a number of previous papers
(e.g. Collin et al. 2003; Rózanska et al. 2002; Dumont et al. 2002), is based
on subsequent (instead of simultaneous) determination of the
density and the temperature. As explained in Paper I, this solution
is smooth (without a sharp discontinuity), but it does not accurately satisfy
the condition of the constant pressure. Such a solution corresponds
to the default, automatic mode of the TITAN code and it will be
denoted here as an ``intermediate state'' (I).
- (ii)
- The second mode, described in Paper I, precisely preserves the
condition of the constant pressure and it locally finds multiple
solutions to the temperature and the density in the instability
zone. In this mode, the code is set to choose either the
hottest solution, or the coldest solution, as the actual solution.
In this way two separate solutions are obtained: one
with the most extended hot zone at the cloud irradiated surface
(state H) and the second one with the narrowest hot zone (state C).
- (iii)
- A third operation mode has been developed for the
purpose of this work. In such a mode,
multiple local solutions for the density and the temperature are
searched for, but the actual discontinuous transition from the
hot to the cold branch is set in advance, somewhere within the
instability strip, at an arbitrarily imposed value of the
column density,
.
2.1 Conduction flux and the location of the transition zone
We used this third operational mode of the code to
calculate a series of models for an ionized cloud
with a given set of parameters, only varying the position
x of the temperature discontinuity. We denote such states as
NH(x), with the parameter x as the column density in
1023 units, i.e.
cm-2. For each of
the models NH(x), we have calculated the value of the
integral A (right hand of Eq. (A.3), see Appendix A)
with the purpose of verifying whether the position of the
discontinuity is consistent with the criterion for the discontinuity
position provided by the conduction flux (A = 0.; see
Eq. (A.4)). In this way we estimated
the position x for which the value of the integral A is the
closest to zero. The corresponding model was taken as the correct
solution to the radiative transfer with the conduction term. We
then analyzed the properties of this model, and compared
them to the properties of the model obtained from the TITAN
computations in its automatic mode (state I). These results were
used to quantitatively estimate the errors made in the transition
spectra when the faster, automatic mode is used. This is important
since the method used in the present paper to account for
the electronic conduction is
highly time-consuming, although physically correct, and it would be
unrealistic to use it as such
for interpretating observational data.
3 Results
3.1 Reference cloud
Our reference cloud was modeled under
constant total pressure equilibrium. It is characterized
by an ionizing parameter erg s-1 cm, a column density
cm-2, and a number density
cm-3, on the illuminated side of the cloud.
Such cloud parameters fall within the range of conditions valid
for the warm absorber in active galactic nuclei.
Abundances are cosmic (Allen 1973). This is a simple choice since the present work is aimed at testing the sensitivity of the cloud structure to the adopted description of the thermal instability. The actual abundances in a given AGN are unknown. Warm absorber material is flowing out from the nuclear region of an AGN, so this material is likely to have a long and complex history, because it originates as an interstellar material of a host galaxy, later reaching the inner region where the intense starburst activity takes place (such activity frequently - or perhaps always - accompanies the nuclear activity), the remnants of this starburst activity flow toward the black hole, thereby forming an accretion disk. Finally, this material flows out in the form of a disk wind before actually reaching a black hole. Therefore, in modeling a given active nucleus, we actually should determine the chemical composition from the fits to the data. In most cases, when such an attempt has been made (mostly for iron, occasionally for lighter elements) the composition was roughly solar within a factor of a few. Therefore, the choice of standard composition is natural for the test purposes.
We assumed this cloud is normally illuminated by a
power-law incident flux with slope
from 10 eV to 100 keV.
Such a cloud is highly stratified in temperature, because it is subject
to thermal instability. To model our cloud, the medium
was divided into a few hundred layers, using an adaptive
grid. The first layer is hot (
K) and shows no
thermal instability. The last layer is cold (
K) and
it shows no thermal instability either. In the innermost regions
of the ionized cloud, two instability zones exist (see Paper I), which
may overlap partially. The whole thermal instability zone spans the
range
cm-2.
The extension of the instability zones and consequently the limits to the exact location of the temperature discontinuity is provided by the two extreme models, H and C, corresponding to the lowest and the highest values of the hydrogen column density measured at the discontinuity.
For our reference cloud, the values of the hydrogen column are
cm-2 for state H, and
cm-2 for state C. In our notations,
those values correspond to the NH(1.92) and NH(1.54) states, respectively.
The temperature profiles with the sharp drops for the states C and H
are shown in Fig. 1.
Since the intermediate
state is closer to the H state than to the C state, we calculate
the sequence of the intermediate states from NH(1.78) to NH(1.89).
For all these states we calculated the value of the integral A from Eq. (A.4) in the Appendix A. This requires the construction of the heating/cooling curves as a function of the temperature in various layers. The curves have a complex shape, since multiple instability zones actually exist within the cloud. The integral gains a positive contribution from the beginning of the instability zone as the heating/cooling is positive, whereas the contribution is negative at the end of the transition region. This result is depicted in Fig. 2.
![]() |
Figure 1: Temperature profile versus the column density for our reference cloud in the extreme C and H states, in the intermediate I state resulting from TITAN in its automatic operational mode, and in the state NH(1.86) with the position of the temperature drop consistent with the conduction based criterion. Inserted box: zoom of the temperature profile on the thermal instability zone. |
Open with DEXTER |
The best agreement with the condition A = 0 is achieved for
a given state NH(x) with
1.85 < x < 1.87, consistent
with the real temperature drop being located somewhere between
cm-2 and
cm-2.
Because the NH(x) model computations are extremely time-consuming,
we did not calculate the
1.85 < x < 1.87 states very densely.
Among the computed models, the NH(1.86) state offers the best
agreement with the condition A = 0.
![]() |
Figure 2: Variation in integral A (units erg2 s-2 cm-4) as a function of the column density (units cm-2) for several intermediate states labeled for the adopted position of the discontinuity. The best agreement between the adopted position and the position indicated by A=0 is obtained for model NH(1.86). |
Open with DEXTER |
Defining our equilibrium condition using
Eq. (A.4) given in the Appendix, we
assumed the left-hand part of this equation to be zero, i.e. A=0.
This requires that either
or else that both
and
are null (or at least negligible
as compared to the integral, when it is non null). We tested
that hypothesis a posteriori: calculating the
left-hand side of Eq. (A.3), we found it
to be negligible.
Finally, it is interesting to note that the models where the adopted temperature drop coincides roughly with the estimated location of the temperature drop from the conduction criterion are numerically stable (see Appendix B).
![]() |
Figure 3:
Transmitted continuum flux
|
Open with DEXTER |
3.2 Transmission spectra
In Paper I we have shown that the transmission spectra for the H and C states are considerably different, thus demonstrating the importance of a sufficiently accurate measure of the temperature drop location. Thanks to the method developed in this paper, we were able to identify the (approximate) location that is consistent with the condition resulting from the introduction of the electron conduction term. We are now capable of comparing the transmission spectrum obtained for state NH(1.86) with the spectra resulting from states H and C. These were also compared with the standard I state resulting from the default TITAN mode (faster and numerically more stable than the other modes).
Figure 3 shows a pure transmitted spectra
(i.e. covering factor )
at resolution
,
and Table 1 the comparison of the
strongest lines present in H, C, I, and NH(1.86) models.
The continuum in NH(1.86) shows somewhat less absorption than
in I, but the maximum difference, located around 1 keV,
is never more than 27%. The maximum relative difference between
the measured lines in NH(1.86) and I models is less than 24%
for the 200 strongest features. Most of the strongest lines
show a difference less than 1%, because they are saturated;
some lines show a 20% difference (e.g. Fe XII and
Fe XIII) that is related to a change of ionization fraction,
due to a different ionic column density (see Paper I).
For the outward and reflected spectra (emitted in all directions
but the normal), the difference between models I and
NH(1.86) is again weak, i.e. less than 25%. We therefore
conclude that the I model can be used safely to describe such a
medium, which falls within the range of conditions valid
for the warm absorber in AGN.
Table 1: Table highlighting the differences between the transmitted spectra resulting from I and NH(1.86) models, with lines given in order of their strength.
3.3 The structure of the conductive layer
Since we determined the position of the transition layer, we can calculate a posteriori the structure of this sharp transition zone in a semi-analytical way, including the conduction flux. For that purpose we first obtained the heating/cooling function as a function of the temperature for the total pressure characterizing the model NH(1.86) in the transition zone. We stress that the heating/cooling is obtained for a constant presure, not constant density, as this assumption esentially influences its shape, as nicely explained in Goncalves et al. (2007). The resulting heating/cooling curve is shown in Fig. 4.
![]() |
Figure 4: The heating/cooling curve for model NH(1.86) at the position of the conductive discontinuity as a function of the temperature, obtained under constant pressure condition. The curve crosses zero (radiative balance) for three values of the temperature. Compton heating as well as all atomic heating/cooling processes are included. |
Open with DEXTER |
Having this curve and assuming constant pressure in the sharp transition zone, we can now obtain the conduction flux as a function of the temperature from the integral form of the conduction equation (see Eq. (4) in the Appendix) by replacing
with a set of intermediate value of T. Next, returning to the expression for conduction flux, we can obtain the spacial temperature gradient in the constant pressure zone. The resulting flux and the temperature profile are shown in Fig. 5. The hydrogen column within this zone is equal to
cm so the zone unresolved in radiative transfer computations should not contribute strongly to the spectrum although it is not totally negligible.
![]() |
Figure 5:
The spacially resolved conduction flux ( upper panel) and temperature ( lower panel) profiles in the transition zone obtained from semi-analytical approximation (see text). The transition is sharp in comparison with extension of a whole reference cloud
|
Open with DEXTER |
In the transition zone the conductive flux is still orders of magnitude smaller than the radiative flux but nevertheless it contributes significantly to the total energy balace in this zone since the derivative of the conduction flux is large. The conductive flux balances the departure of the heating/cooling curve from zero and since the single heating or cooling mechanism gives values of order of 10-24 erg s-1 cm3, the flux derivative is about the same order of magnitude as the heat deposit or loss, just due to the geometrical narrowness of the zone. Therefore, comparing the two fluxes instead of their derivatives, as done e.g. by Sage & Seaton (1973) in their analysis, may be misleading.
The overall significance of the conduction-dominated narrow zone itself is of course affected by its narrowness, and from this point of view the simple argument of Sage & Seaton (1973) is correct and the zone does not contribute strongly to the final transmission spectra so that we neglected this zone completely in our numerical computations of the radiative transfer is fully justified.
The important effect of electron conduction is, however, in determining of the position of the discontinuity, and this is the real effect we search for, as explained in Sect. 2. With these semi-analytical results, we can now check that the conduction is effective enough to move the position of the transition zone.
If the initial position of the transition zone is not consistent with the conduction criterion (see Appendix) for some reason, the transition zone will move with the velocity roughly given by the ratio of the conduction flux to the local energy density, i.e.
km s-1. This speed is strongly subsonic. The whole size of the reference cloud where three solutions for the radiative equilibrium are possible covers the range of
cm-2 in the hydrogen column and the spatial distance of
cm, i.e. the conduction flux will adjust the front position in about 3 years. A slowly moving (v < 500 km s-1) distant (above
cm) cloud in the radiation field of an active nucleus has time to reach such equilibrium since the motion of the front due to systematically decreasing flux is slower than the effect of conduction estimated above. That the irradiation in an AGN is time-dependent complicates the problem, as discussed by Chevallier et al. (2007).
A question remains whether such a narrow zone can exist and be stable. The first issue is related to the mean free path of electrons. In the likely presence of a weak magnetic field, this seems to be possible. The mean free path of an electron for a micro-Gauss magnetic field is a fraction of a cm, much smaller than the size of the transition region. The stability issue is another problem. A study e.g. by Rózanska (1999) in the context of disk/corona border indicates marginal stability; however, a magnetic field may also reduce heat conduction, and in this case we may strongly overestimate the conduction flux. The issue was discussed in other astrophysical contexts like wind collisions, but the theory is rather complex and uncertain. Recent Chandra observations of a planetary nebula PN BD+303639 by Nordon et al. (2009) suggest complex structure of such a zone, with ions penetrating the contact discontinuity boundary. The count rate in AGN observation is much lower, so at present the search for such traces in AGN warm absorber is likely to be premature.
4 Discussion
This study represents a strong contribution to the understanding of the technical aspects behind radiative transfer computations in stratified media under total constant pressure. It is well known that an irradiated gas in constant pressure will develop thermal instabilities, if the electron conduction is neglected (e.g. Krolik et al. 1981). These instabilities appear as technical problems in radiative transfer and photoionization codes; therefore, various codes use different approaches to address this problem.
In this paper we tested the different numerical approaches used by the TITAN code. We constructed a subset of radiative transfer solutions, all for the same column density of the cloud and the same incident radiation flux, and identified the one that is close to satisfying the equilibrium condition derived from the electron conduction theory. We found that the approximate solution obtained by the standard, automatic mode of the TITAN code satisfactorily representats the transmitted spectra. We find that most line intensities resulting from this automatic mode agree with the values derived from the electron conduction theory, reaching an accuracy on the order of 1%, with only a few exceptions (mostly iron lines in an intermediate ionization state) where the discrepancy between the solutions reached 20%.
We tested our methodology and hypothesis on a given cloud, described by a single set of parameters; these were chosen within the range of conditions usually met in the modeling of the X-ray ionized gas near the central regions of AGN, and they are fully compatible with a warm absorber medium. One issue remains, that of computations based on the conduction criterion being much more time-consuming (at least by a factor 20) than the standard radiative transfer computations, which are already quite long. Therefore, it is not feasible to apply this computation method to directly obtaining solutions for a grid of models, for instance. The results of this work provide, nevertheless, important evidence that the standard, automatic computing mode of the TITAN code is not strongly affected by the systematical errors caused by neglecting to include the electron conduction. This is important, since constant pressure models are a much better approximation than constant density models, with applications in several physical circumstances, e.g. the warm absorber clouds, or the disk/corona transition zone in an accretion disk onto a black hole.
On the other hand, our analysis does not address the evolutionary aspect of the presence of thermal instability. In this paper, as in other papers devoted to the modeling of the warm absorber, we assumed that the very existence of the clouds (instead of a continuous wind) is a result of this instability, but we did not address the timescales for the cloud formation and its survival. A first attempt at this exercise can be found in Paper I and in Chevallier et al. (2007). Additional discussions of these aspects of the thermal instability have been reported in other contexts i.e.: the hot medium in galaxy clusters (Kim & Narayan 2003) and the multiphase interstellar medium (Inoue et al. 2006). In the case of AGN, the issue of the medium fragmentation is likewise made more complex by the fact that the irradiation flux itself varies significantly in time, so the cloud can never achieve a full thermal equilibrium (e.g. Chevallier et al. 2007). Those issues are clearly beyond of the scope of the present paper.
Acknowledgements
We thank S. Collin and V. Karas for helpful discussions, and we are grateful for the hospitality provided by the Astronomical Institute, Academy of Sciences of the Czech Republic, Prague, where part of this project was completed. A. C. Gonçalves acknowledges support from the Fundação para a Ciência e a Tecnologia, Portugal, under grant No. BPD/11641/2002. Part of this work was carried out within the framework of the European Associated Laboratory ``Astrophysics Poland-France'', and was supported by grant 1P03D00829 of the Ministry of Science and Education, and the Polish Astroparticle Network 621/E-78/BWSN-0068/2008.
Appendix A
We recall the main relations for energy balance with the thermal
conduction effect included (see Rózanska & Czerny 2000a), which are
used and eventually modified for the purpose of this study.
The energy balance at z (i.e. which permits computing the
temperature at this position) across the boundary between the hot and
cold media in the static case is given by
where Q+ and






where


Under the condition of constant gaseous pressure (i.e. perfect gas),
we can express Eq. (A.1) by performing an integral
(Eq. (8) in Rózanska & Czerny 2000a) over the temperature in the
thermal instability region from the
cold solution
to the hot solution
.
Under the
constant total (gaseous + radiative) equilibrium of our clouds
,
both gaseous and radiative pressures
vary in the thermal instability region, in a non-local dependence of
the temperature. This assumption cannot be made on theoretical
grounds; however, we observe that this variation is weak, so the
non-local dependence is not drastic, and use the local value of the
gaseous and radiative pressure is consistent with the way we
computed the heating and cooling curves for our code. Equation (8) from
Rózanska & Czerny (2000a) becomes, when inserting the factor
,
where
erg K-1 is the Boltzmann constant,
where the conductive heat flux at temperatures




Appendix B
![]() |
Figure B.1: The convergence properties of our reference cloud as a function of the iteration number up to iteration 50, for the states H (no change up to iteration 65) C, and NH(1.86). |
Open with DEXTER |
It is interesting to note that finding the exact location of the discontinuity in the model also provides a solution to some numerical problems that are met in TITAN used in the second (nonstandard) operational mode, used in Paper I. If we attempt to determine the type H state of the cloud (i.e. we always choose the hottest solution in the multiple solution region; see Sect. 2), we actually face the problem of the small regular oscillations in the state properties, which translate into the small but finite oscillations in the total energy balance of the cloud. To illustrate this we present in Fig. B.1 the dependence of our estimator of convergence on the number of iterations for different states. The oscillations in total energy balance correspond to the oscillations in the location of the main temperature drop (see Fig. 1).
The oscillations cease to exist if we impose the position of the temperature discontinuity quite close to the position required by the conduction condition A = 0. Specifically, oscillations are not present in all models from NH(1.82) to NH(1.89).
References
- Allen, C. W. 1973, Astrophysical quantities (London: University of London, Athlone Press), 3rd edn. (In the text)
- Begelman, M. C., & McKee, C. F. 1990, ApJ, 358, 375 [NASA ADS] [CrossRef]
- Chelouche, D., & Netzer, H. 2005, ApJ, 625, 95 [NASA ADS] [CrossRef]
- Chevallier, L., Czerny, B., Rózanska, A., & Gonçalves, A. C. 2007, A&A, 467, 971 [NASA ADS] [CrossRef] [EDP Sciences]
- Collin, S., Coupé, S., Dumont, A.-M., Petrucci, P.-O., & Rózanska, A. 2003, A&A, 400, 437 [NASA ADS] [CrossRef] [EDP Sciences]
- Collin, S., Dumont, A.-M., & Godet, O. 2004, A&A, 419, 877 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Dullemond, C. P. 1999, A&A, 341, 936 [NASA ADS]
- Dumont, A.-M., Abrassart, A., & Collin, S. 2000, A&A, 357, 823 [NASA ADS] (In the text)
- Dumont, A.-M., Czerny, B., Collin, S., & Zycki, P. T. 2002, A&A, 387, 63 [NASA ADS] [CrossRef] [EDP Sciences]
- Field, G. B. 1965, ApJ, 142, 531 [NASA ADS] [CrossRef] (In the text)
- Gonçalves, A. C., Collin, S., Dumont, A.-M., et al. 2006, A&A, 451, L23 [NASA ADS] [CrossRef] [EDP Sciences]
- Gonçalves, A. C., Collin, S., Dumont, A.-M., & Chevallier, L. 2007, A&A, 465, 9 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Inoue, T., Inutsuka, S.-I., & Koyama, H. 2006, ApJ, 652, 1331 [NASA ADS] [CrossRef] (In the text)
- Kaastra, J. S., Steenbrugge, K. C., Raassen, A. J. J., et al. 2002, A&A, 386, 427 [NASA ADS] [CrossRef] [EDP Sciences]
- Kaspi, S., Brandt, W. N., Netzer, H., et al. 2001, ApJ, 554, 216 [NASA ADS] [CrossRef]
- Kim, W.-T., & Narayan, R. 2003, ApJ, 596, 889 [NASA ADS] [CrossRef] (In the text)
- Krolik, J. H. 1998, ApJ, 498, L13 [NASA ADS] [CrossRef] (In the text)
- Krolik, J. H., & Kriss, G. A. 2001, ApJ, 561, 684 [NASA ADS] [CrossRef]
- Krolik, J. H., McKee, C. F., & Tarter, C. B. 1981, ApJ, 249, 422 [NASA ADS] [CrossRef] (In the text)
- Liu, B. F., Yuan, W., Meyer, F., Meyer-Hofmeister, E., & Xie, G. Z. 1999, ApJ, 527, L17 [NASA ADS] [CrossRef]
- Maciolek-Niedzwiecki, A., Krolik, J. H., & Zdziarski, A. A. 1997, ApJ, 483, 111 [NASA ADS] [CrossRef]
- McKee, C. F., & Begelman, M. C. 1990, ApJ, 358, 392 [NASA ADS] [CrossRef]
- Meyer, F., & Meyer-Hofmeister, E. 1994, A&A, 288, 175 [NASA ADS]
- Nayakshin, S., Kazanas, D., & Kallman, T. R. 2000, ApJ, 537, 833 [NASA ADS] [CrossRef] (In the text)
- Netzer, H. 1993, ApJ, 411, 594 [NASA ADS] [CrossRef]
- Nordon, R., Behar, E., Soker, N., Kastner, J. H., & Yu, Y. S. 2009, arXiv e-prints (In the text)
- Raymond, J. C. 1993, ApJ, 412, 267 [NASA ADS] [CrossRef]
- Rózanska, A. 1999, MNRAS, 308, 751 [NASA ADS] [CrossRef]
- Rózanska, A., & Czerny, B. 1996, Acta Astron., 46, 233 [NASA ADS]
- Rózanska, A., & Czerny, B. 2000a, MNRAS, 316, 473 [NASA ADS] [CrossRef] (In the text)
- Rózanska, A., & Czerny, B. 2000b, A&A, 360, 1170 [NASA ADS]
- Rózanska, A., Dumont, A.-M., Czerny, B., & Collin, S. 2002, MNRAS, 332, 799 [NASA ADS] [CrossRef] (In the text)
- Rózanska, A., Goosmann, R., Dumont, A.-M., & Czerny, B. 2006, A&A, 452, 1 [NASA ADS] [CrossRef] [EDP Sciences]
- Sage, G., & Seaton, M. J. 1973, in Les Nébuleuses Planétaires, 241 (In the text)
- Smith, M. D., & Raine, D. J. 1988, MNRAS, 234, 297 [NASA ADS]
- Steenbrugge, K. C., Kaastra, J. S., Crenshaw, D. M., et al. 2005, A&A, 434, 569 [NASA ADS] [CrossRef] [EDP Sciences]
- Torricelli-Ciamponi, G., & Courvoisier, T. J.-L. 1998, A&A, 335, 881 [NASA ADS]
All Tables
Table 1: Table highlighting the differences between the transmitted spectra resulting from I and NH(1.86) models, with lines given in order of their strength.
All Figures
![]() |
Figure 1: Temperature profile versus the column density for our reference cloud in the extreme C and H states, in the intermediate I state resulting from TITAN in its automatic operational mode, and in the state NH(1.86) with the position of the temperature drop consistent with the conduction based criterion. Inserted box: zoom of the temperature profile on the thermal instability zone. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Variation in integral A (units erg2 s-2 cm-4) as a function of the column density (units cm-2) for several intermediate states labeled for the adopted position of the discontinuity. The best agreement between the adopted position and the position indicated by A=0 is obtained for model NH(1.86). |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Transmitted continuum flux
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: The heating/cooling curve for model NH(1.86) at the position of the conductive discontinuity as a function of the temperature, obtained under constant pressure condition. The curve crosses zero (radiative balance) for three values of the temperature. Compton heating as well as all atomic heating/cooling processes are included. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
The spacially resolved conduction flux ( upper panel) and temperature ( lower panel) profiles in the transition zone obtained from semi-analytical approximation (see text). The transition is sharp in comparison with extension of a whole reference cloud
|
Open with DEXTER | |
In the text |
![]() |
Figure B.1: The convergence properties of our reference cloud as a function of the iteration number up to iteration 50, for the states H (no change up to iteration 65) C, and NH(1.86). |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.