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.

Understanding how the release of stored magnetic energy contributes to the generation of non-thermal high energy particles during solar flares is an important open problem in solar physics. Magnetic reconnection plays a fundamental role in the energy release and conversion processes taking place during flares. A common approach for investigating particle acceleration is to use test particles in fields derived from magnetohydrodynamic (MHD) simulations of reconnection. These MHD simulations use anomalous resistivities that are much larger than the Spitzer resistivity based on Coulomb collisions. The processes leading to enhanced resistivity should also affect the test particles. We explore the link between resistivity and particle orbits building on a previous study using a 2D MHD simulation of magnetic reconnection. This paper extends the previous investigation to a 3D magnetic reconnection configuration and to study the effect on test particle orbits. We carried out orbit calculations using a 3D MHD simulation of separator reconnection. We use the relativistic guiding centre approximation including stochastic pitch angle scattering. The effects of varying the resistivity and the models for pitch angle scattering on particle orbit trajectories, final positions, energy spectra, final pitch angle distribution, and orbit duration are all studied in detail. Pitch angle scattering widens collimated beams of orbit trajectories, allowing orbits to access previously unaccessible field lines; this causes final positions to spread to topological structures that were previously inaccessible. Scattered orbit energy spectra are found to be predominantly affected by the level of anomalous resistivity, with the pitch angle scattering model only playing a role in isolated cases. Scattering is found to play a crucial role in determining the pitch angle and orbit duration distributions.

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 $\epsilon=\langle\delta 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 $\tau_{\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\'{e}n point, where the solar wind becomes super-Alfv\'{e}nic. By comparing PSP observations and Monte Carlo simulations, we predict relative density fluctuations $\epsilon$ 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 $\epsilon$ 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 ($\epsilon=0.06-0.07$) obtained by remote measurements of radio sources generated at this radial distance.

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.

Drift pairs are an unusual type of fine structure sometimes observed in dynamic spectra of solar radio emission. They appear as two identical short narrowband drifting stripes separated in time; both positive and negative frequency drifts are observed. Using the Low Frequency Array (LOFAR), we report unique observations of a cluster of drift pair bursts in the frequency range of 30-70 MHz made on 12 July 2017. Spectral imaging capabilities of the instrument have allowed us for the first time to resolve the temporal and frequency evolution of the source locations and sizes at a fixed frequency and along the drifting pair components. Sources of two components of a drift pair have been imaged and found to propagate in the same direction along nearly the same trajectories. Motion of the second component source is delayed in time with respect to that of the first one. The source trajectories can be complicated and non-radial; positive and negative frequency drifts correspond to opposite propagation directions. The drift pair bursts with positive and negative frequency drifts, as well as the associated broadband type-III-like bursts, are produced in the same regions. The visible source velocities are variable from zero to a few $10^4$ (up to ${\sim 10^5}$) km/s, which often exceeds the velocities inferred from the drift rate ($\sim 10^4$ km/s). The visible source sizes are of about $10'-18'$; they are more compact than typical type III sources at the same frequencies. The existing models of drift pair bursts cannot adequately explain the observed features. We discuss the key issues that need to be addressed, and in particular the anisotropic scattering of the radio waves. The broadband bursts observed simultaneously with the drift pairs differ in some aspects from common type III bursts and may represent a separate type of emission.

Solar type III radio bursts are excited by electron beams propagating outward from the Sun. The flux of type III radio burst has a time profile of rising and decay phase at a given frequency, which has been actively studied since 1970s. Several factors that may influence the duration of a type III radio burst has been proposed. However, the major cause of the duration is still an open question. In this work, to study the dominant cause of the duration, we investigate the source positions of the front edge, the peak, and the tail edge in the dynamic spectrum of a single and clear type III radio burst. The duration of this type III burst at a given frequency is about 3 second for decameter wave. The beam-formed observations by the LOw-Frequency ARray (LOFAR) are used, which can provide the radio source positions and the dynamic spectra at the same time. We find that, for this burst, the source positions of the front edge, the peak, and the tail edge split with each other spatially. The radial speed of the front edge, the peak, and the tail edge is 0.42 c, 0.25 c, and 0.16 c, respectively. We estimate the influences of the corona density fluctuation and the electron-velocity dispersion on the duration, and the scattering effect by comparison with a few short-duration bursts from the same region. The analysis yields that, in the frequency range of 30 - 41 MHz, the electron-velocity dispersion is the dominant factor that determines the time duration of type III radio bursts with long duration, while scattering may play important role in the duration of short bursts.

The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is presented, based on the Fokker-Planck and Langevin equations of radio-wave transport in a medium containing anisotropic electron density fluctuations. Using a numerical treatment based on this model, we investigate the characteristic source sizes and burst decay times for Type III solar radio bursts. Comparison of the simulations with the observations of solar radio bursts shows that predominantly perpendicular density fluctuations in the solar corona are required, with an anisotropy factor $\sim 0.3$ for sources observed at around 30~MHz. The simulations also demonstrate that the photons are isotropized near the region of primary emission, but the waves are then focused by large-scale refraction, leading to plasma radio emission directivity that is characterized by a half-width-half-maximum of about 40~degrees near 30~MHz. The results are applicable to various solar radio bursts produced via plasma emission.

One of the key problems in solar flare physics is the determination of the low-energy cut-off; the value that determines the energy of nonthermal electrons and hence flare energetics. We discuss different approaches to determine the low-energy cut-off in the spectrum of accelerated electrons: (i) the total electron number model, (ii) the time-of-flight model (based on the equivalence of the time-of-flight and the collisional deflection time); (iii) the warm target model of Kontar et al.~(2015), and (iv) the model of the spectral cross-over between thermal and nonthermal components. We find that the first three models are consistent with a low-energy cutoff with a mean value of $\approx 10$ keV, while the cross-over model provides an upper limit for the low-energy cutoff with a mean value of $ \approx 21$ keV. Combining the first three models we find that the ratio of the nonthermal energy to the dissipated magnetic energy in solar flares has a mean value of $q_E=0.57\pm0.08$, which is consistent with an earlier study based on the simplified approximation of the warm target model alone ($q_E=0.51\pm0.17$). This study corroborates the self-consistency between three different low-energy cutoff models in the calculation of nonthermal flare energies.

The solar atmosphere is dominated by loops of magnetic flux which connect the multi-million-degree corona to the much cooler chromosphere. The temperature and density structure of quasi-static loops is determined by the continuous flow of energy from the hot corona to the lower solar atmosphere. Loop scaling laws provide relationships between global properties of the loop (such as peak temperature, pressure, and length); they follow from the physical variable dependencies of various terms in the energy equation, and hence the form of the loop scaling law provides insight into the key physics that controls the loop structure. Traditionally, scaling laws have been derived under the assumption of collision-dominated thermal conduction. Here we examine the impact of different regimes of thermal conduction -- collision-dominated, turbulence-dominated, and free-streaming -- on the form of the scaling laws relating the loop temperature and heating rate to its pressure and half-length. We show that the scaling laws for turbulence-dominated conduction are fundamentally different than those for collision-dominated and free-streaming conduction, inasmuch as the form of the scaling laws now depend primarily on conditions at the low-temperature, rather than high-temperature, part of the loop. We also establish regimes in temperature and density space in which each of the applicable scaling laws prevail.

During solar flares, a large fraction of the released magnetic energy is carried by energetic electrons that transfer and deposit energy in the Sun's atmosphere. Electron transport is often approximated by a cold thick-target model (CTTM), assuming that electron energy is much larger than the temperature of the ambient plasma, and electron energy evolution is modeled as a systematic loss. Using kinetic modeling of electrons, we re-evaluate the transport and deposition of flare energy. Using a full collisional warm-target model (WTM), we account for electron thermalization and for the properties of the ambient coronal plasma such as its number density, temperature and spatial extent. We show that the deposition of non-thermal electron energy in the lower atmosphere is highly dependent on the properties of the flaring coronal plasma. In general, thermalization and a reduced WTM energy loss rate leads to an increase of non-thermal energy transferred to the chromosphere, and the deposition of non-thermal energy at greater depths. The simulations show that energy is deposited in the lower atmosphere initially by high energy non-thermal electrons, and later by lower energy non-thermal electrons that partially or fully thermalize in the corona, over timescales of seconds, unaccounted for in previous studies. This delayed heating may act as a diagnostic of both the injected non-thermal electron distribution and the coronal plasma, vital for constraining flare energetics.

The problem of beam propagation in a plasma with small scale and low intensity inhomogeneities is investigated. It is shown that the electron beam propagates in a plasma as a beam-plasma structure and is a source of Langmuir waves. The plasma inhomogeneity changes the spatial distribution of the waves. The spatial distribution of the waves is fully determined by the distribution of plasma inhomogeneities. The possible applications to the theory of radio emission associated with electron beams are discussed.

Dynamics of an spatially limited electron beam in the inhomogeneous solar corona plasma is considered in the framework of weak turbulence theory when the temperature of the beam significantly exceeds that of surrounding plasma. The numerical solution of kinetic equations manifests that generally the beam accompanied by Langmuir waves propagates as a beam-plasma structure with a decreasing velocity. Unlike the uniform plasma case the structure propagates with the energy losses in the form of Langmuir waves. The results obtained are compared with the results of observations of type III bursts. It is shown that the deceleration of type III sources can be explained by the corona inhomogeneity. The frequency drift rates of the type III sources are found in a good agreement with the numerical results of beam dynamics.

The dynamics of a hot electron cloud in the solar corona-like plasma based on the numerical solution of kinetic equations of weak turbulence theory is considered. Different finite difference schemes are examined to fit the exact analytical solutions of quasilinear equations in hydrodynamic limit (gas-dynamic solution). It is shown that the scheme suggested demonstrates correct asymptotic behavior and can be employed to solve initial value problems for an arbitrary initial electron distribution function.

The self-consistent description of Langmuir wave and ion-sound wave turbulence in the presence of an electron beam is presented for inhomogeneous non-isothermal plasmas. Full numerical solutions of the complete set of kinetic equations for electrons, Langmuir waves, and ion-sound waves are obtained for a inhomogeneous unmagnetized plasma. The result show that the presence of inhomogeneity significantly changes the overall evolution of the system. The inhomogeneity is effective in shifting the wavenumbers of the Langmuir waves, and can thus switch between different process governing the weakly turbulent state. The results can be applied to a variety of plasma conditions, where we choose solar coronal parameters as an illustration, when performing the numerical analysis.

Electromagnetic wave scattering off density inhomogeneities in the solar corona is an important process which determines both the apparent source size and the time profile of radio bursts observed at 1 AU. Here we model the scattering process using a Fokker-Planck equation and apply this formalism to several regimes of interest. In the first regime the density fluctuations are considered quasi-static and diffusion in wavevector space is dominated by angular diffusion on the surface of a constant energy sphere. In the small-angle ("pencil beam") approximation, this diffusion further occurs over a small solid angle in wavevector space. The second regime corresponds to a much later time, by which scattering has rendered the photon distribution near-isotropic resulting in a spatial diffusion of the radiation. The third regime involves time-dependent fluctuations and, therefore, Fermi acceleration of photons. Combined, these results provide a comprehensive theoretical framework within which to understand several important features of propagation of radio burst waves in the solar corona: emitted photons are accelerated in a relatively small inner region and then diffuse outwards to larger distances. En route, angular diffusion results both in source sizes which are substantially larger than the intrinsic source, and in observed intensity-versus-time profiles that are asymmetric, with a sharp rise and an exponential decay. Both of these features are consistent with observations of solar radio bursts.

We present observations of electron energization in magnetic reconnection outflows during the pre-impulsive phase of solar flare SOL2012-07-19T05:58. During a time-interval of about 20 minutes, starting 40 minutes before the onset of the impulsive phase, two X-ray sources were observed in the corona, one above the presumed reconnection region and one below. For both of these sources, the mean electron distribution function as a function of time is determined over an energy range from 0.1~keV up to several tens of keV, for the first time. This is done by simultaneous forward fitting of X-ray and EUV data. Imaging spectroscopy with RHESSI provides information on the high-energy tail of the electron distribution in these sources while EUV images from SDO/AIA are used to constrain the low specific electron energies. The measured electron distribution spectrum in the magnetic reconnection outflows is consistent with a time-evolving kappa-distribution with $\kappa =3.5-5.5$. The spectral evolution suggests that electrons are accelerated to progressively higher energies in the source above the reconnection region, while in the source below, the spectral shape does not change but an overall increase of the emission measure is observed, suggesting density increase due to evaporation. The main mechanisms by which energy is transported away from the source regions are conduction and free-streaming electrons. The latter dominates by more than one order of magnitude and is comparable to typical non-thermal energies during the hard X-ray peak of solar flares, suggesting efficient acceleration even during this early phase of the event.

Solar flare hard X-ray spectroscopy serves as a key diagnostic of the accelerated electron spectrum. However, the standard approach using the collisional cold thick-target model poorly constrains the lower-energy part of the accelerated electron spectrum, and hence the overall energetics of the accelerated electrons are typically constrained only to within one or two orders of magnitude. Here we develop and apply a physically self-consistent warm-target approach which involves the use of both hard X-ray spectroscopy and imaging data. The approach allows an accurate determination of the electron distribution low-energy cutoff, and hence the electron acceleration rate and the contribution of accelerated electrons to the total energy released, by constraining the coronal plasma parameters. Using a solar flare observed in X-rays by the {\em RHESSI} spacecraft, we demonstrate that using the standard cold-target methodology, the low-energy cutoff (and hence the energy content in electrons) is essentially undetermined. However, the warm-target methodology can determine the low-energy electron cutoff with $\sim$7\% uncertainty at the $3\sigma$ level and hence permits an accurate quantitative study of the importance of accelerated electrons in solar flare energetics.

Coronal Mass Ejections (CMEs) are believed to be effective in producing shocks in the solar corona and the interplanetary space. One of the important signatures of shocks and shock acceleration are Type II solar radio bursts that drift with the shock speed and produce bands of fundamental and higher harmonic plasma radio emission. An intriguing aspect of Type II radio bursts is the occasional split of a harmonic band into thinner lanes, known as band-splitting. Here, we report a detailed imaging and spectroscopic observation of a CME-driven shock producing band-splitting in a Type II burst. Using the Low Frequency Array (LOFAR), we examine the spatial and temporal relation of the Type II burst to the associated CME event, use source imaging to calculate the apparent coronal density, and demonstrate how source imaging can be used to estimate projection effects. We consider two widely accepted band-splitting models that make opposing predictions regarding the locations of the true emission sources with respect to the shock front. Our observations suggest that the locations of the upper and lower sub-band sources are spatially separated by $\sim 0.2 \pm 0.05 \, \mathrm{R_\odot}$. However, we quantitatively show, for the first time, that such separation is consistent with radio-wave scattering of plasma radio emission from a single region, implying that the split-band Type II sources could originate from nearly co-spatial locations. Considering the effects of scattering, the observations provide supporting evidence for the model that interprets the band-splitting as emission originating in the upstream and downstream regions of the shock front, two virtually co-spatial areas.

The Square Kilometre Array (SKA) will be the largest radio telescope ever built, aiming to provide collecting area larger than 1 km$^2$. The SKA will have two independent instruments, SKA-LOW comprising of dipoles organized as aperture arrays in Australia and SKA-MID comprising of dishes in South Africa. Currently the phase-1 of SKA, referred to as SKA1, is in its late design stage and construction is expected to start in 2020. Both SKA1-LOW (frequency range of 50-350 MHz) and SKA1-MID Bands 1, 2, and 5 (frequency ranges of 350-1050, 950-1760, and 4600-15300 MHz, respectively) are important for solar observations. In this paper we present SKA's unique capabilities in terms of spatial, spectral, and temporal resolution, as well as sensitivity and show that they have the potential to provide major new insights in solar physics topics of capital importance including (i) the structure and evolution of the solar corona, (ii) coronal heating, (iii) solar flare dynamics including particle acceleration and transport, (iv) the dynamics and structure of coronal mass ejections, and (v) the solar aspects of space weather. Observations of the Sun jointly with the new generation of ground-based and space-borne instruments promise unprecedented discoveries.

Observations of solar flares at sub-THz frequencies (mm and sub-mm wavelengths) over the last two decades often show a spectral component rising with frequency. Unlike a typical gyrosynchrotron spectrum decreasing with frequency, or a weak thermal component from hot coronal plasma, the observations can demonstrate a high flux level (up to ~10^4 s.f.u. at 0.4 THz) and fast variability on sub-second time scales. Although, many models has been put forward to explain the puzzling observations, none of them have clear observational support. Here we propose a scenario to explain the intriguing sub-THz observations. We show that the model, based on free-free emission from the plasma of flare ribbons at temperatures 10^4-10^6 K, is consistent with all existing observations of frequency-rising sub-THz flare emission. The model provides a temperature diagnostic of the flaring chromosphere and suggests fast heating and cooling of the dense transition region plasma.

A component of space weather, electron beams are routinely accelerated in the solar atmosphere and propagate through interplanetary space. Electron beams interact with Langmuir waves resulting in type III radio bursts. Electron beams expand along the trajectory, and using kinetic simulations, we explore the expansion as the electrons propagate away from the Sun. Specifically, we investigate the front, peak and back of the electron beam in space from derived radio brightness temperatures of fundamental type III emission. The front of the electron beams travelled at speeds from 0.2c--0.7c, significantly faster than the back of the beam that travelled between 0.12c--0.35c. The difference in speed between the front and the back elongates the electron beams in time. The rate of beam elongation has a 0.98 correlation coefficient with the peak velocity; in-line with predictions from type III observations. The inferred speeds of electron beams initially increase close to the acceleration region and then decrease through the solar corona. Larger starting densities and harder initial spectral indices result in longer and faster type III sources. Faster electron beams have higher beam energy densities, produce type IIIs with higher peak brightness temperatures and shorter FWHM durations. Higher background plasma temperatures also increase speeds, particularly at the back of the beam. We show how our predictions of electron beam evolution influences type III bandwidth and drift-rates. Our radial predictions of electron beam speed and expansion can be tested by the upcoming in situ electron beam measurements made by Solar Orbiter and Parker Solar Probe.

A considerable fraction of the energy in a solar flare is released as suprathermal electrons; such electrons play a major role in energy deposition in the ambient atmosphere and hence the atmospheric response to flare heating. Historically the transport of these particles has been approximated through a deterministic approach in which first-order secular energy loss to electrons in the ambient target is treated as the dominant effect, with second-order diffusive terms (in both energy and angle) being generally either treated as a small correction or neglected. However, it has recently been pointed out that while neglect of diffusion in energy may indeed be negligible, diffusion in angle is of the same order as deterministic scattering and hence must be included. Here we therefore investigate the effect of angular scattering on the energy deposition profile in the flaring atmosphere. A relatively simple compact expression for the spatial distribution of energy deposition into the ambient plasma is presented and compared with the corresponding deterministic result. For unidirectional injection there is a significant shift in heating from the lower corona to the upper corona; this shift is much smaller for isotropic injection. We also compare the heating profiles due to return current Ohmic heating in the diffusional and deterministic models.

The phenomena observed at the Sun have a variety of unique radio signatures that can be used to diagnose the processes in the solar atmosphere. The insights provided by radio obervations are further enhanced when they are combined with observations from space-based telescopes. This special issue demonstrates the power of combination methodology at work and provides new results on i) type I solar radio bursts and thermal emission to study active regions; ii) type II and IV bursts to better understand the structure of coronal mass ejections; iii)~non-thermal gyro-synchrotron and/or type III bursts to improve characterization of particle acceleration in solar flares. The ongoing improvements in time, frequency, and spatial resolutions of ground-based telescopes reveal new levels of solar phenomena complexity and pose new questions.

