Source author record

Eduard P. Kontar

Eduard P. Kontar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2025arXiv

Energy-Containing Electrons in Solar Flares: Improving Hard X-Ray and EUV Diagnostics

Solar flares effectively accelerate particles to non-thermal energies. These accelerated electrons are responsible for energy transport and subsequent emissions in HXR, radio, and UV/EUV radiation. Due to the steeply decreasing electron spectrum, the electron population and consequently the overall flare energetics, are predominantly influenced by low-energy non-thermal electrons. However, deducing the electron distribution in this energy-containing range remains a significant challenge. In this study, we apply the warm-target HXR emission model with kappa-form injected electrons to two well-observed GOES M-class flares. Moreover, we utilize EUV observations to constrain the flaring plasma properties, which enables us to determine the characteristics of accelerated electrons across a range from a few keV to tens of keV. We demonstrate that the warm-target model reliably constrains the properties of flare-associated electrons, even accounting for the uncertainties that had previously been unaddressed. The application of a kappa distribution for the accelerated electrons allows for meaningful comparisons with electron distributions inferred from EUV observations, specifically for energy ranges below the detection threshold of RHESSI. Our results indicate that the accelerated electrons constitute only a small fraction of the total electron population within the flaring region. Moreover, the physical parameters, such as electron escape time and acceleration time scale, inferred from both the warm-target model and the EUV observations further support the scenario in which electrons undergo thermalization within the corona. This study highlights the effectiveness of integrating the warm-target model with EUV observations to accurately characterize energy-containing electrons and their associated acceleration and transport processes.

preprint2021arXiv

Spectral Analysis of Solar Radio Type III Bursts from 20 kHz to 410 MHz

We present the statistical analysis of the spectral response of solar radio type III bursts over the wide frequency range between 20 kHz and 410 MHz. For this purpose, we have used observations that were carried out using both spaced-based (Wind/Waves) and ground-based (Nançay Decameter Array and Nançay Radioheliograph) facilities. In order to compare the flux densities observed by the different instruments, we have carefully calibrated the data and displayed them in Solar Flux Units. The main result of our study is that type III bursts, in the metric to hectometric wavelength range, statistically exhibit a clear maximum of their median radio flux density around 2 MHz. Although this result was already reported by inspecting the spectral profiles of type III bursts in the frequency range 20 kHz - 20 MHz, our study extends such analysis for the first time to metric radio frequencies (i.e., from 20 kHz to 410 MHz) and confirms the maximum spectral response around 2 MHz. In addition, using a simple empirical model we show that the median radio flux $S$ of the studied dataset obeys the polynomial form $Y = 0.04 X^3 - 1.63 X^2 + 16.30 X -41.24$, with $X=\ln{(F_\text{MHz})}$ and with $Y=\ln{(S_\text{SFU})}$. Using the Sittler and Guhathakurtha model for coronal streamers \citep{Sit1999}, we have found that maximum of radio power falls therefore in the range 4 to 10 $R_{\odot}$, depending on whether the type III emissions are assumed to be at the fundamental or the harmonic.

preprint2020arXiv

Density Fluctuations in the Solar Wind Based on Type III Radio Bursts Observed by Parker Solar Probe

Radio waves are strongly scattered in the solar wind, so that their apparent sources seem to be considerably larger and shifted than the actual ones. Since the scattering depends on the spectrum of density turbulence, better understanding of the radio wave propagation provides indirect information on the relative density fluctuations $ε=\langleδn\rangle/\langle n\rangle$ at the effective turbulence scale length. Here, we have analyzed 30 type III bursts detected by Parker Solar Probe (PSP). For the first time, we have retrieved type III burst decay times $τ_{\rm{d}}$ between 1 MHz and 10 MHz thanks to an unparalleled temporal resolution of PSP. We observed a significant deviation in a power-law slope for frequencies above 1 MHz when compared to previous measurements below 1 MHz by the twin-spacecraft Solar TErrestrial RElations Observatory (STEREO) mission. We note that altitudes of radio bursts generated at 1 MHz roughly coincide with an expected location of the Alfvén point, where the solar wind becomes super-Alfvénic. By comparing PSP observations and Monte Carlo simulations, we predict relative density fluctuations $ε$ at the effective turbulence scale length at radial distances between 2.5$R_\odot$ and 14$R_\odot$ to range from $0.22$ and $0.09$. Finally, we calculated relative density fluctuations $ε$ measured in situ by PSP at a radial distance from the Sun of $35.7$~$R_\odot$ during the perihelion \#1, and the perihelion \#2 to be $0.07$ and $0.06$, respectively. It is in a very good agreement with previous STEREO predictions ($ε=0.06-0.07$) obtained by remote measurements of radio sources generated at this radial distance.

preprint2020arXiv

First Observation of a Type II Solar Radio Burst Transitioning Between a Stationary and Drifting State

Standing shocks are believed to be responsible for stationary Type II solar radio bursts, whereas drifting Type II bursts are excited by moving shocks often related to coronal mass ejections (CMEs). Observations of either stationary or drifting Type II bursts are common, but a transition between the two states has not yet been reported. Here, we present a Type II burst which shows a clear, continuous transition from a stationary to a drifting state, the first observation of its kind. Moreover, band splitting is observed in the stationary parts of the burst, as well as intriguing negative and positive frequency-drift fine structures within the stationary emissions. The relation of the radio emissions to an observed jet and a narrow CME was investigated across multiple wavelengths, and the mechanisms leading to the transitioning Type II burst were determined. We find that a jet eruption generates a streamer-puff CME and that the interplay between the CME-driven shock and the streamer is likely to be responsible for the observed radio emissions.

preprint2020arXiv

Radio Echo in the Turbulent Corona and Simulations of Solar Drift-Pair Radio Bursts

