Draft version October 15, 2024 Typeset using LT E X A preprint style in AASTeX631

## Flare-accelerated Electrons in the Kappa Distribution from X-Ray Spectra with the Warm-Target Model

Yingjie Luo ( 骆 英 杰 ), 1 Eduard P. Kontar, 1 and Debesh Bhattacharjee 1

1 School of Physics &amp; Astronomy, University of Glasgow, G12 8QQ, Glasgow, UK

(Published 2024 October 8)

## Submitted to APJ

## ABSTRACT

X-ray observations provide important and valuable insights into the acceleration and propagation of nonthermal electrons during solar flares. Improved X-ray spectral analysis requires a deeper understanding of the dynamics of energetic electrons. Previous studies have demonstrated that the dynamics of accelerated electrons with a few thermal speeds are more complex than those with significantly higher speeds. To better describe the energetic electrons after injection, a model considering energy diffusion and thermalization effects in flare conditions (the warm-target model) has recently been developed for spectral analysis of hard X-rays. This model has demonstrated how the low-energy cutoff, which can hardly be constrained in cold-target modeling, can be determined. However, the power-law form may not be the most suitable representation of injected electrons. The kappa distribution, which is proposed as a physical consequence of electron acceleration, has been applied successfully in RHESSI spectral analysis. In this study, we employ the kappa-form injected electrons in the warm-target model to analyze two M-class flares, observed by RHESSI and STIX, respectively. The best-fit results show that the kappa-form energetic electron spectrum generates lower nonthermal energy than the power-law form when producing a similar photon spectrum in the fit range. We also demonstrated that the fit parameters associated with the kappa-form electron spectrum can be well determined with small uncertainty. Further, the kappa distribution, which covers the entire electron energy range, enables the determination of key electron properties such as total electron number density and average energy at the flare site, providing valuable information on electron acceleration processes.

Keywords: Solar flares (1496); Solar physics (1476); Solar activity (1475); Solar x-ray flares (1816); Active solar corona (1988)

## 1. INTRODUCTION

Corresponding author: Yingjie Luo

Yingjie.Luo@glasgow.ac.uk

Solar flares are events with an explosive release of energy and emissions that span a wide range of electromagnetic wavelengths from radio to /u1D6FE -rays. In order to understand the acceleration mechanism andthetransport of energetic electrons, it is crucial to know where and how the energy is released (Benz 2017). Energetic electrons from solar flares can be indirectly observed through X-ray and radio emissions. Because of its unique advantages, hard X-ray (HXR) diagnostics is essential for determining the electron properties (see Holman et al. (2011); Kontar et al. (2011) for reviews). HXR emission observed from solar flares is dominantly due to electron-ion bremsstrahlung. The intensity of HXR emission is linearly proportional to the nonthermal electron density, so the observed HXR photon spectrum provides a straightforward relationship with the energetic electrons. Furthermore, HXR emission is less affected by propagation effects than radio emissions, resulting in minimal modification of the original electron properties during observations.

The warm-target model has been utilized to analyze RHESSI flares to more precisely fit the characteristics of nonthermal electrons (Kontar et al. 2019; Aschwanden et al. 2016, 2017). One significant chal-

Over the past two decades, advanced X-ray instruments, such as the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI; Lin et al. 2002) and the Spectrometer/ Telescope for Imaging Xrays (STIX; Krucker et al. 2020), have facilitated the collection of high-quality spectra and spectroscopic imaging. The frequent observation of X-ray sources in the corona and footpoints in the chromosphere (Sui et al. 2002; Emslie et al. 2003; Battaglia &amp; Benz 2006; Kontar et al. 2008; Shao &amp; Huang 2009) strongly supports the 'thick-target' interpretation (Brown 1971; Syrovatskii &amp; Shmeleva 1972), where flare-accelerated electrons propagate down to the dense chromosphere and deposit all their energy within the chromosphere. According to the thick-target model, the observed X-ray photon spectrum is determined by the injected electron spectrum, the bremsstrahlung cross section, and the energy loss rate within the target (Brown et al. 2003). When the energy of electrons is high compared to that of the target, i.e. /u1D438 ≫ /u1D458 B /u1D447 ( /u1D458 B is the Boltzmann constant, /u1D447 is the temperature), the model is referred to as a cold target. In cold-target conditions, the kinetic energy loss is predominantly due to Coulomb collisions (Brown 1971; Lin &amp; Hudson 1976), while other effects, such as return current ohmic losses (Knight &amp; Sturrock 1977; Emslie 1980; Zharkova &amp; Gordovskyy 2006; Holman 2012; Alaoui &amp; Holman 2017), could affect the results. While such a cold-target model is widely utilized as the basis of HXR spectral analysis, observations suggest that the coronal plasma during flares can reach high temperatures above 10 7 K, highlighting the need to improve the model. Studies have demonstrated that the behavior of electrons in the warm-target condition ( /u1D438 ∼ /u1D458 B /u1D447 ) is more complex and significantly different from what is predicted by the cold-target model. Emslie et al. (2003) highlighted a rapid drop in the energy loss rate in warm-target conditions, indicating a substantial overestimation of electron flux by the cold-target model in generating observed photon spectra. Additionally, collisional energy diffusion and thermalization (Galloway et al. 2005) have a significant impact on the dynamics of injected electrons. Furthermore, a physically self-consistent model requires consideration of both energy diffusion and spatial diffusion during electron transportation (Jeffrey et al. 2014). Taking into account these advances, Kontar et al. (2015, 2019) developed a warm-target model that provides a physically sound treatment of electron dynamics in the warm-target condition (referred to as the warm-target model hereafter). This model includes collisional energy diffusion, suggesting that injected electrons thermalize within the target rather than losing all their energy and disappearing. The HXR spectrum generated by the warm-target model includes both the contributions from the electrons in the coldtarget chromosphere and the the component from the corona.

lenge in inferring nonthermal electron properties is that the nonthermal emission is typically masked by the thermal emissions below 20-30 keV (Holman et al. 2011). This makes it difficult to determine the low-energy cutoff, /u1D438 /u1D450 , which is important in nonthermal electron energetics (see Aschwanden et al. 2019, for a review). This parameter and the closely related nonthermal power are the key parameters to characterize the acceleration mechanism. The cold-target condition /u1D438 ≫ /u1D458 B /u1D447 is typically applicable in the range above 20 keV, providing solely an upper limit for the low-energy cutoff /u1D438 /u1D450 . On the other hand, the warm-target model considers the thermalization of injected electrons, which enables the electron properties to be effectively constrained, because the spectra of thermalized electrons are proportional to the population of the injected electrons. Kontar et al. (2019) have demonstrated that using the warm-target model with a power-law distribution can significantly enhance the accuracy of determining the low-energy cutoff /u1D438 /u1D450 compared to the cold-target model.

In addition to not requiring an arbitrary energy cutoff /u1D438 /u1D450 , the adoption of a kappa-form injection offers several advantages over the power-law distribution. First, the kappa-form electron spectrum is believed to result from acceleration out of thermal equilibrium, and its parameters can provide crucial information on electron acceleration. Various studies have delved into the origin of the kappa distribution from different perspectives, including collisional relaxation along with wave-particle interaction (Ma &amp; Summers 1998; Vocks &amp; Mann 2003) and beam-plasma interactions (Yoon et al. 2006; Ryu et al. 2007; Rhee et al. 2006). Moreover, originating from the Fokker-Planck equation, Bian et al. 2014 have pointed out that the kappa electron distribution results from stochastic acceleration in the presence of Coulomb collisions and velocity diffusion, specifying insight into electron kinetics during solar flares. Similarly, simulations of electron acceleration during magnetic reconnection lead to a distribution resembling a kappa distribution (Arnold et al. 2021). Another advantage of the kappa distribution of the electron spectrum over the power-law distribution is its entire coverage of electron kinetic energies. The kappa distribution can provide electron information at low energy (a few keV), which is not sensitive to X-

However, it is important to recognize that the power-law form of the injected electron spectrum is not unique and not always the most suitable representation in many cases (Brown et al. 2006). The power-law distribution requires a low-energy cutoff ( /u1D438 /u1D450 ) to prevent an infinite number of electrons. The presence of the low-energy cutoff in the injected electron spectrum, while suggested to have physical significance in terms of the critical velocity of runaway and free acceleration (Holman &amp; Benka 1992), results in a nonphysical discontinuous electron distribution of a Maxwellian plus the power law. Around the low-energy cutoff /u1D438 /u1D450 , where the slope of the distribution is positive ( /u1D451/u1D453 ∕/u1D451/u1D463 &gt; 0 ), Langmuir waves are expected to generate and to grow (Vedenov &amp; Velikhov 1963; Drummond &amp; Pines 1964; Emslie &amp; Smith 1984; Hannah et al. 2009). The interaction between the waves will result in a flattened distribution around the low-energy cutoff (Krall &amp; Trivelpiece 1973), suggesting the injected (accelerated) electron distribution should be a nonincreasing function for all energies. On the other hand, the kappa distribution, consisting of a Maxwellian core and smoothly merged power-law tail, has been studied to represent accelerated electrons. The kappa distribution is often used to model the entire energetic electron spectrum from in-situ observations (Vandas 1989; Mace &amp; Hellberg 1995; Collier et al. 1996; Maksimovic et al. 2005; Imada et al. 2011) and has emerged as a promising alternative for describing the accelerated electrons in HXR spectral fitting (Kašparová &amp; Karlický 2009; Oka et al. 2013, 2015; Battaglia et al. 2015; Effenberger et al. 2017). Previous studies have shown that while spatially integrated spectra cannot be well fitted only with the electrons obeying a kappa distribution, kappa-form accelerated electrons with a prominent thermal component yield promising results.

ray instruments but can be observed through other diagnostics such as EUV observations. This allows for a more comprehensive study of the electron distribution through multiwavelength observations, providing additional constraints (Oka et al. 2015; Battaglia et al. 2015). Furthermore, due to its whole electron energy coverage and convergent nature, the kappa distribution allows for the determination of parameters such as electron number density and average electron energy, which is not achievable with the power-law distribution. The kappa distribution is closely related to the thermal properties of the ambient plasma. Parameters such as electron number density and average electron energy can be directly compared with the corresponding thermal properties of the preaccelerated plasma, serving as an important constraint and cross-check.

The paper investigates the spectrum of accelerated (injected) electrons in solar flares. Section 2 discusses the relationship between the kappa-form injected electron spectrum and the mean electron flux based on the warm-target model. We utilize the warm-target model with kappa-form injected electrons to analyze solar flares observed by RHESSI on 2011 February 24 and by STIX on 2022 March 28. The best-fit results and a comparison between the kappa and power-law forms of injection is provided in Section 3. Section 4 discusses the uncertainty of fitting parameters related to the kappa distribution with different approaches. The resulting characteristics such as electron number density and average electron energy are also shown. The findings and the physical significance of the derived kappa parameters are summarized and discussed in Section 5.

## 2. SOURCE-INTEGRATED SPECTRUM FROM KAPPA-FORM INJECTED ELECTRONS

Brown et al. (2003) highlighted the concept of the source-integrated density-weighted mean electron flux ⟨/u1D45B/u1D449 /u1D439 ⟩(/u1D438 ) [electrons cm -2 s -1 keV -1 ] and expressed the observed HXR spectrum /u1D43C (/u1D700 ) [photons cm -2 s -1 keV -1 ] as:

