A&A 403, 787791 (2003)
DOI: 10.1051/00046361:20030463
A. A. Mihajlov^{1, 2}  D. Jevremovic ^{3, 4}  P. Hauschildt^{5, 6}  M. S. Dimitrijevic ^{2, 4}  Lj. M. Ignjatovic ^{1, 2}  F. Alard^{7}
1  Institute of Physics, PO Box 57, 11001 Belgrade,
Serbia and Montenegro
2  Institute Isaac Newton of Chile, Yugoslavia Branch, Belgrade,
Serbia and Montenegro
3 
APS Division, Phys. Depart., Queen's University Belfast,
BT7 1NN, Northern Ireland, UK
4 
Astronomical Observatory, Volgina 7, 11160 Belgrade, Serbia
and Montenegro
5 
Depart. of Phys. and Astronomy and Center for Simulat. Physics,
University of Georgia, Athens, GA 30602, USA
6 
Hamburger Sternwarte, Universitaet Hamburg, Gojenbergsweg 112, 21029
Hamburg, Germany
7 
CRAL/ENS, 46 Allée d'Italie, 69634 Lyon Cedex 07, France
Received 22 July 2002 / Accepted 4 February 2003
Abstract
We study the influence of a group of chemiionization and chemirecombination
processes on the populations of higher states of hydrogen in the
layers of a stellar atmosphere.
The group of processes includes
ionization:
H^{*}(n) + H(1s)
H_{2 +} + e ,
H^{*}(n) + H(1s)
H(1s) + H^{+} + e,
and inverse recombination:
H_{2 +} + e
H^{*}(n) + H(1s),
H(1s) + H^{+} + e
H^{*}(n) + H(1s),
where
H^{*}(n) is the hydrogen atom in a state with the principal quantum
number n
1, and
H_{2 +} is the hydrogen molecular ion in a weakly
bound rhovibrational state of the ground state. These processes have
been treated within the framework
of the semiclassical approximation, developed in several previous
papers, and have been included in the general stellar atmosphere
code PHOENIX. We present results for an M dwarf
atmosphere with
T_{eff} = 3800 K and find that the
inclusion of chemiionization and chemirecombination processes is
significant in the low temperature parts of the
atmosphere.
Key words: atomic processes  molecular processes  stars: late type
For several years we have studied the chemiionization atomatom collisional processes and inverse chemirecombination processes that may influence the populations of the hydrogen atom Rydberg states in the weakly ionized layers of stellar atmospheres (Mihajlov et al. 1996a,1997). In the case of hydrogen, these processes are:


where H^{*}(n) is a hydrogen atom in a state with the principal quantum number n 1, and H_{2 +} is a hydrogen molecular ion in a weakly bound rhovibrational state of the ground electronic state. In the case of ionization (1a,b), the lower limit of the Rydberg state range that was considered was set to n = 4 (Mihajlov et al. 1996a,1997). The block of rhovibrational states with a dissociation energy less than 1 eV has an essential role in the case of formation of the molecular ion (processes (1a) and (2a)) and within the range of quantum numbers that was considered. In the case of cool stellar atmosphere layers (3000 K T 5000 K), the block of such weakly bound rhovibrational states can be considered to be in equilibrium with the collisional H(1s) + H^{+} complexes, without the requirement that the lower rhovibrational states of molecular ion are in equilibrium (Mihajlov & Ljepojevic 1982; Mihajlov et al. 1992). This fact makes the description of recombination processes (2a) and (2b) easier, since it enables us to avoid the direct determination of the density of molecular ions in weakly ionized rhovibrational states.
In our previous work (Mihajlov et al. 1997, 1998) we estimated the influence of the processes (1a,b; 2a,b) for the low temperature layers of the Solar atmosphere. We compared the ionization and recombination fluxes with and without inclusion of the above processes and found a significant difference in the population of the H^{*}(n) atomic states with 4 < n < 8.
Here we discuss the importance of these processes in stellar atmospheres of lower temperature. This paper is organized as follows: in Sect. 2 we summarize the theoretical framework and in Sect. 3 we describe the changes in the stellar atmosphere code PHOENIX and parameters of the test model atmosphere. We present our results and discussion in Sect. 4 and draw our conclusions in Sect. 5.
Let I_{i}^{(a, b)}(n) and I_{r}^{(a, b)}(n) denote ionization and recombination fluxes affecting population and depopulation of excited hydrogen atom atomic levels due to reactions. By definition, ionization fluxes I_{i}^{(a, b)}(n) determining the rate of depopulation of H^{*}(n) states are given by