Drift-pair bursts are an unusual type of solar low-frequency radio emission, which appear in the dynamic spectra as two parallel drifting bright stripes separated in time. Recent imaging spectroscopy observations allowed for the quantitative characterization of the drifting pairs in terms of source size, position, and evolution. Here, the drift-pair parameters are qualitatively analyzed and compared with the newly-developed Monte Carlo ray-tracing technique simulating radio-wave propagation in the inhomogeneous anisotropic turbulent solar corona. The results suggest that the drift-pair bursts can be formed due to a combination of the refraction and scattering processes, with the trailing component being the result of turbulent reflection (turbulent radio echo). The formation of drift-pair bursts requires an anisotropic scattering with the level of plasma density fluctuations comparable to that in type III bursts, but with a stronger anisotropy at the inner turbulence scale. The anisotropic radio-wave scattering model can quantitatively reproduce the key properties of drift-pair bursts: the apparent source size and its increase with time at a given frequency, the parallel motion of the source centroid positions, and the delay between the burst components. The trailing component is found to be virtually co-spatial and following the main component. The simulations suggest that the drift-pair bursts are likely to be observed closer to the disk center and below 100 MHz due to the effects of free-free absorption and scattering. The exciter of drift-pairs is consistent with propagating packets of whistlers, allowing for a fascinating way to diagnose the plasma turbulence and the radio emission mechanism.

preprint2020arXiv

Spatio-temporal energy partitioning in a non-thermally dominated two-loop solar flare

Solar flares show remarkable variety of the energy partitioning between thermal and nonthermal components. Those with a prominent nonthermal component but only a modest thermal one are particularly well suited to study the direct effect of the nonthermal electrons on plasma heating. Here, we analyze such a well observed, impulsive single-spike nonthermal event, a SOL2013-11-05T035054 solar flare, where the plasma heating can be entirely attributed to the energy losses of these impulsively accelerated electrons. Evolution of the energy budget of thermal and nonthermal components during the flare is analysed using X-ray, microwave, and EUV observations and three-dimensional modeling. The results suggest that (i) the flare geometry is consistent with a two-loop morphology and the magnetic energy is likely released due to interaction between these two loops; (ii) the released magnetic energy converted to the nonthermal energy of accelerated electrons only, which is subsequently converted to the thermal energy of the plasma; (iii) the energy is partitioned in these two flaring loops in comparable amounts; (iv) one of these flaring loops remained relatively tenuous but rather hot, while the other remained relatively cool but denser than the first one. Therefore, this solar flare demonstrates an extreme efficiency of conversion of the free magnetic energy to the nonthermal energy of particle acceleration and the energy flow into two loops from the non-thermal to thermal one with a negligible direct heating.

preprint2016arXiv

Anomalous Cooling of Coronal Loops with Turbulent Suppression of Thermal Conduction

We investigate the impact of turbulent suppression of parallel heat conduction on the cooling of post-flare coronal loops. Depending on the value of the mean free path $λ_T$ associated with the turbulent scattering process, we identify four main cooling scenarios. The overall temperature evolution, from an initial temperature in excess of $10^7$~K, is modeled in each case, highlighting the evolution of the dominant cooling mechanism throughout the cooling process. Comparison with observed cooling times allows the value of $λ_T$ to be constrained, and interestingly this range corresponds to situations where collision-dominated conduction plays a very limited role, or even no role at all, in the cooling of post-flare coronal loops.

preprint2016arXiv

Diagnosing the Source Region of a Solar Burst on 26 September 2011 by Microwave Type III Pairs

This work reports a peculiar and interesting train of microwave type III pair bursts in the impulsive rising phase of a solar flare on 2011 September 26. The observations include radio spectrometers at frequency of 0.80 - 2.00 GHz, hard X-ray (RHESSI and FERMI), EUV images of SWAP/PROBA-2 and magnetogram of HMI/SDO. By using a recently developed method (Tan et al. 2016a), we diagnosed the plasma density, temperature, plasma beta, magnetic field near the source region, the energy of energetic electrons and the distance between the acceleration region and the emission start sites of type III bursts. From the diagnostics, we find that: (1) The plasma density, temperature, magnetic field, and the distance between the acceleration region and the emission start sites almost have no obvious variations during the period of type III pair trains, while the energy of electrons has an obvious peak value which is consistent to the hard X-ray emission. (2) The plasma beta is much higher than an unity showing a highly dynamic process near the emission start site of type III bursts. (3) Although the reversed-slope type III branches drift slower at one order of magnitude than that of the normal type III branches, the related downgoing and upgoing electrons still could have same order of magnitude of energy. These facts indicate that both of the upgoing and downgoing electrons are possibly accelerated by similar mechanism and in a small source region. This diagnostics can help us to understand the microphysics in the source region of solar bursts.

preprint2016arXiv

Narrowband Gyrosynchrotron Bursts: Probing Electron Acceleration in Solar Flares

Recently, in a few case studies we demonstrated that gyrosynchrotron microwave emission can be detected directly from the acceleration region when the trapped electron component is insignificant. For the statistical study reported here, we have identified events with steep (narrowband) microwave spectra that do not show a significant trapped component and at the same time show evidence of source uniformity, which simplifies the data analysis greatly. Initially, we identified a subset of more than 20 radio bursts with such narrow spectra, having low- and high-frequency spectral indices larger than 3 in absolute value. A steep low-frequency spectrum implies that the emission is nonthermal (for optically-thick thermal emission, the spectral index cannot be steeper than 2), and the source is reasonably dense and uniform. A steep high-frequency spectrum implies that no significant electron trapping occurs; otherwise a progressive spectral flattening would be observed. Roughly half of these radio bursts have RHESSI data, which allows for detailed, joint diagnostics of the source parameters and evolution. Based on an analysis of radio-to-X-ray spatial relationships, timing, and spectral fits, we conclude that the microwave emission in these narrowband bursts originates directly from the acceleration regions, which have relatively strong magnetic field, high density, and low temperature. In contrast, the thermal X-ray emission comes from a distinct loop with smaller magnetic field, lower density, but higher temperature. Therefore, these flares occurred likely due to interaction between two (or more) magnetic loops.

