A&A 471, 993-997 (2007)
DOI: 10.1051/0004-6361:20066809

Particle acceleration in solar flares: linking magnetic energy release with the acceleration process

C. Dauphin

LUTh, Observatoire de Paris, 92195 Meudon Cedex, France

Received 24 November 2006 / Accepted 21 May 2007

Abstract
Context. Particles accelerated in solar flares draw their energy from the magnetic field. The dissipative scale of the magnetic energy in the solar corona implies that the magnetic energy is transmitted to the particles through a large number of dissipative regions (DR).
Aims. To finally compute the particle energy distribution, we present an approach to linking the magnetic energy release to the acceleration process that occurs in each DR .
Methods. Although not directly observed, the magnetic energy release process is assumed to evolve in an SOC state that leads to the power-law behaviour of the probability distribution function for the energy released in each DR. We consider that most of the accelerated particles leave the acceleration region after one interaction, but we introduce a distribution of the particle acceleration lengths inside each acceleration region.
Results. We calculate the kinetic energy distribution of the accelerated particles. Our model allows us to identify the physical process that determines each part of the particle energy distribution.
Conclusions. As a conclusion, we discuss the limitations of the approach presented in this paper.

Key words: Sun: flares - acceleration of particles

1 Introduction

A solar flare is defined as a sudden transfer of magnetic energy to the ambient coronal plasma that leads to heating, flows, and particle acceleration. Observations and analyses of hard X-ray emissions are the most direct diagnostic of these non thermal particles. The last decades of observations have shown that a large number of particles are accelerated ($\sim$1038 electrons s-1) to energy going up to several MeV, which implies a high efficiency for the acceleration mechanism.

One of the most relevant pieces of information that allows us to constrain the physical process that occurs during particle acceleration is certainly provided by the particle energy spectrum. Generally, hard X-ray spectra are reproduced by a simple or a double power-law electron spectrum (Holman et al. 2003; Conway et al. 2003; Piana et al. 2003). These electron spectra show a steep break generally between 10 keV and 100 keV, and the spectral index below the break is usually lower than the spectral index at higher energy (Lin et al. 1982).

In this paper, an attempt is made to identify the physical processes that determines each part of the particle energy distribution. The basic idea is to take the particle acceleration process and the magnetic energy release process into account.

Particle acceleration mechanisms can be divided into three different classes: shock acceleration, turbulent acceleration, and direct electric field (further discussion can be found in the review of Miller et al. 1997). The energy of all aspects of the flare process, including particle acceleration, must ultimately come from the release of the stored magnetic energy. One possible way to release the magnetic energy stored in the corona is the formation of small-scale magnetic discontinuities where the frozen-in law is violated and where instabilities and anomalous resistivity can easily rise. Here, we consider that the particles are accelerated by the direct electric field that appears in such a region.

A large number of works have been devoted to analytical and numerical analyses of particle acceleration by a direct electric field in X-point (e.g. Hamilton et al. 2005), by a reconnecting current sheet (Litvinenko 1996; Zharkova & Gordovskyy 2004; Wood & Neukirch 2005), and in a spin and fan geometry (Litvinenko 2006; Dalla & Browning 2005, 2006). These different studies usually found a power-law particle-energy distribution due to the configuration of the magnetic field that generates a chaotic behaviour of the particle trajectories.

The complex magnetic environment of an active region (Arzner & Vlahos 2006), as well as the turbulent nature of photospheric motions that externally and continuously drives the system, will add new stresses to the existing large-scale topologies and lead to the formation of short-lived, small-scale magnetic discontinuities (see e.g., Vlahos et al. 2004; Fragos et al. 2004). Moreover, the properties of coronal turbulence implies that the magnetic energy is dissipated and transmitted to the particles on very small lengthscales (Einaudi & Velli 1994). Thus, every model of particle acceleration during a solar flare has to take into account the large-scale physical process that leads to the appearance of several dissipative regions where the particle acceleration process can takes place.