$$I ( \varepsilon ) = \frac { 1 } { 4 \pi R ^ { 2 } } \int _ { \varepsilon } ^ { \infty } Q ( \varepsilon, E ) \langle n V F \rangle ( E ) \, d E,$$

where /u1D445 is the distance between the observer and the Sun, and /u1D444(/u1D700, /u1D438 ) is angle-averaged bremsstrahlung cross section (see discussion in Kontar et al. 2011 for RHESSI). Determination of the mean electron flux from the observed photon spectrum requires only the well-studied bremsstrahlung cross section (Haug 1998). To obtain the injected electron rate spectrum ̇ /u1D441 (/u1D438 ) [electrons s -1 keV -1 ] from the obtained mean electron flux, a model for the dynamics of injected electrons is necessary. The commonly used coldtarget model establishes the relationship between mean electron flux and injected electron rate as

$$\langle n V F \rangle ( E ) = \frac { E } { K } \int _ { E } ^ { \infty } \dot { N } ( E _ { 0 } ) \, d E _ { 0 } \,,$$

were /u1D43E = 2/u1D70B/u1D452 4 ln(Λ) is the collision parameter, /u1D452 [esu] is the elementary charge and /u1D459 /u1D45B (Λ) is the Coulomb logarithm (Spitzer 1962). To account for the electron dynamics in the warm-target corona and cold chromosphere, Kontar et al. (2015) included the energy diffusion, transport, and thermalization of the accelerated electrons, rewriting Equation (2) as

$$\langle n V F \rangle ( E ) = \frac { 1 } { 2 K } E \, e ^ { - E / k _ { B } T } \, \int _ { E _ { m i n } } ^ { E } \frac { e ^ { E ^ { \prime } / k _ { B } T } } { E ^ { \prime } G ( \sqrt { E ^ { \prime } / k _ { B } T } ) } \, d E ^ { \prime } \times \int _ { E ^ { \prime } } ^ { \infty } \dot { N } ( E _ { 0 } ) \, d E _ { 0 } \,.$$

Kontar et al. (2019) further simplify equation (3) and rewrite it as:

$$\intertext { a l e a l, \, \infty \, \infty \, \infty \, \infty \, \infty \, \ e m p n y \, c q u a t o n \, \vee \, \infty \, \ e w i n e \, \ e a s. } \langle n V F \rangle ( E ) & \approx \Delta E M \sqrt { \frac { 8 } { \pi m _ { e } } } \frac { E } { ( k _ { B } T ) ^ { 3 / 2 } } e ^ { - E / k _ { B } T } + \frac { E } { K } \int _ { E } ^ { \infty } \grave { N } ( E _ { 0 } ) \, d E _ { 0 } \,, \\ \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

where ∆/u1D438/u1D440 , referred to as the 'thermalized emission measure', is due to the thermalization of accelerated electrons and can be explicitly written as

$$\text{com in and can } \in \text{expify when as} \\ \Delta & \text{EM} \approx \frac { \pi } { K } \sqrt { \frac { m _ { e } } { 8 } } ( k _ { B } T ) ^ { 2 } \frac { \dot { N } _ { 0 } } { E _ { \min } ^ { 1 / 2 } } \,, \quad \text{where} \quad E _ { \min } \approx 3 k _ { B } T ( \frac { 5 \lambda } { L } ) ^ { 4 } \,, \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

and ̇ /u1D441 0 [electrons s -1 ] = ∫ ∞ ̇ (/u1D438) /u1D451/u1D438 0 /u1D441 is the total electron injection (acceleration) rate and /u1D706 = (/u1D458 B /u1D447 ) 2 ∕2/u1D43E/u1D45B is the collisional mean free path. The half-loop length /u1D43F , number density /u1D45B , and temperature /u1D447 are the parameters associated with the coronal part of the loop. Equation 4 consists of two parts, referred to as the thermalized and nonthermal components, corresponding to the warm-target (coronal) and cold-target conditions, respectively. It is worth noting that the thermalized component is proportional to the total injected electron rate ̇ /u1D441 0 , which is why this model is effective in constraining the nonthermal component when fitting the observed X-ray spectrum. This warm-target model 'f\_thick\_warm' , which utilizes Equations 4 and 5, is implemented in the Solar Software (SSW) and OSPEX package (Schwartz et al. 2002; Tolbert &amp; Schwartz 2020). The function 'f\_thick\_warm' utilizes the commonly employed broken power-law form of injected electrons, and it is worth noting that the warm-target model (Equations 4 and 5) can accept any form of electron injection spectrum. As mentioned earlier, the use of a kappa distribution for injection is a promising alternative and has been applied to fit the RHESSI spectra.

$$f _ { k } ( v ) = \frac { n _ { k } } { \pi ^ { 3 / 2 } v _ { \text{te} } ^ { 3 } \kappa ^ { 3 / 2 } } \frac { \Gamma ( \kappa ) } { \Gamma ( \kappa - 3 / 2 ) } ( 1 + \frac { v ^ { 2 } } { \kappa v _ { \text{te} } ^ { 2 } } ) ^ { - \kappa } \,, \quad \quad \quad$$

Bian et al. 2014 derive the kappa distribution as a result of stochastic acceleration in collisional plasma:

where /u1D463 te is the thermal speed ( 1∕2/u1D45A /u1D452 /u1D463 2 te = /u1D458 B /u1D447 /u1D705 ), and /u1D705 is the kappa index, which importantly represents the ratio between the timescales of acceleration /u1D70F acc (/u1D463) and collisional deceleration /u1D70F /u1D450 (/u1D463) , /u1D705 = /u1D70F acc (/u1D463) 2/u1D70F /u1D450 (/u1D463) . While different formats of the kappa distribution exist (see, e.g. Kašparová &amp; Karlický 2009; Oka et al. 2013), they are identical in essence and differ only in their format and use of parameters. Furthermore, the physical interpretation of /u1D705 varies when additional factors are taken into account. In this study, we have chosen to adopt Equation 6 for the kappa distribution.

Therelationship between the injected electron rate spectrum and the electron distribution for velocity is ̇ /u1D441 (/u1D438 )/u1D451/u1D438 = /u1D434/u1D463/u1D453 /u1D458 (/u1D463)/u1D451 3 /u1D463 , where /u1D434 is the injection area. For an isotropic electron distribution ( /u1D451 3 /u1D463 = 4/u1D70B/u1D463 2 /u1D451/u1D463 ), the injected electron rate spectrum becomes (see red curves in Figure 1):

√

