The X-ray emission of coronal flare sources can be explained considering the kappa electron distribution. Motivated by this fact, we study the problem of how hot plasma with the kappa distribution of electrons is confined in these sources. For comparison, we analyze the same problem but with the Maxwellian distribution. We use a 3-D particle-in-cell code, which is large in one direction and thus effectively only one-dimensional, but describing all electromagnetic effects. In the case with the Maxwellian distribution, and in agreement with the previous studies, we show a formation of the double layer at the hot-cold transition region that suppresses the flux of hot electrons from hot plasma into the cold one. In the case with the kappa distribution, contrary to the Maxwellian case, we found that there are several fronts with the double layers in the hot-cold transition region. It is caused by a more extended tail in the kappa case than in the Maxwellian one. The electrons from the extended tail freely escape from the hot plasma into a cold one. They form a beam which generates the return current and also Langmuir turbulence, where at some locations Langmuir waves are accumulated. At these locations owing to the ponderomotive force, Langmuir waves generate density depressions, where the double layers with the thermal fronts, suppressing the hot electron flux, are formed. We also show how protons are accelerated in these processes. Finally, we compared the kappa and Maxwellian cases and discussed how these processes could be observed.

Zebras were observed not only in the solar radio emission but also in radio emissions of Jupiter and Crab Nebula pulsar. In their models, growth rates of the electrostatic waves play an important role. Considering the plasma composed from the thermal background plasma and hot and rare component with the Dory-Guest-Harris distribution, we compute the growth rates $\gamma$ and dispersion branches of the electrostatic waves in the $\omega - k_\perp$ domain. We show complexity of the electrostatic wave branches in the upper-hybrid band. In order to compare the results, which we obtained using the kinetic theory and Particle-in-cell (PIC) simulations, we define and compute the integrated growth rate $\Gamma$, where the "characteristic width" of dispersion branches was considered. We found a very good agreement between the integrated growth rates and those from PIC simulations. For maximal and minimal $\Gamma$ we showed locations of dispersion branches in the $\omega - k_\perp$ domain. We found that $\Gamma$ has a maximum when the dispersion branches not only cross the region with high growth rates $\gamma$, but when the dispersion branches in this region are sufficiently long and wide. We also mentioned the effects of changes in the background plasma and hot component temperatures.

Fine structures of radio bursts play an important role in diagnostics of the solar flare plasma. Among them the zebras, which are prevalently assumed to be generated by the double plasma resonance instability, belong to the most important. In this paper we compute the growth rate of this instability for two types of the electron distribution: a) for the power-law distribution and b) for the kappa distribution, in the both cases with the loss-cone type anisotropy. We found that the growth rate of the upper-hybrid waves for the power-law momentum distribution strongly depends on the pitch-angle boundary. The maximum growth rate was found for the pitch-angle $\theta_\mathrm{c} \approx$ 50$^\circ$. For small angles the growth rate profile is very flat and for high pitch-angles the wave absorption occurs. Furthermore, analyzing the growth rate of the upper hybrid waves for the kappa momentum distribution we found that a decrease of the characteristic momentum $p_\kappa$ shifts the maximum of the growth rate to lower values of the ratio of the electron-plasma and electron-cyclotron frequencies, and the frequency widths of the growth rate peaks are very broad. But, if we consider the kappa distribution which is isotropic up to some large momentum $p_m$ and anisotropic with loss-cone above this momentum then distinct peaks of the growth rate appear and thus distinct zebra stripes can be generated. It means that the restriction for small momenta for the anisotropic part of distributions is of principal importance for the zebra stripes generation. Finally, for the 1 August 2010 zebra stripes, the growth rates in dependence on radio frequency were computed. It was shown that in this case the growth rate peaks are more distinct than in usually presented dependencies of growth rates on the ratio of the plasma and cyclotron frequencies.

We study processes occurring after a sudden heating of the chromosphere at the flare arcade footpoints which is assumed to be caused by particle beams. For the numerical simulations we adopt a 2-D magnetohydrodynamic (MHD) model, in which we solve a full set of the time-dependent MHD equations by means of the FLASH code, using the Adaptive Mesh Refinement (AMR) method. In the initial state we consider a model of the solar atmosphere with densities according to the VAL-C model and the magnetic field arcade having the X-point structure above, where the magnetic reconnection is assumed. We found that the sudden pulse-beam heating of the chromosphere at the flare arcade footpoints generates magnetohydrodynamic shocks, one propagating upwards and the second one propagating downwards in the solar atmosphere. The downward moving shock is reflected at deep and dense atmospheric layers and triggers oscillations of these layers. These oscillations generate the upwards moving magnetohydrodynamic waves that can influence the above located magnetic reconnection in a quasi-periodic way. Because these processes require a sudden heating in very localized regions in the chromosphere therefore they can be also associated with seismic waves.

