A&A 437, 1023-1027 (2005)
DOI: 10.1051/0004-6361:20052694

Processes of (n - n')-mixing in collisions of Rydberg H*(n) atoms with H(1s) in the Solar atmosphere

A. A. Mihajlov1,2 - Lj. M. Ignjatovic1,2 - M. S. Dimitrijevic2,3


1 - Institute of Physics, PO Box 57, 11001 Belgrade, Serbia and Montenegro
2 - Institute Isaac Newton of Chile, Yugoslavia Branch, and Universidad Diego Portales, Chile
3 - Astronomical Observatory, Volgina 7, 11160 Belgrade 74, Serbia and Montenegro

Received 13 January 2005 / Accepted 17 March 2005

Abstract
The significance of (n-n')-mixing processes in H*(n)+H(1s) collisions, for the principal quantum number $n \ge
4$, in the Solar photosphere and lower chromosphere has been investigated. These processes have been treated by the mechanism of resonant energy exchange within the electron component of the considered collision system. These processes must have significant influence in comparison with corresponding electron-atom collision processes on the populations of hydrogen Rydberg atoms in weakly ionized layers of the Solar atmosphere (ionization degree of the order of 10-4). From the results obtained it follows that the examined (n-n')-mixing processes have to be included in any modelling and investigation of Solar plasma, especially in the region of the temperature minimum in the Solar photosphere.

Key words: atomic processes - molecular processes - Sun: atmosphere

1 Introduction

Mihajlov et al. (1997, 1999, 2003) drew attention to the importance of inelastic processes in slow collisions of Rydberg state atoms with ground state H atoms in solar and stellar atmospheres. These papers were devoted to an investigation of the processes of chemi-ionization

 \begin{displaymath}{\rm H} + {\rm H}^{*}(n) \rightarrow \vec{e} + \left\{
\begi...
...
\displaystyle{ {\rm H} + {\rm H}^{+}},
\end{array} \right.
\end{displaymath} (1)

and the inverse process of chemi-recombination

 \begin{displaymath}\vec{e} + {\rm H}^{+} + {\rm H} \rightarrow \left\{
\begin{a...
...displaystyle{ {\rm H}^{*}(n) + {\rm H} },
\end{array} \right.
\end{displaymath} (2)

where H designates a ground state $H(1\rm s)$ atom, H*(n) a Rydberg state hydrogen atom with principal quantum number $n \ge
4$ and $\vec{e}$ is a free electron. The processes (1) and (2) were treated using the resonant energy exchange mechanism (Smirnov & Mihajlov 1971; Janev & Mihajlov 1979) between the electrons of the ${\rm H+H}^{*}(n)$system, to which we subsequently refer to as the resonant mechanism. It was shown that the processes (1) and (2) can be of great importance in weakly ionized layers of many stellar atmospheres, in comparison with other possible ionization and recombination processes involving electrons or photons, namely direct ionization by electron impact

 \begin{displaymath}
\vec{e}+{\rm H}^{*}(n) \rightarrow {\rm H}^{+}+\vec{e}+\vec{e},
\end{displaymath} (3)

the inverse 3-body electron recombination

 \begin{displaymath}
\vec{e}+\vec{e}+{\rm H}^{+} \rightarrow \vec{e}+{\rm H}^{*}(n),
\end{displaymath} (4)

and the radiative recombination

 \begin{displaymath}
\vec{e}+{\rm H}^{+} \rightarrow {\rm H}^{*}(n)+\epsilon_{\rm ph},
\end{displaymath} (5)

where $\epsilon_{\rm ph}$ is the energy of the emitted photon.

However, Mihajlov et al. (1997, 1999, 2003) adopted the usual assumption that the distribution of Rydberg states was determined by electron collisions

 \begin{displaymath}
\vec{e}+{\rm H}^{*}(n)\rightarrow \vec{e}+{\rm H}^{*}(n').
\end{displaymath} (6)

More recently, Mihajlov et al. (2004) showed that the Rydberg state distribution in a weakly ionized hydrogen plasma could also be strongly influenced by (n-n')-mixing processes in ${\rm H}^{*}(n) + {\rm H}$collisions. These conclusions were based on the results of an investigation of the excitation processes

 \begin{displaymath}{\rm H}^{*}(n) + {\rm H} \rightarrow \left\{
\begin{array}{l...
...}(n'=n+p)}
\end{array} \right., \qquad n\ge 4, \quad p \ge 1,
\end{displaymath} (7)

and the inverse de-excitation process

 \begin{displaymath}{\rm H}^{*}(n) + {\rm H} \rightarrow \left\{
\begin{array}{l...
... {\rm H}^{*}(n'=n-p)}
\end{array},
\right. \qquad n-p \ge 4,
\end{displaymath} (8)

which were caused by the same resonant mechanism as the processes of chemi-ionization and chemi-recombination. The results obtained impose the necessity of an investigation of the processes (7) and (8) in weakly ionized layers of stellar atmospheres and particularly in the Solar atmosphere. The confirmation of the importance of these processes in the Solar photosphere would support the necessity of including of the whole group of non-elastic processes in ${\rm H}^{*}(n) + {\rm H}$ collisions, i.e. Eqs. (1), (2), (7) and (8), in models of weakly ionized layers of Solar and stellar atmospheres. That should be confirmed by one of the known and widely used non-LTE models of the Solar atmosphere.

For the reason mentioned we shall examine here in more detail the significance of (n-n')-mixing processes (7) and (8) in the Solar photosphere and lower chromosphere. For this purpose we shall calculate the rate coefficients of these processes for different n and p using a semiclassical method developed by Mihajlov et al. (2004). The rate coefficients obtained will be compared with the rate coefficients of electron-atom mixing processes (6) for electron and hydrogen atom densities and temperatures from the model C of the Solar photosphere and the lower chromosphere of Vernazza et al. (1981).

2 Theoretical comments

The resonant mechanism. Since the resonant mechanism for processes (1)-(8) has been discussed in several previous papers, and especially has been described in detail in Mihajlov et al. (2004), only the basic facts will be given here. On the basis of the resonant mechanism the ${\rm H}^{*}(n) + {\rm H}$ collision system is treated within the domain

 \begin{displaymath}R \ll r_{n},
\end{displaymath} (9)

where R is the internuclear distance and $r_{n} \sim n^{2}$ is the average radius of the Rydberg atom ${\rm H}^*(n)$. Within this domain the ${\rm H}^{*}(n) + {\rm H}$ system is treated in the form: ${\rm e} +({\rm H}^{+}+{\rm H})$, where e is the outer electron of the ${\rm H}^*(n)$Rydberg atom. The electronic states of the subsystem $({\rm H}^{+}+{\rm H})$are described using the adiabatic electronic ground state $\vert 1{\rm s}
\sigma_{g};R\rangle$, or of the first excited state $\vert 2{\rm p} \sigma_{u};R\rangle$of the molecular ion ${\rm H}_{2}^{+}$. The relative internuclear motion is described in the approximation of two classical trajectories: the first one, which corresponds to the adiabatic $1{\rm s}
\sigma_{g}$-term, and the other, which corresponds to the $2{\rm p}
\sigma_{u}$-term. It is assumed that each of these trajectories is realized with the same probability $p_{\rm in}$ which is equal to 1/2.

The processes (7) and (8) of (n-n')-mixing, as well as the processes of chemi-ionization/recombination (1) and (2), are treated as a result of the resonant energy exchange between the outer electron e and the electronic component of the ${\rm H}^{+}+{\rm H}$ subsystem. This means that the transition of the outer electron from the initial energetic state to the upper ones occurs simultaneously with the transition of the ${\rm H}^{+}+{\rm H}$ subsystem from the electronic excited state $\vert 2{\rm p} \sigma_{u};R\rangle$ to the ground state $\vert 1{\rm s}
\sigma_{g};R\rangle$, and the transition of the outer electron to the lower energetic state occurs simultaneously with the transition of the $\rm H^{+}+H$subsystem from the ground electronic state $\vert 1{\rm s}
\sigma_{g};R\rangle$ to the excited state $\vert 2{\rm p} \sigma_{u};R\rangle$ as illustrated in Fig. 1. Within the considered mechanism one assumes that all mentioned transitions are caused by the interaction of the outer electron with the dipole momentum of the ${\rm H}^{+}+{\rm H}$subsystem.

  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{2694fig1.eps}
\end{figure} Figure 1: The resonant mechanism in chemi-ionization/recombination channels (dashed arrows), and (n-n')-mixing channels (full arrows).
Open with DEXTER

In Janev & Mihajlov (1979, 1980) the resonant mechanism described was applied in the range n > 4 where there are no crossings of the terms of H $^{*}(n)+{\rm H}(1\rm s)$ systems with the term of ${\rm H}^{+}+{\rm H}^{-}$ system. However, following our later papers, we will take n=4 as a lower boundary of the considered principal quantum number region. In accordance with Sidis et al. (1983) the corresponding crossing point is placed at such a large internuclear distance that it does not influence the considered processes.

The resonant mechanism was also treated within the framework of theories related to a wider class of collision atomic systems A*(n)+B, where in the general case A and B denote different atoms (Flannery 1980; Lebedev 1991b; Lebedev & Fabrikant 1996). However, in the symmetrical case these theories give the same results as the theory presented in Janev & Mihajlov (1979, 1980) and Mihajlov & Janev (1981). Besides, the processes of (n-n')-mixing in the Li *(n)  + Li collisions were experimentally studied in Dubreuil (1983). The experimental results obtained were compared with the results from Janev & Mihajlov (1979). It was demonstrated that the agreement of these results is good enough to confirm the validity of the described resonant mechanism.

The processes (1)-(8), caused by the resonant mechanism, should not be confused with the processes in H*(n)+A collisions caused by the direct interaction of the electron of atom H* with atom A whose electronic state does not change during the collision (Matsuzawa 1974; Gerstein 1976; Olson 1977; Percival 1978; Lebedev 1991a; Lebedev & Fabrikant 1996). The influence of such processes can be neglected in comparison with the influence of processes (7) and (8) at thermal collision energies ($\la$$ 1~\rm eV$), as was discussed in Janev & Mihajlov (1979).

  \begin{figure}
\par\includegraphics[width=7cm,clip]{2694fig2.eps}
\end{figure} Figure 2: The partial "smearing'' of the discrete Coulomb spectrum to a continuum.
Open with DEXTER

The rate coefficients of processes (7) and (8) for fixed n and $n'=n~ \pm~ p$ have been determined semi-classically, describing the internuclear motion by using the trajectories determined in the potential corresponding to the excited electronic state $\vert 2{\rm p} \sigma_{u};R\rangle$ of the ${\rm H}^{+}+{\rm H}$ subsystem. For the calculations of the considered rate coefficient an approximation was used where a block of Rydberg states with the principal quantum numbers from n + p1 to n + p2, where $p_2 \ge p_1 \ge 1$, is "smeared'' to a continuum within the effective principal quantum number range from $n + p_1
- \Delta p$ to $n + p_2 + 1 - \Delta p$, where $0 < \Delta p < 1$, as is shown in Fig. 2. Such "smearing'' is performed under the conditions described in Mihajlov et al. (2004).

The described procedure for determination of the rate coefficients of the examined (n-n')-mixing processes means that the hydrogen (nl)-states with fixed n are populated according to their statistical weights. This assumption is based on the properties of the existing models of the photosphere and lower chromosphere of the Sun and similar stars, of some cooler stars (as for example M and K type stars), and of some other stars (like some of DB white dwarfs for example). It is well known that these models assume only limited deviations from the LTE. For example in Vernazza's model C of the Solar atmosphere used here it is assumed that the distribution of the total populations of the atomic shells with different principal quantum number n can deviate from the corresponding Boltzman distribution, but that there are no such deviations within particular shells. There are physical reasons for such an approach. Namely, it is well known that the processes of (l-l')-mixing for fixed n are practically resonant and are very effective over the distances proportional to an average Rydberg atom radius (Olson 1977; Percival 1978; Lebedev 1991a). Consequently such processes are much faster that the (n-n')-mixing processes, especially within the region where the (n-n') transition energy is of the order of kT. A similar assumption for the population distribution of (nl) - states with fixed n was used and discussed in many previous articles, including Janev & Mihajlov (1979, 1980).

The procedure of "smearing'' a discrete spectrum has already been used. For example in Sobel'man (1979) this has been used for description of the electron - ion recombination with transition of the free electron in a block of Rydberg states. In Duman & Shmatov (1980), "smearing'' the discrete spectrum of the molecular ion A2+ has been used to describe the chemi - ionization processes in A*(n) + A collisions. In these cases it was used as physically acceptable approximation without special justifications.

In our case however, well-determined physical reasons exist for the approximation of "smearing''. Concerning the "smearing'' of the hydrogen Rydberg spectrum, some results from Janev & Mihajlov (1979) have been taken into account. These are the results related to the behavior of the (n - n') transition probability due to $H^*(n) + H(1\rm s)$ collisions as a function of collision parameters for fixed impact energies. On the basis of the results obtained one can consider that the (n - n') transition occurs within some interval of internuclear distances whose center is a point of resonance (where the observed transition may be treated as a result of the strictly resonant energy exchange within the electronic component of the considered system) and that such intervals overlap significantly for n' and $n'~\pm~ 1$. These results are related to the case when the atom - atom impact energy is at least several times greater than energies for the considered (n-n') transitions, which justified the use of the straight line trajectory approximation in the mentioned article.

This allows the considered process to be treated as the continuous decay of the initial electronic state (of the total system), during the internuclear motion within the reaction zone (the internuclear distance range being less or close to the resonant distance for the transition $n \rightarrow n+1$).

The rate coefficients. Since the procedure of obtaining the rate coefficients for the processes (7) and (8) is described in detail in Mihajlov et al. (2004), only the corresponding final expressions are presented here. The rate coefficients Kn;n+p(T) for the excitation processes (7), for given n, p and the temperature T, are obtained in the form

 \begin{displaymath}\begin{array}{l}
{\displaystyle K_{n;n+p}(T)=\frac{2
\pi}{3\s...
...}{kT}\right]\frac{R^{4}\cdot
{\rm d}R}{a_{0}^{5}},}
\end{array}\end{displaymath} (10)

where e and a0 are the electron charge and the atomic unit of length, gn;n+p is the Gaunt factor defined in Johnson (1972), $R_{\rm min}(n,n+p)$ and $R_{\rm max}(n,n+p)$ are roots of equations

 \begin{displaymath}\begin{array}{l}
{\displaystyle \epsilon(R_{\rm min})=\epsilo...
...R_{\rm max})=\epsilon_{n+p-\Delta p}-\epsilon_{n},}
\end{array}\end{displaymath} (11)

with $\epsilon(R)=U_{2}(R)-U_{1}(R)$ and $\Delta p = 0.380$. Here, U1(R) and U2(R) are adiabatic terms corresponding to the electronic states $\vert 1{\rm s}
\sigma_{g};R\rangle$ and $\vert 2{\rm p} \sigma_{u};R\rangle$, respectively, and ${\epsilon }_n = -Ry/n^2$. X(R) denotes the function

 \begin{displaymath}X(R)=\frac{\Gamma\left(3/2;\frac{\vert U_{1}(R)\vert}{kT}\right)}{\Gamma(3/2)},
\end{displaymath} (12)

where $\Gamma(3/2;x)$ and $\Gamma(3/2)$ are the corresponding values of the incomplete and complete Gamma functions. The factor $p_{\rm in}=1/2$ in Eq. (10) is the mentioned above probability for the realization of a trajectory which describes the internuclear motion in the potential U2(R), while the factor 2 reflects the fact that the system ${\rm H}^{*}(n) + {\rm H}$ crosses the transition zone twice: once on entry and once on leaving.

The rate coefficients Kn;n-p(T) for de-excitation processes (8) are obtained from the thermodynamic balance principle in the form

 \begin{displaymath}K_{n;n-p}(T)=K_{n-p;n}(T)\cdot \frac{(n-p)^{2}}{n^{2}}\cdot
\exp{\left(\frac{\epsilon_{n-p;n}}{kT}\right)},
\end{displaymath} (13)

where $\epsilon_{n-p;n} = \epsilon_n - \epsilon_{n-p}$.

3 Results and discussion

As in the previous article related to the chemi-ionization/recombination processes (1) and (2), we will treat the Solar photosphere and lower chromosphere using the standard C model of Vernazza et al. (1981). For our calculations we need to take from this model the electron density $N_{\rm e}$, the density of hydrogen atoms in the ground state (n = 1), denoted here by N(1), and the temperature T. The behavior of $N_{\rm e}$ and N(1) in the Solar atmosphere as a function of height (h) is shown in Fig. 3, while Fig. 4 shows the behavior of T and the parameter $\eta $ defined by the relation

  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{2694fig3.eps}
\end{figure} Figure 3: The behavior of $N_{\rm e}$ and N(1) in [cm-3] as functions of height h within the considered part of the Solar atmosphere from model C of Vernazza et al. (1981)
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{2694fig4.eps}
\end{figure} Figure 4: The behavior of the temperature Tand the parameter  $\eta =N(1)/N_{\rm e}$ as functions of height hwithin the considered part of the Solar atmosphere from model C of Vernazza et al. (1981)
Open with DEXTER


 \begin{displaymath}\eta_{\rm ea}=\frac{N(1)}{N_{\rm e}} \cdot
\end{displaymath} (14)

In our case the parameter $\eta $ is close to the inverse value of the considered plasma's degree of ionization. Also, we will use the data about the hydrogen atom excited state population ($n \ge
2$) from Vernazza et al. (1981). In order to determine the relative influence of processes (7) and (8) in the considered parts of Solar atmosphere, we will compare them with electron-atom collision processes (6). Following the mentioned above papers related to chemi-ionization/recombination processes, we compare the studied collision processes only with other relevant collision processes, omitting the radiative decay of hydrogen Rydberg states. The radiative decay, which causes the deviation from LTE in the considered part of the Solar atmosphere, is already included in Vernazza's model, so that the "optical part'' of this model does not need to undergo any changes. We show that the "collision part'' of Vernazza's model has to be changed by adding atomic collision processes (1)-(8).

The relative efficiency of processes (7) and (8) in comparison to process (6) is characterized by the parameter $F_{n;n\pm p}(T)$

 \begin{displaymath}\displaystyle{F_{n;n \pm p}(T)=\frac{K_{n;n \pm p}(T)N(n)N(1)...
...{\frac{K_{n;n \pm
p}(T)\eta_{\rm ea}}{\alpha_{n;n \pm p}(T)},}
\end{displaymath} (15)

where N(n) is the excited atom states population for given n, the parameter $\eta_{\rm ea}$ is defined by Eq. (14), and $\alpha_{n;n \pm p}(T_{\rm e})$ is the rate coefficient for the electron-atom process (6) taken from Vriens & Smeets (1980).

We have calculated the values of the parameter $F_{n;n\pm p}$ for $n \ge
4$ and $p \ge 1$, within the range of h corresponding to the Vernazza et al. (1981) model C of the Solar photosphere. Our results are illustrated in Figs. 5-7. In these figures the behavior of the parameter Fn;n+p is shown, for $4 \le n \le 8$ and $1 \le p \le 5$, as a function of h within the range $ -50 ~{\rm km} \le h \le 1000$ km. In the largest part of the photosphere the (n-n')-mixing processes (7) and (8) are totally dominant in comparison to the electron-atom processes (6) for n=4, 5, 6 and practically any $p \ge 1$, and have similar intensities for n=7, 8. Comparing Figs. 5-7 with Fig. 4 one can see that the region of the maximal $F_{n;n\pm p}$ values corresponds to an h region near the temperature minimum, where the parameter $\eta $, defined by Eq. (14), is also maximal.

  \begin{figure}
\par\includegraphics[width=7.2cm,clip]{2694fig5.eps}
\end{figure} Figure 5: The behavior of the parameter $F_{n,n+p}=(K_{n,n+p}\cdot \eta )/\alpha _{n,n+p}$, characterizing the relative efficiency of atom-atom and electron-atom (n-n')- mixing processes, as a function of height h for n=4 and p=1-5. The atom-atom rate coefficient Kn,n+p is given by Eqs. (10)-(12) and the parameter $\eta $- by Eq. (14). The electron-atom rate coefficient  $\alpha _{n,n+p}$ is taken from Vriens & Smeets (1980).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.2cm,clip]{2694fig6.eps}
\end{figure} Figure 6: The same as in Fig. 5, but for n=5.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.8cm,clip]{2694fig7.eps}
\end{figure} Figure 7: The same as in Fig. 5, but for n=6,7,8.
Open with DEXTER

To demonstrate the significance of the obtained results we introduce the quantities $n=n_{\rm min}^{(V)}$ and $n=n_{\rm min}^{(\rm eq)}$dependent on h. The first value of n is defined by the condition:

 \begin{displaymath}N^{(V)}\left(n=n^{(V)}_{\rm min}\right) = \min_{n\ge 1} N^{(V)}(n),
\end{displaymath} (16)

where N(V)(n) are excited atom populations from model C for a given h, and the second one is the value of n defined by the condition:

\begin{displaymath}N^{\rm (eq)}(n=n^{\rm (eq)}_{\rm min};T) = \min_{n\ge 1} N^{(\rm eq)}(n;T)
\end{displaymath} (17)

where $N^{\rm (eq)}(n;T)$ are equilibrium (Boltzmanian) excited atom populations for T from model C for the same h as in the previous case. It could be shown that in the considered temperature range, $n^{\rm (eq)}_{\rm (min)}$ is close to the value of (Ry/kT)(1/2). In Table 1 the values of $n=n_{\rm min}^{(V)}$ and $n=n_{\rm min}^{\rm (eq)}$ are presented for a series of h values in -50 km to 1065 km range. Within the considered domain of h the values of $n=n_{\rm min}^{(V)}$ and $n=n_{\rm min}^{(\rm eq)}$ are 4, 5 and 6 for which domination of the processes (7) and (8) in comparison with the processes (6) is demonstrated.

Table 1: The values of $n^{\rm (eq)}_{\rm min}$ and $n^{(V)}_{\rm min}$ for h and T from Vernazza et al. (1981) model C.

For plasmas similar to the considered photospheric one, the exchange rate between both hydrogen basic components, atoms $\rm H(1s)$ and ions H+, is limited by the group of excited atom states with n close to the $n^{(V)}_{\rm min}$ or $n^{(\rm eq)}_{\rm min}$. However, one can see that the processes (7) and (8) play the dominant role in this group of excited states. This means that the exchange between the excited state atom populations within the lower part of the Rydberg region of n is determined by the processes (7) and (8). If we add to this the fact that the exchange between the mentioned group of Rydberg states and the continuum is determined by the chemi-ionization/recombination processes (1) and (2), one can conclude that the processes (1), (2), (7) and (8), closely connected by the resonant mechanism, completely dominate for the excited states with n close to $n^{(V)}_{\rm min}$ or $n^{\rm (eq)}_{\rm min}$.

We gave here all relevant results demonstrating the significance of the processes (7) and (8) for a large part of the Solar photosphere and lower chromosphere, particularly in the neighborhood of the temperature minimum. On the basis of presented results we can conclude that the influence of the (n-n')-mixing processes (7) and (8), as well as the chemi-ionization/recombination processes (1) and (2), must be taken into account for any modelling of the Solar photosphere and lower chromosphere. Our findings are also related to similar stellar atmospheres as well as to stellar atmospheres with lower effective temperatures but where the atomic component is still dominant over the molecular one.

Acknowledgements

The authors wish to thank Prof. E. V. Kononovich for useful suggestions regarding the Solar atmosphere. This work is a part of the "Radiation and transport properties of the non-ideal laboratory and ionospheric plasma'' project (Project number 1466), and "Influence of collision processes on the spectra of astrophysical plasmas'' (Project number GA 1195) and was supported by the Ministry of Science, Technology and Development of the Republic of Serbia.

References

 

Copyright ESO 2005