$$\dot { N } ( E ) = A \frac { n _ { k } \Gamma ( \kappa ) } { \Gamma ( \kappa - 3 / 2 ) \kappa ^ { 3 / 2 } } \sqrt { \frac { 8 } { \pi m _ { e } k _ { B } T _ { \kappa } } } \frac { E / k _ { B } T _ { \kappa } } { ( 1 + E / \kappa k _ { B } T _ { \kappa } ) ^ { \k$$

where /u1D45B /u1D458 [cm -3 ] = ∫ /u1D453 /u1D458 (/u1D463) /u1D451 3 /u1D463 is the total electron number density and Γ(/u1D465) is the gamma function. The total injection rate becomes

$$\alpha \text{ mppound raw events} \\ \dot { N } _ { 0 } = \int _ { 0 } ^ { \infty } \dot { N } ( E ) \, d E = 2 A n _ { k } \sqrt { \frac { 2 k _ { B } T _ { \kappa } } { m _ { e } } } \frac { \kappa ^ { 1 / 2 } } { ( \kappa - 2 ) B ( \kappa - 3 / 2, 1 / 2 ) },$$

Figure 1. Top Panel: the electron injection rate spectrum from distributions of the power-law form and kappa form. For the case shown on the left, the power-law form of injection has parameters of ̇ /u1D441 0 = 1 × 10 10 electrons s -1 , /u1D438 /u1D450 = 20 keV , and /u1D6FF = 4 , and the kappa form of injection is set to ̇ /u1D441 0 = 6.8 × 10 10 electrons s -1 , /u1D458 B /u1D447 /u1D705 = 1.5 keV , and /u1D705 = /u1D6FF +1 = 5 . The selection is for satisfying ̇ /u1D441 /u1D705 (/u1D438 1 ) = /u1D441 ̇ /u1D440 /u1D44E/u1D465 (/u1D438 1 ) , when /u1D438 1 = 100 keV . The case shown in the right panel is with the power-law parameters: ̇ /u1D441 0 = 1×10 10 electrons s -1 , /u1D438 /u1D450 = 10 keV , /u1D6FF = 3 , condition of a high proportion of nonthermal electrons. The same prescription is used to obtain the kappa-form electron spectrum: ̇ /u1D441 0 = 1.17 × 10 10 electrons s -1 , /u1D458 B /u1D447 /u1D705 = 1.5 keV , and /u1D705 = 4 . Bottom Panels: the photon spectrum, based on the warm-target model (WTM), generated by the electron injection in the top panel with a coronal condition of /u1D43F loop = 30Mm /u1D45B , loop = 5 × 10 10 cm -3 , and /u1D458 B /u1D447 loop = 1.5 keV . Photon spectra from powerlaw and kappa form of injection are plotted in black and red, respectively. Solid curves represent the total photon spectrum, while the dashed and dotted-dashed curves represent the 'nonthermal' and 'thermalized' components, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-c5bc82b6-e630-4880-8280-1acd142069c7.jpg)

where /u1D435(/u1D465, /u1D466 ) is the beta function. Equation 7 can be expressed as

$$\dot { N } ( E ) = \frac { \dot { N } _ { 0 } } { k _ { B } T _ { \kappa } } \frac { ( \kappa - 1 ) ( \kappa - 2 ) } { \kappa ^ { 2 } } \frac { E / k _ { B } T _ { \kappa } } { ( 1 + E / \kappa k _ { B } T _ { \kappa } ) ^ { \kappa } }.$$

According to Equation 4, the mean electron flux from a kappa-form injected electrons is given by (Equation 17 in the Appendix of Kontar et al. 2019):

√

$$& \langle n V F \rangle ( E ) \approx \Delta E M \sqrt { \frac { 8 } { \pi m _ { e } } } \frac { E } { ( k _ { B } T ) ^ { 3 / 2 } } e ^ { - E / k _ { B } T } + \dot { N } _ { 0 } \frac { E } { K } \frac { 1 + ( 1 - 1 / \kappa ) E / k _ { B } T _ { \kappa } } { ( 1 + E$$

For high-energy electrons, the kappa distribution can be well approximated by a power-law distribution with a spectral index /u1D705 -1 . Figure 1 shows that the two distributions have an identical injection rate near 100 keV. However, in the energetically important deka-keV range, the kappa distribution is distinct from a power-law distribution, and thus a different photon spectrum emerges. Figure 1 illustrates the electron injection rate (top panels, red and black for kappa and power-law injection, respectively)

and the resulting photon spectra (bottom panels) generated by the warm-target model in a typical flare coronal condition ( /u1D43F loop = 30 Mm /u1D45B , loop = 5 × 10 10 cm -3 , and /u1D458 B /u1D447 loop = 1.5 keV ).

## 3. APPLICATION OF THE WARM-TARGET MODEL TO OBSERVATIONS

In this section, we apply the warm-target model to analyze two well-observed flare events: one on 2011 February 24 observed by RHESSI, and one on 2022 March 28 observed by STIX. Both events are situated in the limb and provide spectroscopic images with two footpoint and loop-top sources. The warm-target model with kappa-form electron injection (function f\_thick\_warm\_kappa , based on Equation 10) encompasses parameters of the kappa distribution (total electron injection rate ̇ /u1D441 0 , kappa temperature /u1D447 /u1D705 , and kappa index /u1D705 ) and thermal properties of the target coronal loop (the half-loop length /u1D43F , number density /u1D45B loop , and temperature /u1D447 loop ).

In this study, the thermal properties of the target plasma are determined by the HXR observations of the thermal component. Following the approach outlined in Kontar et al. (2019), the thermal component of the X-ray spectra prior to the HXR burst can provide the plasma temperature ( /u1D447 ) and emission measure ( /u1D438 /u1D440 ). Moreover, the HXR spectroscopic imaging can estimate the coronal loop length ( /u1D43F ) and the size of the coronal X-ray loop source for determining the electron number density ( /u1D45B th ). While this method is straightforward, we recognize that other methods, such as differential emission measure analysis (DEM; see algorithms from Craig &amp; Brown 1976; Hannah &amp; Kontar 2012; Cheung et al. 2015) from the EUV flare context or combining X-ray data with EUV data (Inglis &amp; Christe 2014; Battaglia et al. 2015), can yield thermal properties (Battaglia &amp; Kontar 2013; Battaglia et al. 2015).

## 3.1. 2011 February 24 RHESSI flare

Figure 2. Left panel: RHESSI light curves detailing a GOES-class M3.5 solar flare on 2011 February 24. The blue, black, and red shaded regions denote the chosen time frames for background, preburst, and HXR burst, respectively, for the subsequent spectral fit. Right panel: RHESSI contours (50, 70, 90% of the maximum; the time range used for imaging is the red shaded region in the left panel) overlaid on the SDO/AIA 193 Å flare context (inverted grayscale). A cyan-colored semicircular loop (half-loop length /u1D43F ) that traverses through two footpoints and the loop-top source is plotted. The 50% contour of the loop-top source also provides the characteristic lengths ( /u1D459 = 2/u1D45F and /u1D451 ) necessary for calculating the size of the loop-top source.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-84d13457-54ba-4cc6-bb41-afc18fc59766.jpg)

The first flare we study is a GOES-class M3.5 flare located on the eastern limb from the Earth's perspective. The flare occurred on 2011 February 24 and the GOES soft X-ray flux peaked at ∼ 07:35 UT

Figure 3. Left panel: fit photon spectrum of the 2011 February 24 flare. The fit was done twice, using functions f\_vth (isothermal, red dashed curve) and f\_thick\_warm f\_thick\_warm\_kappa / (warm-target model with power-law and kappa forms of electron injection, green and blue dashed curves, respectively). The observed background-subtracted spectrum is shown in black (times of burst and background are shown as red and blue shaded regions in Figure 2 left panel), and the fit range spans from 6 to 100 keV. Normalized residuals, which are the difference between observed (black curve) and fit (green and blue solid curves for different fits) photons divided by error, are also displayed in the bottom. Right panel: the photon spectra of the isothermal (red) and warm-target models (power-law and kappa forms of injection are in green and blue, respectively). The nonthermal and thermal components (see Equation 4) of the warm-target model are presented as solid and dashed curves, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-fd3dad47-bfb3-4cf9-b0e2-2ec37fe34f1f.jpg)

(see Battaglia &amp; Kontar 2011, 2012; Stackhouse &amp; Kontar 2018, for previous studies of this event). Our attention is drawn to the first HXR peak at approximately 07:30 UT (red shaded time range in Figure 2 left panel). The HXR spectroscopic imaging of the burst displays a 'classical' configuration of the loop-top source and two HXR footpoints (Figure 2, right panel). The X-ray observations from RHESSI of the selected HXR burst benefit from the low/negligible pileup effect. Furthermore, RHESSI was in optimal working condition during the flare, providing more functional detectors.

For spectral analysis, we first gather the thermal properties of the target coronal loops according to the images of the burst (07:30:00 to 07:30:44 UT, Figure 2, left panel red shaded time range) and spectra before the interval analyzed (short time interval just before the burst analysed, refer as 'preburst' hereafter, in this case 07:29:52 to 07:30:00 UT, Figure 2 left panel, black shaded time range). In the analysis of the preburst X-ray spectra, we utilize functions f\_vth (isothermal) and f\_thick2 (cold thick target). At this stage, we focus solely on the thermal plasma parameters and can overlook the constraints associated with using the cold-target model. f\_vth from the best-fit result of the preburst spectra provides us with the plasma temperature ( /u1D458 B /u1D447 0 = 1.30 keV ) and emission measure ( /u1D438 /u1D440 0 = 0.127 × 10 49 cm -3 ) before the injection/acceleration of nonthermal electrons. Additionally, we use X-ray spectroscopic imaging (Figure 2 right panel) to estimate the size of the X-ray source and target coronal loop. We determine the source size in various dimensions (see pink and red solid lines in Figure 2 right panel, /u1D45F = /u1D459∕2 = 3.7 Mm

Figure 4. Left panel: the injected electron rate spectrum ̇ /u1D441 (/u1D438 ) inferred from the warm-target fit results (see Table 1 for fitted parameters). The green and blue curves correspond to the power-law ( /u1D438 /u1D450 = 7.13 keV ) and kappa forms of electron injection, respectively. Right panel: the distribution of electrons for the energy /u1D453 /u1D458 (/u1D438 ) ( /u1D453 /u1D458 (/u1D438 )/u1D451/u1D438 = /u1D453 /u1D458 (/u1D463)/u1D451 3 /u1D463 ) of the fitted kappa distribution (blue curve, total electron number density /u1D45B /u1D458 = 1.67 × 10 9 cm -3 , from Equation 8) and preburst thermal Maxwellian distribution (red curve, total electron number density /u1D45B /u1D458 = /u1D45B /u1D461 ℎ = 4.4 × 10 10 cm -3 ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-e3ab82b8-2366-46f7-b063-7b7b30477733.jpg)

Table 1. The parameters and fit results of 2011 February 24 M3.5 flare using the warm-target model. The isothermal component and thermal properties of the warm-target model are fixed, and properties associated with the accelerated electron spectrum are set free. The reduced /u1D712 2 and nonthermal electron power ( /u1D443 ) are also listed. For the convenience of comparison, the nonthermal power produced by the kappa-form electron spectrum above the power-law low-energy cutoff /u1D438 /u1D450 is also shown.

| Parameters                            | vth+WTM(power-law)   | vth+WTM(kappa)                             |       |
|---------------------------------------|----------------------|--------------------------------------------|-------|
| EM[ 10 49 cm -3 ]                     | 0.127                | 0.127                                      | Fixed |
| /u1D458 B /u1D447 [keV]               | 1.30                 | 1.30                                       | Fixed |
| /u1D45B loop [ 10 10 cm -3 ]          | 4.4                  | 4.4                                        | Fixed |
| /u1D458 B /u1D447 loop [keV]          | 1.30                 | 1.30                                       | Fixed |
| /u1D43F [Mm]                          | 15.8                 | 15.8                                       | Fixed |
| ̇ /u1D441 0 [ 10 35 electrons sec -1 ] | 20.3 ± 0.8           | 22.5 ± 0.8                                 | Free  |
| /u1D6FF                               | 3.77 ± 0.02          | N/A                                        | Free  |
| /u1D705                               | N/A                  | 5.14 ± 0.04                                | Free  |
| /u1D438 /u1D450 [keV]                 | 7.13 ± 0.19          | N/A                                        | Free  |
| /u1D458 B /u1D447 /u1D705 [keV]       | N/A                  | 1.29 ± 0.04                                | Free  |
| /u1D712 2                             | 1.93                 | 2.11                                       | N/A   |
| Power [ 10 28 erg∕s ]                 | 3.62                 | 2.23 for total; 1.41 above /u1D438 /u1D450 | N/A   |

and /u1D451 = 15.3 Mm , respectively) by considering the 50% contour of X-ray loop-top sources. The crossarea can be calculated as /u1D434 = /u1D70B/u1D45F 2 = 4.21 × 10 10 cm 2 , and the source volume can be obtained through the cylindrical formula: /u1D449 = /u1D434 ×/u1D451 = /u1D70B/u1D45F 2 /u1D451 = 6.46 × 10 26 cm 3 . Thus the thermal electron number density can be estimated using /u1D45B loop = √ /u1D438 /u1D440 0 ∕/u1D449 = 4.4 × 10 10 cm -3 . To obtain the coronal loop length, we construct a semicircular loop (see cyan dashed line in Figure 2 right panel) that connects two footpoints

and passes through the loop-top source. Consequently, we determine the thermal properties: loop temperature ( /u1D458 B /u1D447 loop = /u1D458 B /u1D447 0 = 1.30 keV ), loop number density ( /u1D45B loop = √ /u1D438 /u1D440 0 ∕/u1D449 = 4.4 × 10 10 cm -3 ), and the the half-loop length ( /u1D43F = 15.8 Mm ).

Thefits using the warm-target model, with kappa and power-law forms of electron injection, both produce favorable fit results with reasonable parameters. With the systematic uncertainty of 0.02, which is the default for RHESSI in OSPEX , the reduced /u1D712 2 is 1.93 and 2.11 for the fits with kappa and power-law forms of electron injection, respectively. From Figure 3 bottom left panel, the significant normalized residuals remain below 12 keV, where thermal bremsstrahlung dominates. This may arise from the use of the isothermal assumption f\_vth for the thermal emission (Jeffrey et al. 2015). From the photon spectrum made by an individual component (Figure 3 right panel), we found that for this HXR burst, the emission measure contributed by the thermalized injected electrons ( ∆/u1D438/u1D440 , dashed curves) is of comparable magnitude to the isothermal component ( /u1D438 /u1D440 , red solid curve). Here, we note that although the ratio of ∆/u1D438/u1D440∕/u1D438/u1D440 is affected by the preburst timing chosen, all reasonable selections of preburst timing can yield acceptable fit results. The different ratios only arise from the identification of preaccelerated plasma.

Once the thermal properties are determined, we proceed to fit the spectra of the burst time interval (07:30:00 to 07:30:44 UT) at 6-100 keV using functions f\_vth and f\_thick\_warm\_kappa . In the function f\_vth , the temperature and emission measure are fixed according to the preburst spectra ( /u1D447 = /u1D447 0 , /u1D438 /u1D440 = /u1D438/u1D440 0 ). The warm-target algorithm proposes that the emission measure during the burst interval is a combination of the emission measure from the preburst and the emission measure of thermalized injected electrons ( /u1D438 /u1D440 = /u1D438/u1D440 0 +∆/u1D438/u1D440 ; the photon spectra associated with /u1D438 /u1D440 0 and ∆/u1D438/u1D440 are plotted as the red solid curve and the dashed curves in Figure 3 right panel). Since the plasma is slowly changing during the 40 s interval (electron diffusion time along the loop ∼ 56 s, use Equation 11 of Kontar et al. 2019), the parameters in f\_thick\_warm\_kappa related to thermal properties of the loop are fixed in the fit (the half-loop length /u1D43F , number density /u1D45B loop , temperature /u1D447 loop , and the relative elemental abundances, which are 1 by default; Landi et al. 2013). The parameters of the kappa-form injection spectrum (total electron injection rate ̇ /u1D441 0 , kappa temperature /u1D447 /u1D705 , and kappa index /u1D705 , except for the high-energy cutoff /u1D438 high fixed to 10 4 keV ) are determined from the fit. Moreover, a fit with a single power-law electron injection in the warm-target model 'f\_thick\_warm' was performed for comparison: the total injection rate ̇ /u1D441 0 , power-law index /u1D6FF , and low-energy cutoff /u1D438 /u1D450 are free, but the break energy /u1D438 break and high-energy cutoff /u1D438 high are fixed to 10 4 keV . The best-fit parameters (see Table 1) and generated photon spectra are plotted in Figure 3 (blue and green curves for warm-target model with kappa and power-law forms of injection, respectively, same color scheme for the normalized residuals in the bottom panel).

In Figure 4, left panel, we plot the injected electron rate spectrum (blue and green curves represent the best-fit kappa and power-law forms of injection, respectively). We noticed that the power-law and kappaformsofelectron injection exhibit similarly between 40 and 100 keV but behave differently below approximately 30 keV. At the power-law low-energy cutoff /u1D438 /u1D450 , the electron injection rate ̇ /u1D441 (/u1D438 /u1D450 ) of the kappa distribution is almost one order of magnitude smaller ( 1.11 vs. 7.88 ×10 35 electrons sec -1 keV -1 for the kappa and power-law forms of injection, respectively). Additionally, the fitted injected electron spectrum in the kappa distribution decreases more rapidly ( /u1D705 &gt; /u1D6FF + 1 ) at higher energy. Although the best-fit electron spectra in kappa and power-law forms have similar total injection rates ( ̇ /u1D441 0 = 22.5 vs. 20.3 ×10 35 electrons sec -1 for kappa and power-law forms of injection, respectively), the power-

law electron spectrum produces a more substantial nonthermal power than the kappa-form electron spectrum (see Table 1, 3.62 vs 2.23 ×10 28 erg∕s , /u1D443 /u1D705 /u1D443 /u1D45D/u1D45C/u1D464 /u1D452/u1D45F-/u1D459 /u1D44E/u1D464 ≈ 0.61 ). Moreover, for kappa form electron spectrum, the nonthermal power produced by the electrons above the power-law low-energy cutoff /u1D438 /u1D450 is 1.41 ×10 28 erg∕s . The equations used for calculating the power are listed below. For power-law distribution:

$$P = \int _ { E _ { c } } ^ { \infty } E ^ { \prime } N ( E ^ { \prime } ) \, d E ^ { \prime } = \frac { \delta - 1 } { \delta - 2 } E _ { C } \dot { N } _ { 0 }.$$

For the kappa distribution:

$$P & = \int _ { E _ { 1 } } ^ { E _ { 2 } } E ^ { \prime } N ( E ^ { \prime } ) \, d E ^ { \prime } = \frac { \dot { N } _ { 0 } } { ( k _ { B } T _ { \kappa } \kappa ) ^ { 2 } ( \kappa - 3 ) } ( k _ { B } T _ { \kappa } \kappa + E ^ { \prime \prime } ) \times \\ & \quad ( \frac { E ^ { \prime \prime } } { k _ { B } T _ { \kappa } \kappa } + 1 ) ^ { - \kappa } [ 2 ( k _ { B } T _ { \kappa } \kappa ) ^ { 2 } + 2 k _ { B } T _ { \kappa } ( \kappa - 1 ) \kappa E ^ { \prime \prime } + ( \kappa ^ { 2 } - 3 \kappa + 2 ) E ^ { \prime \prime 2 } ] \Big | _ { E _ { 1 } } ^ { E _ { 2 } }. \\ \tilde { \omega } \dots$$

Thus

$$P = \int _ { 0 } ^ { \infty } E ^ { \prime } N ( E ^ { \prime } ) \, d E ^ { \prime } = \frac { 2 \dot { N } _ { 0 } k _ { \text{B} } T _ { \text{x} } \kappa } { \kappa - 3 } ;$$

$$P & = \int _ { E _ { c } } ^ { \infty } E ^ { \prime } N ( E ^ { \prime } ) \, d E ^ { \prime } = \frac { \dot { N } _ { 0 } } { ( k _ { B } T _ { \kappa } \kappa ) ^ { 2 } ( \kappa - 3 ) } ( k _ { B } T _ { \kappa } \kappa + E _ { c } ) \times \\ & \quad ( \frac { E _ { c } } { k _ { B } T _ { \kappa } \kappa } + 1 ) ^ { - \kappa } [ 2 ( k _ { B } T _ { \kappa } \kappa ) ^ { 2 } + 2 k _ { B } T _ { \kappa } ( \kappa - 1 ) \kappa E _ { c } + ( \kappa ^ { 2 } - 3 \kappa + 2 ) E _ { c } ^ { 2 } ]. \\ \intertext { m a t i n d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d }$$

We also include a distribution plot ( /u1D453 /u1D458 (/u1D438 ) , /u1D453 /u1D458 (/u1D438 )/u1D451/u1D438 = ⟨/u1D453 /u1D458 (/u1D463)⟩/u1D451 3 /u1D463 ) of the accelerated kappa-form electrons (blue curve in Figure 4 right panel, /u1D45B /u1D458 = 1.67 × 10 9 cm -3 , calculated from Equation 8; here /u1D434 = 4.21 × 10 10 cm 2 ). The Maxwellian thermal distribution, which is taken from the isothermal component f\_vth (total electron number density /u1D45B /u1D458 = /u1D45B th = 4.4 × 10 10 cm -3 , /u1D458 B /u1D447 = 1.30 keV), is also shown as a red curve for comparison. We will provide a more detailed discussion of the parameters obtained in the next section.

## 3.2. 2022 March 28 STIX flare

We conducted the same spectral analysis on the other well-observed GOES-class M4.0 flare on 2022 March 28. GOES soft X-ray flux peaks at ∼ 11:29 UT; the flare-associated HXR bursts were detected by STIX on board the Solar Orbiter (Sol-O; Müller et al. 2020). The STIX temporal profile is shown in Figure 5, left panel; with the time difference of 335.7 seconds between Sol-O and the Earth has been applied to the STIX data. The flare site is situated in the eastern limb from Sol-O's perspective. The ExtremeUltraviolet Imager (EUI; Rochus et al. 2020)/Full Sun Imager (FSI) provides EUV context from the same viewpoint (Figure 5, right panel). This flare comprises several HXR bursts, and our focus is on the burst occurring at approximately 11:21 UT (red shaded time range in Figure 5, left panel).

The selected burst and preburst time ranges for this flare are 11:21:24-11:21:44 UT and 11:21:1211:21:17 UT (Figure 5 left panel, red and black shaded regions), respectively. For this particular event,

Figure 5. Left panel: STIX light curves (time has been shifted to the estimated timing observed on Earth) detailing a GOES-class M4.0 solar flare peaking at ∼ 11:29 UT on 2022 March 28. The blue, black, and red shaded regions denote the chosen time frames for background, preburst, and HXR burst we study, respectively, for the subsequent spectral fit. Right panel: STIX contours (50, 70, 90% of the maximum; the time range used for imaging is the red shaded region in the left panel) overlaid on the EUI/FSI 174 Å flare context (inverted grayscale, using the image at 11:16 UT instead of 11:26 UT to avoid the overexposure of the flare site). A cyan-colored semielliptical loop (half-loop length /u1D43F ) that traverses through two footpoints and the loop-top source is plotted. The 50% contour of the two loop-top sources provides the characteristic lengths ( /u1D459 = 2/u1D45F and /u1D451 , northern and southern loop-top sources marked with '1' and '2', respectively) necessary for calculating the size of the loop-top source.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-142007ee-73ed-41fb-8830-12dc34e46e83.jpg)

