A&A 472, L55-L58 (2007)
DOI: 10.1051/0004-6361:20078170
LETTER TO THE EDITOR
J. Stepán^{1,2} - J. Kasparová^{1} - M. Karlický^{1} - P. Heinzel^{1}
1 - Astronomical Institute, Academy of Sciences of the Czech Republic, v.v.i.,
Fricova 298, 251 65 Ondrejov, Czech Republic
2 - LERMA, Observatoire de Paris-Meudon, CNRS UMR 8112, 5
place Jules Janssen, 92195 Meudon Cedex, France
Received 27 June 2007 / Accepted 30 July 2007
Abstract
Aims. We investigate the effect of the electric return currents in solar flares on the profiles of hydrogen Balmer lines. We consider the monoenergetic approximation for the primary beam and runaway model of the neutralizing return current.
Methods. Propagation of the 10 keV electron beam from a coronal reconnection site is considered for the semiempirical chromosphere model F1. We estimate the local number density of return current using two approximations for beam energy fluxes between 4
10^{11} and 1
.
Inelastic collisions of beam and return-current electrons with hydrogen are included according to their energy distributions, and the hydrogen Balmer line intensities are computed using an NLTE radiative transfer approach.
Results. In comparison to traditional NLTE models of solar flares that neglect the return-current effects, we found a significant increase emission in the Balmer line cores due to nonthermal excitation by return current. Contrary to the model without return current, the line shapes are sensitive to a beam flux. It is the result of variation in the return-current energy that is close to the hydrogen excitation thresholds and the density of return-current electrons.
Key words: Sun: flares - plasmas - line: formation - atomic processes
The ongoing study of nonthermal excitation of the flaring chromospheric plasmas has been mainly concentrated on the effect of particle beams coming from the coronal reconnection site (Canfield et al. 1984; Fang et al. 1993; Kasparová & Heinzel 2002; Stepán et al. 2007; Hawley & Fisher 1994, and references therein). However, Karlický & Hénoux (2002) and Karlický et al. (2004) recently suggested that the role of neutralizing return currents can be as important as the role of the primary beam itself, both for intensity and linear polarization profiles. Karlický et al. (2004) proposed a simple model of return current formed by the runaway electrons and compared the rates of atomic transitions due to collisions both with the thermal electrons and with the electrons of the primary beam, and due to collisions with the return current formed by the runaway electrons. They showed that the rates due to the return current would dominate the collisional processes in the atmospheric region of Balmer line formation. However, no calculations of theoretical spectral line profiles were presented.
The aim of this paper is to take a first step towards self-consistent modeling of the Balmer line formation with return-current effects taken into account. We use a semi-empirical model of the flaring atmosphere as a basis for our NLTE radiative transfer model. Then we use a standard model for electron-beam deceleration due to Coulomb collisions with the ambient atmosphere and combine it with the two different physical models of the return-current generation. We incorporate the relevant processes that enter the atomic statistical equilibrium equations and solve them with the non-local equations of radiation transfer. At the end, we discuss the results and validity of our models.
We assume an electron beam that is accelerated in a coronal reconnection site and injected into the cold chromosphere along the magnetic field lines. During its propagation, the beam evolves under the influence of several processes (Karlický 1997): (a) the beam generates the return current that decelerates the beam in the return-current electric field, (b) the beam generates the plasma waves causing the quasi-linear relaxation of the beam, and (c) the beam electrons are decelerated and scattered due to collisions with the background plasma particles. In the following model, we neglect the plasma wave processes and the return current is taken in the form of runaway electrons (Rowland & Vlahos 1985; van den Oord 1990; Norman & Smith 1978). In this form, the return-current losses are strongly reduced (Rowland & Vlahos 1985; Karlický et al. 2004). Thus, only collisional losses, as described by Emslie (1978), decelerate the electron beam in our case.
Let
be the particle flux of the monoenergetic beam of the
energy
,
where
is the density of the beam
electrons,
their velocity, and
the mass of the electron. According
to Norman & Smith (1978) and Karlický et al. (2004), a fraction of background electrons
forms the current that moves in the opposite direction in
order to neutralize the electric current
associated with the primary
beam. We use
for the number density of the return-current electrons and
for the number density of the background electrons. The neutralization
condition can be expressed as
Figure 1: From upper panel: H, H, and Hdisk-center line profiles for (dotted line, black), 10^{11} (thick solid line, blue), 6 10^{11} (thick dash-dotted line, green), 8 10^{11} (thin solid line, red), and 1.0 10^{12} (thin dash-dotted line, violet) . No return-current effects are taken into account. Note that nonthermal profiles almost overlap in this range of beam fluxes. | |
Open with DEXTER |
Figure 2: Same as Fig. 1 plus including collisions with return-current electrons. The value of is given by Eq. (2). | |
Open with DEXTER |
We used two models for estimating .
First, following
Norman & Smith (1978), the number of runaway electrons can be estimated as
The normalized energy distribution of electrons can be locally expressed in
the form
To take the effect of the beam/return current into account, we have to
calculate all the electron-hydrogen excitation rates,
the rates of ionization by electron impacts, and the rates of the
inverse processes. We use the
data for total collisional cross-sections for the bound-bound and bound-free
transitions by Janev & Smith (1993), retrieved through the
GENIE database (http://www-amdis.iaea.org/GENIE/).
We do not consider any atomic polarization or angular dependence of the
collisional processes in this work.
The excitation or deexcitation rate of the transition
between the two shells can be calculated using the common formula
(5) |
(6) |
Figure 3: Same as Fig. 1 plus including collisions with return-current electrons. The value of is constant along whole beam trajectory. | |
Open with DEXTER |
The rates of the inverse process of ionization by an electron impact,
the three-body recombination
,
must be treated separately
due to the nonthermal nature of the problem. Using the arguments of detailed
balance (Jefferies 1968; Fowler 1955) in the quasi-classical approximation of
the electron-hydrogen collisions, one can derive a formula for the recombination
rate. Let us consider the ionization of level n with the ionization energy E_{n} by an electron with energy E_{i}. Once cross-section
of
the encounter after which we find the two electrons with energies E_{a} and
E_{b}=E_{i}-E_{n}-E_{a} is known, one can write the total recombination rate as
We calculated the NLTE radiative transfer for a 5-level plus continuum hydrogen using the semiempirical 1D plane-parallel flare model F1 (Machado et al. 1980) in which the temperature structure was kept fixed; in this way, we found the differential effects on the H, H, and H lines (cf. Kasparová & Heinzel 2002). We used the preconditioned equations of statistical equilibrium (Rybicki & Hummer 1991) and solved the coupled system of NLTE equations by the accelerated lambda iteration (ALI) method. For further details, see Heinzel (1995). We found the equilibrium state for several beam fluxes with or without the return current. The initial beam fluxes chosen in our calculations were: 10^{11}, 6 10^{11}, 8 10^{11}, and 1 . If the fluxes were lower, the number of runaway electrons would decrease fast and their energy would exceed the beam energy in most depths. This gives a limit for using of this simple model. Higher fluxes, on the other hand, would be unrealistic. The initial energy of beam electrons was set to E_{0}=10 keV.
In Fig. 1, there are first three Balmer line profiles that result from the nonthermal bombardment by the primary beam. The effects of return current were completely ignored. In this sense, these calculations are similar to the ones of Fang et al. (1993) and Kasparová & Heinzel (2002). In spite of their probably limited physical relevance, these profiles are useful for demonstrating the effects of return currents in the more appropriate models that follow. Figure 2 shows the situation where is calculated using Eq. (2); i.e., the relative number of runaway return-current electrons is calculated at each depth in the atmosphere. Finally, in Fig. 3, there are profiles for the model with constant along the atmosphere. In the layers of Balmer line formation, remains approximately constant and its mean values are shown in Table 1. Comparing Fig. 1 with Figs. 2 and 3, one can see that the effect of return current is very significant: All three lines show a prominent increase emission in the line center.
Table 1: The properties of the return currents. stands for the initial flux of the 10 keV beam, is the mean relative number of runaway return-current electrons in the Balmer lines formation region, and stands for the typical energy of the return-current electrons in these layers. Both and are average quantities which can roughly characterize the return-current properties in the region of interest.
In the region of Balmer line formation, the energy of the return current can be close to the excitation threshold of levels of hydrogen, and it remains approximately constant along an extended trajectory. The beam is finally stopped on a very short path. The overall path of the beam is, however, sensitive to the initial energy of the beam. In order to model the beam propagation and line formation accurately, one has to interpolate the original F1 model of Machado et al. (1980) by a number of grid points in the layers of the Balmer line formation. Since the return-current energy and density are sensitive to the beam flux (see Table 1), the resulting variation of the nonthermal collisional rates leads to a significant variation in line profiles. For both models under consideration, a maximum emission is found for the beam flux of 6 , although the resulting profiles from these models differ slightly from each other. In contrast to the case of 4 , for which we found the least emission among the studied flux intervals, the return-current density is higher by a factor of 5. It leads to a significant increase in nonthermal excitation rates. The disagreement of the profiles at fluxes below 6 is the result of the significant dependency of on the beam flux and atmospheric depth. The values of at low fluxes (shown in the Table 1) are only a rough approximation for the Balmer line formation layers, and the model with seems to be less accurate than the one given by Eq. (2). On the other hand, good correspondence between the models is found for high fluxes. In this case, the variation in is less sensitive to the beam flux and does not strongly vary with atmospheric depth. Then, the model seems to give the appropriate results. The reason the higher beam fluxes lead to a lower emission in the lines is that the energy of the return current is not sufficient to excite hydrogen atoms as can be seen in Table 1.
In this paper we used a simple model of the 10 keV electron beam propagating in the chromosphere. We used two different models of the return-current formation and calculated the differential effect on the profiles of the hydrogen H, H, and H lines of the semiempirical F1 model.
The return-current flux only depends on the total flux of the beam. For a realistic power-law distribution of the beam, the main part of the beam flux is given by electrons with energy close to the low-energy cutoff of this distribution. Therefore, for simplicity we used the monoenergetic beam in our model. Moreover, the excitation and ionization cross-sections oflow-energy return-current electrons are larger than those for beam electrons, which makes the return-current effects on line core formation stronger. Taking high-energy beam electrons into account (i.e. using power-law distribution) would lead to increased emission in the line wings due to penetration of those electrons into the deeper atmospheric layers. However, the total flux of the beam in these layers is significantly lower than the initial beam flux, and the return current and corresponding effects are also strongly reduced.
Even our simple model shows that the effect of return current is very important for future study of the hydrogen lines formation since the energy of the return current can be expected to be on the order of the excitation threshold energies of upper hydrogen levels, for which the excitation cross-sections are high. Moreover, the fluxes are high enough to excite a sufficient number of atoms. As shown by Karlický et al. (2004), the collisional rates from the nonthermal collisions can dominate the collisional rates in the Balmer line formation regions. The two models used in this work lead to similar results for higher energy fluxes, but the result differs for lower fluxes. The excitation threshold effects seem to play an important role for higher fluxes, but they are very likely only a consequence of the monoenergetic model we used. The difference between the two models shows the used approximations to be incompatible for lower fluxes, where the the approximation of constant is not applicable. A detailed description of the energy distribution of the return-current electrons would lead to more realistic line intensities. This complex issue will be subject of a forthcoming paper, which will also study the impact polarization of the Balmer lines.
Acknowledgements
This work was partially supported by the grant 205/06/P135 of the Grant Agency of the Czech Republic, partially by the grant IAA300030701 of the Grant Agency of the Academy of Sciences of the Czech Republic, and partially by the project LC06014 Center for Theoretical Astrophysics.