The double plasma resonance (DPR) instability plays a basic role in the generation of solar radio zebras. In the plasma, consisting of the loss-cone type distribution of hot electrons and much denser and colder background plasma, this instability generates the upper-hybrid waves, which are then transformed into the electromagnetic waves and observed as radio zebras. In the present paper we numerically study the double plasma resonance instability from the point of view of the zebra interpretation. We use a 3-dimensional electromagnetic particle-in-cell (3-D PIC) relativistic model. First using the multi-mode model, we study details of the double plasma resonance instability. We show how the distribution function of hot electrons changes during this instability. Then we show that there is a very good agreement between results obtained by the multi-mode and specific-mode models, which is caused by a dominance of the wave with the maximal growth rate. Therefore, for computations in a broad range of model parameters, we use the specific-mode model. We compute the maximal growth rates of the double plasma resonance instability. The results are compared with the analytical ones. We find a very good agreement between numerical and analytical growth rates. We also compute saturation energies of the upper-hybrid waves in a very broad range of parameters. We find that the saturation energies of the upper-hybrid waves show maxima and minima at almost the same values of $\omega_\mathrm{UH}/\omega_\mathrm{ce}$ as the growth rates. Furthermore, we find that the saturation energy of the upper-hybrid waves is proportional to the density of hot electrons. The maximum saturated energy can be up to one percent of the kinetic energy of hot electrons. All these findings can be used in the interpretation of solar radio zebras.

We estimated the brightness temperature of radio zebras (zebra pattern -- ZP), considering that ZPs are generated in loops having an exponential density profile in their cross-section. We took into account that when in plasma there is a source emitting in all directions, then in the escape process from the plasma the emission obtains a directional character nearly perpendicular to the constant-density profile. Owing to the high directivity of the plasma emission the region from which the emission escapes can be very small. We estimated the brightness temperature of three observed ZPs for two values of the density scale height (1 and 0.21 Mm) and two values of the loop width (1 and 2 arcsec). In all cases high brightness temperatures were obtained. For the higher value of the density scale height, the brightness temperature was estimated as 1.1 $\times$ 10$^{15}$ - 1.3 $\times$ 10$^{17}$ K, and for the lower value as 4.7 $\times$ 10$^{13}$ - 5.6 $\times$ 10$^{15}$ K. We also computed the saturation energy density of the upper-hybrid waves (which according to the double plasma resonance model are generated in the zebra source) using a 3D particle-in-cell model with the loss-cone type of distribution of hot electrons. We found that this saturated energy is proportional to the ratio of hot electron and background plasma densities. Thus, comparing the growth rate and collisional damping of the upper-hybrid waves, we estimated minimal densities of hot electrons as well as the minimal value of the saturation energy density of the upper-hybrid waves. Finally, we compared the computed energy density of the upper-hybrid waves with the energy density of the electromagnetic waves in the zebra source and thus estimated the efficiency of the wave transformation.

Using the FLASH code, which solves the full set of the two-dimensional (2-D) non-ideal (resistive) time-dependent magnetohydrodynamic (MHD) equations, we study processes during the magnetic reconnection in a vertical gravitationally stratified current sheet. We show that during these processes, which correspond to processes in solar flares, plasmoids are formed due to the tearing mode instability of the current sheet. These plasmoids move upwards or downwards along the vertical current sheet, and some of them merge into larger plasmoids. We study the density and temperature structure of these plasmoids and their time evolution in details. We found that during the merging of two plasmoids the resulting larger plasmoid starts to oscillate; in our model with a $\sim 25~\mathrm{s}$ period. On the other hand, the plasmoid moving downwards merges with the underlying flare arcade which also starts to oscillate during this process; in our model with a $\sim 35~\mathrm{s}$ period. It is shown that the merging process of plasmoid with the flare arcade is a complex process as presented by complex density and temperature structures of the oscillating arcade. Moreover, all these processes are associated with magnetoacoustic waves produced by the motion and merging of plasmoids.