Solar radio emission features a large number of fine structures demonstrating great variability in frequency and time. We present spatially resolved spectral radio observations of type IIIb bursts in the $30-80$ MHz range made by the Low Frequency Array (LOFAR). The bursts show well-defined fine frequency structuring called "stria" bursts. The spatial characteristics of the stria sources are determined by the propagation effects of radio waves; their movement and expansion speeds are in the range of 0.1-0.6c. Analysis of the dynamic spectra reveals that both the spectral bandwidth and the frequency drift rate of the striae increase with an increase of their central frequency; the striae bandwidths are in the range of ~20-100 kHz and the striae drift rates vary from zero to ~0.3 MHz s^-1. The observed spectral characteristics of the stria bursts are consistent with the model involving modulation of the type III burst emission mechanism by small-amplitude fluctuations of the plasma density along the electron beam path. We estimate that the relative amplitude of the density fluctuations is of dn/n~10^-3, their characteristic length scale is less than 1000 km, and the characteristic propagation speed is in the range of 400-800 km/s. These parameters indicate that the observed fine spectral structures could be produced by propagating magnetohydrodynamic waves.

Observational detection of quasi-periodic drifting fine structures in a type III radio burst associated with a solar flare SOL2015-04-16T11:22, with Low Frequency Array, is presented. Although similar modulations of the type III emission have been observed before and were associated with the plasma density fluctuations, the origin of those fluctuations was unknown. Analysis of the striae of the intensity variation in the dynamic spectrum allowed us to reveal two quasi-oscillatory components. The shorter component has the apparent wavelength of $\sim2$ Mm, phase speed of $\sim657$ km s$^{-1}$, which gives the oscillation period of $\sim3$ s, and the relative amplitude of $\sim0.35$%. The longer component has the wavelength of $\sim12$ Mm, and relative amplitude of $\sim5.1$%. The short frequency range of the detection does not allow us to estimate its phase speed. However, the properties of the shorter oscillatory component allowed us to interpret it as a fast magnetoacoustic wave guided by a plasma non-uniformity along the magnetic field outwards from the Sun. The assumption that the intensity of the radio emission is proportional to the amount of plasma in the emitting volume allowed us to show that the superposition of the plasma density modulation by a fast wave and a longer-wavelength oscillation of an unspecified nature could readily reproduce the fine structure of the observed dynamic spectrum. The observed parameters of the fast wave give the absolute value of the magnetic field in the emitting plasma of $\sim1.1$ G which is consistent with the radial magnetic field model.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

Solar radio observations provide a unique diagnostic of the outer solar atmosphere. However, the inhomogeneous turbulent corona strongly affects the propagation of the emitted radio waves, so decoupling the intrinsic properties of the emitting source from the effects of radio-wave propagation has long been a major challenge in solar physics. Here we report quantitative spatial and frequency characterization of solar radio burst fine structures observed with the LOw Frequency Array (LOFAR), an instrument with high time resolution that also permits imaging at scales much shorter than those corresponding to radio-wave propagation in the corona. The observations demonstrate that radio-wave propagation effects, and not the properties of the intrinsic emission source, dominate the observed spatial characteristics of radio burst images. These results permit more accurate estimates of source brightness temperatures, and open opportunities for quantitative study of the mechanisms that create the turbulent coronal medium through which the emitted radiation propagates.

Solar type III radio bursts contain a wealth of information about the dynamics of electron beams in the solar corona and the inner heliosphere; currently unobtainable through other means. However, the motion of different regions of an electron beam (front, middle and back) have never been systematically analysed before. We characterise the type III burst frequency-time evolution using the enhanced resolution of LOFAR in the frequency range 30 to 70 MHz and use this to probe electron beam dynamics. Methods. The rise, peak and decay times with a 0.2 MHz spectral resolution were defined for a collection of 31 type III bursts. The frequency evolution is used to ascertain the apparent velocities of the front, middle and back of the type III sources and the trends are interpreted using theoretical and numerical treatments. The type III time profile was better approximated by an asymmetric Gaussian profile, not an exponential as previously used. Rise and decay times increased with decreasing frequency and showed a strong correlation. Durations were smaller than previously observed. Drift rates from the rise times were faster than from the decay times, corresponding to inferred mean electron beam speeds for the front, middle and back of 0.2, 0.17, 0, 15 c, respectively. Faster beam speeds correlate with smaller type III durations. We also find type III frequency bandwidth decreases as frequency decreases. The different speeds naturally explain the elongation of an electron beam in space as it propagates through the heliosphere. The rate of expansion is proportional to the mean speed of the exciter; faster beams expand faster. Beam speeds are attributed to varying ensembles of electron energies at the front, middle and back of the beam.