bottom pixels have noticeable larger counts than the top pixels. Thus, the spectrum for fitting and spectroscopic imaging relies solely on the data from the bottom pixels. The fit range is set to be 6 50 keV. The preburst spectra yield values of temperature and emission measure of /u1D458 B /u1D447 0 = 1.67 keV and /u1D438 /u1D440 0 = 0.0532 × 10 49 cm -3 . The right panel of Figure 5 displays X-ray spectroscopic images using the CLEAN algorithm (Högbom 1974), which are overlaid on EUI/FUI 174 Å flare context. Two looptop sources (red contours, 8-12 keV), of which the southern one existed before the flare (See Figure 2 in Purkhart et al. 2023, for detailed flare evolution), and two footpoint sources (blue contours, 20-50 keV) are observed. According to the 50% contour of the two loop-top sources, we obtain the source size in different dimensions ( /u1D45F 1 = /u1D459 1 ∕2 = 3.4 Mm , /u1D45F 2 = /u1D459 2 ∕2 = 3.2 Mm , pink lines, and /u1D451 1 = 12.1 Mm , /u1D451 1 = 9.6 Mm , red lines). The cross-area (using the northern source, /u1D434 = /u1D70B/u1D45F 2 1 = 3.64 × 10 10 cm 2 ) and source volume ( /u1D449 = /u1D70B/u1D45F 2 1 /u1D451 1 +/u1D70B/u1D45F 2 2 /u1D451 2 = 7.59 × 10 26 cm 3 ) can be calculated accordingly. The constructed semi-elliptical coronal loop (cyan dashed line in Figure 5 right panel), which is rooted at two HXR footpoint sources and passes through the loop-top sources, features a half-loop length of /u1D43F = 18.4 Mm . As a result, we can determine the target loop's thermal properties ( /u1D458 B /u1D447 loop = /u1D458 B /u1D447 0 = 1.67 keV and /u1D45B loop = √ /u1D438 /u1D440 0 ∕/u1D449 = 2.4 × 10 10 cm -3 ).

The best-fit results using f\_vth and f\_thick\_warm\_kappa f\_thick\_warm / are shown in Figure 6 and Table 2. To account for STIX calibration uncertainty, which is still unknown, 6% systematic error was added in OSPEX following Jeffrey et al. (2024). From the left panel of Figure 7, a distinct difference between the kappa and power-law forms of the injection electron spectrum was noted below 30 keV, consistent with the findings from the 2011 February 24 flare. For this event, the low-energy cutoff /u1D438 /u1D450 obtained from f\_thick\_warm is 12.46 keV. For comparison, the injection rate ̇ /u1D441 (/u1D438 = /u1D438 /u1D450 ) is 0.641 vs. 4.21 ×10 35 electrons sec -1 keV -1 for kappa and power-law forms of injection, respectively. The

Figure 6. Left panel: fit photon spectrum of the 2022 March 28 M4.0 flare. As in Figure 3, the fit was done twice, using functions f\_vth (isothermal, red dashed curve) and f\_thick\_warm f\_thick\_warm\_kappa / (warmtarget model with power-law and kappa forms of electron injection, green and blue dashed curves, respectively). The observed background-subtracted spectrum is shown in black (times of burst and background are shown as red and blue shaded regions in Figure 5 left panel), and the fit range spans from 6 to 50 keV. Normalized residuals are displayed at the bottom. Right panel: the photon spectrum of the isothermal (red) and warm-target models (power-law and kappa injection, in green and blue, respectively). The nonthermal and thermal components of the warm-target model are presented as solid and dashed curves, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-a3fb6de4-cb78-4af1-8b0d-1e115ac0aeab.jpg)