preprint2016arXiv

On the origin of 140 GHz emission from the 4 July 2012 solar flare

The sub-THz event observed on the 4 July 2012 with the Bauman Moscow State Technical University Radio Telescope RT-7.5 at 93 and 140~GHz as well as Kislovodsk and Metsähovi radio telescopes, Radio Solar Telescope Network (RSTN), GOES, RHESSI, and SDO orbital stations is analyzed. The spectral flux between 93 and 140 GHz has been observed increasing with frequency. On the basis of the SDO/AIA data the differential emission measure has been calculated. It is shown that the thermal coronal plasma with the temperature above 0.5~MK cannot be responsible for the observed sub-THz flare emission. The non-thermal gyrosynchrotron mechanism can be responsible for the microwave emission near $10$~GHz but the observed millimeter spectral characteristics are likely to be produced by the thermal bremsstrahlung emission from plasma with a temperature of about 0.1~MK.

preprint2016arXiv

Suppression of parallel transport in turbulent magnetized plasmas and its impact on non-thermal and thermal aspects of solar flares

The transport of the energy contained in electrons, both thermal and suprathermal, in solar flares plays a key role in our understanding of many aspects of the flare phenomenon, from the spatial distribution of hard X-ray emission to global energetics. Motivated by recent {\em RHESSI} observations that point to the existence of a mechanism that confines electrons to the coronal parts of flare loops more effectively than Coulomb collisions, we here consider the impact of pitch-angle scattering off turbulent magnetic fluctuations on the parallel transport of electrons in flaring coronal loops. It is shown that the presence of such a scattering mechanism in addition to Coulomb collisional scattering can significantly reduce the parallel thermal and electrical conductivities relative to their collisional values. We provide illustrative expressions for the resulting thermoelectric coefficients that relate the thermal flux and electrical current density to the temperature gradient and the applied electric field. We then evaluate the effect of these modified transport coefficients on the flare coronal temperature that can be attained, on the post-impulsive-phase cooling of heated coronal plasma, and on the importance of the beam-neutralizing return current on both ambient heating and the energy loss rate of accelerated electrons. We also discuss the possible ways in which anomalous transport processes have an impact on the required overall energy associated with accelerated electrons in solar flares.

preprint2015arXiv

3D Radio and X-Ray Modeling and Data Analysis Software: Revealing Flare Complexity

We have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The object-based architecture provides an interactive graphical user interface that allows the user to import photospheric magnetic field maps and perform magnetic field extrapolations to almost instantly generate 3D magnetic field models, to investigate the magnetic topology of these models by interactively creating magnetic field lines and associated magnetic flux tubes, to populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; to investigate the spatial and spectral properties of radio and X-ray emission calculated from the model, and to compare the model-derived images and spectra with observational data. The application integrates shared-object libraries containing fast gyrosynchrotron emission codes developed in FORTRAN and C++, soft and hard X-ray codes developed in IDL, a FORTRAN-based potential-field extrapolation routine and an IDL-based linear force free field extrapolation routine. The interactive interface allows users to add any user-defined radiation code that adheres to our interface standards, as well as user-defined magnetic field extrapolation routines. Here we use this tool to analyze a simple single-loop flare and use the model to constrain the 3D structure of the magnetic flaring loop and 3D spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic fluxtubes obtained from an extrapolation of a magnetogram taken prior to the flare, and compare them with imaging data obtained by SDO, NoRH, and RHESSI instruments. We use this event to illustrate use of the tool for general interpretation of solar flares to address disparate problems in solar physics.

preprint2015arXiv

Multi-thermal representation of the kappa-distribution of solar flare electrons and application to simultaneous X-ray and EUV observations

Acceleration of particles and plasma heating is one of the fundamental problems in solar flare physics. An accurate determination of the spectrum of flare energized electrons over a broad energy range is crucial for our understanding of aspects such as the acceleration mechanism and the total flare energy. Recent years have seen a growing interest in the kappa-distribution as representation of the total spectrum of flare accelerated electrons. In this work we present the kappa-distribution as a differential emission measure. This allows for inferring the electron distribution from X-ray observations and EUV observations by simultaneously fitting the proposed function to RHESSI and SDO/AIA data. This yields the spatially integrated electron spectra of a coronal source between less than 0.1 keV up to several tens of keV. The method is applied to a single-loop GOES C4.1 flare. The results show that the total energy can only be determined accurately by combining RHESSI and AIA observations. Simultaneously fitting the proposed representation of the kappa-distribution reduces the electron number density in the analyzed flare by a factor of ~30 and the total flare energy by a factor of ~5 compared with the commonly used fitting of RHESSI spectra. The spatially integrated electron spectrum of the investigated flare between 0.043 keV and 24 keV is consistent with the combination of a low-temperature (~2 MK) component and a hot (~11 MK) kappa-like component with spectral index 4, reminiscent of solar wind distributions.

preprint2015arXiv

On the speed and acceleration of electron beams triggering interplanetary type III radio bursts

Type III radio bursts are intense radio emissions triggered by beams of energetic electrons often associated with solar flares. These exciter beams propagate outwards from the Sun along an open magnetic field line in the corona and in the interplanetary (IP) medium. We performed a statistical survey of 29 simple and isolated IP type III bursts observed by STEREO/Waves instruments between January 2013 and September 2014. We investigated their time-frequency profiles in order to derive the speed and acceleration of exciter electron beams. We show these beams noticeably decelerate in the IP medium. Obtained speeds range from $\sim$ 0.02c up to $\sim$ 0.35c depending on initial assumptions. It corresponds to electron energies between tens of eV and hundreds of keV, and in order to explain the characteristic energies or speeds of type III electrons ($\sim 0.1$c) observed simultaneously with Langmuir waves at 1 au, the emission of type III bursts near the peak should be predominately at double plasma frequency. Derived properties of electron beams can be used as input parameters for computer simulations of interactions between the beam and the plasma in the IP medium.

