A&A 464, 631-634 (2007)
DOI: 10.1051/0004-6361:20066564
Y. Zhang^{1,2} - B. Ercolano^{3} - X.-W. Liu^{1}
1 - Department of Astronomy, Peking University,
Beijing 100871, PR China
2 -
Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong, PR China
3 -
Harvard-Smithsonian Centre for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
Received 13 October 2006 / Accepted 22 November 2006
Abstract
Aims. We investigate temperature fluctuations in H II regions in terms of a two-phase model, which assumes that the nebular gas consists of a hot and a cold phase.
Methods. We derive general formulae for T([O III]), the [O III] forbidden line temperature, and T(H I), the hydrogen Balmer jump temperature, in terms of the temperatures of the hot and cold phases,
and .
Results. For large temperature differences, the values of t^{2} required to account for the observed difference between T([O III]) and T(H I) are much lower than those deduced using the classical formulae that assume random and small amplitude temperature fluctuations. One should therefore be cautious when using a two-phase model to account for empirically derived t^{2} values. We present a correction of a recent work by Giammanco & Beckman, who use a two-phase model to estimate the ionization rate of H II regions by cosmic rays. We show that a very small amount of cold gas is sufficient to account for t^{2} values typically inferred for H II regions.
Key words: ISM: general - ISM: H II regions
Temperature fluctuations in H II regions are a much-discussed problem. Peimbert (1967) investigated for the first time the effects of such fluctuations on temperatures empirically derived from spectroscopic observations and found that they may lead to higher electron temperatures being derived from the collisionally excited [O III] nebular-to-auroral forbidden-line ratio, T([O III]), than those derived from the Balmer jump of the H I recombination spectrum, T(H I). If significant temperature fluctuations exist in H II regions and yet are ignored in the analysis, they may lead to underestimating ionic abundances calculated from collisionally excited lines (CELs) (e.g. Esteban et al. 2002). The parameter t^{2} (Peimbert 1967) was introduced to quantitatively characterize temperature fluctuations. Since then the parameter has been extensively used in nebular studies. Under the conditions that , where T(r) is the local electron temperature and T_{0} is the average value weighted by the square of density, the value of t^{2} can be determined either by comparing T([O III]) and T(H I) or by comparing ionic abundances derived from CELs and from recombination lines (RLs) (see Peimbert et al. 2004, for further details).
Two-phase models, which approximate a nebula by two components of different physical conditions, represent an over-simplified, yet frequently used method of studying nebular physics (e.g. Viegas & Clegg 1994; Zhang et al. 2005). Using a two-phase model, Stasinska (2002) points out that the classical picture of temperature fluctuations may be misleading under certain conditions, and three parameters are needed to characterize temperature inhomogeneities. One of the open questions in the study of H II regions is that values of t^{2} derived from observations are consistently higher than those predicted by photoionization models (Stasinska 2000). Recently, Giammanco & Beckman (2005, GB05 thereafter) constructed a two-temperature-phase model capable of explaining t^{2} values deduced for a number of H II regions by means of incorporating a component of cool ionic gas ionized by cosmic rays. However, t^{2} values deduced from observations cannot be applied directly to two-phase models. This is because empirical values of t^{2} deduced from observations were calculated from formulae derived by assuming random and small-amplitude temperature fluctuations, assumptions that are apparently broken for a two-phase model.
The purpose of the current work is to quantitatively study the relationship between the t^{2} values predicted by two-phase models and those measured by observations. We show that, in the case of a very cold ionic gas component embedded in a "normal'' H II region, t^{2} deduced from observations using the empirical method may have significantly overestimated the real values.
According to Peimbert (1967), for a given ionic species of number
density ,
the thermal structure of an H II region can be
characterized by an average temperature T_{0} and a mean square temperature
fluctuation parameter t^{2}, defined as
(1) |
(2) |
In the framework of the two-phase model, which assumes that the electron
temperature structure of an H II region consists of a hot and a
cold phase, we follow GB05 and assume equal densities of the two phases and an
ionization fraction of unity for the hot gas. Electron temperatures are
designated as
and
for the hot and cold phases,
respectively. The intensity of an [O III] forbidden line transition of
wavelength
is given by
(6) |
Similarly, from the flux of the Balmer jump,
(8) |
(9) |
In Fig. 1 for given values of and , we compare t^{2} as a function of derived from the empirical method using Eqs. (3) and (4), and that derived in the scenario of two-phase model using Eqs. (7), (10), and (13). The plots show that, depending on above a critical value of , the empirical method significantly overestimates t^{2}, particularly for the case of small . The amount of deviation is insensitive to the value adopted for . In addition, we find that as the temperature difference between the two phases is larger than a critical value (typically 6000 K), empirical t^{2}deduced from observations can no longer be used at their face values to constrain two-phase models, a point overlooked by GB05 as discussed in the following section.
Table 1: Estimated values of t^{2}, T_{0}, and for , 1000, and 4000 K for a sample of H II regions. and are taken from Esteban et al. (2002). The numbers in parentheses are values derived by GB05 and included here for comparison.
Figure 1: t^{2} versus (T([O III])-T(H I)), deduced from the empirical method (Eqs. (3) and (4), dotted lines) and for two-phase model (Eqs. (7), (10), (12), and (13), solid lines). Left panels: = 0.01; right panels: = 0.1; upper panels: = 10 000 K; lower panels: = 15 000 K. | |
Open with DEXTER |
GB05 showed that the ionization of cold neutral gas by cosmic rays may significantly contribute to temperature fluctuations. They used a two-phase model to explain t^{2} values obtained by Esteban et al. (2002) for a number of H II regions. However, the high temperature difference between the two phases in GB05 model (see their Table 1) suggests that t^{2} values obtained by Esteban et al. cannot be applied directly to two-phase models. Values of derived by GB05 need to be re-considered.
We re-estimate
values for the sample of H II regions of
Esteban et al. (2002). Following GB05, three values of
are
assumed, 100, 1000, and 4000 K. Under these conditions, temperature in the
cold gas is too low to collisionally excite the [O III] lines, and
consequently Eq. (7) can be simplified to
In Table 1 we compare our values to those of GB05 (given in parentheses); for low values of differences of up to a factor of a hundred are found. It can easily be seen that the discrepancies increase with decreasing temperature of the cold gas, as suggested by Fig. 1. For K, our derived values of are very low, suggesting that the values of t^{2} reported by Esteban et al. (2002) can be explained by the existence of a very small amount of cold gas.
Table 1 also gives t^{2} and T_{0} values derived from the Eqs. (12) and (13). As the Table shows, the real t^{2} values are lower when than those derived from the empirical method. As a result, values of the cosmic ray ionization rate, , derived by GB05 have been grossly overestimated (cf. their Eq. (17)). Our conclusion is consistent with the range of values inferred for the Orion nebula from Gamma ray observations.
We have studied the relationship between values of t^{2} predicted by a two-phase model and those derived empirically from observations (empirical method). Our results show that the existence of extremely cold gas within H II regions may lead to overestimated t^{2} calculated from empirically determined T([O III]) and T(H I). We stress that care should be taken when using the two-phase model to study large temperature fluctuations of H II regions. In this model, CELs are hardly produced by the cold gas, which on the other hand makes a large contribution to the flux at the Balmer jump, due to the dependence of . Accordingly, the existence of a very small amount of cold material may lead to a large discrepancy between T([O III]) and T(H I). In other words, in spite of its small mass, the existence of extremely cold material can reproduce apparently large t^{2} (as derived from the empirical method), much larger than the actual value (as defined by Eq. (2)).
Finally, we revisited the GB05 study of cosmic ray ionization as a mechanism for creating temperature fluctuations in H II regions. While this provides a potential mechanism for creating cold ionized plasma in H II regions, we show that the values of required to produce the ionization have been overestimated in their treatment, due to the t^{2} discrepancy discussed above. Based on the formulae presented here, we re-estimated their model parameters. The corrections are apparent, particularly in cases where temperature of the cold gas component is low, resulting in lower values of that agree better with the estimates for the Orion nebula published in the literature.
Acknowledgements
We thank the referee, Dr. C. Morisset, for helpful comments that improved clarity of the paper. We would also like to thank Dr. Morisset for computing new values of t^{2} and T_{0}, based on our estimates of . Those values are now tabulated in Table 1. Y.Z. and X.W.L. acknowledge support by NSFC grant #10325312.