The idea of multiple dissipation sites in flares has been modelled by Lu & Hamilton (1991) by using the self-organised criticality (SOC) concept. Evolution of this model followed by satisfying the Maxwell's equations (e.g. Vlahos et al. 1995; Vassiliadis et al. 1998; Isliker et al. 2000). These SOC and associated cellular automata models identify dissipation occurring in many spatially separated regions. Several authors have investigated the problem of particle acceleration by using a model in which the accelerated particles move between the dissipation regions (Vlahos et al. 2004; Anastasiadis et al. 2004; Dauphin et al. 2007). Particles gain or lose energy in each dissipative region. The rate of energy change is proportional to the electric field, which is a function of the magnetic energy.

Another approach is to calculate the acceleration of particles in an electromagnetic field taken from MHD simulation or random-phase turbulence proxies or p-model (Decamp & Malara 2006). To study particle acceleration, the electromagnetic field is frozen in time, and test particles are injected into the simulation box. The energy gain of the test particles is calculated by solving the relativistic equations of motion. Such an approach has been used by Dmitruk et al. (2004, 2005) in direct simulations of homogeneous 3D MHD turbulence and by Arzner & Vlahos (2004) and Arzner et al. (2006), who have worked with random-phase turbulence proxies. Recently Onofri et al. (2006) addressed the problem of stochastic acceleration in turbulent electric field generated by 3D reconnection and Turkmani et al. (2005, 2006) addressed the issue of particle acceleration in a 3D coronal MHD model. Some particles undergo rapid and systematic energy gain and leave the computational box with relastivistic energy, whereas other particles do not intereact with the current sheets and gain no energy or little.

These different studies cannot simultaneously take into account the process that leads to the appearance of the acceleration regions and the particle acceleration process inside each of these regions. In this paper, we consider these processes by introducing two probability distribution functions that represent the spatial distribution of acceleration regions and the distribution of the particle acceleration lengths inside each acceleration region.

2 Model

2.1 Particle energy distribution from one acceleration region


  \begin{figure}
\includegraphics[width=11.5cm,height=8cm]{6809fig1.eps}\end{figure} Figure 1: Schematic view of a global acceleration region composed of multiple acceleration regions. The size of one AR is proportional to the value of the magnetic energy coming into the acceleration region. The small grey arrows represent the direction of accelerated electrons in each AR. The grey and white arrows outside the global acceleration region represent, respectively, the accelerated electrons and ions.
Open with DEXTER

We consider an isolated acceleration region (AR). Particles are accelerated by the electric field into the AR and gain energy from the magnetic energy. We equate the magnetic energy flux to the particle energy flux in order to determine the electric field into the acceleration region. The magnetic energy flux is given by the equation:

\begin{displaymath}P_{\rm mag}=\frac{1}{\pi}v_{\rm in}B_{\rm free}^{2}S_{\rm AR},
\end{displaymath} (1)

where $v_{\rm in}$ is the velocity of the plasma entering into the AR, $B_{\rm free}^{2}$ is the magnetic free energy, and $S_{\rm AR}$ is the surface of the AR. The particle energy flux is proportional to the electric field $E_{\rm AR}$ into the AR and to the average of the electron acceleration length $\Delta$l$_{\rm e}$ and proton acceleration length $\Delta$l$_{\rm p}$ over the total number of electrons and protons accelerated. Thus, the particle energy flux can be written as

\begin{displaymath}P_{\rm par}=\dot{N}_{\rm e}eE_{\rm AR}\left<\Delta l_{\rm e}\...
...}_{\rm p}eE_{\rm AR}\left<\Delta l_{\rm p}\right>_{n_{\rm p}},
\end{displaymath} (2)

where the number of particle entering into the AR per unit of time can be expressed by

\begin{displaymath}\dot{N}_{\rm e,p}=4S_{\rm AR}v_{\rm in}n_{\rm e,p},
\end{displaymath} (3)

where n is the particle density for protons and electrons. In this expression, we neglect the flux and energy gain of ions heavier than protons since they are less abundant than protons. Equating Eqs. (1) and (2) with the assumption $P_{\rm mag}=\alpha$ $P_{\rm par}$, we find the expression of the electric field,