Solar radio type III bursts are believed to be the most sensitive signature of near-relativistic electron beam propagation in the corona. A solar radio type IIIb-III pair burst with fine frequency structures, observed by the Low Frequency Array (LOFAR) with high temporal ($\sim10$ ms) and spectral (12.5 kHz) resolutions at 30 - 80 MHz, is presented. The observations show that the type III burst consists of many striae, which have a frequency scale of about 0.1 MHz in both the fundamental (plasma) and the harmonic (double plasma) emission. We investigate the effects of background density fluctuations based on the observation of striae structure to estimate the density perturbation in solar corona. It is found that the spectral index of the density fluctuation spectrum is about $-1.7$, and the characteristic spatial scale of the density perturbation is around $700$ km. This spectral index is very close to a Kolmogorov turbulence spectral index of $-5/3$, consistent with a turbulent cascade. This fact indicates that the coronal turbulence may play the important role of modulating the time structures of solar radio type III bursts, and the fine structure of radio type III bursts could provide a useful and unique tool to diagnose the turbulence in the solar corona.

The imaging spectroscopy capabilities of the Reuven Ramaty high energy solar spectroscopic imager (RHESSI) enable the examination of the accelerated electron distribution throughout a solar flare region. In particular, it has been revealed that the energisation of these particles takes place over a region of finite size, sometimes resolved by RHESSI observations. In this paper, we present, for the first time, a spatially distributed acceleration model and investigate the role of inhomogeneous acceleration on the observed X-ray emission properties. We have modelled transport explicitly examining scatter-free and diffusive transport within the acceleration region and compare with the analytic leaky-box solution. The results show the importance of including this spatial variation when modelling electron acceleration in solar flares. The presence of an inhomogeneous, extended acceleration region produces a spectral index that is, in most cases, different from the simple leaky-box prediction. In particular, it results in a generally softer spectral index than predicted by the leaky-box solution, for both scatter-free and diffusive transport, and thus should be taken into account when modelling stochastic acceleration in solar flares.

Using the "enthalpy-based thermal evolution of loops" (EBTEL) model, we investigate the hydrodynamics of the plasma in a flaring coronal loop in which heat conduction is limited by turbulent scattering of the electrons that transport the thermal heat flux. The EBTEL equations are solved analytically in each of the two (conduction-dominated and radiation-dominated) cooling phases. Comparison of the results with typical observed cooling times in solar flares shows that the turbulent mean free-path $\lambda_T$ lies in a range corresponding to a regime in which classical (collision-dominated) conduction plays at most a limited role. We also consider the magnitude and duration of the heat input that is necessary to account for the enhanced values of temperature and density at the beginning of the cooling phase and for the observed cooling times. We find through numerical modeling that in order to produce a peak temperature $\simeq 1.5 \times 10^7$~K and a 200~s cooling time consistent with observations, the flare heating profile must extend over a significant period of time; in particular, its lingering role must be taken into consideration in any description of the cooling phase. Comparison with observationally-inferred values of post-flare loop temperatures, densities, and cooling times thus leads to useful constraints on both the magnitude and duration of the magnetic energy release in the loop, as well as on the value of the turbulent mean free-path $\lambda_T$.

Imaging spectroscopy in X-rays with RHESSI provide the possibility to investigate the spatial evolution of the X-ray emitting electron distribution and therefore to study the transport effects on energetic electrons during solar flares. We study the energy dependence of the energetic electron scattering mean free path in the solar corona. We use the imaging spectroscopy technique with RHESSI to study the evolution of energetic electrons distribution in different part of the magnetic loop during the 2004 May 21 flare. These observations are compared with the radio observations of the gyrosynchrotron radiation of the same flare by Kuznetsov and Kontar (2015), and with the predictions of the diffusive transport model described by Kontar et al. (2014). The X-ray analysis shows a trapping of energetic electrons in the corona and a spectral hardening of the energetic electron distribution between the top of the loop and the footpoints. Coronal trapping of electrons is stronger for the radio-emitting electrons than for the X-ray-emitting electrons. These observations can be explained by the diffusive transport model derived by Kontar et al. (2014). We show that the combination of X-ray and radio diagnostics is a powerful tool to study electron transport in the solar corona in different energy domains. We show that the diffusive transport model can explain our observations; and in the range 25-500 keV, the electron scattering mean free path decreases with electron energy. We can estimate for the first time the scattering mean free path dependence on energy in the corona.

Radio U-bursts and J-bursts are signatures of electron beams propagating along magnetic loops confined to the corona. The more commonly observed type III radio bursts are signatures of electron beams propagating along magnetic loops that extend into interplanetary space. Given the prevalence of solar magnetic flux to be closed in the corona, it is an outstanding question why type III bursts are more frequently observed than U-bursts or J-bursts. We use LOFAR imaging spectroscopy between 30-80 MHz of low-frequency U-bursts and J-bursts, for the first time, to understand why electron beams travelling along coronal loops produce radio emission less often. The different radio source positions were used to model the spatial structure of the guiding magnetic flux tube and then deduce the energy range of the exciting electron beams without the assumption of a standard density model. The radio sources infer a magnetic loop 1 solar radius in altitude, with the highest frequency sources starting around 0.6 solar radii. Electron velocities were found between 0.13 c and 0.24 c, with the front of the electron beam travelling faster than the back of the electron beam. The velocities correspond to energy ranges within the beam from 0.7-11 keV to 0.7-43 keV. The density along the loop is higher than typical coronal density models and the density gradient is smaller. We found that a more restrictive range of accelerated beam and background plasma parameters can result in U-bursts or J-bursts, causing type III bursts to be more frequently observed. The large instability distances required before Langmuir waves are produced by some electron beams, and the small magnitude of the background density gradients make closed loops less facilitating for radio emission than loops that extend into interplanetary space.

During the next solar maximum, two upcoming space-borne X-ray missions, STIX on board Solar Orbiter and MiSolFA, will perform stereoscopic X-ray observations of solar flares at two different locations: STIX at 0.28 AU (at perihelion) and up to inclinations of $\sim25^{\circ}$, and MiSolFA in a low-Earth orbit. The combined observations from these cross-calibrated detectors will allow us to infer the electron anisotropy of individual flares confidently for the first time. We simulated both instrumental and physical effects for STIX and MiSolFA including thermal shielding, background and X-ray Compton backscattering (albedo effect) in the solar photosphere. We predict the expected number of observable flares available for stereoscopic measurements during the next solar maximum. We also discuss the range of useful spacecraft observation angles for the challenging case of close-to-isotropic flare anisotropy. The simulated results show that STIX and MiSolFA will be capable of detecting low levels of flare anisotropy, for M1-class or stronger flares, even with a relatively small spacecraft angular separation of 20-30{\deg}. Both instruments will directly measure the flare X-ray anisotropy of about 40 M- and X-class solar flares during the next solar maximum. Near-future stereoscopic observations with Solar Orbiter/STIX and MiSolFA will help distinguishing between competing flare-acceleration mechanisms, and provide essential constraints regarding collisional and non-collisional transport processes occurring in the flaring atmosphere for individual solar flares.

The conversion of magnetic energy into other forms during solar flares is one of the outstanding open problems in solar physics. It is generally accepted that magnetic reconnection plays a crucial role in these conversion processes. To achieve the rapid energy release required in solar flares, an anomalous resistivity, orders of magnitude higher than the Spitzer resistivity, is often used in MHD simulations of reconnection. Spitzer resistivity is based on Coulomb scattering, which becomes negligible at the high energies achieved by accelerated particles. As a result, simulations of particle acceleration in reconnection events are often performed in the absence of any interaction between accelerated particles and any background plasma. This need not be the case for scattering associated with anomalous resistivity caused by turbulence within solar flares, as the higher resistivity implies an elevated scattering rate. We present results of test particle calculations, with and without pitch angle scattering, subject to fields derived from MHD simulations of two-dimensional (2D) X-point reconnection. Scattering rates proportional to the ratio of the anomalous resistivity to the local Spitzer resistivity, as well as at fixed values, are considered. Pitch angle scattering, which is independent of the anomalous resistivity, causes higher maximum energies in comparison to those obtained without scattering. Scattering rates which are dependent on the local anomalous resistivity tend to produce fewer highly energised particles due to weaker scattering in the separatrices, even though scattering in the current sheet may be stronger when compared to resistivity-independent scattering. Strong scattering also causes an increase in the number of particles exiting the computational box in the reconnection outflow region, as opposed to along the separatrices as is the case in the absence of scattering.

We study the influence of kinetic-scale Alfv\'enic turbulence on the generation of plasma radio emission in the solar coronal regions where the plasma/magnetic pressure ratio $\beta $ is smaller than the electron/ion mass ratio $m_{e}/m_{i}$. The present study is motivated by the phenomenon of solar type I radio storms associated with the strong magnetic field of active regions. The measured brightness temperature of the type I storms can be up to $10^{10}$ K for continuum emission, and can exceed $10^{11}$ K for type I bursts. At present, there is no generally accepted theory explaining such high brightness temperatures and some other properties of the type I storms. We propose the model with the imbalanced turbulence of kinetic-scale Alfv\'en waves producing an asymmetric quasilinear plateau on the upward half of the electron velocity distribution. The Landau damping of resonant Langmuir waves is suppressed and their amplitudes grow spontaneously above the thermal level. The estimated saturation level of Langmuir waves is high enough to generate observed type I radio emission at the fundamental plasma frequency. Harmonic emission does not appear in our model because the backward-propagating Langmuir waves undergo a strong Landau damping. Our model predicts $100\%$ polarization in the sense of the ordinary (o-) mode of type I emission.

We study the energetics of nonthermal electrons produced in small acceleration episodes in the solar corona. We carried out an extensive survey spanning 2004--2015 and shortlisted 6 impulsive electron events detected at 1 AU that was not associated with large solar flares(GOES soft X-ray class $>$ C1) or with coronal mass ejections. Each of these events had weak, but detectable hard Xray (HXR) emission near the west limb, and were associated with interplanetary type III bursts. In some respects, these events seem like weak counterparts of "cold/tenuous" flares. The energy carried by the HXR producing electron population was $\approx 10^{23}$ -- $10^{25}$ erg, while that in the corresponding population detected at 1 AU was $\approx 10^{24}$--$10^{25}$ erg. The number of electrons that escape the coronal acceleration site and reach 1 AU constitute 6 % to 148 % of those that precipitate downwards to produce thick target HXR emission.

Magnetic reconnection and particle acceleration due to the kink instability in twisted coronal loops can be a viable scenario for confined solar flares. Detailed investigation of this phenomenon requires reliable methods for observational detection of magnetic twist in solar flares, which may not be possible solely through extreme UV and soft X-ray thermal emission. Polarisation of microwave emission in flaring loops can be used as one of the detection criteria. The aim of this study is to investigate the effect of magnetic twist in flaring coronal loops on the polarisation of gyro-synchrotron microwave (GSMW) emission, and determine whether it could provide a means for magnetic twist detection. We use time-dependent magnetohydrodynamic and test-particle models developed using LARE3D and GCA codes to investigate twisted coronal loops relaxing following the kink-instability. Synthetic GSMW emission maps (I and V Stokes components) are calculated using GX simulator. It is found that flaring twisted coronal loops produce GSMW radiation with a gradient of circular polarisation across the loop. However, these patterns may be visible only for a relatively short period of time due to fast magnetic reconfiguration after the instability. Their visibility also depends on the orientation and position of the loop on solar disk. Typically, it would be difficult to see these characteristic polarisation pattern in a twisted loop seen from the top (close to the centre of the solar disk), but easier in a twisted loop seen from the side (i.e. observed very close to the limb).