Context. The zebra patterns observed in solar radio emission are very important for are plasma diagnostics. The most promising model of these patterns is based on double plasma resonance instability, which generates upper-hybrid waves, which can be then transformed into the zebra emission. Aims. We aim to study in detail the double plasma resonance instability of hot electrons, together with a much denser thermal background plasma. In particular, we analyse how the growth rate of the instability depends on the temperature of both the hot plasma and background plasma components. Methods. We numerically integrated the analysed model equations, using Python and Wolfram Mathematica. Results. We found that the growth-rate maxima of the upper-hybrid waves for non-zero temperatures of both the hot and background plasma are shifted towards lower frequencies comparing to the zero temperature case. This shift increases with an increase of the harmonic number s of the electron cyclotron frequency and temperatures of both hot and background plasma components. We show how this shift changes values of the magnetic eld strength estimated from observed zebras. We confirmed that for a relatively low hot electron temperature, the dependence of growth rate vs. both the ratio of the electron plasma and electron cyclotron frequencies expresse distinct peaks, and by increasing this temperature these peaks become smoothed. We found that in some cases, the values of wave number vector components for the upper-hybrid wave for the maximal growth rate strongly deviate from their analytical estimations. We confirmed the validity of the assumptions used when deriving model equations.

A good observation of preflare activities is important for us to understand the origin and triggering mechanism of solar flares, and to predict the occurrence of solar flares. This work presents the characteristics of microwave spectral fine structures as preflare activities of four solar flares observed by Ond\v{r}ejov radio spectrograph in the frequency range of 0.8--2.0 GHz. We found that these microwave bursts which occurred 1--4 minutes before the onset of flares have spectral fine structures with relatively weak intensities and very short timescales. They include microwave quasi-periodic pulsations (QPP) with very short period of 0.1-0.3 s and dot bursts with millisecond timescales and narrow frequency bandwidths. Accompanying these microwave bursts, there are filament motions, plasma ejection or loop brightening on the EUV imaging observations and non-thermal hard X-ray emission enhancements observed by RHESSI. These facts may reveal certain independent non-thermal energy releasing processes and particle acceleration before the onset of solar flares. They may be conducive to understand the nature of solar flares and predict their occurrence.

We address a possibility of the flare process initiation and further maintenance of its energy release due to a transformation of sunspot longitudinal waves into transverse magnetic loop oscillations with initiation of reconnection. This leads to heating maintaining after the energy release peak and formation of a flat stage on the X-ray profile. We applied the time-distance plots and pixel wavelet filtration (PWF) methods to obtain spatio-temporal distribution of wave power variations in SDO/AIA data. To find magnetic waveguides, we used magnetic field extrapolation of SDO/HMI magnetograms. The propagation velocity of wave fronts was measured from their spatial locations at specific times. In correlation curves of the 17 GHz (NoRH) radio emission we found a monotonous energy amplification of 3-min waves in the sunspot umbra before the 2012 June 7 flare. This dynamics agrees with an increase in the wave-train length in coronal loops (SDO/AIA, 171 {\AA}) reaching the maximum 30 minutes prior to the flare onset. A peculiarity of this flare time profile in soft X-rays (RHESSI, 3-25 keV) is maintaining the constant level of the flare emission for 10 minutes after the short impulse phase, which indicates at the energy release continuation. Throughout this time, we found 30-sec period transverse oscillations of the flare loop in the radio-frequency range (NoRH, 17 GHz). This periodicity is apparently related to the transformation of propagating longitudinal 3-min waves from the sunspot into the loop transverse oscillations. The magnetic field extrapolation showed the existence of the magnetic waveguide (loop) connecting the sunspot with the energy release region. A flare loop heating can be caused by the interaction (reconnections) of this transversally oscillating waveguide with the underlying twisted loops.

We present SDO/AIA observations of an eruptive X-class flare of July 12, 2012, and compare its evolution with the predictions of a 3D numerical simulation. We focus on the dynamics of flare loops that are seen to undergo slipping reconnection during the flare. In the AIA 131A observations, lower parts of 10 MK flare loops exhibit an apparent motion with velocities of several tens of km/s along the developing flare ribbons. In the early stages of the flare, flare ribbons consist of compact, localized bright transition-region emission from the footpoints of the flare loops. A DEM analysis shows that the flare loops have temperatures up to the formation of Fe XXIV. A series of very long, S-shaped loops erupt, leading to a CME observed by STEREO. The observed dynamics are compared with the evolution of magnetic structures in the "standard solar flare model in 3D". This model matches the observations well, reproducing both the apparently slipping flare loops, S-shaped erupting loops, and the evolution of flare ribbons. All of these processes are explained via 3D reconnection mechanisms resulting from the expansion of a torus-unstable flux rope. The AIA observations and the numerical model are complemented by radio observations showing a noise storm in the metric range. Dm-drifting pulsation structures occurring during the eruption indicate plasmoid ejection and enhancement of reconnection rate. The bursty nature of radio emission shows that the slipping reconnection is still intermittent, although it is observed to persist for more than an hour.