preprint2015arXiv

Resonance broadening due to particle scattering and mode-coupling in the quasi-linear relaxation of electron beams

Of particular interest for radio and hard X-ray diagnostics of accelerated electrons during solar flares is the understanding of the basic non-linear mechanisms regulating the relaxation of electron beams propagating in turbulent plasmas. In this work, it is shown that in addition to scattering of beam electrons, scattering of the beam-generated Langmuir waves via for instance mode-coupling, can also result in broadening of the wave-particle resonance. We obtain a resonance-broadened version of weak-turbulence theory with mode-coupling to ion-sound modes. Resonance broadening is presented here as a unified framework which can quantitatively account for the reduction and possible suppression of the beam instability due to background scattering of the beam electrons themselves or due to scattering of the beam-generated Langmuir waves in fluctuating plasmas. Resonance broadening being essentially equivalent to smoothing of the electron phase-space distribution, it is used to construct an intuitive physical picture for the stability of inverted populations of fast electrons that are commonly observed \emph{in-situ} to propagate in the solar-wind.

preprint2015arXiv

Solar and Heliospheric Physics with the Square Kilometre Array

The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.

preprint2015arXiv

Stopping Frequency of Type III Solar Radio Bursts in Expanding Magnetic Flux Tubes

Understanding the properties of type III radio bursts in the solar corona and interplanetary space is one of the best ways to remotely deduce the characteristics of solar accelerated electron beams and the solar wind plasma. One feature of all type III bursts is the lowest frequency they reach (or stopping frequency). This feature reflects the distance from the Sun that an electron beam can drive the observable plasma emission mechanism. The stopping frequency has never been systematically studied before from a theoretical perspective. Using numerical kinetic simulations, we explore the different parameters that dictate how far an electron beam can travel before it stops inducing a significant level of Langmuir waves, responsible for plasma radio emission. We use the quasilinear approach to model self-consistently the resonant interaction between electrons and Langmuir waves in inhomogeneous plasma, and take into consideration the expansion of the guiding magnetic flux tube and the turbulent density of the interplanetary medium. We find that the rate of radial expansion has a significant effect on the distance an electron beam travels before enhanced leves of Langmuir waves, and hence radio waves, cease. Radial expansion of the guiding magnetic flux tube rarefies the electron stream to the extent that the density of non-thermal electrons is too low to drive Langmuir wave production. The initial conditions of the electron beam have a significant effect, where decreasing the beam density or increasing the spectral index of injected electrons would cause higher type III stopping frequencies. We also demonstrate how the intensity of large-scale density fluctuations increases the highest frequency that Langmuir waves can be driven by the beam and how the magnetic field geometry can be the cause of type III bursts only observed at high coronal frequencies.

preprint2014arXiv

On the variation of solar flare coronal x-ray source sizes with energy

Observations with {\em RHESSI} have enabled the detailed study of the structure of dense hard X-ray coronal sources in solar flares. The variation of source extent with electron energy has been discussed in the context of streaming of non-thermal particles in a one-dimensional cold-target model, and the results used to constrain both the physical extent of, and density within, the electron acceleration region. Here we extend this investigation to a more physically realistic model of electron transport that takes into account the finite temperature of the ambient plasma, the initial pitch-angle distribution of the accelerated electrons, and the effects of collisional pitch-angle scattering. The finite temperature results in the thermal diffusion of electrons, that leads to the observationally-inferred value of the acceleration region volume being an overestimate of its true value. The different directions of the electron trajectories, a consequence of both the non-zero injection pitch-angle and scattering within the target, cause the projected propagation distance parallel to the guiding magnetic field to be reduced, so that a one-dimensional interpretation can overestimate the actual density by a factor of up to $\sim 6$. The implications of these results for the determination of acceleration region properties (specific acceleration rate, filling factor, etc.) are discussed.

preprint2014arXiv

The Formation of Kappa-Distribution Accelerated Electron Populations in Solar Flares

Driven by recent RHESSI observations of confined loop-top hard X-ray sources in solar flares, we consider stochastic acceleration of electrons in the presence of Coulomb collisions. If electron escape from the acceleration region can be neglected, the electron distribution function is determined by a balance between diffusive acceleration and collisions. Such a scenario admits a stationary solution for the electron distribution function that takes the form of a kappa distribution. We show that the evolution toward this kappa distribution involves a "wave front" propagating forwards in velocity space, so that electrons of higher energy are accelerated later; the acceleration time scales with energy according to $τ_{\rm acc} \sim E^{3/2}$. At sufficiently high energies escape from the finite-length acceleration region will eventually dominate. For such energies, the electron velocity distribution function is obtained by solving a time-dependent Fokker-Planck equation in the "leaky-box" approximation. Solutions are obtained in the limit of a small escape rate from an acceleration region that can effectively be considered a thick target.

preprint2014arXiv

The Low-High-Low Trend of Type III Radio Burst Starting Frequencies and Solar Flare Hard X-rays

Using simultaneous X-ray and radio observations from solar flares, we investigate the link between the type III radio burst starting frequency and hard X-ray spectral index. For a proportion of events the relation derived between the starting height (frequency) of type III radio bursts and the electron beam velocity spectral index (deduced from X-rays) is used to infer the spatial properties (height and size) of the electron beam acceleration region. Both quantities can be related to the distance travelled before an electron beam becomes unstable to Langmuir waves. To obtain a list of suitable events we considered the RHESSI catalogue of X-ray flares and the Phoenix 2 catalogue of type III radio bursts. From the 200 events that showed both type III and X-ray signatures, we selected 30 events which had simultaneous emission in both wavelengths, good signal to noise in the X-ray domain and > 20 seconds duration. We find that > 50 % of the selected events show a good correlation between the starting frequencies of the groups of type III bursts and the hard X-ray spectral indices. A low-high-low trend for the starting frequency of type III bursts is frequently observed. Assuming a background electron density model and the thick target approximation for X-ray observations, this leads to a correlation between starting heights of the type III emission and the beam electron spectral index. Using this correlation we infer the altitude and vertical extents of the flare acceleration regions. We find heights from 183 Mm down to 25 Mm while the sizes range from 13 Mm to 2 Mm. These values agree with previous work that places an extended flare acceleration region high in the corona. We analyse the assumptions required and explore possible extensions to our assumed model. We discuss these results with respect to the acceleration heights and sizes derived from X-ray observations alone.