The energy released in solar flares derives from a reconfiguration of magnetic fields to a lower energy state, and is manifested in several forms, including bulk kinetic energy of the coronal mass ejection, acceleration of electrons and ions, and enhanced thermal energy that is ultimately radiated away across the electromagnetic spectrum from optical to X-rays. Using an unprecedented set of coordinated observations, from a suite of instruments, we here report on a hitherto largely overlooked energy component -- the kinetic energy associated with small-scale turbulent mass motions. We show that the spatial location of, and timing of the peak in, turbulent kinetic energy together provide persuasive evidence that turbulent energy may play a key role in the transfer of energy in solar flares. Although the kinetic energy of turbulent motions accounts, at any given time, for only $\sim (0.5-1)$\% of the energy released, its relatively rapid ($\sim$$1-10$~s) energization and dissipation causes the associated throughput of energy (i.e., power) to rival that of major components of the released energy in solar flares, and thus presumably in other astrophysical acceleration sites.

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

The transport of the energy contained in suprathermal electrons in solar flares plays a key role in our understanding of many aspects of flare physics, from the spatial distributions of hard X-ray emission and energy deposition in the ambient atmosphere to global energetics. Historically the transport of these particles has been largely treated through a deterministic approach, in which first-order secular energy loss to electrons in the ambient target is treated as the dominant effect, with second-order diffusive terms (in both energy and angle) being generally either treated as a small correction or even neglected. We here critically analyze this approach, and we show that spatial diffusion through pitch-angle scattering necessarily plays a very significant role in the transport of electrons. We further show that a satisfactory treatment of the diffusion process requires consideration of non-local effects, so that the electron flux depends not just on the local gradient of the electron distribution function but on the value of this gradient within an extended region encompassing a significant fraction of a mean free path. Our analysis applies generally to pitch-angle scattering by a variety of mechanisms, from Coulomb collisions to turbulent scattering. We further show that the spatial transport of electrons along the magnetic field of a flaring loop can be modeled rather effectively as a Continuous Time Random Walk with velocity-dependent probability distribution functions of jump sizes and occurrences, both of which can be expressed in terms of the scattering mean free path.

Solar electron beams responsible for type III radio emission generate Langmuir waves as they propagate out from the Sun. The Langmuir waves are observed via in-situ electric field measurements. These Langmuir waves are not smoothly distributed but occur in discrete clumps, commonly attributed to the turbulent nature of the solar wind electron density. Exactly how the density turbulence modulates the Langmuir wave electric fields is understood only qualitatively. Using weak turbulence simulations, we investigate how solar wind density turbulence changes the probability distribution functions, mean value and variance of the beam-driven electric field distributions. Simulations show rather complicated forms of the distribution that are dependent upon how the electric fields are sampled. Generally the higher magnitude of density fluctuations reduce the mean and increase the variance of the distribution in a consistent manor to the predictions from resonance broadening by density fluctuations. We also demonstrate how the properties of the electric field distribution should vary radially from the Sun to the Earth and provide a numerical prediction for the in-situ measurements of the upcoming Solar Orbiter and Solar Probe Plus spacecraft.

Flare-driven coronal rain can manifest from rapidly cooled plasma condensations near coronal loop-tops in thermally unstable post-flare arcades. We detect 5 phases that characterise the post-flare decay: heating, evaporation, conductive cooling dominance for ~120 s, radiative / enthalpy cooling dominance for ~4700 s and finally catastrophic cooling occurring within 35-124 s leading to rain strands with s periodicity of 55-70 s. We find an excellent agreement between the observations and model predictions of the dominant cooling timescales and the onset of catastrophic cooling. At the rain formation site we detect co-moving, multi-thermal rain clumps that undergo catastrophic cooling from ~1 MK to ~22000 K. During catastrophic cooling the plasma cools at a maximum rate of 22700 K s-1 in multiple loop-top sources. We calculated the density of the EUV plasma from the DEM of the multi-thermal source employing regularised inversion. Assuming a pressure balance, we estimate the density of the chromospheric component of rain to be 9.21x10^11 +-1.76x10^11 cm-3 which is comparable with quiescent coronal rain densities. With up to 8 parallel strands in the EUV loop cross section, we calculate the mass loss rate from the post-flare arcade to be as much as 1.98x10^12 +/-4.95x10^11 g s-1. Finally, we reveal a close proximity between the model predictions of 10^5.8 K and the observed properties between 10^5.9 K and 10^6.2 K, that defines the temperature onset of catastrophic cooling. The close correspondence between the observations and numerical models suggests that indeed acoustic waves (with a sound travel time of 68 s) could play an important role in redistributing energy and sustaining the enthalpy-based radiative cooling.

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 $\lambda_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 $\lambda_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.