\begin{displaymath}E_{\rm AR}=\alpha\frac{B_{}^{2}}{4\pi e(\left<\Delta l_{\rm e...
..._{\rm e}+\left<\Delta l_{\rm p}\right>_{n_{\rm p}}n_{\rm p})},
\end{displaymath} (4)

where $\alpha$, taken randomly between 0 and 1, represents the amount of magnetic energy transferred to the particles. This expression allows us to determine the energy gained for one particle (electron or proton) by the expression

\begin{displaymath}\varepsilon_{\rm e,p}=eE_{\rm AR}\left<\Delta l_{\rm e,p}\right>P(\Delta l),
\end{displaymath} (5)

where $\Delta$l is the acceleration length for one particle, index e or p represents the electron or proton energy gained, and $P(\Delta l)$ is the probability distribution function of the particle acceleration length. By inserting Eq. (4), we obtain the electron energy gained:

\begin{displaymath}\varepsilon_{\rm e}=\alpha\frac{B_{\rm free}^{2}}{2\pi n(1+\beta)}P(\Delta l).
\end{displaymath} (6)

Here $n=n_{\rm e}/2=n_{\rm p}/2$ is the plasma density, and $\beta$ is defined as the ratio of the average acceleration lenghts $\frac{<\Delta l_{\rm p}>_{n_{\rm p}}}{<\Delta l_{\rm e}>_{n_{\rm e}}}\cdot$ The complete determination of the particle energy gained requires knowing $\beta$ and $P(\Delta l)$. Below, we assume $\beta=1$ for the salve of simplicity and thus $\varepsilon_{\rm e}=\varepsilon_{\rm p}=\varepsilon$.

We assume that the probability distribution function of the acceleration length is given by

\begin{displaymath}P(\Delta l)=k_1(\Delta l)^{-\delta},
\end{displaymath} (7)

where k1 is constant and $\delta$ the spectral index of the power law distribution. The particle energy distribution is given by

\begin{displaymath}N(\varepsilon)=nP(\Delta l)\frac{{\rm d}(\Delta l)}{{\rm d} \varepsilon}\cdot
\end{displaymath} (8)

By inserting the different expressions, we find

\begin{displaymath}N(\varepsilon)=k_1n(eE_{\rm AR})^{\delta-1}\varepsilon^{-\delta}
\end{displaymath} (9)

for the particle energy distribution. For one active region, the spectral index of the particle energy distribution is the one for the particle acceleration length distribution.

2.2 Particle energy distribution in a "gas'' of acceleration regions

We now consider several DR leading to a "gas'' of acceleration regions with different values of magnetic energy incoming into each AR. Figure 1 shows a schematic view of a global acceleration region (GAR) composed of multiple acceleration regions. Particles (electrons and ions) are accelerated in a different direction into each AR. Thus, electrons and ions are accelerated in the same direction outside of the global acceleration region. The total particle energy distribution outside of the global acceleration region corresponds to the sum of the particle energy distribution in each AR. In other words, every particle interacts only with one acceleration region.

We consider that the probability distribution function of the magnetic energy entering into the different AR is given by a power law (Vlahos et al. 2004, and references therein):

\begin{displaymath}P(B_{\rm free}^{2})=k_2(B_{\rm free}^{2})^{-\zeta}.
\end{displaymath} (10)

As mentioned in the introduction, such a distribution can result from the evolution of the magnetic field in a SOC state. Thus, the electric field distribution is given by:

\begin{displaymath}P(E_{\rm AR})=k_3(E_{\rm AR})^{-\zeta},
\end{displaymath} (11)

where k2 and k3 are constants. The total particle energy distribution can be written as

\begin{displaymath}N(\varepsilon)_{\rm total}=\sum_{i=1}^{n_{\rm AR}}N(\varepsilon)_{{\rm AR},i},
\end{displaymath} (12)

where $n_{\rm AR}$ is the total number of AR and N( $\varepsilon$) is given by Eq. (9).

3 Results

We compute the total particle energy distribution given by Eq. (12) for different values of $n_{\rm AR}$. For each AR, we randomly select $\alpha$ between 0 and 1. Particles are injected into the acceleration volume with an initial Maxwellian distribution with a temperature of 106 K and with initial velocities in the range $0<v<8v_{\rm th}$ where $v_{\rm th}$ is the thermal velocity. The particle energy distribution for electrons and ions is calculated with a density of $n=10^{10}~{\rm cm}^{-3}$. We consider only super Dreicer events, thereby normalizing the minimal value of the electric field distribution to the Driecer electric field equal to $5.5 \times 10^{-4}~{\rm V~cm}^{-1}$.

3.1 Particle energy distribution in the case of one interaction for each particle

We consider an example of particle acceleration from 10 000 AR. The electric field value of each AR is simulated between ED and 105 ED with a spectral index equal to 1.6 (Vlahos et al. 1995). As an example, the particle's acceleration length values are included between $\sim$0.5 km and $\sim$16 km with a spectral index equal to 3.


  \begin{figure}
\par\includegraphics[width=7cm,height=6.5cm]{6809fig2.eps}\end{figure} Figure 2: Kinetic particle energy distribution obtained for $n_{\rm AR}=10~000$, $\zeta =5/3$ and with a distribution of particle acceleration length given by a power law of spectral index $\delta =3$. The temperature of the injected Maxwellian is 106 K.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=14cm,height=7cm]{6809fig3.eps}\end{figure} Figure 3: Kinetic particle energy distributions obtained for $n_{\rm AR}=10~000$, $\zeta =1.6$, and $\delta =3$. In the left panel, 10% of the total number of particles undergo between 1 and 1000 interactions (dashed line). In the right panel, 10% of the total number of particles undergo between 1 and 10 000 interactions (dashed line). The dash-dotted line represents the thermal part of the spectra.
Open with DEXTER