We study two modifications of the collisional thick-target model (CTTM) based on the global and local re-acceleration of non-thermal electrons by static and stochastic electric fields during their transport from the coronal acceleration site to the thick-target region in the chromosphere. We concentrate on a comparison of the non-thermal electron distribution functions, chromospheric energy deposits, and HXR spectra obtained for both considered modifications with the CTTM itself. The results were obtained using a relativistic test-particle approach. We simulated the transport of non-thermal electrons with a power-law spectrum including the influence of scattering, energy losses, magnetic mirroring, and also the effects of the electric fields corresponding to both modifications of the CTTM. We show that both modifications of the CTTM change the outcome of the chromospheric bombardment in several aspects. The modifications lead to an increase in chromospheric energy deposit, change of its spatial distribution, and a substantial increase in the corresponding HXR spectrum intensity.

We present a model for dm-fiber bursts that is based on assuming fast sausage magnetoacoustic wave trains that propagate along a dense vertical filament or current sheet. Eight groups of dm-fiber bursts that were observed during solar flares were selected and analyzed by the wavelet analysis method. To model these fiber bursts we built a semi-empirical model. We also did magnetohydrodynamic simulations of a propagation of the magnetoacoustic wave train in a vertical and gravitationally stratified current sheet. In the wavelet spectra of the fiber bursts computed at different radio frequencies we found the wavelet tadpoles, whose head maxima have the same frequency drift as the drift of fiber bursts. It indicates that the drift of these fiber bursts can be explained by the propagating fast sausage magnetoacoustic wave train. Using new semi-empirical and magnetohydrodynamic models with a simple radio emission model we generated the artificial radio spectra of the fiber bursts, which are similar to the observed ones.

Magnetic reconnection is commonly considered as a mechanism of solar (eruptive) flares. A deeper study of this scenario reveals, however, a number of open issues. Among them is the fundamental question, how the magnetic energy is transferred from large, accumulation scales to plasma scales where its actual dissipation takes place. In order to investigate this transfer over a broad range of scales we address this question by means of a high-resolution MHD simulation. The simulation results indicate that the magnetic-energy transfer to small scales is realized via a cascade of consecutive smaller and smaller flux-ropes (plasmoids), in analogy with the vortex-tube cascade in (incompressible) fluid dynamics. Both tearing and (driven) "fragmenting coalescence" processes are equally important for the consecutive fragmentation of the magnetic field (and associated current density) to smaller elements. At the later stages a dynamic balance between tearing and coalescence processes reveals a steady (power-law) scaling typical for cascading processes. It is shown that cascading reconnection also addresses other open issues in solar flare research such as the duality between the regular large-scale picture of (eruptive) flares and the observed signatures of fragmented (chaotic) energy release, as well as the huge number of accelerated particles. Indeed, spontaneous current-layer fragmentation and formation of multiple channelised dissipative/acceleration regions embedded in the current layer appears to be intrinsic to the cascading process. The multiple small-scale current sheets may also facilitate the acceleration of a large number of particles. The structure, distribution and dynamics of the embedded potential acceleration regions in a current layer fragmented by cascading reconnection are studied and discussed.

Based on our recent MHD simulations, first, a concept of the successive merging of plasmoids and fragmentation in the current sheet in the standard flare model is presented. Then, using a 2.5-D electromagnetic particle-in-cell model with free boundary conditions, these processes were modelled on the kinetic level of plasma description. We recognized the plasmoids which mutually interacted and finally merged into one large plasmoid. Between interacting plasmoids further plasmoids and current sheets on smaller and smaller spatial scales were formed in agreement with the fragmentation found in MHD simulations. During interactions (merging - coalescences) of the plasmoids the electrons were very efficiently accelerated and heated. We found that after a series of such merging processes the electrons in some regions reached the energies relevant for the emission in the hard X-ray range. Considering these energetic electrons and assuming the plasma density 10^9-10^10 cm^{-3} and the source volume as in the December 31, 2007 flare (Krucker at al. 2010, ApJ 714, 1108), we computed the X-ray spectra as produced by the bremsstrahlung emission process. Comparing these spectra with observations, we think that these processes can explain the observed above-the-loop-top hard X-ray sources. Furthermore, we show that the process of a fragmentation between two merging plasmoids can generate the narrowband dm-spikes. Formulas for schematic fractal reconnection structures were derived. Finally, the results were discussed.