preprint2013arXiv

Electron Distribution Functions in Solar Flares from combined X-ray and EUV Observations

Simultaneous solar flare observations with SDO and RHESSI provide spatially resolved information about hot plasma and energetic particles in flares. RHESSI allows the properties of both hot (> 8 MK) thermal plasma and nonthermal electron distributions to be inferred, while SDO/AIA is more sensitive to lower temperatures. We present and implement a new method to reconstruct electron distribution functions from SDO/AIA data. The combined analysis of RHESSI and AIA data allows the electron distribution function to be inferred over the broad energy range from ~0.1 keV up to a few tens of keV. The analysis of two well observed flares suggests that the distributions in general agree to within a factor of three when the RHESSI values are extrapolated into the intermediate range 1-3 keV, with AIA systematically predicting lower electron distributions. Possible instrumental and numerical effects, as well as potential physical origins for this discrepancy are discussed. The inferred electron distribution functions in general show one or two nearly Maxwellian components at energies below ~ 15 keV and a non-thermal tail above.

preprint2013arXiv

Implications for electron acceleration and transport from non-thermal electron rates at looptop and footpoint sources in solar flares

The interrelation of hard X-ray (HXR) emitting sources and the underlying physics of electron acceleration and transport presents one of the major questions in high-energy solar flare physics. Spatially resolved observations of solar flares often demonstrate the presence of well-separated sources of bremsstrahlung emission, so-called coronal and footpoint sources. Using spatially resolved X-ray observations by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and recently improved imaging techniques, we investigate in detail the spatially resolved electron distributions in a few well-observed solar flares. The selected flares can be interpreted as having a standard geometry with chromospheric HXR footpoint sources related to thick-target X-ray emission and the coronal sources characterised by a combination of thermal and thin-target bremsstrahlung. Using imaging spectroscopy techniques, we deduce the characteristic electron rates and spectral indices required to explain the coronal and footpoint X-ray sources. We found that, during the impulsive phase, the electron rate at the looptop is several times (a factor of 1.7-8) higher than at the footpoints. The results suggest that a sufficient number of electrons accelerated in the looptop explain the precipitation into the footpoints and imply that electrons accumulate in the looptop. We discuss these results in terms of magnetic trapping, pitch-angle scattering, and injection properties. Our conclusion is that the accelerated electrons must be subject to magnetic trapping and/or pitch-angle scattering, keeping a fraction of the population trapped inside the coronal loops. These findings put strong constraints on the particle transport in the coronal source and provide quantitative limits on deka-keV electron trapping/scattering in the coronal source.

preprint2013arXiv

Probing Dynamics of Electron Acceleration with Radio and X-ray Spectroscopy, Imaging, and Timing in the 2002 Apr 11 Solar Flare

Based on detailed analysis of radio and X-ray observations of a flare on 2002 April 11 augmented by realistic 3D modeling, we have identified a radio emission component produced directly at the flare acceleration region. This acceleration region radio component has distinctly different (i) spectrum, (ii) light curves, (iii) spatial location, and, thus, (iv) physical parameters from those of the separately identified, trapped or precipitating electron components. To derive evolution of physical parameters of the radio sources we apply forward fitting of the radio spectrum time sequence with the gyrosynchrotron source function with 5 to 6 free parameters. At the stage when the contribution from the acceleration region dominates the radio spectrum, the X-ray- and radio- derived electron energy spectral indices agree well with each other. During this time the maximum energy of the accelerated electron spectrum displays a monotonic increase with time from ~ 300 keV to ~ 2 MeV over roughly one minute duration indicative of an acceleration process in the form of growth of the power-law tail; the fast electron residence time in the acceleration region is about 2 - 4 s, which is much longer than the time of flight and so requires a strong diffusion mode there to inhibit free-streaming propagation. The acceleration region has a relatively strong magnetic field, B ~ 120 G, and a low thermal density, less than 2x10^9 cm^-3. These acceleration region properties are consistent with a stochastic acceleration mechanism.

preprint2013arXiv

Temporal variations of X-ray solar flare loops: length, corpulence, position, temperature, plasma pressure and spectra

The spatial and spectral properties of three solar flare coronal X-ray loops are studied before, during and after the peak X-ray emission. Using observations from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI), we deduce the temporal changes in emitting X-ray length, corpulence, volume, position, number density and thermal pressure. We observe a decrease in the loop length, width and volume before the X-ray peak, and an increasing number density and thermal pressure. After the X-ray peak, volume increases and loop corpulence grows due to an increasing width. The volume variations are more pronounced than the position variations, often known as magnetic line contraction. We believe this is the first dedicated study of the temporal evolution of X-ray loop lengths and widths. Collectively, the observations also show for the first time three temporal phases given by peaks in temperature, X-ray emission and thermal pressure, with minimum volume coinciding with the X-ray peak. Although the volume of the flaring plasma decreases before the peak in X-ray emission, the relationship between temperature and volume does not support simple compressive heating in a collapsing magnetic trap model. Within a low beta plasma, shrinking loop widths perpendicular to the guiding field can be explained by squeezing the magnetic field threading the region. Plasma heating leads to chromospheric evaporation and growing number density, producing increasing thermal pressure and decreasing loop lengths as electrons interact at shorter distances and we believe after the X-ray peak, the increasing loop corpulence.