Figure 2 shows the particle energy distribution obtained from these 10 000 AR. Several results can be drawn from this figure:

1.
The distribution of the particle energy gain can be divided in three parts. The first one corresponds to the thermal plasma injected with a temperature of 106 K leading to a peak in the particle energy distribution at 43 eV. The second part of the spectrum is a power law of spectral index 1.6 corresponding to the spectral index of the magnetic energy release process. The last part of the spectrum is the particle energy distribution accelerated from the AR of the larger electric field. This part corresponds here to a power law of spectral index equal to 3.
2.
The break in the kinetic energy distribution occurs at $\varepsilon_{\rm break}$, which corresponds to the energy gain of the particles of lower acceleration length in the AR with a higher electric field value. Thus, $\varepsilon_{\rm break}$ is given by

\begin{displaymath}\varepsilon_{\rm break}=eE_{\rm max}\Delta l_{\rm min}.
\end{displaymath} (13)

In this case, the maximum value of electric field (105 ED) is $55 ~ {\rm V~cm}^{-1}$ and the minimum acceleration length is 500 cm. The resulting energy break occurs at 27.5 keV. Moreover, the break in the kinetic energy distribution can be express as a function of the Alfven velocity. Indeed, $\varepsilon_{\rm break}$ is proportional to $V^{2}_{\rm A_{\rm max}} \Delta l_{\rm min}/\left<\Delta l\right>$.
The establishment of the double power law spectra is a function of the number of acceleration regions. For a low number of acceleration regions ( $n_{\rm AR}=10$), the probability of high electric field values is very low. Thus, particles are not efficiently accelerated, and the resulting particles energy distribution is the spectrum of the thermal plasma. By increasing the number of acceleration regions ( $n_{\rm AR}$ = 100), some of the acceleration regions have a large electric field and particles are then strongly accelerated outside the thermal part of the spectrum. This effect creates bumps in the spectrum at high energy. From $n_{\rm AR}=1000$, electric field values for the different AR are regularly distributed between $E_{\rm min}$ and $E_{\rm max}$. In this case, a double power law appears clearly outside the thermal part of the spectra.

Finally, we have to underline, in the case of one interaction for each particle, that the values of the spectral index below and above the break corresponds, respectively, to the spectral index value chosen for the probability distribution of the magnetic energy release and particle acceleration length.

The parameters $\beta$ and $P(\Delta l)$ determine the acceleration process of electrons and ions inside each acceleration region. Thus, the total ion energy distribution depends on the values of the two parameters. Although the exact determination of $\beta$ and $P(\Delta l)$ for electrons and ions is beyond the scope of this paper, we can discuss the effect of these two parameters on the ion energy distribution. The parameter $\beta$ models the possibility that electrons and ions are accelerated at different energies inside each acceleration region. Compared to the kinetic energy distribution of electrons, this parameter changes the value of the maximum energy gained by the ions and changes the value of the break in the ion energy distribution. The probability distribution function $P(\Delta l)$ changes the particle energy distribution above the break.

