Issue 
A&A
Volume 533, September 2011



Article Number  A81  
Number of page(s)  16  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201116741  
Published online  30 August 2011 
Diagnostics of nonthermal distributions in solar flare spectra observed by RESIK and RHESSI^{⋆}
^{1}
Astronomical Institute of the Academy of Sciences of the Czech Republic, v. v. i., Fričova 298, 251 65 Ondřejov, Czech Republic
email: kulinova@asu.cas.cz; kasparov@asu.cas.cz; elena@asu.cas.cz; karlicky@asu.cas.cz
^{2}
Department of Astronomy, Physics of the Earth and Meteorology, Faculty of Mathematics, Physics and Informatics, Comenius University, Mlynská dolina, 842 48 Bratislava, Slovakia
email: kulinova@fmph.uniba.sk
^{3}
Space Research Centre, Polish Academy of Sciences, 51622, Kopernika 11, Wrocław, Poland
email: js@cbk.pan.wroc.pl; bs@cbk.pan.wroc.pl
Received: 17 February 2011
Accepted: 7 July 2011
Context. During solar flares an enormous amount of energy is released, and the charged particles, like electrons, are accelerated. These nonthermal electrons interact with the plasma in various parts of solar flares, where the distribution function of electrons can therefore be nonMaxwellian.
Aims. We focus on the nonthermal components of the electron distribution in the keV range and analyse highenergy resolution Xray spectra detected by RESIK and RHESSI for three solar flares.
Methods. In the 2–4 keV range we assume that the electron distribution can be modelled by an ndistribution. Using a method of lineintensity ratios, we analyse allowed and satellite lines of Si observed by RESIK and estimate the parameters of this ndistribution. At higher energies we explore RHESSI bremsstrahlung spectra. Adopting a forwardfitting approach and thicktarget approximation, we determine the characteristics of injected electron beams.
Results. RHESSI nonthermal component associated with the electron beam is correlated well with presence of the nonthermal ndistribution obtained from the RESIK spectra. In addition, such an ndistribution occurs during radio bursts observed in the 0.61–15.4 GHz range. Furthermore, we show that the ndistribution could also explain RHESSI emission below ~5 keV. Therefore, two independent diagnostics methods indicate the flare plasma being affected by the electron beam can have a nonthermal component in the ~2–5 keV range, which is described by the ndistribution well. Finally, spectral line analysis reveals that the ndistribution does not occupy the same location as the thermal component detected by RHESSI at ~10 keV.
Key words: Sun: flares / Sun: Xrays, gamma rays / Sun: radio radiation / radiation mechanisms: nonthermal / methods: data analysis
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2011
1. Introduction
The energy release during the solar flares is accompanied by many phenomena, i.e. magnetic field reconnection and reconfiguration, particle acceleration, plasma heating, and electromagnetic radiation across the whole spectral range.
Xray spectra in deka keV range allow us to diagnose a thermal part of the flare plasma, as well as its nonthermal, highenergy tail. Recent highenergy resolution, hard Xray spectra obtained by Reuven Ramaty HighEnergy Solar Spectroscopic Imager, RHESSI (Lin et al. 2002) enable us to infer the corresponding electron spectrum. This can be achieved with a modelindependent approach, such as regularised techniques (Piana et al. 2003; Kontar et al. 2004). Also, the forwardfitting method is often used after assuming that the spectrum of the electrons emitting in this energy range is the sum of an isothermal component and a nonthermal powerlaw distribution: e.g., Brown (1971) and Holman et al. (2003). This approach is adopted here.
In space plasma, the nonthermal electron distributions are commonly modelled by κdistributions, based on in situ measurements, e.g., Vandas (1989); Maksimovic et al. (1997). Tsallis (1988) and Leubner (2002) have shown that the κdistribution is a consequence of the generalised entropy favoured by nonextensive statistics, therefore having a deeper physical meaning. Kašparová & Karlický (2009) interpreted RHESSI Xray spectra of several flares using κdistributions. They found that some looptop Xray sources can be described by such a distribution.
Spectral lines used for diagnostics of the nonthermal distribution from RESIK spectrometer.
Seely et al. (1987) analysed SOLFLEX line spectra of Fe xxv and Fe xxiv during several flares. They found that enhanced intensities of the satellite lines can be explained by the presence of nonthermal energy distributions that have higher and narrower peaks than the Maxwell distribution. This type of distribution functions has been used in the analysis of nonthermal distribution in laboratory plasma (Hares et al. 1979) as well. In this paper we refer to them as ndistributions. Recently, Dzifčáková & Karlický (2008) have introduced a simple model of a monoenergetic electron beam penetrating the solar atmosphere and creating a return current. They showed that adding a drift velocity corresponding to the electric current results in a nonthermal electron distribution that can be approximated by the ndistribution. Both the electron beam and the ndistribution influence the line intensities but each in a different way. The latter enhances the intensities of satellite lines, while the electron beam affects the intensities of allowed lines dominantly. Dzifčáková et al. (2008) analysed spectra of a flare using this kind of nonthermal distribution. They showed that the spectra, calculated using the ndistribution with the parameters estimated from observations, could mimic the observed relative intensities of allowed and satellite lines very well. Additionally, they showed that the high intensities of the satellite lines observed by RESIK cannot be reproduced by a multithermal distribution, which further supports the preferred choice of the ndistribution.
In this paper we apply two independent diagnostic techniques to analyse soft and hard Xray flare spectra. In the soft Xray range (~2 keV), where the bulk of the electron distribution is important, we consider the ndistribution. In the hard Xray range (above 6 keV), we assume a single powerlaw distribution describing the highenergy tail. A Maxwell distribution accounts for a possible thermal component at keV energies. We analyse flare spectra detected by two spectrometers and determine the parameters of these model distribution functions.
2. Model distribution functions
In the present work we consider three types of distributions: a Maxwell, a single power law, and an ndistribution, and for their determination we use different diagnostic methods. The thermal component is generally present and dominates the Xray emission around 10 keV. We suppose it is an isothermal emission so we approximate the corresponding electron distribution function as Maxwell.
We assume that the hard Xray emission is caused by the thicktarget bremsstrahlung (Brown 1971) of an electron beam injected into the plasma with a single powerlaw spectrum: (1)where δ is the spectral index, E_{C} the lowenergy cutoff, F_{T} represents the total rate of injected electrons of E ≥ E_{C}, and C is a constant. Symbol F(E) in Eq. (1) is the electron rate spectrum (Brown 2005) in units [s^{1} eV^{1}].
Fig. 1
ndistribution functions for different values of the parameter n = 1, 3, 5, 9, 13, 19. n = 1 represents the Maxwell distribution. 
At the lower energies of Xray flare spectra, i.e. 2–2.5 keV, we use a diagnostics based on soft Xray line spectra without continuum. We follow the diagnostic technique developed by Dzifčáková et al. (2008), which uses the normalized ndistribution (Fig. 1) to describe the electron distribution: (2)where is the normalization constant, E the energy of free electrons, k Boltzmann constant, and n and T are the free parameters of the distribution. The parameter T in this distribution is not the thermodynamic temperature, although it is also given in Kelvin. The parameter n is a dimensionless number describes the degree of deviation from the Maxwell distribution function. If n = 1, the distribution function becomes Maxwell. The mean energy of the ndistribution, Eq. (3), depends on two parameters n and T. To compare the nonthermal component with the thermal one, we use the pseudotemperature, τ, which relates the mean energy of the ndistribution to the mean energy of the Maxwell distribution (Dzifčáková 1998; Dzifčáková & Kulinová 2001): (3)where τ = (n + 2)T/3.
2.1. Diagnostics of the ndistribution
The line ratios of the three spectral lines listed in Table 1 were used for the diagnostics of the ndistribution.The intensities of two allowed lines (Si xiv and Si xiii) depend on excitation rates that are integrals of a collisional cross section, velocity, and an electron distribution function. The intensities of the allowed lines are thus sensitive to the highenergy tail part of the electron distribution function (Fig. 2).
Fig. 2
Schematic figure demonstrates which parts of a distribution function influence the intensities of spectral lines (Table 1) used for diagnostics of the ndistribution. The hatched regions of the distribution influence the intensities of the allowed lines of ions Si xiv and Si xiii. The discrete energy of the doubly excited state of Si xiid is marked by a thick vertical line. 
Fig. 3
Synthetic line ratios are functions of both n and τ. The plot shows the synthetic line ratios for parameter n = 1, 3, 5, 7, 9, 11, and the curves of constant log τ (grey numbers). The curve n = 1 represents the line ratios for the Maxwell distribution. 
Besides the excitation equilibrium, the intensity ratio of the allowed lines is affected by the ionization equilibrium as well. The electron ionization rate of the Si xiii is also an integral quantity that depends on the shape of the electron distribution function above the ionization energy of Si xiii. However, the radiative recombination rate depends on the shape of the whole electron distribution function. Consequently, the ratio of the allowed lines is mainly sensitive to the pseudotemperature τ but only weakly depends on the parameter n, see Fig. 3. The Si xiv/Si xiii ratio is very sensitive to T or τ, and it varies from ~0.1 for log τ = 6.8 to ~2 for log τ = 7.2. The magnitude of changes in this range of τ is smaller for the higher values of parameter n. This behaviour corresponds to the changes in the ionization and excitation equilibrium because of the ndistribution (Dzifčáková 1998; Dzifčáková & Tóthová 2007). The satellite lines of Si xii (further as Si xiid) show different behaviours. Their intensities depend on atomic parameters and the number of electrons with the energy corresponding to the energy of the doubly excited state of ion Si xii above its ionization limit (Fig. 2). Thus, the ratio Si xiii/Si xiid is extremely sensitive to the shape of the electron distribution function. For detailed information see Dzifčáková et al. (2011a), who investigate the influence of the nonthermal plasma bulk and the highenergy tail on the Si ionization and its line spectra. The observed enhanced intensities of satellite lines (Dzifčáková et al. 2008) require an increased number of particles at Si xiid energy in comparison to the Maxwell distribution (Figs. 1–3). The ndistribution with higher n satisfies this condition.
Although the electron distribution is nonthermal, this does not necessarily imply that the ionization is out of equilibrium. The ionization rates of Si XIV are of the order of 10^{10} cm^{3} s^{1}, which results in the ionization time of about 1 s for the electron density 10^{10} cm^{3}. As the results show (this paper and Dzifčáková et al. 2008), the typical time scale of changes in the derived plasma parameters is longer, several tens of seconds. Thus, the ionization equilibrium can be achieved if the conditions for formation of the nonthermal distributions change on a longer time scale than the ionization time. That seems to be the case for the present analysis, therefore the assumption of the ionization equilibrium is valid here.
In our model the ndistribution represents the bulk of the plasma electron distribution. It does not include the enhanced highenergy tail, representing a beam of accelerated electrons, usually observed at hard Xray energies.
We would like to point that, although the κdistribution also provides the enhanced number of particles in its highenergy tail, its shape in the region of 1–4 keV is not able to explain the enhanced intensities of the satellite lines. This can occur only when the gradient of the electron distribution function is high enough in this energy region, and this is not fulfilled either by κdistribution or the Maxwell one.
3. Instruments and data analysis
We analysed flare spectra observed by two Xray spectrometers: REntgenovsky Spektrometr s Izognutymi KristalamiRESIK (Sylwester et al. 2005) and RHESSI.
3.1. RESIK
The first spectrometer, RESIK, is an uncollimated spectrometer with two bent quartz crystals that covers the 2.0–3.7 keV energy range. It has four spectral channels with nominal wavelength ranges: 3.4–3.8 Å, 3.83–4.27 Å, 4.35–4.86 Å, and 5.00–6.05 Å. In the present work we analysed spectra within the 5.15–5.95 Å spectral range coming from the fourth channel where the allowed lines of ions S xv, Si xiv, Si xiii, and the satellite lines of ion Si xii dominate. The raw RESIK spectra were reduced to the absolute flux units by incorporating all the known instrumental factors (Sylwester et al. 2005). Then, the reduced spectra were timeweighted according to the exposition times of individual spectra. To improve the signaltonoise ratio, we averaged the spectra in nearly equidistant time intervals, see Tables A.1, A.3, and A.5. The diagnostic method of RESIK data is based on ratios of line fluxes, see Sect. 2.1, therefore we subtracted a linear continuum from the timeweighted spectra. However, RESIK spectra are contaminated by an artificial continuum flux that are like mounds. The exact origin of these “mounds” is still unclear, but we removed them using the same method as described in Dzifčáková et al. (2008). The imperfectly removed “mounds” can lead to relatively high errors of Si xiv line intensity, in particular, thus causing higher uncertainty of n for the spectra registered during a decay phase when this line intensity is low. Finally, the line spectra were fitted by Gaussian profiles with a constant full width at half maximum (FWHM). The errors of the line fluxes from the count statistics and fitting of the spectral lines reach 35%.
Diagnostics of the ndistribution from the RESIK spectra makes use of the precalculated grid of synthetic Si line ratios (Dzifčáková et al. 2008). These synthetic line ratios are functions of both parameters n and τ, see Fig. 3. Applying the spline interpolation to the synthetic line ratios, we determined the values of n and log τ with their errors from the measured line ratios.
Our line ratio diagnostic method can reliably detect only n ≤ 11. For higher values of n, the curves of constant n start to overlap and lead to ambiguity in n. Therefore for all the measured line ratios lying below the curve n = 11 (Fig. 3), the parameter n was set artificially to 11 and the upper limit of n to 12. Such an upper limit was chosen only for display purposes, and physically it means that n ≥ 11. In contrast to the previous paper by Dzifčáková et al. (2008), log τ was determined using both line ratios (Table 1). There is a lower limit for log τ: Fig. 3 shows that only log τ ≳ 6.8 can be reliably detected. In the case of n = 11, log τ and its upper limit were determined as the values corresponding to n = 11 and 12, respectively.
Obviously, the ndistribution with any parameter n > 1 means a deviation from the Maxwell distribution. However, the errors from the count statistics and from removing the “mounds” mainly affect the intensities of weak lines, and the higher uncertainties can usually lead to overestimated values of n during the decay phase. Therefore, in this study we concentrate on the time intervals of n ≥ 5.
Fig. 4
RESIK and RHESSI analysis results of the 2003 February 22 flare. Top left: the synthetic line ratios for n = 1, 3, 5, 7, 9, 11 (solid, dotted, dashed, dotdashed, dotdotdot dashed, and longdashed lines, respectively), and the measured line ratios (stars) with their error bars. The arrows follow the time evolution of the RESIK spectra (0 marks the first observed spectrum, see Table A.1). Top right: the time evolution of the parameter n (with error bars, upward arrows indicate n > 12) of the ndistribution and the total flux (dashed line) from all four RESIK channels. Bottom left: the time evolution of log τ (diamonds) from RESIK and log T (triangles) from the RHESSI thermal component. Bottom right: the time evolution of the spectral index δ (triangles) in comparison with the parameter n (diamonds and dashed line). This panel zooms in on the time evolution of n depicted in the panel above. The points in all the plots are depicted for the middle time of each spectrum. 
3.2. RHESSI
The second spectrometer, RHESSI, was designed to investigate the flare spectra in the energy range from 3 keV up to 7 MeV. We were interested in nonthermal bremsstrahlung radiation produced by an electron beam injected to ambient plasma. In all the analysed flares the enhanced hard Xray emission was observed approximately up to 100 keV. Among the RHESSI flare spectra, we studied in detail those time intervals during which we could observe both thermal and nonthermal components.
In the hard Xrays, the nonthermal components of the flares lasted from two to about four minutes. At later times the hard Xray emission of the flares could be fitted by one or two thermal components. We assumed the emission is composed of continuum emission caused by electronproton bremsstrahlung plus contributions of Fe (~7 keV) and Fe/Ni (~8 keV) line complexes (Phillips 2004).
We fitted the spectra in the Object SPectral EXecutive (OSPEX) environment using the model for the isothermal emission and lines below about 10 keV (vth function with the default set of the solar coronal abundances, Feldman et al. 1992) and a thicktarget nonthermal emission at higher energies (photon_thick function, Brown et al. 2008) assuming a single powerlaw electron beam, Eq. (1). If necessary, the pileup_mod pseudofunction^{1} with fixed default parameters was used. For spectra without the attenuator, state A0, the drm_mod pseudofunction^{2} was applied. We analysed the output from a single RHESSI detector 4F only because of its fine energy resolution of 0.98 keV (Smith et al. 2002). Also, it produces a photon spectrum close to the mean spectrum of all front detectors usable for spectroscopy (Phillips et al. 2006). Moreover, the drm_mod function can only be used to fit singledetector spectra. The lowenergy boundary for the fits was set according to the attenuator state: ~6 keV in A0, ~8 keV in A1, and ~12 keV in A3. The upper boundary of the energy fitting ranges was set to the values where the flare data were nearly equal to the background. In this way we obtained the time evolution of temperature of the thermal component T and the parameters of the injected electron beam, δ, E_{C}, and F_{T}.
3.3. Radio data
As complimentary data we used radio observations provided by USAF Radio Solar Telescope Network (RSTN) and Ondřejov radiospectrographs (Jiřička et al. 1993). RSTN network consists of four solar radio observatories that monitor radio emission from the Sun on eight discrete frequencies (in GHz): 0.245, 0.410, 0.610, 1.415, 2.695, 4.995, 8.800, and 15.400. The solar radio flux at these frequencies is recorded each second and is given in SFU units. Two radiospectrographs RT4 (2.0–4.5 GHz) and RT5 (0.8–2.0 GHz) at Ondřejov, Czech Republic, observed the Sun with 0.1 s time resolution and sample the frequency band in 256 frequency channels.
4. Flare events
RESIK detected several flares between August 2001 and April 2003. From its catalogue of flares we picked up three flares that were simultaneously observed by RHESSI and RSTN. One of the flares was observed by radiospectrographs at Ondřejov Observatory.
4.1. 2003 February 22 flare
The first flare was observed on February 22, 2003 in the active region NOAA 10290. According to Solar Event Reports it was of GOES class C5.8 with the soft Xray beginning at 09:24 UT, and with its maximum at 09:29 UT and end at 09:31 UT.
Fig. 5
Radio bursts observed on February 22, 2003 at six frequencies 15.40, 8.80, 4.99, 2.69, 1.41, and 0.61 GHz by RSTN. The dashed vertical line at 09:29 UT marks the time of the GOES soft Xray flare maximum. The heavy black line denotes the time interval with the parameter n ≥ 5 determined from RESIK and the heavy grey line represents time intervals with n < 5. The short black line denotes the time interval when groups of dmtype III bursts were observed, see Fig. 6. 
Fig. 6
Composite radio spectrum at 0.8–4.5 GHz observed on February 22, 2003 by the Ondřejov radiospectrographs showing a group of type III bursts in the time interval 09:28:16–09:28:33 UT. 
The results from the RESIK line spectra and RHESSI analysis are shown in Fig. 4. The top lefthand panel shows the comparison of the Si xiv/Si xiii vs. Si xiii/Si xiid synthetic line ratios with the measured ones. The significant deviation of electron distribution from the thermal one occurred slightly before the GOES soft Xray maximum and during the time of increasing RESIK Xray flux (09:27:43–09:29:03 UT), see Fig. 4 (right column). Later on, n ~ 3, which suggests slow thermalization of the flare plasma.
Figure 4 (bottom left) shows that the log τ decreased monotonically from its highest value, log τ ~ 7.3 (1.5 keV) observed in the rise phase of the flare. The log T from RHESSI is higher by ~0.1 dex than log τ. Tables A.1 and A.2 list the determined RESIK and RHESSI parameters, respectively.
The RSTN radio flux at several frequencies is plotted in Fig. 5. The peaks of radio bursts are observed within one minute before the GOES maximum, i.e. at ~09:28 UT. The flare shows a single peak at 0.61–2.69 GHz, whereas there are multiple peaks at 4.99–15.40 GHz frequencies. The observed radio bursts correlate well with n ≥ 5 detected from RESIK and with low δ derived from RHESSI. Furthermore, the bursts occurred simultaneously with a group of dmtype III bursts observed in the time interval 09:28:16–09:28:33 UT, as shown in the radio spectrum in Fig. 6. Some of these bursts have positive frequency drifts and other negative ones. This indicates the electron beams moving in both directions: downwards and upwards in the solar atmosphere.
4.2. 2003 January 7 flare
This M4.9 flare occurred in the active region NOAA 10251, starting in GOES soft Xrays at 23:25 UT, reaching its maximum at 23:33 UT and ending at 23:40 UT. During this flare, see Fig. 7, parameter n started to rise about five minutes before the GOES soft Xray maximum and remained ≥ 5 until the end of the soft Xray flare, except for a single drop at ~23:29 UT. After 23:42 UT it decreased below five. The parameter n generally follows the total RESIK flux plotted in the top right of Fig. 7. This is different from the previous flare where the maximum of the total RESIK flux occurred about two minutes after the maximum values of n. The spectral index δ shows a distinct time evolution from two other analysed flares. Its value stays almost constant, ~7, with two noticeable drops, i.e. local minima, at 23:30 UT and 23:32 UT.
On the other hand, the radio data suggest correlations with the time behaviour of δ and n. The radio burst observed by RSTN (Fig. 8) displays two peaks at frequencies 2.69–15.40 GHz appearing before the soft Xray maximum. The times of these radio peaks, i.e. ~23:30 and ~23:32 UT, correspond well to the two drops in the spectral index δ. Furthermore, the parameter n shows an increase up to 11 at the time of the second RSTN peak and remains at that value for several minutes (see RESIK spectra 5 and 6 in Table A.3 and Fig. 7 top right). Additionally, the enhanced radio flux is still visible up to 23:40 UT. During that time, RESIK spectra still indicate deviations from the Maxwell distribution, exhibiting n ≥ 5 (Fig. 8). In contrast to it, RHESSI spectra suggest thermalization of the hard Xray emitting plasma sooner, at ~23:33 UT, when δ rises to ~9. Since then, no significant nonthermal component is present in the RHESSI spectra.
Time evolutions of log τ and log T during the flare are plotted in the bottom left of the Fig. 7. The determined values of log τ span the range of 6.81–7.29 (0.6–1.7 keV). After reaching its maximum at ~23:31 UT, i.e. two minutes before the soft Xray maximum, τ starts to decrease until the end of the flare. Again, log T obtained from RHESSI thermal component are higher than log τ with values between 7.29–7.44 (1.7–2.4 keV). In contrast to τ, T has a rather spiky character. Its maximum was observed before the soft Xray maximum, too, at 23:32 UT. Tables A.3 and A.4 list the determined RESIK and RHESSI parameters, respectively.
4.3. 2002 October 4 flare
The third flare was of GOES class M4.0 and occurred in the active region NOAA 10137. It started at 05:34 UT, reached the GOES maximum at 05:38 UT, and ended at 05:41 UT. For this flare the parameter n ≥ 5 was observed between ~05:35–05:42 UT (Fig. 9, top right). After ~05:39 UT n started to decrease. Like the 2003 January 7 flare, the global behaviour of n follows the RESIK total flux. On the other hand, the time evolution of the spectral index δ (Fig. 9, bottom right) has a similar shape to the 2003 February 22 flare (Fig. 4, bottom right), but in this case the RHESSI spectra are much steeper, therefore δ has higher values. Fast variations of n and its large uncertainties during the rise and the maximum phase are not a real feature, but are probably due to low photon flux and related high statistical errors. The time evolution of log T obtained from RHESSI is nearly constant, and its values are again higher than log τ. Tables A.5 and A.6 list the determined RESIK and RHESSI parameters, respectively.
Radio burst observed by RSTN (Fig. 10) occurred before the flare’s soft Xray maximum and started at ~05:36 UT. The shape of the radio flux shows a double peak at frequencies 4.99 GHz and 8.80 GHz with the first peak occurring at about 05:36 UT and the second one about 05:37 UT. At the highest RSTN frequency, 15.40 GHz corresponding to the lower chromosphere, the burst shows a triple peak between 05:36–05:37 UT. In contrast to this, the 2.69 GHz flux shows a single broad peak with a maximum at about 05:37 UT. There radio bursts correlate well with the presence of the RHESSI nonthermal component observed at 05:35:28–05:37:32 UT. Furthermore, enhanced radio flux at all four frequencies, as well as ndistribution of n ≥ 5, were still visible after the GOES maximum.
5. RHESSI and ndistribution
If there is any nonthermal ndistribution really present in flaring plasma, its electrons should also contribute to the production of bremsstrahlung radiation. Therefore, we included the thintarget model of bremsstrahlung from ndistribution into fitting of RHESSI spectra and tried to test the parameters of ndistribution independently from RESIK. We performed such an analysis only for the first two flares, 2003 February 22 and 2003 January 7, as we needed to use the RHESSI spectra observed in A0 attenuator state for this purpose. This allowed us to examine the lowest part of the RHESSI energy range that is closest to the energies observed by RESIK. Unlike analysis of the RESIK Xray line spectra, we suppose that RHESSI will give us the information about continuum radiation produced by the ndistribution.
We therefore revisited A0 RHESSI fits, as seen in the first row in Tables A.2 and A.4 and Fig. A.1, where the fits are shown with their normalized residuals. Both fits started at 6 keV and they had ( means the reduced χ^{2}). Then, using the bestfit parameters as a starting point, we pushed the lowenergy boundary down to 4 keV and let the parameters of the isothermal and the powerlaw component vary. (When either pileup_mod or drm_mod pseudofunctions were used, the values of their parameters were taken from the previous fits and kept fixed.) However, no satisfactory fits could be obtained in this way. In this way was found for both flares, and the residuals were unacceptably large, showing a systematic pattern below 10 keV (Fig. A.2).
Therefore, in the last set of fits we added the thintarget bremsstrahlung radiation from the ndistribution. This component was calculated by a new fitting function, thin_ndistr, which has the four parameters to fit: n, τ, “emission measure” (i.e. ambient proton density × volume × electron density in the ndistribution), and a highenergy cutoff. (The function, thin_ndistr, was incorporated into SSW tree and can be accessed from OSPEX.) The highenergy cutoff for the ndistribution is introduced for a numerical purpose only. It was set to 1 MeV and fixed for all measurements. The initial values of other three parameters of the thin_ndistr function were set to those obtained from the RESIK diagnostics.
These fits (Fig. 11) improved significantly both in terms of residuals and . The component due to the ndistribution describes the lowest energies of the fits, whereas the spectrum at higher energies is not affected, and its parameters are similar to those from the original fits starting from 6 keV. Tables A.7 and A.8 summarize the parameters for individual flares and for all three fits: the original fit, the fit with the enlarged fitting range, and the one combined with the ndistribution bremsstrahlung. To complete the derived parameters, we also list RESIK results for the corresponding time intervals in those tables. Tables A.7 and A.8 and Fig. 11 show the fits for a specific case of the parameter n being fixed at a value derived from RESIK. The reason is that the studied Si lines provide the best possibility for determining the parameter n because the line ratios with the satellite lines are very sensitive to the shape of the distribution considered. In addition to the fits with fixed n, we also analysed fits setting n as a free parameter. Both sets of fits were similar in terms of χ^{2} and residuals. However, in the latter case the best fit was not well defined since a wide range of ndistribution parameters resulted in a similar value for . Such behaviour suggested that the parameter uncertainties could be large.
In order to estimate the uncertainties of the fitted parameters, we examined the χ^{2} values by keeping a selected parameter fixed within a range of values and by fitting the other parameters (see also chi2_map method in OSPEX). Having done so, we found out that χ^{2} does not exhibit well localized minima for ndistribution parameters. Moreover, in some cases only their lower or upper limits could be estimated. This is a situation similar to the lowenergy cutoff (and consequently the total flux) for which it is known that in some cases the spectral fits could only provide an upper limit estimate. A cutoff below the upper limit still fits the data because at low energies the nonthermal flux is usually well below the thermal component (Sui et al. 2005). On the other hand, this χ^{2} mapping technique allowed us to obtain thermal parameters and spectral index δ uncertainties, i.e. as values corresponding to . However, we note that this approach does not provide joint uncertainties of fitted parameters, but instead corresponds to a single parameter uncertainty independent of any knowledge about the others. To obtain a joint estimate of several parameters at a 68% (“1σ”) confidence level, a much larger region, , would have to be investigated (Lampton et al. 1976; Nakamura et al. 2010). But this is beyond the scope of the paper.
The estimations of the parameter uncertainties are given in Tables A.7 and A.8. For the case of the 2003 February 22 flare, the fit with n fixed at the RESIK value, i.e. n = 11, gives ndistribution parameters that are slightly above the estimated 1σ range, i.e. n = 2 − 9, kτ = (1.0 − 1.5) keV (but well within 2σ). However, due to the weak dependence of χ^{2} on the ndistribution parameters, such a fit could be considered acceptable. Although the uncertainties of n and τ obtained either from RHESSI or RESIK data are large for both events, their ranges overlap.
We also performed a test to check the reality of the excess emission detected in the lowest RHESSI energies. A spectrum from a different detector (3F) was analysed and resulting fitting parameters agreed well within the uncertainties with those discussed above. Therefore, the excess emission, fitted here by the ndistribution, is not an instrumental artefact but a real feature.
Using the reconstructed images, we estimated the electron densities in plasma for different electron distributions. From the RHESSI fits we have the information about the emission measure for the thermal component (EM_{th}) and for the ndistribution (EM_{n}), as well as the total rate of injected electrons (F_{T}) for the power law. The ndistribution appears at the lowest part of the RHESSI spectrum. We therefore reconstructed the images in the 3–6 keV range for it, the 6–12 keV channel was used for the thermal component, and the 25–50 keV channel was dedicated to the hard Xray emission. From the images we determined the area S within the 70% contour of the maximum and supposed that the particular volume of emitting plasma has a spherical shape with a circular cross section S. Both for the ndistribution and the Maxwell distribution we supposed that N_{p} = N_{e} and and , where index “th” is for the thermal component and “n” corresponds to the ndistribution. The 70% contours in the 3–6 keV and 6–12 keV channels were almost the same (Fig. 12), so we used only the 6–12 keV channel for the electron density estimation, assuming V_{th} = V_{n}. The electron density for the powerlaw component was derived using the relation , where S_{hard} is the area in the 25–50 keV hard Xray channel.
In both flares the electron densities for ndistribution are three to four times higher than the electron densities of the thermal component. Moreover, they are about three to four orders of magnitude higher than the electron densities that correspond to the powerlaw distribution. Furthermore, the thermal and ndistribution electron densities are similar in both events, while N_{e} for the powerlaw distribution differs nearly by one order. The 2003 January 7 flare was a limb flare (Fig. 12, right), so we register only its projection. It is likely that the flare footpoints, from where the dominant contribution to the 25–50 keV channel comes, are the most affected by that projection.
6. Discussion
Basically, we suppose that a nonthermality in soft Xray range manifested by the ndistribution is connected with a nonthermality in the hard Xray region. Then at a first sight, one could expect to see an “Xpattern” in n and δ plots, i.e. when δ is the lowest, n reaches the highest values. But, in reality this connection is not so simple and does not depend on relation between n and δ only. It very likely also depends on electron flux density (Dzifčáková & Karlický 2008) if we assume that the ndistribution is produced because of the return current. None of the three analysed flares exhibits a clear “Xpattern”. There is a hint of such tendency, i.e. a smooth decrease in n versus a rise in δ, for the simple compact flare from February 22, 2003 of the GOES class C (see Fig. 4, bottom right), however the large uncertainties of n do not allow us to consider this as significant. The other two M class flares were more complex. The 2002 October 4 flare was a two ribbon flare occurring near the limb centre. The 2003 January 7 flare was a limb flare exhibiting a rather complex magnetic configuration, possibly with a plasmoid structure (Fig. 12 and Dzifčáková et al. 2008, Fig. 2). Additionally, both these flares show longer time occurrence of the nonthermal ndistribution than the powerlaw component observed in RHESSI spectra.
In all three flares the time interval with n ≥ 5 corresponds well with radio bursts observed at frequencies 4.99–8.80 GHz. Also, we found a good correlation of radio emission with hard Xrays related to electron beams. This behaviour is known and has been studied, e.g., by Aschwanden et al. (1995), Benz et al. (2005), and Fárník & Karlický (2007). For the 2003 February 22 flare, no enhanced radio emission was observed after the burst. On the other hand, the long duration enhanced radio emission was presented, together with n ≥ 5 for both M class flares, but a correlation with δ shows a shorter duration.
Derived electron densities for Maxwell , power law , and ndistributions .
For the 2003 February 22 flare, the type III bursts with positive and negative frequency drifts were observed within the time interval with n ≥ 5. These bursts indicate the electron beams are moving downwards and upwards in the solar atmosphere. Because each electron beam needs to be accompanied also by the return current electrons (Karlický 2009), such a radio spectrum can be considered as evidence of the nonMaxwell electron distribution function representing an electron beam with a return current in the region of the radio source. Assuming radio emission generated on the plasma frequency (800–4500 GHz, see Fig. 6), these type III bursts thus show that deviations from the Maxwell distribution are present in the atmospheric layers with the plasma density N_{e} from 7.9 × 10^{9} cm^{3} to 2.5 × 10^{11} cm^{3}, i.e. in the chromosphere through the transition layer to the low corona. The electron density of the ndistribution determined from the RHESSI spectral fits and images, see Table 2, agrees well with the above estimate from the radio emission. This further supports the idea of electron beams and of an accompanying return current that causes the ndistribution. Additionally, assuming that the ratio (Melrose 1980, pp. 203–204), where is the plasma density of the electron beam generating type III bursts, then we can see that is much lower, ≈ 7.9 × 10^{3} − 2.5 × 10^{5} cm^{3}, than the plasma densities of the ndistributions presented in Table 2. Thus a contribution of the electron beam electrons producing type III bursts to the ndistribution is very small.
Furthermore, we made the following comparison of parameters of flares in February 22, 2003 and January 7, 2003. In the model in which the ndistribution is formed by electrons of the return current induced by the electron beam (Dzifčáková & Karlický 2008), we can write where is the electron density corresponding to the ndistribution, N_{b} the density of the powerlaw tail (electron beam), the mean return current velocity in the proton coordinate system, and v_{b} the mean velocity of the powerlaw tail (beam). The ndistribution densities in both the flares are nearly the same (Table 2). Therefore, if we assume the same velocity , i.e. the same ndistributions in both flares, then the mean velocity of the powerlaw tail (electron beam) v_{b} in the 2003 February 22 flare should be greater than that in the 2003 January 7 flare. However, estimation of the mean velocities significantly depends on the lowenergy cutoff that was determined with large uncertainty. Nevertheless, the powerlaw index δ ~ 3 in the 2003 February 22 flare is much smaller than that (δ ~ 6.5) in the 2003 January 7 flare, which means that the mean velocity in the 2003 February 22 flare is greater than in 2003 January 7, as estimated. That also supports the model of the ndistribution based on the return current.
We expected that the bremsstrahlung from ndistribution would very likely manifest itself at low energies of the RHESSI spectral window. This was confirmed by the RHESSI fits of ndistribution discussed in the previous section. We note that in both cases the ndistribution significantly contributes to the Xray spectrum in a rather narrow energy range, therefore the determined parameters of the ndistribution have large uncertainties. In all spectra the thermal component had to be considered to account for the emission in the ~6–10 keV energy range. Forcing the ndistribution component to higher energies instead of using the thermal component did not work: fits consisting of an ndistribution, two Gaussian lines (as approximations of the line complexes at ~ 6.7 and 8 keV), and a powerlaw highenergy tail did not lead to acceptable results. Finally, we also fitted the RHESSI spectra using a double power law and an isothermal component, i.e. without the ndistribution, to examine a possibility that the lowenergy emission is due to a flat powerlaw component extending down to keV energies. Although an acceptable fit of was obtained, the resulting lowenergy powerlaw part was very steep; i.e., δ_{L} ~ 11 with δ_{L} > δ, thus resembling the ndistribution component.
Neither RHESSI nor RESIK spectra are spatially resolved. The radiation produced by the thermal component and the ndistribution comes from the coronal plasma, while hard Xray bremsstrahlung origins primarily at the footpoints (Fig. 12). The 25–50 keV contour is the largest one with 3–6 keV and 6–12 keV countours in the middle (Fig. 12, right) or in its upper part (Fig. 12, left). Even though the contours of the 3–6 keV and 6–12 keV channels overlay each other, we suppose that the ndistribution and the thermal one are localized in different substructures of flaring plasma. There are several reasons that led to this conclusion.
Fig. 11
RHESSI fits (from 4 keV) of the single power law (thicktarget model), the thermal component, and the ndistribution (thintarget model) in the low energy range, together with residuals: the 2003 February 22 flare (left) and the 2003 January 7 flare (right). 
Fig. 12
RHESSI images in the 12–25 keV channel with contours of 70% from channels: 3–6 keV (full line), 6–12 keV (dashdot line), and 25–50 keV (long dash line): the 2003 February 22 flare (left) and the 2003 January 7 flare (right). The thin solid line marks the solar limb for the 2003 January 7 flare. 
First, T determined from the RHESSI thermal component is always higher than the values of τ from the ndistribution. This difference between log T and log τ is usually 0.1–0.2 dex, and it cannot be explained by the measured uncertainties (Figs. 4, 7, 9, bottom left). The enhanced intensities of Si xiid lines observed by RESIK are the second argument for a different spatial location of the Maxwell and the ndistribution. The high ratio of Si xiid to Si xiii lines needs a strong gradient of particles in the energy range above approximately 1.5 keV (Figs. 1, 2). The high number of electrons with the energy corresponding to the excitation energy of doubly excited state results in the high intensities of the satellite lines. On the contrary, the low number of electrons with energies higher than 2 keV results in the lower excitation rates, hence the lower intensities of the allowed lines. As a result, the observed line ratio of high n can be obtained. However, the relation between the intensities of the satellite and the allowed lines is not so straightforward because electrons with the energy above 2 keV influence also the ionization equilibrium and control the Si xiii/Si xiv ratio. The ndistribution alone fulfils the condition of a strong gradient of the electron distribution function in its highenergy range.
The electron distribution composed of the ndistribution and the Maxwell one does not have such a strong gradient in its highenergy part because it has a much higher number of electrons there. Such a composed distribution has the same number of electrons with the energy corresponding the excitation energy of doubly excited state and therefore the same intensities of the satellite lines. But, it also has much higher intensities of the allowed lines due to the high number of particles with the energy higher than 2 keV. This must result in a much higher Si xiii to Si xiid ratio than is observed. The line emission from different locations in the plasma is proportional to their emission measures. The RHESSI spectra indicated that the emission measure of the ndistribution is much higher than for the Maxwell one (Tables A.7, A.8); so the line spectrum of the ndistribution dominates over the Maxwell spectrum, which agrees with the observations (Dzifčáková et al. 2011b). Finally, the electron densities of the powerlaw distribution, as well as of the electrons producing the type III bursts, are too low compared to the other two assumed distributions to have any significant effects on the studied RESIK Si line spectrum; see also the recent results on a combined distribution composed of the ndistribution and the powerlaw tail (Dzifčáková et al. 2011a).
In summary, those facts point to different spatial origins for the ndistribution and the thermal emission.
7. Conclusions
In this analysis we have determined that during flares the plasma affected by an injected beam of accelerated electrons can have the nonthermal component in the 2–6 keV region and that it can be described by the ndistribution. We derived the parameters n and τ of this ndistribution using the flare spectra from two satellite spectrometers: RESIK and RHESSI. The spectra were analysed by two distinct diagnostic methods. The soft Xray lines of Si xiv, xiii, xiid observed by RESIK provided us with the best possibility of analysing the behaviour of parameter n. The RHESSI analysis of the ndistribution is based on studying its continuous, bremsstrahlung radiation. Although the parameters n and τ obtained from the RHESSI fits have large uncertainties, their values still overlap with the RESIK results quite well. The RESIK and RHESSI results indicate that departure from the Maxwell distribution at 2–6 keV energies can really occur. The nonthermal ndistribution and the powerlaw component are observed in good temporal correlation with the radio bursts detected in the 0.61–15.40 GHz range by RSTN. The thermal component is still present in RHESSI spectra. RHESSI images show that the emission from both thermal and nonthermal ndistributions comes from nearly the same area. Despite this we argue that in reality each distribution occupies different substructures within the depicted area. Such a conclusion is supported by the character of the observed RESIK line spectra with high fluxes of satellite lines. Finally, our joint analysis of radio and Xray emission also supports the model of formation of the ndistribution because of the returncurrent.
Acknowledgments
This work has been supported by the Scientific Grant Agency VEGA, Slovakia, grant No. 1/0240/11 and Scientific Grant Agency of Czech Republic GAČR, grants Nos. 205/09/1705 and P209/10/1680, grant IAA300030701 (GA AS CR), and the research project AV0Z10030501. We are very thankful for the open data policy of RESIK, RHESSI and GOES. B.S. and J.S. acknowledge support from the Polish Ministry of Science grant N N203 381736, as well as from the European Commission’s Seventh Framework Programme (FP7/20072013) under grant agreement No. 218816 (SOTERIA project, http://www.soteriaspace.eu).
References
 Aschwanden, M. J., Benz, A. O., Dennis, B. R., & Schwartz, R. A. 1995, ApJ, 455, 347 [NASA ADS] [CrossRef] [Google Scholar]
 Benz, A. O., Grigis, P. C., Csillaghy, A., & SaintHilarie, P. 2005, Sol. Phys., 226, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Fárník, F., & Karlický, M. 2007, Sol. Phys., 240, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, J. C. 1971, Sol. Phys., 18, 489 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, J. C. 2005, in Sol. Mag. Phenomena Proc. of the 3rd Summerschool and Workshop held at the Solar Observatory Kanzelhöhe, Kärneten, Austria, August 25–September 5, 2003, ed. A. Hanslmeier, A. Veronig, & M. Messerotti (The Netherlands: Springer) [Google Scholar]
 Brown, J. C., Kašparová, J., Massone, A. M., & Piana, M. 2008, A&A, 486, 1023 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dzifčáková, E. 1998, Sol. Phys., 178, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Dzifčáková, E. 2005, in Proc. of the 11th European Solar Physics Meeting, The Dynamic Sun: Challenges for Theory and Observations, 11–16 September, Leuven, Belgium, ed. D. Danesy, S. Poedts, A. De Groof, & J. Andries, Published on CDROM, ESA SP600, 119 [Google Scholar]
 Dzifčáková, E. 2006a, in Proc. of SoHO17: 10 Years of SoHO and Beyond, 7–12 May, Giardini Naxos, Sicily, Italy, ed. H. Lacoste, & L. Ouwehand, ESA SP617, 89.1 [Google Scholar]
 Dzifčáková, E. 2006b, Sol. Phys., 234, 243 [NASA ADS] [CrossRef] [Google Scholar]
 Dzifčáková, E., & Karlický M. 2008, Sol. Phys., 250, 329 [NASA ADS] [CrossRef] [Google Scholar]
 Dzifčáková, E., & Kulinová, A. 2001, Sol. Phys., 203, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Dzifčáková, E., & Kulinová, A. 2006, in Proc. of SoHO17: 10 Years of SoHO and Beyond, 7–12 May, Giardini Naxos, Sicily, Italy, ed. H. Lacoste, & L. Ouwehand, ESA SP617, 88.1 [Google Scholar]
 Dzifčáková, E., & Tóthová, D. 2007, Sol. Phys., 240, 211 [NASA ADS] [CrossRef] [Google Scholar]
 Dzifčáková, E., Kulinová, A., Chifor, C., et al. 2008, A&A, 488, 311 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dzifčáková, E., Homola, M., & Dudík, J. 2011a, A&A, 531, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dzifčáková, E., Kulinová, A., & Kašparová, J. 2011b, Cool Stars 16, ASP Conf. Ser., submitted [Google Scholar]
 Feldman, U., Mandelbaum, P., Seely, J. L., Doschek, G. A., & Gursky, H. 1992, ApJS, 81, 387 [NASA ADS] [CrossRef] [Google Scholar]
 Hares, D. J., Kilkenny, J. D., Key, M. H., & Lunney, J. G. 1979, Phys. Rev. Lett., 42, 1216 [NASA ADS] [CrossRef] [Google Scholar]
 Holman, G. D., Sui, L., Schwartz, R. A., & Emslie, A. G. 2003, ApJ, 595, L97 [NASA ADS] [CrossRef] [Google Scholar]
 Jiřička, K., Karlický, M., Kepka, O., & Tlamicha, A. 1993, Sol. Phys., 147, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Karlický, M. 2009, ApJ, 690, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Kašparová, J., & Karlický, M. 2009, A&A, 497, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kontar, E. P., Piana, M., Massone, A. M., Emslie, A. G., & Brown, J. C. 2004, Sol. Phys., 225, 293 [NASA ADS] [CrossRef] [Google Scholar]
 Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177 [NASA ADS] [CrossRef] [Google Scholar]
 Landi, E., Del Zanna, G., Young, P. R., et al. 2006, ApJS, 162, 261 [NASA ADS] [CrossRef] [Google Scholar]
 Leubner, M. P. 2002, Ap&SS, 282, 573 [NASA ADS] [CrossRef] [Google Scholar]
 Lin R. P., Dennis, B. R., & Hurford, G. J. 2002, Sol. Phys., 210, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Maksimovic, M., Pierrard, V., & Lemaire, J. F. 1997, A&A, 324, 725 [NASA ADS] [Google Scholar]
 Melrose, D. B. 1980, Plasma Astrophysics (New York: Gordon and Breach Science Publ.), Vol. II, 203 [Google Scholar]
 Nakamura, K., & Particle Data Group 2010, J. Phys. G Nucl. Phys., 37, 075021 [NASA ADS] [CrossRef] [Google Scholar]
 Phillips, K. J. H. 2004, ApJ, 605, 921 [NASA ADS] [CrossRef] [Google Scholar]
 Phillips, K. J. H., Chifor, C., & Dennis, B. R. 2006, ApJ, 647, 1480 [NASA ADS] [CrossRef] [Google Scholar]
 Piana, M., Massone, A. M., Kontar, E. P., et al. 2003, ApJ, 595, L127 [NASA ADS] [CrossRef] [Google Scholar]
 Seely, J. F., Feldman, U., & Doschek, G. A. 1987, ApJ, 319, 541 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, D. M., Lin, R. P., & Turin, P. 2002, Sol. Phys., 210, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Sui, L., Holman, G. D., & Dennis, B. R. 2005, ApJ, 626, 1102 [NASA ADS] [CrossRef] [Google Scholar]
 Sylwester, J., Gaicki, I., & Kordylewski, Z. 2005, Sol. Phys., 226, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Tsallis, C. 1988, J. Stat. Phys., 52, 476 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Vandas, M. 1989, Bull. Astron. Inst. Czechosl., 40, 175 [NASA ADS] [Google Scholar]
Online material
Appendix A: Complementary figures and tables
Timeweighted RESIK spectra for the C5.8 flare on February 22, 2003 and the determined ranges of the parameters n and log τ.
Results from RHESSI analysis of the 2003 February 22 flare.
Details of the RHESSI spectrum fitting (2003 February 22, 09:28:12–09:28:20 UT) and comparison with the corresponding RESIK results (09:27:43–09:28:31 UT).
Fig. A.1
RHESSI fits (from 6 keV) of the single power law (thicktarget model) and the thermal component together with residuals: the 2003 February 22 flare (left) and the 2003 January 7 flare (right). 
All Tables
Spectral lines used for diagnostics of the nonthermal distribution from RESIK spectrometer.
Timeweighted RESIK spectra for the C5.8 flare on February 22, 2003 and the determined ranges of the parameters n and log τ.
Details of the RHESSI spectrum fitting (2003 February 22, 09:28:12–09:28:20 UT) and comparison with the corresponding RESIK results (09:27:43–09:28:31 UT).
All Figures
Fig. 1
ndistribution functions for different values of the parameter n = 1, 3, 5, 9, 13, 19. n = 1 represents the Maxwell distribution. 

In the text 
Fig. 2
Schematic figure demonstrates which parts of a distribution function influence the intensities of spectral lines (Table 1) used for diagnostics of the ndistribution. The hatched regions of the distribution influence the intensities of the allowed lines of ions Si xiv and Si xiii. The discrete energy of the doubly excited state of Si xiid is marked by a thick vertical line. 

In the text 
Fig. 3
Synthetic line ratios are functions of both n and τ. The plot shows the synthetic line ratios for parameter n = 1, 3, 5, 7, 9, 11, and the curves of constant log τ (grey numbers). The curve n = 1 represents the line ratios for the Maxwell distribution. 

In the text 
Fig. 4
RESIK and RHESSI analysis results of the 2003 February 22 flare. Top left: the synthetic line ratios for n = 1, 3, 5, 7, 9, 11 (solid, dotted, dashed, dotdashed, dotdotdot dashed, and longdashed lines, respectively), and the measured line ratios (stars) with their error bars. The arrows follow the time evolution of the RESIK spectra (0 marks the first observed spectrum, see Table A.1). Top right: the time evolution of the parameter n (with error bars, upward arrows indicate n > 12) of the ndistribution and the total flux (dashed line) from all four RESIK channels. Bottom left: the time evolution of log τ (diamonds) from RESIK and log T (triangles) from the RHESSI thermal component. Bottom right: the time evolution of the spectral index δ (triangles) in comparison with the parameter n (diamonds and dashed line). This panel zooms in on the time evolution of n depicted in the panel above. The points in all the plots are depicted for the middle time of each spectrum. 

In the text 
Fig. 5
Radio bursts observed on February 22, 2003 at six frequencies 15.40, 8.80, 4.99, 2.69, 1.41, and 0.61 GHz by RSTN. The dashed vertical line at 09:29 UT marks the time of the GOES soft Xray flare maximum. The heavy black line denotes the time interval with the parameter n ≥ 5 determined from RESIK and the heavy grey line represents time intervals with n < 5. The short black line denotes the time interval when groups of dmtype III bursts were observed, see Fig. 6. 

In the text 
Fig. 6
Composite radio spectrum at 0.8–4.5 GHz observed on February 22, 2003 by the Ondřejov radiospectrographs showing a group of type III bursts in the time interval 09:28:16–09:28:33 UT. 

In the text 
Fig. 7
2003 January 7 flare. Notation is the same as in Fig. 4. 

In the text 
Fig. 8
2003 January 7 flare. Notation is the same as in Fig. 5. 

In the text 
Fig. 9
2002 October 4 flare. Notation is the same as in Fig. 4. 

In the text 
Fig. 10
2002 October 4 flare. Notation is the same as in Fig. 5. 

In the text 
Fig. 11
RHESSI fits (from 4 keV) of the single power law (thicktarget model), the thermal component, and the ndistribution (thintarget model) in the low energy range, together with residuals: the 2003 February 22 flare (left) and the 2003 January 7 flare (right). 

In the text 
Fig. 12
RHESSI images in the 12–25 keV channel with contours of 70% from channels: 3–6 keV (full line), 6–12 keV (dashdot line), and 25–50 keV (long dash line): the 2003 February 22 flare (left) and the 2003 January 7 flare (right). The thin solid line marks the solar limb for the 2003 January 7 flare. 

In the text 
Fig. A.1
RHESSI fits (from 6 keV) of the single power law (thicktarget model) and the thermal component together with residuals: the 2003 February 22 flare (left) and the 2003 January 7 flare (right). 

In the text 
Fig. A.2
Widened RHESSI fits (from 4 keV). Notation is the same as in Fig. A.1. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.