preprint2012arXiv

A Classification Scheme For Turbulent Acceleration Processes In Solar Flares

We establish a classification scheme for stochastic acceleration models involving low-frequency plasma turbulence in a strongly magnetized plasma. This classification takes into account both the properties of the accelerating electromagnetic field, and the nature of the transport of charged particles in the acceleration region. We group the acceleration processes as either resonant, non-resonant or resonant-broadened, depending on whether the particle motion is free-streaming along the magnetic field, diffusive or a combination of the two. Stochastic acceleration by moving magnetic mirrors and adiabatic compressions are addressed as illustrative examples. We obtain expressions for the momentum-dependent diffusion coefficient $D(p)$, both for general forms of the accelerating force and for the situation when the electromagnetic force is wave-like, with a specified dispersion relation $ω=ω(k)$. Finally, for models considered, we calculate the energy-dependent acceleration time, a quantity that can be directly compared with observations of the time profile of the radiation field produced by the accelerated particles, such as during solar flares.

preprint2012arXiv

Density fluctuations and the acceleration of electrons by beam-generated Langmuir waves in the solar corona

Non-thermal electron populations are observed throughout the heliosphere. The relaxation of an electron beam is known to produce Langmuir waves which, in turn, may substantially modify the electron distribution function. As the Langmuir waves are refracted by background density gradients and as the solar and heliospheric plasma density is naturally perturbed with various levels of inhomogeneity, the interaction of Langmuir waves with non-thermal electrons in inhomogeneous plasmas is an important topic. We investigate the role played by ambient density fluctuations on the beam-plasma relaxation, focusing on the effect of acceleration of beam electrons. The scattering of Langmuir waves off turbulent density fluctuations is modeled as a wavenumber diffusion process which is implemented in numerical simulations of the one-dimensional quasilinear kinetic equations describing the beam relaxation. The results show that a substantial number of beam electrons are accelerated when the diffusive time scale in wavenumber space τ_D is of the order of the quasilinear time scale τ_ql, while when τ_D << τ_ql, the beam relaxation is suppressed. Plasma inhomogeneities are therefore an important means of energy redistribution for waves and hence electrons, and so must be taken into account when interpreting, for example, hard X-ray or Type III emission from flare-accelerated electrons.

preprint2012arXiv

Determination of the Acceleration Region Size in a Loop-structured Solar Flare

In order to study the acceleration and propagation of bremsstrahlung-producing electrons in solar flares, we analyze the evolution of the flare loop size with respect to energy at a variety of times. A GOES M3.7 loop-structured flare starting around 23:55 on 2002 April 14 is studied in detail using \textit{Ramaty High Energy Solar Spectroscopic Imager} (\textit{RHESSI}) observations. We construct photon and mean-electron-flux maps in 2-keV energy bins by processing observationally-deduced photon and electron visibilities, respectively, through several image-processing methods: a visibility-based forward-fit (FWD) algorithm, a maximum entropy (MEM) procedure and the uv-smooth (UVS) approach. We estimate the sizes of elongated flares (i.e., the length and width of flaring loops) by calculating the second normalized moments of the intensity in any given map. Employing a collisional model with an extended acceleration region, we fit the loop lengths as a function of energy in both the photon and electron domains. The resulting fitting parameters allow us to estimate the extent of the acceleration region which is between $\sim 13 \rm{arcsec}$ and $\sim 19 \rm{arcsec}$. Both forward-fit and uv-smooth algorithms provide substantially similar results with a systematically better fit in the electron domain.The consistency of the estimates from these methods provides strong support that the model can reliably determine geometric parameters of the acceleration region. The acceleration region is estimated to be a substantial fraction ($\sim 1/2$) of the loop extent, indicating that this dense flaring loop incorporates both acceleration and transport of electrons, with concurrent thick-target bremsstrahlung emission.

preprint2012arXiv

Evolution of the Solar Flare Energetic Electrons in the Inhomogeneous Inner Heliosphere

Solar flare accelerated electrons escaping into the interplanetary space and seen as type III solar radio bursts are often detected near the Earth. Using numerical simulations we consider the evolution of energetic electron spectrum in the inner heliosphere and near the Earth. The role of Langmuir wave generation, heliospheric plasma density fluctuations, and expansion of magnetic field lines on the electron peak flux and fluence spectra is studied to predict the electron properties as could be observed by Solar Orbiter and Solar Probe Plus. Considering various energy loss mechanisms we show that the substantial part of the initial energetic electron energy is lost via wave-plasma processes due to plasma inhomogeneity. For the parameters adopted, the results show that the electron spectra changes mostly at the distances before $\sim20 R_\odot$. Further into the heliosphere, the electron flux spectra of electrons forms a broken power-law relatively similar to what is observed at 1 AU.

preprint2012arXiv

Numerical simulations of chromospheric hard X-ray source sizes in solar flares

X-ray observations are a powerful diagnostic tool for transport, acceleration, and heating of electrons in solar flares. Height and size measurements of X-ray footpoints sources can be used to determine the chromospheric density and constrain the parameters of magnetic field convergence and electron pitch-angle evolution. We investigate the influence of the chromospheric density, magnetic mirroring and collisional pitch-angle scattering on the size of X-ray sources. The time-independent Fokker-Planck equation for electron transport is solved numerically and analytically to find the electron distribution as a function of height above the photosphere. From this distribution, the expected X-ray flux as a function of height, its peak height and full width at half maximum are calculated and compared with RHESSI observations. A purely instrumental explanation for the observed source size was ruled out by using simulated RHESSI images. We find that magnetic mirroring and collisional pitch-angle scattering tend to change the electron flux such that electrons are stopped higher in the atmosphere compared with the simple case with collisional energy loss only. However, the resulting X-ray flux is dominated by the density structure in the chromosphere and only marginal increases in source width are found. Very high loop densities (>10^{11} cm^{-3}) could explain the observed sizes at higher energies, but are unrealistic and would result in no footpoint emission below about 40 keV, contrary to observations. We conclude that within a monolithic density model the vertical sizes are given mostly by the density scale-height and are predicted smaller than the RHESSI results show.