In the paper by B\'arta et al. (arXive:astro-ph:/1011.4035, 2010) the authors addressed some open questions of the CSHKP scenario of solar flares by means of high-resolution MHD simulations. They focused, in particular, on the problem of energy transfer from large to small scales in decaying flare current sheet (CS). Their calculations suggest, that magnetic flux-ropes (plasmoids) are formed in full range of scales by a cascade of tearing and coalescence processes. Consequently, the initially thick current layer becomes highly fragmented. Thus, the tearing and coalescence cascade can cause an effective energy transfer across the scales. In the current paper we investigate whether this mechanism actually applies in solar flares. We extend the MHD simulation by deriving model-specific features that can be looked for in observations. The results of the underlying MHD model showed that the plasmoid cascade creates a specific hierarchical distribution of non-ideal/acceleration regions embedded in the CS. We therefore focus on the features associated with the fluxes of energetic particles, in particular on the structure and dynamics of emission regions in flare ribbons. We assume that the structure and dynamics of diffusion regions embedded in the CS imprint themselves into structure and dynamics of flare-ribbon kernels by means of magnetic-field mapping. Using the results of the underlying MHD simulation we derive the expected structure of ribbon emission and we extract selected statistical properties of the modelled bright kernels. Comparing the predicted emission and its properties with the observed ones we obtain a good agreement of the two.

Imaging microwave observations of an eruptive, partially occulted solar flare on 18 April 2001 suggest that the global structure of the event can be described by the helical kink instability of a twisted magnetic flux rope. This model is suggested by the inverse gamma shape of the source exhibiting crossing legs of a rising flux loop and by evidence that the legs interact at or near the crossing point. The interaction is reflected by the location of peak brightness near the crossing point and by the formation of superimposed compact nonthermal sources most likely at or near the crossing point. These sources propagate upward along both legs, merge into a single, bright source at the top of the structure, and continue to rise at a velocity $>\!1000$~km\,s$^{-1}$. The compact sources trap accelerated electrons which radiate in the radio and hard X-ray ranges. This suggests that they are plasmoids, although their internal structure is not revealed by the data. They exhibit variations of the radio brightness temperature at a characteristic time scale of $\sim\!40$~s, anti-correlated to their area, which also support their interpretation as plasmoids. Their propagation path differs from the standard scenario of plasmoid formation and propagation in the flare current sheet, suggesting the helical current sheet formed by the instability instead.

Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $\Phi\gtrsim6\pi$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the original rope, it can form a new, low-lying, less twisted flux rope. The new flux rope is pushed downward by the reconnection outflow, which typically forces it to break as well by reconnecting with the ambient field. The top part of the original rope, largely rooted in the sources of the ambient flux after the break-up, can fully erupt or be halted at low heights, producing a "failed eruption." The helical current sheet associated with the instability is squeezed between the approaching legs, temporarily forming a double current sheet. The leg-leg reconnection proceeds at a high rate, producing sufficiently strong electric fields that it would be able to accelerate particles. It may also form plasmoids, or plasmoid-like structures, which trap energetic particles and propagate out of the reconnection region up to the top of the erupting flux rope along the helical current sheet. The kinking of a highly twisted flux rope involving leg-leg reconnection can explain key features of an eruptive but partially occulted solar flare on 18 April 2001, which ejected a relatively compact hard X-ray and microwave source and was associated with a fast coronal mass ejection.

With the invention of new far-infrared (FIR) and radio mm and sub-mm instruments (DESIR on SMESE satellite, ESO-ALMA), there is a growing interest in observations and analysis of solar flares in this so far unexplored wavelength region. Two principal radiation mechanisms play a role: the synchrotron emission due to accelerated particle beams moving in the magnetic field and the thermal emission due to the energy deposit in the lower atmospheric layers. In this contribution we explore the time-dependent effects of beams on thermal FIR and radio continua. We show how and where these continua are formed in the presence of time dependent beam heating and non-thermal excitation/ionisation of the chromospheric hydrogen plasma.

Aims: To calculate the power budget for electron acceleration and the efficiency of the plasma emission mechanism in a post-flare decimetric continuum source. Methods: We have imaged a high brightness temperature ($\sim 10^{9}$K) post-flare source at 1060 MHz with the Giant Metrewave Radio Telescope (GMRT). We use information from these images and the dynamic spectrum from the Hiraiso spectrograph together with the theoretical method described in Subramanian & Becker (2006) to calculate the power input to the electron acceleration process. The method assumes that the electrons are accelerated via a second-order Fermi acceleration mechanism. Results: We find that the power input to the nonthermal electrons is in the range $3\times 10^{25}$--$10^{26}$ erg/s. The efficiency of the overall plasma emission process starting from electron acceleration and culminating in the observed emission could range from $2.87\times 10^{-9}$ to $2.38 \times 10^{-8}$.