Figure 7. Left panel: the injected electron rate spectrum ̇ /u1D441 (/u1D438 ) inferred from warm-target fit results (see Table 2 for fitted parameters). The green and blue curves correspond to the power-law ( /u1D438 /u1D450 = 12.46 keV ) and kappa forms of electron injection spectrum, respectively. Right panel: the distribution of electrons for the energy /u1D453 /u1D458 (/u1D438 ) of the fitted kappa distribution (blue curve, total electron number density /u1D45B /u1D458 = 1.1 × 10 9 cm -3 , from Equation 8) and the ambient thermal plasma Maxwellian distribution (red curve, total electron number density /u1D45B /u1D458 = /u1D45B th = 2.4 × 10 10 cm -3 ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-79d07132-ca11-436c-9237-f253e90fd51e.jpg)

Table 2. The best-fit parameters and results of 2022 March 28 M4.0 flare using the warm-target model. The isothermal component and thermal properties of the warm-target model are fixed, and properties associated with the accelerated electrons are set free. The reduced /u1D712 2 and nonthermal electron power /u1D443 (for kappa-form electrons, the total power and power above power-law low-energy cutoff /u1D438 /u1D450 ) are also listed.

| Parameters                            | vth+WTM(power-law)   | vth+WTM(kappa)                             |       |
|---------------------------------------|----------------------|--------------------------------------------|-------|
| EM[ 10 49 cm -3 ]                     | 0.0532               | 0.0532                                     | Fixed |
| /u1D458 B /u1D447 [keV]               | 1.67                 | 1.67                                       | Fixed |
| /u1D45B loop [ 10 10 cm -3 ]          | 2.4                  | 2.4                                        | Fixed |
| /u1D458 B /u1D447 loop [keV]          | 1.67                 | 1.67                                       | Fixed |
| /u1D43F [Mm]                          | 18.4                 | 18.4                                       | Fixed |
| ̇ /u1D441 0 [ 10 35 electrons sec -1 ] | 13.6 ± 1.3           | 20.5 ± 1.4                                 | Free  |
| /u1D6FF                               | 4.86 ± 0.05          | N/A                                        | Free  |
| /u1D705                               | N/A                  | 8.85 ± 0.41                                | Free  |
| /u1D438 /u1D450 [keV]                 | 12.46 ± 0.45         | N/A                                        | Free  |
| /u1D458 B /u1D447 /u1D705 [keV]       | N/A                  | 2.97 ± 0.17                                | Free  |
| /u1D712 2                             | 0.91                 | 1.06                                       | N/A   |
| Power [ 10 28 erg∕s ]                 | 3.66                 | 2.95 for total; 1.51 above /u1D438 /u1D450 | N/A   |

kappa-form electron distribution, like that for the 2011 February 24 flare, will decrease more rapidly at higher energy ( /u1D705 &gt; /u1D6FF + 1 ). The nonthermal power generated by the injected electron spectrum ( 3.66 vs 2.95 × 10 28 erg∕s , for electrons in power-law and kappa forms of distribution, respectively) also exhibits significant dissimilarities. The kappa-form electrons generate less nonthermal power /u1D443 /u1D705 ∕/u1D443 power-law ≈ 0.81 despite larger total injection rates: 20.5 vs. 13.6 [ 10 35 electrons sec -1 ]. Figure 7 also illustrates a comparison between the distribution of accelerated kappa-form electrons (blue curve, total electron number density /u1D45B /u1D458 = 1.1 × 10 9 cm -3 ) and the ambient thermal Maxwellian distribution (red curve, total electron number density /u1D45B /u1D461 ℎ = 2.4 × 10 10 cm -3 , /u1D458 B /u1D447 = 1.67 keV) in the right panel.

## 4. ANALYSIS OF THE FIT PARAMETERS AND FIT RESULTS

In this section, we will further investigate how well the three free kappa parameters ( ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 ) can be constrained by the fit with the warm-target model. We will provide the uncertainty of the fit parameters using different methods. Additionally, we will compare the total electron density /u1D45B /u1D458 and average electron energy /u1D438 avg derived from the fit kappa-form spectrum with the corresponding thermal properties.

## 4.1. 2011 February 24 RHESSI flare

In the previous section, we show that fitting using the functions f\_vth and f\_thick\_warm\_kappa involves three free parameters: ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 , all associated with the kappa distribution. To estimate the uncertainty of each parameter at different confidence levels, we apply the method used in Kontar et al. (2019). We began by assessing the uncertainty of total injection rate ̇ /u1D441 0 . We carry out spectral fits with ̇ /u1D441 0 set at a fixed value while setting the two other kappa parameters ( /u1D705 and /u1D447 /u1D705 ) free. By doing this, we can obtain the best-fit full /u1D712 2 as well as the corresponding kappa index /u1D705 and kappa temperature /u1D447 /u1D705 at this fixed ̇ /u1D441 0 . Subsequently, we repeat these fits multiple times with different fixed ̇ /u1D441 0 values and record the corresponding full /u1D712 2 values. Further, we define ∆/u1D712 2 = /u1D712 2 -/u1D45A/u1D456/u1D45B(/u1D712 2 ) , where /u1D45A/u1D456 /u1D45B(/u1D712 2 ) is the

Figure 8. (a): ∆/u1D712 2 against the given total injected electron rate ( ̇ /u1D441 0 ) fixed for the fit. The red and green horizontal lines represent the 1 /u1D70E and 3 /u1D70E uncertainty levels, respectively. Red and green vertical lines correspond to the ̇ /u1D441 0 ranges within the 1 /u1D70E and 3 /u1D70E uncertainty levels, 22.5 +0.8 -0.6 ×10 35 electrons sec -1 , and 22.5 +2.4 -1.8 ×10 35 electrons sec -1 , respectively. (b): Corresponding average electron energy /u1D438 avg (black curve) and number density /u1D45B /u1D458 (blue curve) derived from the fit kappa-form electron spectrum with different fixed ̇ /u1D441 0 in (a). The red and green vertical lines represent the ̇ /u1D441 0 range within the 1 /u1D70E and 3 /u1D70E uncertainty levels, same in panel (a). The black horizontal line represents the average electron energy from the thermal loop ( 3 2 /u1D458 B /u1D447 loop ) at 1.95 keV. The blue horizontal line represents the total number density of 4.4 × 10 10 cm -3 , which is the number density of ambient coronal loop plasma and serves as an upper limit of the fit number density. We note that the fit results within the acceptable range satisfy the two constraints from average electron energy ( /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 &gt; /u1D438 /u1D44E/u1D463 /u1D454-/u1D459 /u1D45C/u1D45C/u1D45D ) and number density ( /u1D45B /u1D458 &lt; /u1D45B loop ). (c): Same as panel (a) but the fixed parameter is kappa temperature /u1D458 B /u1D447 /u1D705 . The /u1D458 B /u1D447 /u1D705 ranges within the 1 /u1D70E and 3 /u1D70E uncertainty levels are 1.29 +0.04 -0.04 keV and 1.29 +0.13 -0.12 keV, respectively. (d): Same as panel (a) but the fixed parameter is kappa index /u1D705 . The /u1D705 ranges within the 1 /u1D70E and 3 /u1D70E uncertainty levels are 5.14 +0.04 -0.04 and 5.14 +0.12 -0.12 , respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-473be25a-efe0-45fa-b1d4-4f8a627c9219.jpg)

minimum of all the obtained full /u1D712 2 values. ∆/u1D712 2 reaches 0 at the ̇ /u1D441 0 value of 22.5 ×10 35 electrons sec -1 (Figure 8 (a)), consistent with the best-fit results shown in Table 1. The ∆/u1D712 2 -̇ /u1D441 0 curve (Figure 8 (a)) displays a prominent minimum, and the ∆/u1D712 2 value changes rapidly in the vicinity of this minimum. The ̇ /u1D441 0 within the 1 /u1D70E uncertainty level (corresponding to a 68% confidence level, red horizontal dashed line in Figure 8 (a)) is 22.5 +0.8 -0.6 ×10 35 electrons sec -1 (red vertical dashed lines in Figure 8 (a)), and ̇ /u1D441 0

Figure 9. Corner plot of the posterior probability function for the 2011 February 24 M3.5 flare obtained from the Monte Carlo analysis performed using OSPEX . The diagonal panels display the one-dimensional projected probability function for each parameter ( ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 ), and the remaining panels show the two-dimensional projection. The final fit results of the Monte Carlo analysis and the best-fit results of the forward fit are depicted in red and blue, respectively. Additionally, the plot includes the 1 /u1D70E level uncertainty range for each parameter: ̇ /u1D441 0 = 22.5 +0.6 -0.6 ×10 35 electrons sec -1 , /u1D705 = 5.14 +0.03 -0.03 , and /u1D458 B /u1D447 /u1D705 = 1.29 +0.04 -0.03 keV.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-a953950d-f19c-4eff-a8f9-8e545d350791.jpg)

within the 3 /u1D70E uncertainty level (corresponding to a 99% confidence level, green horizontal dashed line in Figure 8 (a)) is 22.5 +2.4 -1.8 ×10 35 electrons sec -1 (green vertical dashed lines in Figure 8 (a)). The ∆/u1D712 2 -̇ /u1D441 0 curve and the estimated uncertainty suggest that ̇ /u1D441 0 can be effectively constrained by the fit with the warm-target model.

The same approach is used to analyze the two other kappa parameters, /u1D705 and /u1D447 /u1D705 , for uncertainty. Both the ∆/u1D712 2 -/u1D447 /u1D705 (Figure 8 (c)) and ∆/u1D712 2 -/u1D705 (Figure 8 (d)) curves reveal a distinct minimum, similar to the ∆/u1D712 2 -̇ /u1D441 0 curve. The ranges of /u1D458 B /u1D447 /u1D705 for 1 /u1D70E and 3 /u1D70E uncertainty levels (red and green horizontal dashed lines in Figure 8 (c)) are 1.29 +0.04 -0.04 keV and 1.29 +0.13 -0.12 keV (red and green vertical dashed lines in Figure 8 (c)), respectively. The ranges of /u1D705 for 1 /u1D70E and 3 /u1D70E uncertainty levels (red and green horizontal

Figure 10. The upper panel depicts the photon spectrum of the observation (background subtracted, black curve), the best-fit results obtained from the forward fit (blue curve), the final fit result generated by the Monte Carlo analysis (red curve), and the photon spectrum generated by the Monte Carlo sampling parameters (gray curves) of the 2011 February 24 M3.5 flare. The lower panel displays the normalized residuals associated with the corresponding fit spectrum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-3d2dc54b-8a32-4296-9a17-02bf1f3918db.jpg)

dashed lines in Figure 8 (d)) are 5.14 +0.04 -0.04 and 5.14 +0.12 -0.12 (red and green vertical dashed lines in Figure 8 (d)), respectively. The small fit uncertainties for all three kappa parameters indicate that all parameters can be accurately determined using the warm-target model.

$$E _ { \text{avg} } = \frac { \int _ { 0 } ^ { \infty } E f _ { k } ( E ) \, d E } { \int _ { 0 } ^ { \infty } f _ { k } ( E ) \, d E } = \frac { 3 } { 2 } \frac { \kappa } { \kappa - 2. 5 } k _ { \text{B} } T _ { \kappa }.$$