preprint2012arXiv

RHESSI and SDO/AIA observations of the chromospheric and coronal plasma parameters during a solar flare

X-ray and EUV observations are an important diagnostic of various plasma parameters of the solar atmosphere during solar flares. Soft X-ray and EUV observations often show coronal sources near the top of flaring loops, while hard X-ray emission is mostly observed from chromospheric footpoints. Combining RHESSI with simultaneous SDO/AIA observations, it is possible for the first time to determine the density, temperature, and emission profile of the solar atmosphere over a wide range of heights during a flare, using two independent methods. Here we analyze a near limb event during the first of three hard X-ray peaks. The emission measure, temperature, and density of the coronal source is found using soft X-ray RHESSI images while the chromospheric density is determined using RHESSI visibility analysis of the hard X-ray footpoints. A regularized inversion technique is applied to AIA images of the flare to find the differential emission measure (DEM). Using DEM maps we determine the emission and temperature structure of the loop, as well as the density, and compare it with RHESSI results. The soft X-ray and hard X-ray sources are spatially coincident with the top and bottom of the EUV loop, but the bulk of the EUV emission originates from a region without co-spatial RHESSI emission. The temperature analysis along the loop indicates that the hottest plasma is found near the coronal loop top source. The EUV observations suggest that the density in the loop legs increases with increasing height while the temperature remains constant within uncertainties.

preprint2011arXiv

Cold, tenuous solar flare: acceleration without heating

We report the observation of an unusual cold, tenuous solar flare, which reveals itself via numerous and prominent non-thermal manifestations, while lacking any noticeable thermal emission signature. RHESSI hard X-rays and 0.1-18 GHz radio data from OVSA and Phoenix-2 show copious electron acceleration (10^35 electrons per second above 10 keV) typical for GOES M-class flares with electrons energies up to 100 keV, but GOES temperatures not exceeding 6.1 MK. The imaging, temporal, and spectral characteristics of the flare have led us to a firm conclusion that the bulk of the microwave continuum emission from this flare was produced directly in the acceleration region. The implications of this finding for the flaring energy release and particle acceleration are discussed.

preprint2011arXiv

Hard X-ray footpoint sizes and positions as diagnostics of flare accelerated energetic electrons in the low solar atmosphere

The hard X-ray (HXR) emission in solar flares comes almost exclusively from a very small part of the flaring region, the footpoints of magnetic loops. Using RHESSI observations of solar flare footpoints, we determine the radial positions and sizes of footpoints as a function of energy in six near-limb events to investigate the transport of flare accelerated electrons and the properties of the chromosphere. HXR visibility forward fitting allows to find the positions/heights and the sizes of HXR footpoints along and perpendicular to the magnetic field of the flaring loop at different energies in the HXR range. We show that in half of the analyzed events, a clear trend of decreasing height of the sources with energy is found. Assuming collisional thick-target transport, HXR sources are located between 600 and 1200 km above the photosphere for photon energies between 120 and 25 keV respectively. In the other events, the position as a function of energy is constant within the uncertainties. The vertical sizes (along the path of electron propagation) range from 1.3 to 8 arcseconds which is up to a factor 4 larger than predicted by the thick-target model even in events where the positions/heights of HXR sources are consistent with the collisional thick-target model. Magnetic mirroring, collisional pitch angle scattering and X-ray albedo are discussed as potential explanations of the findings.

preprint2011arXiv

Height structure of X-ray, EUV and white-light emission in a solar flare

The bulk of solar flare emission originates from very compact sources located in the lower solar atmosphere and seen in various wavelength ranges: near optical, UV, EUV, soft and hard X-rays, and gamma-ray emission, yet very few spatially resolved imaging observations to determine the structure of these regions exist. We investigate the above-the-photosphere heights of hard X-ray (HXR), EUV and white-light continuum sources in the low atmosphere and the corresponding densities at these heights. Simultaneous EUV/continuum images from SDO and HXR RHESSI images are compared to study a well observed gamma-ray limb flare. Using RHESSI X-ray visibilities we determine the height of the HXR sources as a function of energy above the photosphere. Co-aligning AIA/SDO and HMI/SDO images with RHESSI we infer, for the first time, the heights and characteristic densities of HXR, EUV and continuum sources in a flaring footpoint. 35-100 keV HXR sources are found at heights between 1.7 and 0.8 Mm above the photosphere, below the white-light continuum emission which appears at heights 1.5-3 Mm, and the peak of EUV emission originating near 3 Mm. The EUV emission locations are consistent with energy deposition from low energy electrons of ~12 keV occurring in the top layers of the fully ionized chromosphere/low corona and not by >20 keV electrons that produce HXR footpoints in the lower neutral chromosphere. The maximum of white-light emission appears between the HXR and EUV emission, presumably in the transition between ionized and neutral atmospheres suggesting free-bound and free-free continuum emission. We note that the energy deposited by low energy electrons is sufficient to explain the energetics of optical and UV emissions.

preprint2011arXiv

Implications of X-ray Observations for Electron Acceleration and Propagation in Solar Flares

High-energy X-rays and gamma-rays from solar flares were discovered just over fifty years ago. Since that time, the standard for the interpretation of spatially integrated flare X-ray spectra at energies above several tens of keV has been the collisional thick-target model. After the launch of the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in early 2002, X-ray spectra and images have been of sufficient quality to allow a greater focus on the energetic electrons responsible for the X-ray emission, including their origin and their interactions with the flare plasma and magnetic field. The result has been new insights into the flaring process, as well as more quantitative models for both electron acceleration and propagation, and for the flare environment with which the electrons interact. In this article we review our current understanding of electron acceleration, energy loss, and propagation in flares. Implications of these new results for the collisional thick-target model, for general flare models, and for future flare studies are discussed.

