Determination of the total accelerated electron rate and power using solar flare hard X-ray spectra

Dec 22, 2018
18 pages
Published in:
  • Astrophys.J. 871 (2019) 2, 225
  • Published: Feb 4, 2019
e-Print:

Citations per year

201820202022202420241403
Abstract: (IOP)
Solar flare hard X-ray (HXR) spectroscopy serves as a key diagnostic of the accelerated electron spectrum. However, the standard approach using the collisional cold thick-target model poorly constrains the lower-energy part of the accelerated electron spectrum, hence the overall energetics of the accelerated electrons are typically constrained only to within one or two orders of magnitude. Here, we develop and apply a physically self-consistent, warm-target approach that involves the use of both HXR spectroscopy and imaging data. This approach allows an accurate determination of the electron distribution low-energy cutoff, and hence the electron acceleration rate and the contribution of accelerated electrons to the total energy released, by constraining the coronal plasma parameters. Using a solar flare observed in X-rays by RHESSI, we demonstrate that using the standard cold-target methodology, the low-energy cutoff (hence the energy content in electrons) is essentially undetermined. However, the warm-target methodology can determine the low-energy electron cutoff with ~7% uncertainty at the 3σ level, hence it permits an accurate quantitative study of the importance of accelerated electrons in solar flare energetics.
Note:
  • Accepted for publication in the Astrophysical Journal, 18 pages, 5 figures
  • Sun: activity
  • Sun: flares
  • Sun: X-rays, gamma rays
  • [1]
    First, we evaluate the plasma parameters in the flaring loop. We have chosen to evaluate the thermal parameters by using the short time interval of 02:08:52-02:09:00 UT just before the studied interval of 02:09:00-02:10:00 UT. The HXR spectrum during the 02:08:52-02:09:00 UT interval was fitted with a combined thermal function plus a cold thick-target nonthermal function (e.g., using the SSW/OSPEX functions f_vth + f_thick2). Although both nonthermal and thermal components are included in the modeling, this step is primarily used to estimate, 7: #apjaafad3fn5 the thermal plasma parameters of the coronal loop-top source, namely the temperature T (keV) and the emission measure EM0 (cm-3)
    • [1]
      the nonthermal parameters will be determined more accurately later. The minimum value of the reduced χ2 gives the best estimates of T and EM0. Even though the fitting is performed with an assumed power-law form for the cold thick-target nonthermal component, we have found that using another form (such as a kappa distribution) does not substantially affect the inferred values of the thermal parameters EM0 and T. (For SOL-05-13T02:12, the relative abundance parameter was fixed at the default value of unity, and we determined best-fit values of T = (2.52 ± 0.03) keV and EM0 = (6.51 ± 0.30) × 1048 cm-3 for 02:08:52-02:09:00 UT. Note that throughout this work, 1σ OSPEX uncertainties are used, unless stated otherwise.)
    • [2]
      With the values of EM0 and T now determined, we fit the HXR spectrum during our study time interval of 02:09:00-02:10:00 UT. The HXR spectrum during this time interval is fitted with f_vth and the warm-target fitting function f_thick_warm. Note that the total emission measure EM = EM0 + ΔEM obtained by this fit includes any contribution ΔEM from the thermalization of injected electrons. Any additional parameters/contributions pertinent to the observed flux (e.g., elemental abundances, albedo, pulse pile-up, lines, drm_mod, etc.) should also be included in this fit (Figure, 2: #apjaafad3f2). The function f_thick_warm contains 10 parameters: the first six parameters are exactly the same as for f_thick2 (i.e., the electron acceleration rate \dot{{N}_{0}} (s-1), the spectral index (low) δlow, the break energy EB (keV), the spectral index (high) δhigh, the low-energy cutoff Ec (keV), and the high energy cutoff EH (keV)). The other four parameters relate to the properties of the thermal properties of the corona: the loop temperature Tloop (keV), loop density nloop (cm-3), and half-length of the loop L (cm), and the relative elemental abundances.a. As discussed in Step 1, an isothermal estimate of the loop-top temperature, T, is determined using f_vth during 02:08:50-02:09:00 UT, and we fix {T}_{\mathrm{loop}}\,=T=2.52\,\mathrm{keV}.b. To determine the density {n}_{\mathrm{loop}}, we use the estimate of EM0 obtained in Step 1. The number density can then be estimated using {n}_{\mathrm{loop}}=\sqrt{{\mathrm{EM}}_{0}/V}, where V is the volume of the emitting plasma, found from X-ray imaging (see Figure, 3: #apjaafad3f3). (For SOL-05-13T02:12, the volume was estimated by applying the imaging algorithm Visibility Forward Fitting (VIS_FWDFIT
    • [2]
      2007: #apjaafad3bib42). Using VIS_FWDFIT, the volume of the coronal source varies with energy (see Jeffrey et al
      • Schmahl
    • a. As discussed in Step 1, an isothermal estimate of the loop-top temperature, T, is determined using f_vth during 02:08:50-02:09:00 UT, and we fix {T}_{\mathrm{loop}}\,=T=2.52\,\mathrm{keV}
      • b. To determine the density {n}_{\mathrm{loop}}, we use the estimate of EM0 obtained in Step 1. The number density can then be estimated using {n}_{\mathrm{loop}}=\sqrt{{\mathrm{EM}}_{0}/V}, where V is the volume of the emitting plasma, found from X-ray imaging (see Figure, 3: #apjaafad3f3). (For SOL-05-13T02:12, the volume was estimated by applying the imaging algorithm Visibility Forward Fitting (VIS_FWDFIT
      • 2007: #apjaafad3bib42). Using VIS_FWDFIT, the volume of the coronal source varies with energy (see Jeffrey et al
        • Schmahl
      • c. The half-length L of the loop (the length from the coronal source to the chromospheric footpoint), is estimated directly from the image (see Figure, 3: #apjaafad3f3). (For SOL-05-13T02:12, we estimate from the X-ray image (blue curve in Figure, 3: #apjaafad3f3 left) that L=(24\pm 5)^{\prime\prime} \sim (17\pm 3.5) Mm.)
      • [3]
        Once the plasma parameters {n}_{\mathrm{loop}}, {T}_{\mathrm{loop}}, and L have been determined and fixed, then the remaining nonthermal parameters of the f_thick_warm fit, namely (\dot{{N}_{0}}, {\delta }_{\mathrm{low}}, EB, {\delta }_{\mathrm{high}}, Ec and EH), can be determined during 02:09:00-02:10:00 UT. As usual, we suggest that EH be fixed at the high default value, and that EB and {\delta }_{\mathrm{high}} are left free only if required (otherwise keep them fixed with EB above the highest energy in the fit). (For SOL-05-13T02:12, we find {E}_{c}=(12.1\,\pm 0.3) keV, \dot{{N}_{0}}=(12.8\pm 0.6)\times {10}^{35} electrons s-1 and {\delta }_{\mathrm{low}}\,=\delta =(4.25\pm 0.02).)
      • [4]
        The nonthermal electron power can now be determined as P=(\delta -1)/(\delta -2)\,\dot{{N}_{0}}{E}_{c}, as for the cold thick-target case (Equation (, 4: #apjaafad3eqn4)), as we assume a power-law form for the accelerated electrons. (For SOL-05-13T02:12, we find P\,=(35.8\pm 1.9)\times {10}^{27} erg s-1.)
      • [1]
        As L becomes smaller (decreasing from 40 to 5 Mm), Ec decreases from 20 to 5 keV and the power P correspondingly increases. This is because the smaller column density in the corona means that fewer accelerated electrons thermalize there. The estimated uncertainty in L of ±3.5 Mm leads to an uncertainty of ±2 keV in Ec that is not accounted for in OSPEX
        • [2]
          Decreasing the plasma number density {n}_{\mathrm{loop}} similarly lowers the coronal column density, hence the number of accelerated electrons thermalized there
          • [2]
            a similar variation in Ec hence results. Again, the estimated uncertainty in {n}_{\mathrm{loop}} of ±2.1 × 1010 cm-3 leads to an uncertainty of ±2 keV in Ec that is not accounted for in OSPEX
            • [3]
              Keeping the temperature, T, from f_vth fixed, but varying the f_thick_warm parameter {T}_{\mathrm{loop}} (\ne T) between {T}_{\mathrm{loop}}=4.0\,\mathrm{keV} and {T}_{\mathrm{loop}}=1.5\,\mathrm{keV}, decreases Ec from ~16 keV down to ~7 keV. This is because, once again, a lower value of {T}_{\mathrm{loop}} results in a smaller number of accelerated electrons blending into the thermal background, so that the nonthermal tail now extends to lower values of Ec
              • ADS:
                • Alaoui M.
                  ,
                • Holman G. D.
              • ADS:
                • Aschwanden M. J.
                  ,
                • Caspi A.
                  ,
                • Cohen C.M. S.
              • ADS:
                • Aschwanden M. J.
                  ,
                • Holman G.
                  ,
                • O'Flannagain A.
              • Crossref:
                • Bai T.
                  ,
                • Ramaty R.
              • ADS:
                • Battaglia M.
                  ,
                • Motorina G.
                  ,
                • Kontar E. P.
              • LRSP 5 1, Crossref:
                • Benz A. O.
              • ADS:
                • Bian N. H.
                  ,
                • Emslie A. G.
                  ,
                • Kontar E. P.
              • SoPh 18 489, Crossref:
                • Brown J. C.
              • ADS:
                • Brown J. C.
                  ,
                • Emslie A. G.
                  ,
                • Kontar E. P.
              • AGU Fall Meeting #SH44A-07, ADS:
                • Christe S.
                  ,
                • Shih A. Y.
                  ,
                • Krucker S.