In this study, we utilize the 3 /u1D70E uncertainty level to define the acceptable range of the fit parameters and investigate the physical properties within this acceptable range. One valuable property we examine is the average energy of the electrons /u1D438 avg . Here, we first adopt the acceptable range of ̇ /u1D441 0 , which falls between 20.3 and 24.9 ×10 35 electrons sec -1 . Each given ̇ /u1D441 0 within the acceptable range corresponds to a fit kappa-form electron spectrum. The average electron energy for the kappa-form electron spectrum /u1D438 avg is given by

For comparison, the average electron energy of the Maxwellian electron spectrum is /u1D438 avg = 3 2 /u1D458 B /u1D447 . The obtained /u1D438 avg -̇ /u1D441 0 curve is depicted by the black curve in Figure 8 (b). It is worth noting that, within the acceptable range (green vertical dashed lines in Figure 8 (b)), the average electron energies derived from the kappa-form electron spectrum /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 are consistently higher than the average electron energy of the thermal plasma in the target loop ( /u1D438 /u1D44E/u1D463 /u1D454-/u1D459 /u1D45C/u1D45C/u1D45D = 3 2 /u1D458 B /u1D447 loop = 1.95 keV, black horizontal dashed lines in Figure 8 (b)). According to Equation 8, the total electron number density /u1D45B /u1D458 can be determined once the injection area /u1D434 is known. We use the cross-area of the loop-top X-ray coronal source ( /u1D434 = /u1D70B/u1D45F 2 = 4.21 × 10 10 cm 2 , red outlier contour in Figure 2 right panel) as the injection area /u1D434 to calculate the total electron number density /u1D45B /u1D458 . The /u1D45B /u1D458 -̇ /u1D441 0 curve is shown as the blue curve in Figure 8 (b).

Within the acceptable range of ̇ /u1D441 0 , the total electron number density /u1D45B /u1D458 is found to be lower than the thermal electron density of the target coronal loop /u1D45B loop (blue horizontal dashed lines in Figure 8 (b), 4.4 × 10 10 cm -3 ). While the electron injection area cannot be precisely determined, and it may be smaller than the observed loop cross-area. /u1D45B /u1D458 is unlikely to be higher than the loop density /u1D45B loop , which requires an extremely small injection area /u1D434 &lt; 1.60 × 10 9 cm 2 . The results indicate that the values of /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 and /u1D45B /u1D458 are physically reasonable when compared with the ambient thermal plasma. The fitted kappa-form electron spectrum and associated properties ( /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 and /u1D45B /u1D458 ) within the acceptable range for /u1D447 /u1D705 and /u1D705 closely resemble those obtained within the acceptable range for ̇ /u1D441 0 , and are consistent with the previous findings. Consequently, we have opted not to display them in the plot.

In order to further evaluate how well the fit parameters can be determined and cross-check their uncertainty, we utilized the Monte Carlo analysis provided by OSPEX . Detailed information on this method is available in the 'Fit Parameter Uncertainty Analysis' section of Tolbert &amp; Schwartz (2020). In addition, Ireland et al. (2013) discussed commonly used methods, including the Monte Carlo method, for estimating the constraint of fit parameters. The Monte Carlo analysis requires a predetermined best-fit result, which is available in Table 1. The Monte Carlo sampling explores the multidimensional parameter space ( ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 ) surrounding the best-fit results by assuming a Gaussian probability function for each parameter. Figure 9 shows the posterior probability density function for the three free parameters obtained from the Monte Carlo analysis. The final fit results of the Monte Carlo analysis, as well as the predetermined best-fit results from the forward fit, are also depicted by red and blue dots, respectively. The diagonal panels of the corner plot show the one-dimensional projection of the probability density function. For the 2011 February 24 M3.5 flare, each parameter exhibits a prominent peak and a relatively narrow width, indicating that they are well constrained around the best-fit results. The 1 /u1D70E uncertainty range for each parameter is shown in Figure 9: ̇ /u1D441 0 = 22.5 +0.6 -0.6 ×10 35 electrons sec -1 , /u1D705 = 5.14 +0.03 -0.03 , and /u1D458 B /u1D447 /u1D705 = 1.29 +0.04 -0.03 keV, comparable with the uncertainty obtained from the previous method. The nondiagonal panels display the probability density function projected for different parameter pairs, and illustrate the correlations between different parameters. We observed that /u1D705 and /u1D447 /u1D705 are positively correlated, while ̇ /u1D441 0 is anticorrelated with /u1D705 and /u1D447 /u1D705 . The photon spectrum and normalized residuals generated by Monte Carlo sampling, the Monte Carlo final fit result, and the best forward fit result are shown as gray, red, and blue curves in Figure 10.

## 4.2. 2022 March 28 STIX flare

We conducted the same analysis on the 2022 March 28 M4.0 flare observed by STIX. Similar to the 2011 February 24 M3.5 flare, the ∆/u1D712 2 -̇ /u1D441 0 , ∆/u1D712 2 -/u1D447 /u1D705 , and ∆/u1D712 2 -/u1D705 curves (black curves in Figures 11 (a), (c), and(d), respectively) all exhibit a prominent minimum. The ranges for each parameter at 1 /u1D70E and3 /u1D70E uncertainty levels are as follows: for ̇ /u1D441 0 , 20.5 +1.3 -1.2 ×10 35 electrons sec -1 , and 20.5 +3.9 -3.8 ×10 35 electrons sec -1 ; for /u1D458 B /u1D447 /u1D705 , 2.97 +0.17 -0.16 keV, and 2.97 +0.55 -0.45 keV; for /u1D705 , 8.85 +0.41 -0.38 , and 8.85 +1.45 -1.01 . The spectral fit for the 2022 March 28 M4.0 flare, due to a narrow fit range (6-50 keV) and a 6% systematic error from STIX (2% for RHESSI), results in relatively larger uncertainties, particularly for /u1D447 /u1D705 and /u1D705 . However, these uncertainties still demonstrate that the warm-target fit effectively constrains the fit parameters. We adopt the 3 /u1D70E uncertainty level as the acceptable range, consistent with our approach for the 2011 February 24 M3.5 flare. The derived average electron energy /u1D438 avg and total electron density /u1D45B /u1D458 both satisfy /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 &gt; /u1D438 /u1D44E/u1D463 /u1D454-/u1D459 /u1D45C/u1D45C/u1D45D and /u1D45B /u1D458 &lt; /u1D45B loop within the acceptable range (Figure 11 (b)). Additionally, the results of Monte Carlo analysis are displayed in Figures 12 and 13. The 1 /u1D70E uncertainty level ranges for each parameter are: ̇ /u1D441 0 = 20.5 +0.2 -0.2 ×10 35 electrons sec -1 , /u1D705 = 8.85 +0.17 -0.14 , and /u1D458 B /u1D447 /u1D705 = 2.97 +0.06 -0.05 keV. The correlations

Figure 11. Same as Figure 8 but for the 2022 March 28 M4.0 flare. (a): ∆/u1D712 2 against the different total injected electron rate ( ̇ /u1D441 0 ). The red and green horizontal lines represent the 1 /u1D70E and 3 /u1D70E uncertainty levels, respectively. Red and green vertical lines correspond to the ̇ /u1D441 0 range within the 1 /u1D70E and 3 /u1D70E uncertainty levels, 20.5 +1.3 -1.2 ×10 35 electrons sec -1 and 20.5 +3.9 -3.8 ×10 35 electrons sec -1 , respectively. (b): Corresponding average electron energy /u1D438 avg (black curve) and number density /u1D45B /u1D458 (blue curve) derived from the fit kappa-form electron spectrum for different fixed ̇ /u1D441 0 . The red and green vertical lines represent the ̇ /u1D441 0 range within the 1 /u1D70E and 3 /u1D70E uncertainty levels, the same as in panel (a). The black horizontal line represents the average electron energy from the thermal loop ( 3 2 /u1D458 B /u1D447 loop ) at 2.51 keV. The blue horizontal line represents the electron number density of the ambient thermal plasma, 2.4 × 10 10 cm -3 . (c): Same as panel (a) but the fixed parameter is kappa temperature. The /u1D458 B /u1D447 /u1D705 ranges within the 1 /u1D70E and 3 /u1D70E uncertainty levels are 2.97 +0.17 -0.16 keV and 2.97 +0.55 -0.45 keV, respectively. (d): Same as panel a but the fixed parameter is kappa index. The /u1D705 ranges within the 1 /u1D70E and 3 /u1D70E uncertainty levels are 8.85 +0.41 -0.38 and 8.85 +1.45 -1.01 , respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-25dc1aa8-d476-4dc8-9ec0-4dd00ba24eed.jpg)

between these parameters are consistent with the findings from the 2011 February 24 M3.5 flare (Figure 12 diagonal panels).

## 5. SUMMARY

In this study, we utilized the warm-target model to analyze two GOES M-class limb flares. The M3.5 class flare occurred on 2011 February 24 and was observed by RHESSI, while STIX observed the other M4.0 class flare on 2022 March 28. We used the warm-target model due to its ability to accurately de-

Figure 12. Corner plot of the posterior probability function for the 2022 March 28 M4.0 flare obtained from the Monte Carlo analysis performed using OSPEX . As in Figure 9, the diagonal panels display the one-dimensional projected probability function for each parameter ( ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 ), and the nondiagonal panels show the twodimensional projection. The final fit results of the Monte Carlo analysis and the best-fit results of the forward fit are depicted in red and blue, respectively. The plot also includes the 1 /u1D70E level uncertainty range for each parameter: ̇ /u1D441 0 = 20.5 +0.2 -0.2 ×10 35 electrons sec -1 , /u1D705 = 8.85 +0.17 -0.14 , and /u1D458 B /u1D447 /u1D705 = 2.97 +0.06 -0.05 keV.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-631a8673-e2c8-4010-beb3-b6b0ec8c4190.jpg)

termine nonthermal electron properties, such as nonthermal power, which can hardly be constrained by the commonly used cold-target model. Unlike previous investigations, the kappa distribution is used to characterize the accelerated/injected electrons. Unlike the power-law distribution, the kappa distribution is finite as the electron speed approaches zero, thereby effectively covering the entire electron energy range without the requirement of an arbitrary low-energy cutoff. Consequently, the fitted kappa-form electron spectrum can characterize electrons whose kinetic energy falls below the X-ray instrument's sensitive range. The results of the fitting and subsequent analysis are summarized below.

- · The best-fit results using the warm-target model in kappa-form injected electrons feature reasonable electron spectrum. For comparison, we also employed the warm-target model in the

Figure 13. SameasFigure 10, but for the 2022 March 28 M4.0 flare. The upper panel shows the photon spectrum of the observation (background subtracted, black curve), the best-fit results obtained from the forward fit (blue curve), the final fit result generated by the Monte Carlo analysis (red curve), and the photon spectrum generated by the Monte Carlo sampling parameters (gray curves). The lower panel displays the normalized residuals associated with the corresponding fit spectrum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00213v2-figure-83ae3b47-abf7-4bab-88e6-59e8b0c0c620.jpg)

power-law electron distribution to fit the HXR spectrum. We found that the obtained power-lawform and kappa-form electron spectra generate similar photon spectra within the fitting range. However, the kappa-form electron spectrum produces less nonthermal power despite having a comparable or higher total injection rate ̇ /u1D441 0 . This energy difference is due to the behavior of the electron spectrum. The obtained kappa-form and power-law-form electron spectra exhibited similar behavior at the deka-keV level, while the kappa-form spectrum displayed a significantly lower flux below ∼ 30 keV. Moreover, the kappa-form electron spectrum decreased more rapidly than the power-law spectrum ( /u1D705 &gt; /u1D6FF +1 ) at higher energies (typically above 100 keV, out of the fit range).