preprint2011arXiv

Recent Advances in Understanding Particle Acceleration Processes in Solar Flares

We review basic theoretical concepts in particle acceleration, with particular emphasis on processes likely to occur in regions of magnetic reconnection. Several new developments are discussed, including detailed studies of reconnection in three-dimensional magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell approaches are used and results compared. While these studies show considerable promise in accounting for the various observational manifestations of solar flares, they are limited by a number of factors, mostly relating to available computational power. Not the least of these issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated particles themselves on the environment in which they are accelerated. A brief prognosis for future advancement is offered.

preprint2010arXiv

Relationship between Hard and Soft X-ray Emission Components of a Solar Flare

X-ray observations of solar flares routinely reveal an impulsive high-energy and a gradual low-energy emission component, whose relationship is one of the key issues of solar flare study. The gradual and impulsive emission components are believed to be associated with, respectively, the thermal and nonthermal components identified in spectral fitting. In this paper, a prominent about 50 second hard X-ray (HXR) pulse of a simple GOES class C7.5 flare on 20 February 2002 is used to study the association between high energy, non-thermal and impulsive evolution, and low energy, thermal and gradual evolution. We use regularized methods to obtain time derivatives of photon fluxes to quantify the time evolution as a function of photon energy, obtaining a break energy between impulsive and gradual behavior. These break energies are consistent with a constant value of about 11 keV in agreement with those found spectroscopically between thermal and non-thermal components, but the relative errors of the former are greater than 15% and much greater than the a few percent errors found from the spectral fitting. These errors only weakly depend on assuming an underlying spectral model for the photons, pointing to the current data being inadequate to reduce the uncertainties rather than there being a problem associated with an assumed model. The time derivative method is used to test for the presence of a 'pivot energy' in this flare. Although these pivot energies are marginally consistent with a constant value of about 9 keV, its values in the HXR rise phase appear to be lower than those in the decay phase.

preprint2010arXiv

The influence of albedo on the size of hard X-ray flare sources

Context: Hard X-rays from solar flares are an important diagnostic of particle acceleration and transport in the solar atmosphere. Any observed X-ray flux from on-disc sources is composed of direct emission plus Compton backscattered photons (albedo). This affects both the observed spectra and images as well as the physical quantities derived from them such as the spatial and spectral distributions of accelerated electrons or characteristics of the solar atmosphere. Aims: We propose a new indirect method to measure albedo and to infer the directivity of X-rays in imaging using RHESSI data. Methods: Visibility forward fitting is used to determine the size of a disc event observed by RHESSI as a function of energy. This is compared to the sizes of simulated sources from a Monte Carlo simulation code of photon transport in the chromosphere for different degrees of downward directivity and true source sizes to find limits on the true source size and the directivity. Results: The observed full width half maximum of the source varies in size between 7.4 arcsec and 9.1 arcsec with the maximum between 30 and 40 keV. Such behaviour is expected in the presence of albedo and is found in the simulations. A source size smaller than 6 arcsec is improbable for modest directivities and the true source size is likely to be around 7 arcsec for small directivities. Conclusions: While it is difficult to image the albedo patch directly, the effect of backscattered photons on the observed source size can be estimated. The increase in source size caused by albedo has to be accounted for when computing physical quantities that include the size as a parameter such as flare energetics. At the same time, the study of the albedo signature provides vital information about the directivity of X-rays and related electrons.

preprint2009arXiv

Sub-THz radiation mechanisms in solar flares

Observations in the sub-THz range of large solar flares have revealed a mysterious spectral component increasing with frequency and hence distinct from the microwave component commonly accepted to be produced by gyrosynchrotron (GS) emission from accelerated electrons. Evidently, having a distinct sub-THz component requires either a distinct emission mechanism (compared to the GS one), or different properties of electrons and location, or both. We find, however, that the list of possible emission mechanisms is incomplete. This Letter proposes a more complete list of emission mechanisms, capable of producing a sub-THz component, both well-known and new in this context and calculates a representative set of their spectra produced by a) free-free emission, b) gyrosynchrotron emission, c) synchrotron emission from relativistic positrons/electrons, d) diffusive radiation, and e) Cherenkov emission. We discuss the possible role of the mechanisms in forming the sub-THz emission and emphasize their diagnostics potential for flares.

preprint2005arXiv

Compton backscattered and primary X-rays from solar flares: angle dependent Green's function correction for photospheric albedo

The observed hard X-ray (HXR) flux spectrum $I(ε)$ from solar flares is a combination of primary bremsstrahlung photons $I_P(ε)$ with a spectrally modified component from photospheric Compton backscatter of downward primary emission. The latter can be significant, distorting or hiding the true features of the primary spectrum which are key diagnostics for acceleration and propagation of high energy electrons and of their energy budget. For the first time in solar physics, we use a Green's function approach to the backscatter spectral deconvolution problem, constructing a Green's matrix including photoelectric absorption. This approach allows spectrum-independent extraction of the primary spectrum for several HXR flares observed by the {\it Ramaty High Energy Solar Spectroscopic Imager} (RHESSI). We show that the observed and primary spectra differ very substantially for flares with hard spectra close to the disk centre. We show in particular that the energy dependent photon spectral index $γ(ε)=-d \log I/d \log ε$ is very different for $I_P(ε)$ and for $I(ε)$ and that inferred mean source electron spectra ${\bar F}(E)$ differ greatly. Even for a forward fitting of a parametric ${\bar F}(E)$ to the data, a clear low-energy cutoff required to fit $I(ε)$ essentially disappears when the fit is to $I_P(ε)$ - i.e. when albedo correction is included. The self-consistent correction for backscattered photons is thus shown to be crucial in determining the energy spectra of flare accelerated electrons, and hence their total number and energy.