This study entails the third part of a global flare energetics project, in which Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data of 191 M and X-class flare events from the first 3.5 yrs of the Solar Dynamics Observatory (SDO) mission are analyzed. We fit a thermal and a nonthermal component to RHESSI spectra, yielding the temperature of the differential emission measure (DEM) tail, the nonthermal power law slope and flux, and the thermal/nonthermal cross-over energy $e_{\mathrm{co}}$. From these parameters we calculate the total nonthermal energy $E_{\mathrm{nt}}$ in electrons with two different methods: (i) using the observed cross-over energy $e_{\mathrm{co}}$ as low-energy cutoff, and (ii) using the low-energy cutoff $e_{\mathrm{wt}}$ predicted by the warm thick-target bremsstrahlung model of Kontar et al. {\bf Based on a mean temperature of $T_e=8.6$ MK in active regions we find low-energy cutoff energies of $e_{\mathrm{wt}} =6.2\pm 1.6$ keV for the warm-target model, which is significantly lower than the cross-over energies $e_{\mathrm{co}}=21 \pm 6$ keV. Comparing with the statistics of magnetically dissipated energies $E_{\mathrm{mag}}$ and thermal energies $E_{\mathrm{th}}$ from the two previous studies, we find the following mean (logarithmic) energy ratios with the warm-target model: $E_{\mathrm{nt}} = 0.41 \ E_{\mathrm{mag}}$, $E_{\mathrm{th}} = 0.08 \ E_{\mathrm{mag}}$, and $E_{\mathrm{th}} = 0.15 \ E_{\mathrm{nt}}$. The total dissipated magnetic energy exceeds the thermal energy in 95% and the nonthermal energy in 71% of the flare events, which confirms that magnetic reconnection processes are sufficient to explain flare energies. The nonthermal energy exceeds the thermal energy in 85\% of the events, which largely confirms the warm thick-target model.

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.

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.

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\"ahovi 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.

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.

Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate the thermal structure and flows in flaring atmospheres. This can be used to estimate the partitioning between the thermal and kinetic energies. Our aim is to model large-scale (50-10000km) velocity distributions in order to interpret non-thermal broadening of different spectral EUV lines. The developed models allow us to understand the origin of the observed spectral line shifts and broadening, and link these features to particular physical effects in flaring atmospheres. We use ideal MHD to derive twisted magnetic fluxtube configurations in a stratified atmosphere. The evolution of these fluxtubes is followed using resistive MHD, with anomalous resistivity depending on the local density and temperature. The model also takes into account the thermal conduction and radiative losses. The model allows us to evaluate average velocities and velocity dispersions, which would be interpreted as `non-thermal' velocities in observations, at different temperatures for different parts of the models. Our models show qualitative and quantitative agreement with observations. The line-of-sight (LOS) velocity dispersions demonstrate substantial correlation with the temperature, increasing from about 20-30 km/s around 1 MK to about 200-400 km/s near 10-20 MK. The average LOS velocities also correlate with velocity dispersions, although they demonstrate a very strong scattering, compared to observations. We also note that near foot-points the velocity dispersions across the magnetic field are systematically lower than those along the field. We conclude, that the correlation between the flow velocities, velocity dispersions and temperatures are likely to indicate that the same heating mechanism is responsible for heating the plasma, its turbulisation and expansion/evaporation.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

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.

To solve a number of problems in solar physics related to mechanisms of energy release in solar corona parameters of hot coronal plasma are required, such as energy distribution, emission measure, differential emission measure, and their evolution with time. Of special interest is the distribution of solar plasma by energies, which can evolve from a nearly Maxwellian distribution to a distribution with a more complex structure during a solar flare. The exact form of this distribution for low-energy particles, which receive the bulk of flare energy, is still poorly known; therefore, detailed investigations are required. We present a developed method of simultaneous fitting of data from two spacecrafts Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) and Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), using a differential emission measure and a thin target model for the August 14, 2010 flare event.

The detailed knowledge of plasma heating and acceleration region properties presents a major observational challenge in solar flare physics. Using the Ramaty High Energy Solar Spectroscopic Imager (RHESSI), the high temperature differential emission measure, DEM(T), and the energy-dependent spatial structure of solar flare coronal sources are studied quantitatively. The altitude of the coronal X-ray source is observed to increase with energy by ~+0.2 arcsec/keV between 10 and 25 keV. Although an isothermal model can fit the thermal X-ray spectrum observed by RHESSI, such a model cannot account for the changes in altitude, and multi-thermal coronal sources are required where the temperature increases with altitude. For the first time, we show how RHESSI imaging information can be used to constrain the DEM(T) of a flaring plasma. We develop a thermal bremsstrahlung X-ray emission model with inhomogeneous temperature and density distributions to simultaneously reproduce: i) DEM(T), ii) altitude as a function of energy, and iii) vertical extent of the flaring coronal source versus energy. We find that the temperature-altitude gradient in the region is ~+0.08 keV/arcsec (~1.3 MK/Mm). Similar altitude-energy trends in other flares suggest that the majority of coronal X-ray sources are multi-thermal and have strong vertical temperature and density gradients with a broad DEM(T).

EUV (Extreme-Ultraviolet) waves are globally propagating disturbances that have been observed since the era of the SoHO/EIT instrument. Although the kinematics of the wave front and secondary wave components have been widely studied, there is not much known about the generation and plasma properties of the wave. In this paper we discuss the effect of an EUV wave on the local plasma as it passes through the corona. We studied the EUV wave, generated during the 2011 February 15 X-class flare/CME event, using Differential Emission Measure diagnostics. We analyzed regions on the path of the EUV wave and investigated the local density and temperature changes. From our study we have quantitatively confirmed previous results that during wave passage the plasma visible in the Atmospheric Imaging Assembly (AIA) 171A channel is getting heated to higher temperatures corresponding to AIA 193A and 211A channels. We have calculated an increase of 6 - 9% in density and 5 - 6% in temperature during the passage of the EUV wave. We have compared the variation in temperature with the adiabatic relationship and have quantitatively demonstrated the phenomenon of heating due to adiabatic compression at the wave front. However, the cooling phase does not follow adiabatic relaxation but shows slow decay indicating slow energy release being triggered by the wave passage. We have also identified that heating is taking place at the front of the wave pulse rather than at the rear. Our results provide support for the case that the event under study here is a compressive fast-mode wave or a shock.

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.

X-ray observations are a direct diagnostic of fast electrons produced in solar flares, energized during the energy release process and directed towards the Sun. Since the properties of accelerated electrons can be substantially changed during their transport and interaction with the background plasma, a model must ultimately be applied to X-ray observations in order to understand the mechanism responsible for their acceleration. A cold thick target model is ubiquitously used for this task, since it provides a simple analytic relationship between the accelerated electron spectrum and the emitting electron spectrum in the X-ray source, with the latter quantity readily obtained from X-ray observations. However, such a model is inappropriate for the majority of solar flares in which the electrons propagate in a hot megaKelvin plasma, because it does not take into account the physics of thermalization of fast electrons. The use of a more realistic model, properly accounting for the properties of the background plasma, and the collisional diffusion and thermalization of electrons, can alleviate or even remove many of the traditional problems associated with the cold thick target model and the deduction of the accelerated electron spectrum from X-ray spectroscopy, such as the number problem and the need to impose an ad hoc low energy cut-off.

Extending previous studies of nonthermal electron transport in solar flares which include the effects of collisional energy diffusion and thermalization of fast electrons, we present an analytic method to infer more accurate estimates of the accelerated electron spectrum in solar flares from observations of the hard X-ray spectrum. Unlike for the standard cold-target model, the spatial characteristics of the flaring region, especially the necessity to consider a finite volume of hot plasma in the source, need to be taken into account in order to correctly obtain the injected electron spectrum from the source-integrated electron flux spectrum (a quantity straightforwardly obtained from hard X-ray observations). We show that the effect of electron thermalization can be significant enough to nullify the need to introduce an {\it ad hoc} low-energy cutoff to the injected electron spectrum in order to keep the injected power in non-thermal electrons at a reasonable value. Rather the suppression of the inferred low-energy end of the injected spectrum compared to that deduced from a cold-target analysis allows the inference from hard X-ray observations of a more realistic energy in injected non-thermal electrons in solar flares.

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.

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

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.

We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.

Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray sources. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and the extended coronal source, whose intensity strongly depends on the coronal loop density.

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.

Non-thermal electrons accelerated in the solar corona can produce intense coherent radio emission, known as solar type III radio bursts. This intense radio emission is often observed from hundreds of MHz in the corona down to the tens of kHz range in interplanetary space. It involves a chain of physical processes from the generation of Langmuir waves to nonlinear processes of wave-wave interaction. We develop a self-consistent model to calculate radio emission from a non-thermal electron population over large frequency range, including the effects of electron transport, Langmuir wave-electron interaction, the evolution of Langmuir waves due to non-linear wave-wave interactions, Langmuir wave conversion into electromagnetic emission, and finally escape of the electromagnetic waves. For the first time we simulate escaping radio emission over a broad frequency range from 500~MHz down to a few MHz and infer key properties of the radio emission observed: the onset (starting) frequency, {identification as fundamental or harmonic emission}, peak flux density, instantaneous frequency bandwidth, and timescales for rise and decay. Comparing with the observations, these large scale simulations enable us to identify the processes governing the key type III solar radio burst characteristics.

We have performed a statistical study of $152$ Type III radio bursts observed by Solar TErrestrial RElations Observatory (STEREO)/Waves between May 2007 and February 2013. We have investigated the flux density between $125$kHz and $16$MHz. Both high- and low-frequency cutoffs have been observed in $60\,%$ of events suggesting an important role of propagation. As already reported by previous authors, we observed that the maximum flux density occurs at $1$MHz on both spacecraft. We have developed a simplified analytical model of the flux density as a function of radial distance and compared it to the STEREO/Waves data.

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 $\tau_{\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.

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.

We investigate in detail the 21 May 2004 flare using simultaneous observations of the Nobeyama Radioheliograph, Nobeyama Radiopolarimeters, Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and Solar and Heliospheric Observatory (SOHO). The flare images in different spectral ranges reveal the presence of a well-defined single flaring loop in this event. We have simulated the gyrosynchrotron microwave emission using the recently developed interactive IDL tool GX Simulator. By comparing the simulation results with the observations, we have deduced the spatial and spectral properties of the non-thermal electron distribution. The microwave emission has been found to be produced by the high-energy electrons ($>100$ keV) with a relatively hard spectrum ($\delta\simeq 2$); the electrons were strongly concentrated near the loop top. At the same time, the number of high-energy electrons near the footpoints was too low to be detected in the RHESSI images and spatially unresolved data. The SOHO Extreme-ultraviolet Imaging Telescope images and the low-frequency microwave spectra suggest the presence of an extended "envelope" of the loop with lower magnetic field. Most likely, the energetic electron distribution in the considered flare reflects the localized (near the loop top) particle acceleration (injection) process accompanied by trapping and scattering.

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.

Recent observations from {\em RHESSI} have revealed that the number of non-thermal electrons in the coronal part of a flaring loop can exceed the number of electrons required to explain the hard X-ray-emitting footpoints of the same flaring loop. Such sources cannot, therefore, be interpreted on the basis of the standard collisional transport model, in which electrons stream along the loop while losing their energy through collisions with the ambient plasma; additional physical processes, to either trap or scatter the energetic electrons, are required. Motivated by this and other observations that suggest that high energy electrons are confined to the coronal region of the source, we consider turbulent pitch angle scattering of fast electrons off low frequency magnetic fluctuations as a confinement mechanism, modeled as a spatial diffusion parallel to the mean magnetic field. In general, turbulent scattering leads to a reduction of the collisional stopping distance of non-thermal electrons along the loop and hence to an enhancement of the coronal HXR source relative to the footpoints. The variation of source size $L$ with electron energy $E$ becomes weaker than the quadratic behavior pertinent to collisional transport, with the slope of $L(E)$ depending directly on the mean free path $\lambda$ again pitch angle scattering. Comparing the predictions of the model with observations, we find that $\lambda \sim$$(10^8-10^9)$ cm for $\sim30$ keV, less than the length of a typical flaring loop and smaller than, or comparable to, the size of the electron acceleration region.

The evolution of a solar flare accelerated non-thermal electron population and associated plasma emission is considered in collisional inhomogeneous plasma. Non-thermal electrons collisionally evolve to become unstable and generate Langmuir waves, which may lead to intense radio emission. We self-consistently simulated the collisional relaxation of electrons, wave-particle interactions, and non-linear Langmuir wave evolution in plasma with density fluctuations. Additionally, we simulated the scattering, decay, and coalescence of the Langmuir waves which produce radio emission at the fundamental or the harmonic of the plasma frequency, using an angle-averaged emission model. Long-wavelength density fluctuations, such as are observed in the corona, are seen to strongly suppress the levels of radio emission, meaning that a high level of Langmuir waves can be present without visible radio emission. Additionally, in homogeneous plasma, the emission shows time and frequency variations that could be smoothed out by density inhomogeneities.

The capability of maintaining two satellites in precise relative position, stable in a celestial coordinate system, would enable major advances in a number of scientific disciplines and with a variety of types of instrumentation. The common requirement is for formation flying of two spacecraft with the direction of their vector separation in inertial coordinates precisely controlled and accurately determined as a function of time. We consider here the scientific goals that could be achieved with such technology and review some of the proposals that have been made for specific missions. Types of instrumentation that will benefit from the development of this type of formation flying include 1) imaging systems, in which an optical element on one spacecraft forms a distant image recorded by a detector array on the other spacecraft, including telescopes capable of very high angular resolution; 2) systems in which the front spacecraft of a pair carries an occulting disk, allowing very high dynamic range observations of the solar corona and exoplanets; 3) interferometers, another class of instrument that aims at very high angular resolution and which, though usually requiring more than two spacecraft, demands very much the same developments.

Major solar eruptive events (SEEs), consisting of both a large flare and a near simultaneous large fast coronal mass ejection (CME), are the most powerful explosions and also the most powerful and energetic particle accelerators in the solar system, producing solar energetic particles (SEPs) up to tens of GeV for ions and hundreds of MeV for electrons. The intense fluxes of escaping SEPs are a major hazard for humans in space and for spacecraft. Furthermore, the solar plasma ejected at high speed in the fast CME completely restructures the interplanetary medium (IPM) - major SEEs therefore produce the most extreme space weather in geospace, the interplanetary medium, and at other planets. Thus, understanding the flare/CME energy release process(es) and the related particle acceleration processes are major goals in Heliophysics. To make the next major breakthroughs, we propose a new mission concept, SEE 2020, a single spacecraft with a complement of advanced new instruments that focus directly on the coronal energy release and particle acceleration sites, and provide the detailed diagnostics of the magnetic fields, plasmas, mass motions, and energetic particles required to understand the fundamental physical processes involved.

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.

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.

The acceleration of charged particles in magnetized plasmas is considered during turbulent multi-island magnetic reconnection. The particle acceleration model is constructed for an ensemble of islands which produce adiabatic compression of the particles. The model takes into account the statistical fluctuations in the compression rate experienced by the particles during their transport in the acceleration region. The evolution of the particle distribution function is described as a simultaneous first and second-order Fermi acceleration process. While the efficiency of the first-order process is controlled by the average rate of compression, the second order process involves the variance in the compression rate. Moreover, the acceleration efficiency associated with the second-order process involves both the Eulerian properties of the compression field and the Lagrangian properties of the particles. The stochastic contribution to the acceleration is non-resonant and can dominate the systematic part in the case of a large variance in the compression rate. The model addresses the role of the second-order process, how the latter can be related to the large-scale turbulent transport of particles and explains some features of the numerical simulations of particle acceleration by multi-island contraction during magnetic reconnection.

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.

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.

The aim of this work is to determine the multi-thermal characteristics and plasma energetics of an eruptive plasmoid and occulted flare observed by Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). We study an event from 03-Nov-2010 (peaking at 12:20UT in GOES soft X-rays) of a coronal mass ejection and occulted flare which demonstrates the morphology of a classic erupting flux rope. The high spatial, and time resolution, and six coronal channels, of the SDO/AIA images allows the dynamics of the multi-thermal emission during the initial phases of eruption to be studied in detail. The Differential Emission Measure (DEM) is calculated, using an optimised version of a regularized inversion method (Hannah & Kontar 2012), for each pixel across the six channels at different times, resulting in emission measure maps and movies in a variety of temperature ranges. We find that the core of the erupting plasmoid is hot (8-11, 11-14MK) with a similarly hot filamentary "stem" structure connecting it to the lower atmosphere, which could be interpreted as the current sheet in the flux rope model, though is wider than these models suggest. The velocity of the leading edge of the eruption is 597-664 km s$^{-1}$ in the temperature range $\ge$3-4MK and between 1029-1246 km s$^{-1}$ for $\le$2-3MK. We estimate the density (in 11-14 MK) of the erupting core and stem during the impulsive phase to be about $3\times10^9$ cm$^{-3}$, $6\times10^9$ cm$^{-3}$, $9\times10^8$ cm$^{-3}$ in the plasmoid core, stem and surrounding envelope of material. This gives thermal energy estimates of $5\times10^{29}$ erg, $1\times10^{29}$ erg and $2\times10^{30}$ erg. The kinetic energy for the core and envelope is slightly smaller. The thermal energy of the core and current sheet grows during the eruption, suggesting continuous influx of energy presumably via reconnection.

The aim of this paper is to demonstrate the effect of turbulent background density fluctuations on flare-accelerated electron transport in the solar corona. Using the quasi-linear approximation, we numerically simulated the propagation of a beam of accelerated electrons from the solar corona to the chromosphere, including the self-consistent response of the inhomogeneous background plasma in the form of Langmuir waves. We calculated the X-ray spectrum from these simulations using the bremsstrahlung cross-section and fitted the footpoint spectrum using the collisional "thick-target" model, a standard approach adopted in observational studies. We find that the interaction of the Langmuir waves with the background electron density gradient shifts the waves to a higher phase velocity where they then resonate with higher velocity electrons. The consequence is that some of the electrons are shifted to higher energies, producing more high-energy X-rays than expected if the density inhomogeneity is not considered. We find that the level of energy gain is strongly dependent on the initial electron beam density at higher energy and the magnitude of the density gradient in the background plasma. The most significant gains are for steep (soft) spectra that initially had few electrons at higher energies. If the X-ray spectrum of the simulated footpoint emission are fitted with the standard "thick-target" model (as is routinely done with RHESSI observations) some simulation scenarios produce more than an order-of-magnitude overestimate of the number of electrons $>50$keV in the source coronal distribution.

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 \tau_D is of the order of the quasilinear time scale \tau_ql, while when \tau_D << \tau_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.

The angular distribution of electrons accelerated in solar flares is a key parameter in the understanding of the acceleration and propagation mechanisms that occur there. However, the anisotropy of energetic electrons is still a poorly known quantity, with observational studies producing evidence for an isotropic distribution and theoretical models mainly considering the strongly beamed case. We use the effect of photospheric albedo to infer the pitch angle distribution of X-ray emitting electrons using Hard X-ray data from RHESSI. A bi-directional approximation is applied and a regularized inversion is performed for eight large flare events to deduce the electron spectra in both downward (towards the photosphere) and upward (away from the photosphere) directions. The electron spectra and the electron anisotropy ratios are calculated for broad energy range from about 10 and up to ~ 300 keV near the peak of the flares. The variation of electron anisotropy over short periods of time intervals lasting 4, 8 and 16 seconds near the impulsive peak has been examined. The results show little evidence for strong anisotropy and the mean electron flux spectra are consistent with the isotropic electron distribution. The 3-sigma level uncertainties, although energy and event dependent, are found to suggest that anisotropic distribution with anisotropy larger than ~ 3 are not consistent with the hard X-ray data. At energies above 150-200 keV, the uncertainties are larger and thus the possible electron anisotropies could be larger.

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.

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.

We investigate the effects of acceleration during non-linear electron-beam relaxation in magnetized plasma in the case of electron transport in solar flares. The evolution of electron distribution functions is computed using a three-dimensional particle-in-cell electromagnetic code. Analytical estimations under simplified assumptions are made to provide comparisons. We show that, during the non-linear evolution of the beam-plasma system, the accelerated electron population appears. We found that, although the electron beam loses its energy efficiently to the thermal plasma, a noticeable part of the electron population is accelerated. For model cases with initially monoenergetic beams in uniform plasma, we found that the amount of energy in the accelerated electrons above the injected beam-electron energy varies depending the plasma conditions and could be around 10-30% of the initial beam energy. This type of acceleration could be important for the interpretation of non-thermal electron populations in solar flares. Its neglect could lead to the over-estimation of accelerated electron numbers. The results emphasize that collective plasma effects should not be treated simply as an additional energy-loss mechanism, when hard X-ray emission in solar flares is interpreted, notably in the case of RHESSI data.

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 $\omega=\omega(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.

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.

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.

We develop and apply an enhanced regularization algorithm, used in RHESSI X-ray spectral analysis, to constrain the ill-posed inverse problem that is determining the DEM from solar observations. We demonstrate this computationally fast technique applied to a range of DEM models simulating broadband imaging data from SDO/AIA and high resolution line spectra from Hinode/EIS, as well as actual active region observations with Hinode/EIS and XRT. As this regularization method naturally provides both vertical and horizontal (temperature resolution) error bars we are able to test the role of uncertainties in the data and response functions. The regularization method is able to successfully recover the DEM from simulated data of a variety of model DEMs (single Gaussian, multiple Gaussians and CHIANTI DEM models). It is able to do this, at best, to over four orders of magnitude in DEM space but typically over two orders of magnitude from peak emission. The combination of horizontal and vertical error bars and the regularized solution matrix allows us to easily determine the accuracy and robustness of the regularized DEM. We find that the typical range for the horizontal errors is $\Delta$log$T\approx 0.1 -0.5$ and this is dependent on the observed signal to noise, uncertainty in the response functions as well as the source model and temperature. With Hinode/EIS an uncertainty of 20% greatly broadens the regularized DEMs for both Gaussian and CHIANTI models although information about the underlying DEMs is still recoverable. When applied to real active region observations with Hinode/EIS and XRT the regularization method is able to recover a DEM similar to that found via a MCMC method but in considerably less computational time.

Context. High energy electrons accelerated during solar flare are abundant in the solar corona and in the interplanetary space. Commonly, the number and the energy of non-thermal electrons at the Sun is estimated using hard X-ray (HXR) spectral observations (e.g. RHESSI) and a single-particle collisional approximation. Aims. To investigate the role of the spectrally evolving Langmuir turbulence on the population of energetic electrons in the solar corona. Methods. We numerically simulate the relaxation of a power-law non-thermal electron population in a collisional inhomogeneous plasma including wave-particle, and wave-wave interactions. Results. The numerical simulations show that the long-time evolution of electron population above 20 keV deviates substantially from the collisional approximation when wave-particle interactions in non-uniform plasma are taken into account. The evolution of Langmuir wave spectrum towards smaller wavenumbers, due to large-scale density fluctuations and wave-wave interactions, leads to an effective acceleration of electrons. Furthermore, the time-integrated spectrum of non-thermal electrons, which is normally observed with HXR above 20 keV, is noticeably increased due to acceleration of non-thermal electrons by Langmuir waves. Conclusions. The results show that the observed HXR spectrum, when interpreted in terms of collisional relaxation, can lead to an overestimated number and energy of energetic electrons accelerated in the corona.

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.

This paper aims to study the polarization of hard X-ray (HXR) sources in the solar atmosphere, including Compton backscattering of photons in the photosphere (the albedo effect) and the spatial distribution of polarization across the source. HXR photon polarization and spectra produced via electron-ion bremsstrahlung are calculated from electron distributions typical for solar flares. Compton scattering and photoelectric absorption are then modelled using Monte Carlo simulations of photon transport in the photosphere. Polarization maps across HXR sources (primary and albedo components) for each of the modelled electron distributions are calculated at various source locations from the solar centre to the limb. We show that Compton scattering produces a distinct polarization variation across the albedo patch at peak albedo energies of 20-50 keV for all anisotropies modelled. The results show that there are distinct spatial polarization changes in both the radial and perpendicular to radial directions across the extent of the HXR source at a given disk location. In the radial direction, the polarization magnitude and direction at specific positions along the HXR source will either increase or decrease with increased photon distribution directivity towards the photosphere. We also show how high electron cutoff energies influence the direction of polarization at above ~100 keV. Spatially resolved HXR polarization measurements can provide important information about the directivity and energetics of the electron distribution. Our results indicate the preferred angular resolution of polarization measurements required to distinguish between the scattered and primary components. We also show how spatially resolved polarization measurements could be used to probe the emission pattern of an HXR source, using both the magnitude and the direction of polarization.

A fundamental problem in astrophysics is the interaction between magnetic turbulence and charged particles. It is now possible to use \emph{Ramaty High Energy Solar Spectroscopic Imager (RHESSI)} observations of hard X-rays (HXR) emitted by electrons to identify the presence of turbulence and to estimate the magnitude of the magnetic field line diffusion coefficient at least in dense coronal flaring loops.} {We discuss the various possible regimes of cross-field transport of non-thermal electrons resulting from broadband magnetic turbulence in coronal loops. The importance of the Kubo number $K$ as a governing parameter is emphasized and results applicable in both the large and small Kubo number limits are collected.} {Generic models, based on concepts and insights developed in the statistical theory of transport, are applied to the coronal loops and to the interpretation of hard X-ray imaging data in solar flares. The role of trapping effects, which become important in the non-linear regime of transport, is taken into account in the interpretation of the data.} For this flaring solar loop, we constrain the ranges of parallel and perpendicular correlation lengths of turbulent magnetic fields and possible Kubo numbers. We show that a substantial amount of magnetic fluctuations with energy $\sim 1%$ (or more) of the background field can be inferred from the measurements of the magnetic diffusion coefficient inside thick-target coronal loops.

X-radiation from energetic electrons is the prime diagnostic of flare-accelerated electrons. The observed X-ray flux (and polarization state) is fundamentally a convolution of the cross-section for the hard X-ray emission process(es) in question with the electron distribution function, which is in turn a function of energy, direction, spatial location and time. To address the problems of particle propagation and acceleration one needs to infer as much information as possible on this electron distribution function, through a deconvolution of this fundamental relationship. This review presents recent progress toward this goal using spectroscopic, imaging and polarization measurements, primarily from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Previous conclusions regarding the energy, angular (pitch angle) and spatial distributions of energetic electrons in solar flares are critically reviewed. We discuss the role and the observational evidence of several radiation processes: free-free electron-ion, free-free electron-electron, free-bound electron-ion bremsstrahlung, photoelectric absorption and Compton back-scatter (albedo), using both spectroscopic and imaging techniques. This unprecedented quality of data allows for the first time inference of the angular distributions of the X-ray-emitting electrons using albedo, improved model-independent inference of electron energy spectra and emission measures of thermal plasma. Moreover, imaging spectroscopy has revealed hitherto unknown details of solar flare morphology and detailed spectroscopy of coronal, footpoint and extended sources in flaring regions. Additional attempts to measure hard X-ray polarization were not sufficient to put constraints on the degree of anisotropy of electrons, but point to the importance of obtaining good quality polarization data.

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.

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.

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.

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.

Investigate the spatial and spectral evolution of hard X-ray (HXR) emission from flare accelerated electron beams subject to collisional transport and wave-particle interactions in the solar atmosphere. We numerically follow the propagation of a power-law of accelerated electrons in 1D space and time with the response of the background plasma in the form of Langmuir waves using the quasilinear approximation.}{We find that the addition of wave-particle interactions to collisional transport for a transient initially injected electron beam flattens the spectrum of the footpoint source. The coronal source is unchanged and so the difference in the spectral indices between the coronal and footpoint sources is \Delta \gamma > 2, which is larger than expected from purely collisional transport. A steady-state beam shows little difference between the two cases, as has been previously found, as a transiently injected electron beam is required to produce significant wave growth, especially at higher velocities. With this transiently injected beam the wave-particle interactions dominate in the corona whereas the collisional losses dominate in the chromosphere. The shape of the spectrum is different with increasing electron beam density in the wave-particle interaction case whereas with purely collisional transport only the normalisation is changed. We also find that the starting height of the source electron beam above the photosphere affects the spectral index of the footpoint when Langmuir wave growth is included. This may account for the differing spectral indices found between double footpoints if asymmetrical injection has occurred in the flaring loop.

Aims: We seek to illustrate the analysis problems posed by RHESSI spacecraft motion by studying persistent instrumental oscillations found in the lightcurves measured by RHESSI's X-ray detectors in the 6-12 keV and 12-25 keV energy range during the decay phase of the flares of 2004 November 4 and 6. Methods: The various motions of the RHESSI spacecraft which may contribute to the manifestation of oscillations are studied. The response of each detector in turn is also investigated. Results: We find that on 2004 November 6 the observed oscillations correspond to the nutation period of the RHESSI instrument. These oscillations are also of greatest amplitude for detector 5, while in the lightcurves of many other detectors the oscillations are small or undetectable. We also find that the variation in detector pointing is much larger during this flare than the counterexample of 2004 November 4. Conclusions: Sufficiently large nutation motions of the RHESSI spacecraft lead to clearly observable oscillations in count rates, posing a significant hazard for data analysis. This issue is particularly problematic for detector 5 due to its design characteristics. Dynamic correction of the RHESSI counts, accounting for the livetime, data gaps, and the transmission of the bi-grid collimator of each detector, is required to overcome this issue. These corrections should be applied to all future oscillation studies.

Plasma turbulence is thought to be associated with various physical processes involved in solar flares, including magnetic reconnection, particle acceleration and transport. Using Ramaty High Energy Solar Spectroscopic Imager ({\it RHESSI}) observations and the X-ray visibility analysis, we determine the spatial and spectral distributions of energetic electrons for a flare (GOES M3.7 class, April 14, 2002 23$:$55 UT), which was previously found to be consistent with a reconnection scenario. It is demonstrated that because of the high density plasma in the loop, electrons have to be continuously accelerated about the loop apex of length $\sim 2\times 10^9$cm and width $\sim 7\times 10^8$cm. Energy dependent transport of tens of keV electrons is observed to occur both along and across the guiding magnetic field of the loop. We show that the cross-field transport is consistent with the presence of magnetic turbulence in the loop, where electrons are accelerated, and estimate the magnitude of the field line diffusion coefficient for different phases of the flare. The energy density of magnetic fluctuations is calculated for given magnetic field correlation lengths and is larger than the energy density of the non-thermal electrons. The level of magnetic fluctuations peaks when the largest number of electrons is accelerated and is below detectability or absent at the decay phase. These hard X-ray observations provide the first observational evidence that magnetic turbulence governs the evolution of energetic electrons in a dense flaring loop and is suggestive of their turbulent acceleration.

We investigate the type III radio bursts and X-ray signatures of accelerated electrons in a well observed solar flare in order to find the spatial properties of the acceleration region. Combining simultaneous RHESSI hard X-ray flare data and radio data from Phoenix-2 and the Nan\c{c}ay radioheliograph, the outward transport of flare accelerated electrons is analyzed. The observations show that the starting frequencies of type III bursts are anti-correlated with the HXR spectral index of solar flare accelerated electrons. We demonstrate both analytically and numerically that the type III burst starting location is dependent upon the accelerated electron spectral index and the spatial acceleration region size, but weakly dependent on the density of energetic electrons for relatively intense electron beams. Using this relationship and the observed anti-correlation, we estimate the height and vertical extent of the acceleration region, giving values of around 50 Mm and 10 Mm respectively. The inferred acceleration height and size suggests that electrons are accelerated well above the soft X-ray loop-top, which could be consistent with the electron acceleration between 40 Mm and 60 Mm above the flaring loop.

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.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

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.

Solar flare accelerated electron beams propagating away from the Sun can interact with the turbulent interplanetary media, producing plasma waves and type III radio emission. These electron beams are detected near the Earth with a double power-law energy spectrum. We simulate electron beam propagation from the Sun to the Earth in the weak turbulent regime taking into account the self-consistent generation of plasma waves and subsequent wave interaction with density fluctuations from low frequency MHD turbulence. The rate at which plasma waves are induced by an unstable electron beam is reduced by background density fluctuations, most acutely when fluctuations have large amplitudes or small wavelengths. This suppression of plasma waves alters the wave distribution which changes the electron beam transport. Assuming a 5/3 Kolmogorov-type power density spectrum of fluctuations often observed near the Earth, we investigate the corresponding energy spectrum of the electron beam after it has propagated 1 AU. We find a direct correlation between the spectrum of the double power-law below the break energy and the turbulent intensity of the background plasma. For an initial spectral index of 3.5, we find a range of spectra below the break energy between 1.6-2.1, with higher levels of turbulence corresponding to higher spectral indices.

{This work aims to investigate the spectral structure of the parallel electric field generated by strong anisotropic and balanced Alfvenic turbulence in relation with the problem of electron acceleration from the thermal population in solar flare plasma conditions.} {We consider anisotropic Alfvenic fluctuations in the presence of a strong background magnetic field. Exploiting this anisotropy, a set of reduced equations governing non-linear, two-fluid plasma dynamics is derived. The low-$\beta$ limit of this model is used to follow the turbulent cascade of the energy resulting from the non-linear interaction between kinetic Alfven waves, from the large magnetohydrodynamics (MHD) scales with $k_{\perp}\rho_{s}\ll 1$ down to the small "kinetic" scales with $k_{\perp}\rho_{s} \gg 1$, $\rho_{s}$ being the ion sound gyroradius.} {Scaling relations are obtained for the magnitude of the turbulent electromagnetic fluctuations, as a function of $k_{\perp}$ and $k_{\parallel}$, showing that the electric field develops a component parallel to the magnetic field at large MHD scales.} {The spectrum we derive for the parallel electric field fluctuations can be effectively used to model stochastic resonant acceleration and heating of electrons by Alfven waves in solar flare plasma conditions}

Previous numerical studies have identified "phase mixing" of low-frequency Alfven waves as a mean of parallel electric field amplification and acceleration of electrons in a collisionless plasma. Theoretical explanations are given of how this produces an amplification of the parallel electric field, and as a consequence, also leads to enhanced collisionless damping of the wave by energy transfer to the electrons. Our results are based on the properties of the Alfven waves in a warm plasma which are obtained from drift-kinetic theory, in particular, the rate of their electron Landau damping. Phase mixing in a collisionless low-$\beta$ plasma proceeds in a manner very similar to the visco-resistive case, except for the fact that electron Landau damping is the primary energy dissipation channel. The time and length scales involved are evaluated. We also focus on the evolution of the parallel electric field and calculate its maximum value in the course of its amplification.

Anisotropic Alfv\'{e}nic fluctuations with $k_{\parallel}/k_{\perp}\ll 1$ remain at frequencies much smaller than the ion cyclotron frequency in the presence of a strong background magnetic field. Based on the simplest truncation of the electromagnetic gyrofluid equations in a homogeneous plasma, a model for the energy cascade produced by Alfv\'{e}nic turbulence is constructed, which smoothly connect the large magnetohydrodynamics (MHD) scales and the small "kinetic" scales. Scaling relations are obtained for the electromagnetic fluctuations, as a function of $k_{\perp}$ and $k_{\parallel}$. Moreover, a particular attention is paid to the spectral structure of the parallel electric field which is produced by Alfv\'{e}nic turbulence. The reason is the potential implication of this parallel electric field in turbulent acceleration and transport of particles. For electromagnetic turbulence, this issue was raised some time ago in [A. Hasegawa, K. Mima, J. Geophys. Res. {\bf 83} 1117 (1978)].

The newly developed X-ray visibility forward fitting technique is applied to Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) data of a limb flare to investigate the energy and height dependence on sizes, shapes, and position of hard X-ray chromospheric footpoint sources. This provides information about the electron transport and chromospheric density structure. The spatial distribution of two footpoint X-ray sources is analyzed using PIXON, Maximum Entropy Method, CLEAN and visibility forward fit algorithms at nonthermal energies from $\sim 20$ to $\sim 200$ keV. We report, for the first time, the vertical extents and widths of hard X-ray chromospheric sources measured as a function of energy for a limb event. Our observations suggest that both the vertical and horizontal sizes of footpoints are decreasing with energy. Higher energy emission originates progressively deeper in the chromosphere consistent with downward flare accelerated streaming electrons. The ellipticity of the footpoints grows with energy from $\sim 0.5$ at $ \sim 20$ keV to $\sim 0.9$ at $\sim 150$ keV. The positions of X-ray emission are in agreement with an exponential density profile of scale height $\sim 150$~km. The characteristic size of the hard X-ray footpoint source along the limb is decreasing with energy suggesting a converging magnetic field in the footpoint. The vertical sizes of X-ray sources are inconsistent with simple collisional transport in a single density scale height but can be explained using a multi-threaded density structure in the chromosphere.

We investigate the positions and source sizes of X-ray solar flare sources taking into account Compton backscattering (albedo). Using a Monte Carlo simulation of X-ray photon transport including photo-electric absorption and Compton scattering, we calculate the apparent source sizes and positions of X-ray sources at the solar disk for various source sizes, spectral indices and directivities of the primary source. We show that the albedo effect will alter the true source positions and substantially increase the measured source sizes. The source positions are shifted up to $\sim 0.5"$ radially towards the disk centre and 5 arcsecond source sizes can be two times larger even for an isotropic source (minimum albedo effect) at 1 Mm above the photosphere. X-ray sources therefore should have minimum observed sizes, thus FWHM source size (2.35 times second-moment) will be as large as $\sim 7"$ in the 20-50 keV range for a disk-centered point source at a height of 1 Mm ($\sim 1.4"$) above the photosphere. The source size and position change is the largest for flatter primary X-ray spectra, stronger downward anisotropy, for sources closer to the solar disk centre, and between the energies of 30 and 50 keV. Albedo should be taken into account when X-ray footpoint positions, footpoint motions or source sizes from e.g. RHESSI or Yohkoh data are interpreted, and suggest that footpoint sources should be larger in X-rays than in optical or EUV ranges.

Using the potential of two unprecedented missions, STEREO and RHESSI, we study three well observed fast CMEs that occurred close to the limb together with their associated high energy flare emissions in terms of RHESSI HXR spectra and flux evolution. From STEREO/EUVI and STEREO/COR1 data the full CME kinematics of the impulsive acceleration phase up to 4 Rs is measured with a high time cadence of less equal 2.5 min. For deriving CME velocity and acceleration we apply and test a new algorithm based on regularization methods. The CME maximum acceleration is achieved at heights h < 0.4 Rs, the peak velocity at h < 2.1 Rs (in one case as small as 0.5 Rs). We find that the CME acceleration profile and the flare energy release as evidenced in the RHESSI hard X-ray flux evolve in a synchronized manner. These results support the standard flare/CME model which is characterized by a feed-back relationship between the large-scale CME acceleration process and the energy release in the associated flare.

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.

Solar flare hard X-ray spectra from RHESSI are normally interpreted in terms of purely collisional electron beam propagation, ignoring spatial evolution and collective effects. In this paper we present self-consistent numerical simulations of the spatial and temporal evolution of an electron beam subject to collisional transport and beam-driven Langmuir wave turbulence. These wave-particle interactions represent the background plasma's response to the electron beam propagating from the corona to chromosphere and occur on a far faster timescale than coulomb collisions. From these simulations we derive the mean electron flux spectrum, comparable to such spectra recovered from high resolution hard X-rays observations of solar flares with RHESSI. We find that a negative spectral index (i.e. a spectrum that increases with energy), or local minima when including the expected thermal spectral component at low energies, occurs in the standard thick-target model, when coulomb collisions are only considered. The inclusion of wave-particle interactions does not produce a local minimum, maintaining a positive spectral index. These simulations are a step towards a more complete treatment of electron transport in solar flares and suggest that a flat spectrum (spectral index of 0 to 1) down to thermal energies maybe a better approximation instead of a sharp cut-off in the injected electron spectrum.

The collisional thick target model (CTTM) of solar hard X-ray (HXR) bursts has become an almost 'Standard Model' of flare impulsive phase energy transport and radiation. However, it faces various problems in the light of recent data, particularly the high electron beam density and anisotropy it involves.} {We consider how photon yield per electron can be increased, and hence fast electron beam intensity requirements reduced, by local re-acceleration of fast electrons throughout the HXR source itself, after injection.} {We show parametrically that, if net re-acceleration rates due to e.g. waves or local current sheet electric (${\cal E}$) fields are a significant fraction of collisional loss rates, electron lifetimes, and hence the net radiative HXR output per electron can be substantially increased over the CTTM values. In this local re-acceleration thick target model (LRTTM) fast electron number requirements and anisotropy are thus reduced. One specific possible scenario involving such re-acceleration is discussed, viz, a current sheet cascade (CSC) in a randomly stressed magnetic loop.} {Combined MHD and test particle simulations show that local ${\cal E}$ fields in CSCs can efficiently accelerate electrons in the corona and and re-accelerate them after injection into the chromosphere. In this HXR source scenario, rapid synchronisation and variability of impulsive footpoint emissions can still occur since primary electron acceleration is in the high Alfv\'{e}n speed corona with fast re-acceleration in chromospheric CSCs. It is also consistent with the energy-dependent time-of-flight delays in HXR features.

Impulsive solar energetic electrons are often observed in the interplanetary space near the Earth and have an attractive diagnostic potential for poorly understood solar flare acceleration processes. We investigate the transport of solar flare energetic electrons in the heliospheric plasma to understand the role of transport to the observed onset and spectral properties of the impulsive solar electron events. The propagation of energetic electrons in solar wind plasma is simulated from the acceleration region at the Sun to the Earth, taking into account self-consistent generation and absorption of electrostatic electron plasma (Langmuir) waves, effects of non-uniform plasma, collisions and Landau damping. The simulations suggest that the beam-driven plasma turbulence and the effects of solar wind density inhomogeneity play a crucial role and lead to the appearance of a) spectral break for a single power-law injected electron spectrum, with the spectrum flatter below the break, b) apparent early onset of low-energy electron injection, c) the apparent late maximum of low-energy electron injection. We show that the observed onsets, spectral flattening at low energies, and formation of a break energy at tens of keV is the direct manifestation of wave-particle interactions in non-uniform plasma of a single accelerated electron population with an initial power-law spectrum.

A novel method of using hard X-rays as a diagnostic for chromospheric density and magnetic structures is developed to infer sub-arcsecond vertical variation of magnetic flux tube size and neutral gas density.Using Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) X-ray data and the newly developed X-ray visibilities forward fitting technique we find the FWHM and centroid positions of hard X-ray sources with sub-arcsecond resolution ($\sim 0.2"$) for a solar limb flare. We show that the height variations of the chromospheric density and the magnetic flux densities can be found with unprecedented vertical resolution of $\sim$ 150 km by mapping 18-250 keV X-ray emission of energetic electrons propagating in the loop at chromospheric heights of 400-1500 km. Our observations suggest that the density of the neutral gas is in good agreement with hydrostatic models with a scale height of around $140\pm 30$ km. FWHM sizes of the X-ray sources decrease with energy suggesting the expansion (fanning out) of magnetic flux tube in the chromosphere with height. The magnetic scale height $B(z)(dB/dz)^{-1}$ is found to be of the order of 300 km and strong horizontal magnetic field is associated with noticeable flux tube expansion at a height of $\sim$ 900 km.

The Ramaty High Energy Solar Spectroscopic Imager (RHESSI) X-ray data base (February 2002 -- May 2006) has been searched to find solar flares with weak thermal components and flat photon spectra. Using a regularised inversion technique, we determine the mean electron flux distribution from count spectra of a selection of events with flat photon spectra in the 15--20 keV energy range. Such spectral behaviour is expected for photon spectra either affected by photospheric albedo or produced by electron spectra with an absence of electrons in a given energy range, e.g. a low-energy cutoff in the mean electron spectra of non-themal particles. We have found 18 cases which exhibit a statistically significant local minimum (a dip) in the range of 10--20 keV. The positions and spectral indices of events with low-energy cutoff indicate that such features are likely to be the result of photospheric albedo. It is shown that if the isotropic albedo correction was applied, all low-energy cutoffs in the mean electron spectrum were removed and hence the low energy cutoffs in the mean electron spectrum of solar flares above $\sim$12 keV cannot be viewed as real features in the electron spectrum. If low-energy cutoffs exist in the mean electron spectra, the energy of low energy cutoffs should be less than $\sim$12 keV.

Although both electron-ion and electron-electron bremsstrahlung contribute to the hard X-ray emission from solar flares, the latter is normally ignored. Such an omission is not justified at electron (and photon) energies above $\sim 300$ keV, and inclusion of the additional electron-electron bremsstrahlung in general makes the electron spectrum required to produce a given hard X-ray spectrum steeper at high energies. Unlike electron-ion bremsstrahlung, electron-electron bremsstrahlung cannot produce photons of all energies up to the maximum electron energy involved. The maximum possible photon energy depends on the angle between the direction of the emitting electron and the emitted photon, and this suggests a diagnostic for an upper cutoff energy and/or for the degree of beaming of the accelerated electrons. We analyze the large event of January 17, 2005 observed by RHESSI and show that the upward break around 400 keV in the observed hard X-ray spectrum is naturally accounted for by the inclusion of electron-electron bremsstrahlung. Indeed, the mean source electron spectrum recovered through a regularized inversion of the hard X-ray spectrum, using a cross-section that includes both electron-ion and electron-electron terms, has a relatively constant spectral index $\delta$ over the range from electron kinetic energy $E = 200$ keV to $E = 1$ MeV. However, the level of detail discernible in the recovered electron spectrum is not sufficient to determine whether or not any upper cutoff energy exists.

We investigate the signature of the photospheric albedo contribution in solar flare hard X-ray spectra, the effect of low energy cutoffs in electron spectra, and the directivity of hard X-ray emission. Using Ramaty High Energy Solar Spectroscopic Imager (RHESSI) flare data we perform a statistical analysis of spatially integrated spectra and positions of solar flares. We demonstrate clear centre-to-limb variation of photon spectral indices in the 15-20 keV energy range and a weaker dependency in the 20-50 keV range which is consistent with photospheric albedo as the cause. The results also suggest that low-energy cutoffs sometimes inferred in mean electron spectra are an artefact of albedo. We also derive the anisotropy (ratio of downward/observer directed photons) of hard X-ray emission in the 15-20 keV range for various heliocentric angles.

Hard X-ray (HXR) spectroscopy is the most direct method of diagnosing energetic electrons in solar flares. Here we present a technique which allows us to use a single HXR spectrum to determine an effectively stereoscopic electron energy distribution. Considering the Sun's surface to act as a 'Compton mirror' allows us to look at emitting electrons also from behind the source, providing vital information on downward-propagating particles. Using this technique we determine simultaneously the electron spectra of downward and upward directed electrons for two solar flares observed by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The results reveal surprisingly near-isotropic electron distributions, which contrast strongly with the expectations from the standard model which invokes strong downward beaming, including collisional thick-target model.

Hard X-rays and gamma-rays are the most direct signatures of energetic electrons and ions in the sun's atmosphere which is optically thin at these energies and their radiation involves no coherent processes. Being collisional they are complementary to gyro-radiation in probing atmospheric density as opposed to magnetic field and the electrons are primarily 10--100 keV in energy, complementing the (>100 keV) electrons likely responsible for microwave bursts. The pioneering results of the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) are raising the first new major questions concerning solar energetic particles in many years. Some highlights of these results are discussed -- primarily around RHESSI topics on which the authors have had direct research involvement -- particularly when they are raising the need for re-thinking of entrenched ideas. Results and issues are broadly divided into discoveries in the spatial, temporal and spectral domains, with the main emphasis on flare hard X-rays/fast electrons but touching also on gamma-rays/ions, non-flare emissions, and the relationship to radio bursts.

The observed hard X-ray (HXR) flux spectrum $I(\epsilon)$ from solar flares is a combination of primary bremsstrahlung photons $I_P(\epsilon)$ 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 $\gamma (\epsilon)=-d \log I/d \log \epsilon$ is very different for $I_P(\epsilon)$ and for $I(\epsilon)$ 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(\epsilon)$ essentially disappears when the fit is to $I_P(\epsilon)$ - 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.

A multi-wavelength spatial and temporal analysis of solar high energy electrons is conducted using the August 20, 2002 flare of an unusually flat (gamma=1.8) hard X-ray spectrum. The flare is studied using RHESSI, Halpha, radio, TRACE, and MDI observations with advanced methods and techniques never previously applied in the solar flare context. A new method to account for X-ray Compton backscattering in the photosphere (photospheric albedo) has been used to deduce the primary X-ray flare spectra. The mean electron flux distribution has been analysed using both forward fitting and model independent inversion methods of spectral analysis. We show that the contribution of the photospheric albedo to the photon spectrum modifies the calculated mean electron flux distribution, mainly at energies below 100 keV. The positions of the Halpha emission and hard X-ray sources with respect to the current-free extrapolation of the MDI photospheric magnetic field and the characteristics of the radio emission provide evidence of the closed geometry of the magnetic field structure and the flare process in low altitude magnetic loops. In agreement with the predictions of some solar flare models, the hard X-ray sources are located on the external edges of the Halpha emission and show chromospheric plasma heated by the non-thermal electrons. The fast changes of Halpha intensities are located not only inside the hard X-ray sources, as expected if they are the signatures of the chromospheric response to the electron bombardment, but also away from them.

Recent progress in solar Hard X-ray (HXR) observations with RHESSI data and methods for spectral inversion allow us to study model-independent mean electron flux spectra in solar flares. We report several hard X-ray events observed by RHESSI in which the photon spectra $I(\epsilon)$ are such that the inferred source mean electron spectra are not consistent with the standard model of collisional transport in solar flares. The observed photon spectra are so flat locally that the recovered mean electron flux spectra show a dip around 17-31 keV. While we note that alternative explanations, unrelated to electron transport, have not been ruled out, we focus on the physical implications of this tentative result for the collisional thick-target model.

Recent progress in the diagnosis of flare fast particles is critically discussed with the main emphasis on high resolution Hard X-Ray (HXR) data from RHESSI and coordinated data from other instruments. Spectacular new photon data findings are highlighted as are advances in theoretical aspects of their use as fast particle diagnostics, and some important comparisons made with interplanetary particle data. More specifically the following topics are addressed (a) RHESSI data on HXR (electron) versus gamma-ray line (ion) source locations. (b) RHESSI hard X-ray source spatial structure in relation to theoretical models and loop density structure. (c) Energy budget of flare electrons and the Neupert effect. (d) Spectral deconvolution methods including blind target testing and results for RHESSI HXR spectra, including the reality and implications of dips inferred in electron spectra (e) The relation between flare in-situ and interplanetary particle data.

The deduction from solar flare X-ray photon spectroscopic data of the energy dependent model-independent spectral index is considered as an inverse problem. Using the well developed regularization approach we analyze the energy dependency of spectral index for a high resolution energy spectrum provided by Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The regularization technique produces much smoother derivatives while avoiding additional errors typical of finite differences. It is shown that observations imply a spectral index varying significantly with energy, in a way that also varies with time as the flare progresses. The implications of these findings are discussed in the solar flare context.

Finite thermal velocity modifications to electron slowing-down rates may be important for the deduction of solar flare total electron energy. Here we treat both slowing-down and velocity diffusion of electrons in the corona at flare temperatures, for the case of a simple, spatially homogeneous source. Including velocity diffusion yields a consistent treatment of both `accelerated' and `thermal' electrons. It also emphasises that one may not invoke finite thermal velocity target effects on electron lifetimes without simultaneously treating the contribution to the observed X-ray spectrum from thermal electrons. We present model calculations of the X-ray spectra resulting from injection of a power-law energy distribution of electrons into a source with finite temperature. Reducing the power-law distribution low-energy cutoff to lower and lower energies only increases the relative magnitude of the thermal component of the spectrum, because the lowest energy electrons simply join the background thermal distribution. Acceptable fits to RHESSI flare data are obtained using this model. These also demonstrate, however, that observed spectra may in consequence be acceptably consistent with rather a wide range of injected electron parameters.

Hard X-ray spectra in solar flares provide knowledge of the electron spectrum that results from acceleration and propagation in the solar atmosphere. However, the inference of the electron spectra from solar X-ray spectra is an ill-posed inverse problem. Here we develop and apply an enhanced regularization algorithm for this process making use of physical constraints on the form of the electron spectrum. The algorithm incorporates various features not heretofore employed in the solar flare context: Generalized Singular Value Decomposition (GSVD) to deal with different orders of constraints; rectangular form of the cross-section matrix to extend the solution energy range; regularization with various forms of the smoothing operator; and "preconditioning" of the problem. We show by simulations that this technique yields electron spectra with considerably more information and higher quality than previous algorithms.

Kontar et al. (2004) have shown how to recover mean source electron spectra in solar flares through a physical constraint regularization analysis of the bremsstrahlung photon spectra that they produce. They emphasize the use of non-square inversion techniques, and preconditioning combined with physical properties of the spectra to achieve the most meaningful solution to the problem. Higher-order regularization techniques may be used to generate ${overline F}(E)$ forms with certain desirable properties (e.g., higher order derivatives). They further note that such analyses may be used to infer properties of the electron energy spectra at energies well above the maximum photon energy observed. In this paper we apply these techniques to data from a solar flare observed by RHESSI on 26 February, 2002. Results using different orders of regularization are presented and compared for various time intervals. Clear evidence is presented for a change in the value of the high-energy cutoff in the mean source electron spectrum with time. We also show how the construction of the injected (accelerated) electron spectrum $F_0(E_0)$ (assuming that Coulomb collisions in a cold target dominate the electron energetics) is facilitated by the use of higher-order regularization methods.