- · The fit with the warm-target model in kappa-form electron injection involves three free parameters ( ̇ /u1D441 0 , /u1D705 , and /u1D447 /u1D705 ). The ∆/u1D712 2 -̇ /u1D441 0 curve reveals a prominent minimum, suggesting an accurate determination of ̇ /u1D441 0 within a narrow range. For the M3.5 flare on 2011 February 24, the range of ̇ /u1D441 0 within the 1 /u1D70E uncertainty level is 22.5 +0.8 -0.6 ×10 35 electrons sec -1 , and within the 3 /u1D70E uncertainty level range it is 22.5 +2.4 -1.8 ×10 35 electrons sec -1 . For the M4.0 flare on 2022 March 28, the range of ̇ /u1D441 0 within the 1 /u1D70E uncertainty level is 20.5 +1.3 -1.2 ×10 35 electrons sec -1 , and within the 3 /u1D70E uncertainty level it is 20.5 +3.9 -3.8 ×10 35 electrons sec -1 . We apply the same analyses to the other two fit parameters, /u1D705 and /u1D447 /u1D705 , for their uncertainty. For the M3.5 flare on 2011 February 24, /u1D705 within 1 /u1D70E and 3 /u1D70E uncertainty levels is 5.14 +0.04 -0.04 and 5.14 +0.12 -0.12 , respectively. /u1D458 B /u1D447 /u1D705 within 1 /u1D70E and 3 /u1D70E uncertainty levels is 1.29 +0.04 -0.04 keV and 1.29 +0.13 -0.12 keV, respectively. For the M4.0 flare on March 28, 2022, /u1D705 within 1/u1D70E and 3/u1D70E uncertainty level range is 8.85 +0.41 -0.38 and 8.85 +1.45 -1.01 , respectively. /u1D458 B /u1D447 /u1D705 within 1/u1D70E and 3/u1D70E uncertainty level range is 2.97 +0.17 -0.16 keV and 2.97 +0.55 -0.45 keV, respectively. All the three

parameters, ̇ /u1D441 0 , /u1D447 /u1D705 and /u1D705 are tightly constrained within a narrow range. Additionally, we carry out a Monte Carlo analysis to further confirm the effectiveness in determining the fit parameters with the warm-target model. The resulting posterior probability density function for all three kappa parameters demonstrates distinct peaks and relatively narrow widths, indicating wellconstrained fit parameters. The Monte Carlo analysis also provides a 1 /u1D70E level uncertainty for each parameter: for the 2011 February 24 M3.5 flare, ̇ /u1D441 0 = 22.5 +0.6 -0.6 ×10 35 electrons sec -1 , /u1D705 = 5.14 +0.03 -0.03 , and /u1D458 B /u1D447 /u1D705 = 1.29 +0.04 -0.03 keV; for the 2022 March 28 M4.0 flare, ̇ /u1D441 0 = 20.5 +0.2 -0.2 ×10 35 electrons sec -1 , /u1D705 = 8.85 +0.17 -0.14 , and /u1D458 B /u1D447 /u1D705 = 2.97 +0.06 -0.05 keV.

- · Unlike the power-law distribution, which requires a low-energy cutoff /u1D438 /u1D450 , the kappa-form electron spectrum works for all energies, anchored at zero speed and extending to the speed of light. Therefore, the kappa distribution facilitates the determination of the accelerated electron number density /u1D45B /u1D458 (with the information of the injection area /u1D434 ) and the average electron energy /u1D438 avg . For both flares in this study, we find that the derived total electron energy density is lower than the thermal number density of the target loop ( /u1D45B /u1D458 &lt; /u1D45B loop , 1.7 × 10 9 cm -3 &lt; 4.4 × 10 10 cm -3 and 1.1 × 10 9 cm -3 &lt; 2.4 × 10 10 cm -3 for 2011 February 24 and 2022 March 28 flares, respectively), while the average electron energy surpasses that of the ambient Maxwellian thermal plasma ( /u1D438 /u1D44E/u1D463 /u1D454-/u1D705 &gt; /u1D438 /u1D44E/u1D463 /u1D454-/u1D459 /u1D45C/u1D45C/u1D45D , 3.77 keV &gt; 1.95 keV and 6.21 keV &gt; 2.51 keV for 2011 February 24 and 2022 March 28 flares, respectively).

In this study, we have demonstrated that the fit with the warm-target model provides a convincing kappa-form accelerated electron spectrum. The associated kappa distribution parameters are valuable for understanding electron acceleration and transport. Here, we follow the model proposed by Bian et al. (2014), which uses the Fokker-Planck equation to analyze the evolution of flare-associated electrons. Bian et al. 2014 proposed that the kappa distribution arises from stochastic acceleration in the presence of Coulomb collisions and velocity diffusion. The kappa index /u1D705 is given by /u1D705 = /u1D70F /u1D44E/u1D450 /u1D450 (/u1D463) 2/u1D70F /u1D450 (/u1D463) . The collisional deceleration time can be estimated as: /u1D70F /u1D450 (/u1D463) ≈ /u1D463 3 ∕Γ , where the collisional parameter Γ = 2/u1D43E/u1D45B∕/u1D45A 2 /u1D452 . Here, we take the parameters from the 2011 February 24 M3.5 flare event for study. For the upper limit of electron energy /u1D438 /u1D459 enabling the complete stop within the target, we have /u1D438 /u1D459 = √ 2/u1D43E/u1D45B/u1D43F ≈ 19.0 keV . At /u1D438 /u1D459 , we have the electron speed /u1D463 /u1D459 = 8.0 × 10 9 cm∕s and the collisional deceleration time /u1D70F /u1D450 (/u1D463 /u1D459 ) ≈ 0.71 seconds . Thus the acceleration time can be determined as /u1D70F acc (/u1D463 /u1D459 ) = 2/u1D705/u1D70F /u1D450 (/u1D463 /u1D459 ) ≈ 7.33 seconds . The collisional diffusion time /u1D70F /u1D451 (/u1D463 /u1D459 ) = 2/u1D463 5 /u1D459 Γ/u1D463 2 /u1D461 /u1D452 ≈ 19.77 seconds , where √

The ratio between the density of accelerated nonthermal electrons ( /u1D45B nth ) and the density of ambient thermal electrons ( /u1D45B th ) can provide valuable insights into electron acceleration during solar flares. Previous studies have shown that the ratio of nonthermal electrons to protons ( /u1D45B p = /u1D45B nth + /u1D45B th in the hydrogen plasma) is of the order of approximately 1% (Oka et al. 2013). In this study, we determined the density of accelerated nonthermal electrons in the corona by integrating the kappa electron spectrum over 10 keV( /u1D45B nth [&gt; 10 keV ] ∼ 1.2 × 10 8 cm -3 , and 2.0 × 10 8 cm -3 for the 2011 February 24 M3.5 and 2022 March 28 M4.0 flares, respectively). For the thermal electron density at the flare site, we take the loop electron density obtained from the X-ray thermal component ( /u1D45B th = /u1D45B loop ). The ratio between the density of coronal accelerated nonthermal electrons and ambient protons is found to be less than thermal speed /u1D463 te = 2/u1D458 B /u1D447 /u1D705 ∕/u1D45A /u1D452 ≈ 2.1 × 10 9 cm∕s .

1% ( /u1D45B nth ∕/u1D45B p ∼ 0.003 and 0.008, for the 2011 February 24 M3.5 and 2022 March 28 flares, respectively), consistent with recent studies (Kontar et al. 2023).

Since both RHESSI and Sol-O/STIX are not sensitive to photons below ∼ 3-6 keV, the thermal properties obtained from RHESSI X-ray spectral diagnostics could differ from those obtained from SDO/AIA EUV diagnostics (Battaglia &amp; Kontar 2013). Furthermore, thermal properties inferred from the observed X-ray spectrum are likely to be height-dependent (Jeffrey et al. 2015). The obtained kappa-form energetic electrons allows for a more comprehensive combined study of the X-ray and EUV diagnostics. Potential future applications of the derived kappa distribution also include studying acceleration mechanisms based on the average electron energy and nonthermal power.

Here, we note that the electron number density is calculated using Equation 8 and is dependent on the injection area /u1D434 . We estimated the injection area /u1D434 using 50% of the loop-top X-ray source contour. The resulting total injected electron density is approximately an order of magnitude lower than the estimated coronal loop thermal electron density ( 1.7 × 10 9 cm -3 &lt; 4.4 × 10 10 cm -3 and 1.1 × 10 9 cm -3 &lt; 2.4 × 10 10 cm -3 for the 2011 February 24 and 2022 March 28 flares, respectively). It is important to note that the electron injection in the coronal region may be nonuniform in height. According to Jeffrey et al. (2015), the actual injection might be smaller. In this study, we are utilizing the spatially integrated spectrum for analysis. At present, the signal-to-noise ratio does not sufficiently support precise spatially resolved spectral diagnostics at the arcsecond scale. Once we obtain the spatially resolved spectral information from the improved X-ray observations, the electron density derived from the kappa distribution will be more persuasive.

The authors thank the referee for improving the manuscript. The work was supported via the STFC/UKRI grant ST/T000422/1. 1 2

## REFERENCES

Alaoui, M., &amp; Holman, G. D. 2017, ApJ, 851, 78, doi: 10.3847/1538-4357/aa98de

Arnold, H., Drake, J. F., Swisdak, M., et al. 2021, PhRvL, 126, 135101,

doi: 10.1103/PhysRevLett.126.135101

Aschwanden, M. J., Holman, G., O'Flannagain, A., et al. 2016, ApJ, 832, 27,

doi: 10.3847/0004-637X/832/1/27

Aschwanden, M. J., Kontar, E. P., &amp; Jeffrey, N. L. S.

2019, ApJ, 881, 1, doi: 10.3847/1538-4357/ab2cd4

Aschwanden, M. J., Caspi, A., Cohen, C. M. S., et al. 2017, ApJ, 836, 17,

doi: 10.3847/1538-4357/836/1/17

Battaglia, M., &amp; Benz, A. O. 2006, A&amp;A, 456, 751,

doi: 10.1051/0004-6361:20065233

Battaglia, M., &amp; Kontar, E. P. 2011, A&amp;A, 533, L2,

doi: 10.1051/0004-6361/201117605

-. 2012, ApJ, 760, 142,

doi: 10.1088/0004-637X/760/2/142

-. 2013, ApJ, 779, 107,

doi: 10.1088/0004-637X/779/2/107

Battaglia, M., Motorina, G., &amp; Kontar, E. P. 2015,

ApJ, 815, 73, doi: 10.1088/0004-637X/815/1/73

Benz, A. O. 2017, Living Reviews in Solar Physics, 14, 2, doi: 10.1007/s41116-016-0004-3

Bian, N. H., Emslie, A. G., Stackhouse, D. J., &amp;

Kontar, E. P. 2014, ApJ, 796, 142,

doi: 10.1088/0004-637X/796/2/142

Brown, J. C. 1971, SoPh, 18, 489,

doi: 10.1007/BF00149070

Brown, J. C., Emslie, A. G., Holman, G. D., et al.