3.2 Particle energy distribution in the case of several interactions for each particle

We now consider that some of the accelerated particles undergo multiple interactions. Moreover, we consider that the interaction of an accelerated particle with an acceleration region leads to a loss or a gain in energy for the particle. Thus, the amount of energy transmitted to a given accelerated particle during its interaction with an acceleration region is given by

\begin{displaymath}\varepsilon_{\rm e}=\pm \alpha\frac{B_{\rm free}^{2}}{2\pi n}P(\Delta l).
\end{displaymath} (14)

The plus or minus sign, selected randomly for each particle, represents the gain or the loss in particle energy.

Figure 3 shows the particle energy distribution obtained with 10% of the total number of particles that undergo several interactions. The distribution of the particle energy gain can be divided into four parts. The first one is the thermal plasma. The second part of the spectrum is due to the magnetic energy release process. The third part of the spectrum is due to the particle acceleration process in each acceleration region. The last part of the spectrum is created by the particles that undergo several interactions.Two comments can be drawn from this figure:

1.
The spectrum of the particles that undergo multiple interactions wanders from a power law. With the increase in the number of interactions, each particle has a higher and equal probability of being accelerated by every value of the electric field. This results in a spectrum that is close to an inverse exponential function, since the particle trajectory is formally, equivalent to a Levy walk. In this case, the probability distribution function of energy gained at each interaction is given by a power law.
2.
The increase in the number of interactions leads to the increase in the maximal energy gained by the particles.
Thus, particles accelerated by this stochastic process add a high-energy component in the particle energy distribution, which can leads to the formation of a new break at high energy. Such a break could explain the breaks observed at high energy in several flares. The value of the break mainly depends on the number of interaction.

4 Discussion and conclusions

We have presented an acceleration model that assumes the acceleration of particles drawing their energy from the magnetic energy release in several dissipative regions. We consider that each dissipative region transmits a different amount of magnetic energy to the particles. Into each dissipative region, each particle has a different acceleration length that leads to a different energy gain. Thus, the small-scale physical processes that occur in acceleration regions are taken into account by introducing a distribution of the particles acceleration length inside each acceleration region.

In this model, the probability distribution function of the magnetic energy release is modelled by a power law. Furthermore, numerous acceleration models with various acceleration mechanisms (direct electric field in different magnetic field configurations) simulate the particle acceleration process in a single acceleration region (Dalla & Browning 2005; Wood & Neukirch 2004, 2005; Zharkova & Gordovskyy 2004, 2005a). The results of these models support our assumption of a probability distribution function of the particle acceleration length given by a power law. Moreover, we consider that the probability function of the particle acceleration length and of the magnetic energy release are independent (i.e. the particle acceleraton length is independent of the magnitude of the electric field).

One of the new assumptions introduced in this paper is that most of the particles are accelerated one time. We use this model to compute the total particle energy distribution. Several results can be drawn from this distribution:

1.
Outside of the thermal part, we can distinguish three parts in the particle energy distribution. The shape of the low-energy part is purely due to the magnetic energy release process. The high-energy component is due to the particles that undergo multiple interactions. Between these two parts, the particle energy distribution results from the particle acceleration process in each acceleration region. We can note that the values of the spectral index below and above the break are due to the choice of the free parameters of the model.
2.
The energy of the break in the total particle energy distribution is given by the combination of the higher magnetic energy-release value with the minimum acceleration length.
3.
The high-energy component in the particle energy distribution is due to the particles that undergo several interactions between the different acceleration regions. This high-energycomponent can lead to the formation of a new break at high energy. Such a break could explain the breaks observed at high energy in several flares.

In conclusion, we note that this model allows us to identify the physical processes that lead to the formation of the different parts of the particle energy distribution deduced from the observations. Such a particle acceleration mechanism can explain the main part of the energy distribution of the accelerated particles. However, other particle acceleration mechanisms, as acceleration by waves (Miller & Roberts 1995), are certainly acting during a solar flare. The relative importance of each mechanism as a function of the particle energy has to be investigated further.

Acknowledgements
The author would like to thank L. Vlahos for helpful comments.

References

 

Copyright ESO 2007