where N(H) and N(H^{*}(n)) are H(1s) and H^{*}(n) atom densities, and K_{i}^{a}(n) and K_{i}^{b}(n) are rate coefficients of the reactions (1a) and (1b) respectively. The recombination flux I_{r}^{(a)}(n), determining the velocity of population of H^{*}(n) states due to the reaction (2a), is given by definition as: I_{r}^{(a)} = K_{dr}^{ . }N(H_{2 +})^{ . }N(e), where N(e) and N(H_{2 +}) are the electron density and the density of molecular ion H_{2 +} in the abovementioned weakly bound rhovibrational states, and K_{dr} is the corresponding rate coefficient of the dissociative recombination reaction (2a). The expression for the recombination flux I_{r}^{(a)} may be written by:

where K_{r}^{(a)} = K_{dr}[N(H)^{ . }N(H^{+})/N(H_{2 +})]^{1}. The advantage of expression (4a) for the recombination flux I_{r}^{(a)}(n) is the fact that it introduces a new effective rate coefficient of the reaction (2a). If the ionization rate coefficient K_{i}^{a}(n) is known, the K_{r}^{(a)}(n)may be determined directly, under assumption of the equilibrium between H_{2 +} weakly bound rhovibrational states and the H + H^{+} collisional complex (Mihajlov & Ljepojevic 1982; Mihajlov et al. 1992). This makes the treatment of the problem easier and enables us to avoid the determination of the K_{dr}(n) rate coefficient. The recombination flux I_{r}^{(b)}(n), determining the rate of population of H^{*}(n) states due to the reaction (2b), is given, by definition, by:
I_{r}^{(b)}(n) = K_{r}^{(b)}(n)^{ . }N(H)^{ . }N(H^{+})^{ . }N(e),  (4b) 
where K_{r}^{(b)}(n) is the rate coefficient of the reaction (2b). The influence of the processes (1a,b) and (2a,b) on H^{*}(n) atomic states populations may be estimated using only the N(H), N(H^{*}(n)), N(H^{+}) and N(e) populations.
Rate coefficients for the processes (1, 2) have been determined within the semiclassical theory described in the previous papers (Janev & Mihajlov 1980; Mihajlov & Ljepojevic 1982; Mihajlov et al. 1992). This theory takes into account atomic collisions, where the H^{*}(n) + H system in the processes (1a,b) as well as the e + H^{+} + H system in the processes (2b) go through the phase where they may be treated as a (H^{+} + H) + e system. Here, (H^{+} + H)denotes that the ionatom subsystem may be treated as a quasimolecular ionatom complex under assumption that its electronic states may be described as superpositions of the ground and the first excited states of the H_{2 +} molecular ion. Furthemore, this means that during the quasimolecular phase the average electron distance from the center of mass of the H^{+} + H subsystem is much larger than the internuclear distance between H^{+} and H. Consequently, the dipole interaction between the (H^{+} + H) and an electron (Smirnov & Mihajlov 1971) is the origin of the chemiionization/recombination processes (1a,b) and (2a,b). One can notice that within the frame of this approach, the processes (2a) and (2b) may be described practically in the same way. Expressions for the rate coefficients in the processes (1a,b) and (2a,b) which are used here as well as the method of determination of the rate coefficients are described in details in Mihajlov et al. (1996b ).
We included procesess (1a,b) and (2a,b) into the general stellar atmosphere code PHOENIX described in Hauschildt et al. (1999), Short et al. (1999), Baron & Hauschildt (1998). PHOENIX solves simultaneously the equations of hydrostatic equilibrium, multilevel and multispecies statistical equilibrium and NLTE radiative transfer. The major advantage of PHOENIX is that it uses extensive atomic and molecular databases and can be used for modeling of very cool objects.
The main change in the code is in the addition of a new term in the statistical equilibrium equation which describes additional collisional processes. This change allows us to treat chemiionization/recombination processes consistently with all the other processes. Due to limits of the semiclassical approach, we include additional collisional term only for levels larger than n = 4.
We use for our calculation the planeparallel, static option of PHOENIX which has been adopted for use in "chromospheric'' type of problems. As a test atmosphere we use a "schematic'' chromospheric structure with an effective temperature of 3800 K, logarithm gravity of 5.0, log m_{0} = 2.5 being the column mass at temperature minimum, and log m_{TR} 4.8, m_{TR} standing for the column mass in the transition region. The choosen parameters represent a typical M0 dwarf of weak activity. The procedure of building this model is identical to that described in Jevremovic et al. (2000) and its structure is shown in Fig. 1. As we are interested in behavior of hydrogen populations, we solve the NLTE problem for the 30 lowest levels. The standard collisional data for hydrogen are those from Vriens & Smeets (1980).
In Figs. 24 we present the behavior of population ratios without and with inclusion of chemiionization (1a,b) and chemirecombination (2a,b) processes for the levels 3 to 20. The ratios are close to unity for n = 3, and start to differ for n = 4 and 5 (Fig. 2). For n = 5 the largest difference is around 10%. This is to be expected as we include the procesess for n4, which is based on the validity of the semiclassical approximation. The population of lower levels (n = 2, 3) is influenced by other processes included in the statistical equilibrium equation.
As expected, the additional collisional term in statistical equilibrium brings the solution closer to local thermodynamical equilibrium (LTE). In terms of our population ratios, in the parts of atmosphere where populations, determined in our calculations using PHOENIX, are larger than LTE populations (departure coefficients b_{i} > 1), our parameter is larger than one. This means that in such a case, chemiionization processes depopulate excited states and dominate in comparison with chemirecombination processes. The result is that populations determined with processes (1, 2) included are smaller in comparison with populations calculated without them and consequently is larger than one. In the reverse case (b_{i} < 1), parameter is smaller than one, which means that in this case chemirecombination processes dominate over chemiionization ones and increase the population of excited levels. As a result, in both cases, processes (1, 2) act as important factors closing the difference between the calculated and LTE populations.
As can be seen from Fig. 2, the behaviour of (n) for 6 < n < 9 shows a tendence to converge to a pattern, which may differ up to 40% from unity. For higher levels (Figs. 3 and 4) the convergence to a pattern is complete and it can be understood by assuming that highly excited hydrogen level populations have quasiBoltzman distributions different from distributions for partial LTE. It is not suprising that the population ratios for the levels up to 30 and the proton density show the same behaviour as levels with n > 15. One also has to remember that in the layers with temperatures lower than 4500 K the main electron contributors are metals so, as expected, there is very little change in electron density.
We have demonstrated the importance of the inclusion of chemiionization (1a,b) and chemirecombination (2a,b) processes in modeling of atmospheres of late type stars. In previous work (Mihajlov et al. 1996a, 1997, 1998), it had been assumed that these processes were important up to n = 10. Here we have demonstrated that inclusion of processes (1a,b) and (2a,b) significantly changes the population of all Rydberg levels due to PLTE. However, in some atmospheric layers selfconsistent modeling shows differences of up to 40% between calculations with and without inclusion of these processes.
We plan to continue these investigations for atmospheres with different effective temperatures, including stars with higher T_{eff}, such as the Sun and Sunlike stars. Namely, on the basis of our previous results and those presented here, one could expect that for such stars, processes (1a,b) and (2a,b) should influence not only the population of higher levels of hydrogen, but also the electron density. Furthermore, we plan to include in PHOENIX the chemiionization/recombination processes similar to (1a,b) and (2a,b) in the He(1s^{2}) atom which would be used in modeling of heliumrich stars.
Also, there is a need for a developement of a theory of processes similar to processes (1a,b) and (2a,b) for the lower levels of hydrogen, with the principal quantum numbers n = 2 and 3. This extension would enable more sophisticated modeling of stellar atmospheres as the population of lower levels may significantly alter the electron density and consequently the atmospheric structure.
Acknowledgements
This work is a part of the projects "Radiation and transport properties of the nonideal laboratory and ionospheric plasma'' (Project number 1466) and "Influence of collisional processes on astrophysical plasma lineshapes" (Project number GA 1195) and was supported by the Ministry of Science, Technologies and Development of Serbia. D.J. was supported at QUB through the PPARC grant PPA/G/S/1998/00034. Authors acknowledge the care of reviewer to improve this paper.