2006, ApJ, 643, 523, doi: 10.1086/501497

Brown, J. C., Emslie, A. G., &amp; Kontar, E. P. 2003,

ApJL, 595, L115, doi: 10.1086/378169

Cheung, M. C. M., Boerner, P., Schrijver, C. J., et al. 2015, ApJ, 807, 143,

doi: 10.1088/0004-637X/807/2/143

| Collier, M. R., Hamilton, D. C., Gloeckler, G., Bochsler, P., &Sheldon, R. B. 1996, Geophys. Res. Lett., 23, 1191, doi: 10.1029/96GL00621   |
|---------------------------------------------------------------------------------------------------------------------------------------------|
| Craig, I. J. D., &Brown, J. C. 1976, A&A, 49, 239                                                                                           |
| Drummond, W. E., &Pines, D. 1964, Annals of Physics, 28, 478, doi: 10.1016/0003-4916(64)90205-2                                             |
| Effenberger, F., Rubio da Costa, F., Oka, M., et al. 2017, ApJ, 835, 124, doi: 10.3847/1538-4357/835/2/124                                  |
| Emslie, A. G. 1980, ApJ, 235, 1055, doi: 10.1086/157709                                                                                     |
| Emslie, A. G., Kontar, E. P., Krucker, S., &Lin, R. P. 2003, ApJL, 595, L107, doi: 10.1086/378931                                           |
| Emslie, A. G., &Smith, D. F. 1984, ApJ, 279, 882, doi: 10.1086/161959                                                                       |
| Galloway, R. K., MacKinnon, A. L., Kontar, E. P., & Helander, P. 2005, A&A, 438, 1107, doi: 10.1051/0004-6361:20042137                      |
| Hannah, I. G., &Kontar, E. P. 2012, A&A, 539, A146, doi: 10.1051/0004-6361/201117576                                                        |
| Hannah, I. G., Kontar, E. P., &Sirenko, O. K. 2009, ApJL, 707, L45, doi: 10.1088/0004-637X/707/1/L45                                        |
| Haug, E. 1998, SoPh, 178, 341, doi: 10.1023/A:1005098624121                                                                                 |
| Högbom, J. A. 1974, A&AS, 15, 417 Holman, G. D. 2012, ApJ, 745, 52,                                                                         |
| G. D., &Benka, S. G. 1992, ApJL, 400, L79,                                                                                                  |
| Aurass, H., et al. 2011, SSRv, 159, 107,                                                                                                    |
| 2011, Journal of Geophysical Research (Space Physics), 116, A08217, doi: 10.1029/2011JA016576                                               |
| Ireland, J., Tolbert, A. K., Schwartz, R. A., G. D., &Dennis, B. R. 2013, ApJ, 769, 89,                                                     |
| Holman, doi: 10.1088/0004-637X/769/2/89                                                                                                     |
| Jeffrey, N. L. S., Kontar, E. P., Bian, N. H., &Emslie, A. G. 2014, ApJ, 787, 86,                                                           |
| doi: 10.1088/0004-637X/787/1/86 2015,                                                                                                       |
| A&A, 584, A89,                                                                                                                              |
| doi: 10.1051/0004-6361/201526665                                                                                                            |
| R.                                                                                                                                          |
| Jeffrey, N. L. S., Kontar, E. P., &Dennis, B.                                                                                               |
| A. R., &Christe, S. 2014, ApJ, doi: 10.1088/0004-637X/789/2/116                                                                             |
| doi: 10.1007/s11214-010-9680-9 Imada, S., Hirai, M., Hoshino, M., Inglis,                                                                   |
| Holman, doi: 10.1086/186654                                                                                                                 |
| doi:                                                                                                                                        |
| 10.1088/0004-637X/745/1/52                                                                                                                  |
| Holman, G. D., Aschwanden, M. J.,                                                                                                           |
| &Mukai,                                                                                                                                     |
| T.                                                                                                                                          |
| 789, 116,                                                                                                                                   |

| Jeffrey, N. L. S., Krucker, S., Stores, M., et al. 2024, ApJ, 964, 145, doi: 10.3847/1538-4357/ad236f                                                                                                                                                                           |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Kašparová, J., &Karlický, M. 2009, A&A, 497, L13, doi: 10.1051/0004-6361/200911898                                                                                                                                                                                              |
| Knight, J. W., &Sturrock, P. A. 1977, ApJ, 218, 306, doi: 10.1086/155683                                                                                                                                                                                                        |
| Kontar, E. P., Emslie, A. G., Motorina, G. G., & Dennis, B. R. 2023, ApJL, 947, L13, doi: 10.3847/2041-8213/acc9b7                                                                                                                                                              |
| Kontar, E. P., Hannah, I. G., &MacKinnon, A. L. 2008, A&A, 489, L57, doi: 10.1051/0004-6361:200810719                                                                                                                                                                           |
| Kontar, E. P., Jeffrey, N. L. S., &Emslie, A. G. 2019, ApJ, 871, 225, doi: 10.3847/1538-4357/aafad3                                                                                                                                                                             |
| Kontar, E. P., Jeffrey, N. L. S., Emslie, A. G., &Bian, N. H. 2015, ApJ, 809, 35, doi: 10.1088/0004-637X/809/1/35                                                                                                                                                               |
| Kontar, E. P., Brown, J. C., Emslie, A. G., et al. 2011, SSRv, 159, 301, doi: 10.1007/s11214-011-9804-x Krall, N. A., &Trivelpiece, A. W. 1973, Principles of plasma physics (Tokyo: McGraw-Hill Kogakusha) Krucker, S., Hurford, G. J., Grimm, O., et al. 2020, A&A, 642, A15, |
| doi: 10.1051/0004-6361/201937362 Landi, E., Young, P. R., Dere, K. P., Del Zanna, G., & Mason, H. E. 2013, ApJ, 763, 86, doi: 10.1088/0004-637X/763/2/86                                                                                                                        |
| Lin, R. P., &Hudson, H. S. 1976, SoPh, 50, 153, doi: 10.1007/BF00206199                                                                                                                                                                                                         |
| Lin, R. P., Dennis, B. R., Hurford, G. J., et al. 2002, SoPh, 210, 3, doi: 10.1023/A:1022428818870                                                                                                                                                                              |
| Ma, C.-y., &Summers, D. 1998, Geophys. Res. Lett., 25, 4099, doi: 10.1029/1998GL900108                                                                                                                                                                                          |
| Mace, R. L., &Hellberg, M. A. 1995, Physics of Plasmas, 2, 2098, doi: 10.1063/1.871296                                                                                                                                                                                          |
| Maksimovic, M., Zouganelis, I., Chaufray, J. Y., et al. 2005, Journal of Geophysical Research (Space Physics), 110, A09104, doi: 10.1029/2005JA011119                                                                                                                           |
| Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1, doi: 10.1051/0004-6361/202038467                                                                                                                                                                         |
| Oka, M., Ishikawa, S., Saint-Hilaire, P., Krucker, S.,                                                                                                                                                                                                                          |
| &Lin, R. P. 2013, ApJ, 764, 6, doi: 10.1088/0004-637X/764/1/6                                                                                                                                                                                                                   |
| Oka, M., Krucker, S., Hudson, H. S., &Saint-Hilaire, P. 2015, ApJ, 799, 129, doi: 10.1088/0004-637X/799/2/129                                                                                                                                                                   |
| Purkhart, S., Veronig, A. M., Dickson, E. C. M., et al. 2023, A&A, 679, A99,                                                                                                                                                                                                    |

doi: 10.1051/0004-6361/202346354

| Rhee, T., Ryu, C.-M., &Yoon, P. H. 2006, Journal of Geophysical Research (Space Physics), 111, A09107, doi: 10.1029/2006JA011682   |
|------------------------------------------------------------------------------------------------------------------------------------|
| Rochus, P., Auchère, F., Berghmans, D., et al. 2020, A&A, 642, A8, doi: 10.1051/0004-6361/201936663                                |
| Ryu, C. M., Rhee, T., Umeda, T., Yoon, P. H., & Omura, Y. 2007, Physics of Plasmas, 14, 100701, doi: 10.1063/1.2779282             |
| Schwartz, R. A., Csillaghy, A., Tolbert, A. K., et al. 2002, SoPh, 210, 165, doi: 10.1023/A:1022444531435                          |
| Shao, C., &Huang, G. 2009, ApJ, 691, 299, doi: 10.1088/0004-637X/691/1/299                                                         |
| Spitzer, L. 1962, Physics of Fully Ionized Gases                                                                                   |
| Stackhouse, D. J., &Kontar, E. P. 2018, A&A, 612, A64, doi: 10.1051/0004-6361/201730708                                            |
| Sui, L., Holman, G. D., Dennis, B. R., et al. 2002, SoPh, 210, 245, doi: 10.1023/A:1022417832283                                   |

| 16, 273                                                                                                                                              |
|------------------------------------------------------------------------------------------------------------------------------------------------------|
| Tolbert, K., &Schwartz, R. 2020, OSPEX: Object Spectral Executive, Astrophysics Source Code Library, record ascl:2007.018. https://ascl.net/2007.018 |
| Vandas, M. 1989, Bulletin of the Astronomical Institutes of Czechoslovakia, 40, 175 Vedenov,                                                         |
| A. A., &Velikhov, E. P. 1963, Soviet Journal of Experimental and Theoretical Physics, 16, 682                                                        |
| Vocks, C., &Mann, G. 2003, ApJ, 593, 1134, doi: 10.1086/376682                                                                                       |
| Yoon, P. H., Rhee, T., &Ryu, C.-M. 2006, Journal of Geophysical Research (Space Physics), 111, A09106, doi: 10.1029/2006JA011681                     |
| Zharkova, V. V., &Gordovskyy, M. 2006, ApJ, 651, 553, doi: 10.1086/506423                                                                            |
Yingjie Luo Eduard P. Kontar Debesh Bhattacharjee
Flare-accelerated Electrons in the Kappa Distribution from X-Ray Spectra with the Warm-Target Model
X-ray observations provide important and valuable insights into the
acceleration and propagation of nonthermal electrons during solar flares.
Improved X-ray spectral analysis requires a deeper understanding of the
dynamics of energetic electrons. Previous studies have demonstrated that the
dynamics of accelerated electrons with a few thermal speeds are more complex
than those with significantly higher speeds. To better describe the energetic
electrons after injection, a model considering energy diffusion and
thermalization effects in flare conditions (the warm-target model) has recently
been developed for spectral analysis of hard X-rays. This model has
demonstrated how the low-energy cutoff, which can hardly be constrained in
cold-target modeling, can be determined. However, the power-law form may not be
the most suitable representation of injected electrons. The kappa distribution,
which is proposed as a physical consequence of electron acceleration, has been
applied successfully in RHESSI spectral analysis. In this study, we employ the
kappa-form injected electrons in the warm-target model to analyze two M-class
flares, observed by RHESSI and STIX, respectively. The best-fit results show
that the kappa-form energetic electron spectrum generates lower nonthermal
energy than the power-law form when producing a similar photon spectrum in the
fit range. We also demonstrated that the fit parameters associated with the
kappa-form electron spectrum can be well determined with small uncertainty.
Further, the kappa distribution, which covers the entire electron energy range,
enables the determination of key electron properties such as total electron
number density and average energy at the flare site, providing valuable
information on electron acceleration processes.