Type III radio bursts are radio emissions associated with solar flares. They are considered to be caused by electron beams in the solar corona. Magnetic reconnection is a possible accelerator of electron beams in the course of solar flares which causes unstable distribution functions, and density inhomogeneities. The properties of radio emission by electron beams in such environment are, however, still poorly understood. We capture the non-linear kinetic plasma processes of radio emissions in such plasmas by utilizing fully-kinetic Particle-In-Cell (PIC) code numerical simulations. Our model takes into account initial velocity distribution functions as they are supposed to be created by magnetic reconnection. These velocity distribution functions allow two distinct mechanisms of radio wave emissions: plasma emissions due to plasma-wave interactions and so called electron cyclotron maser emissions (ECME) due to wave-particle interactions. Our most important finding is that the number of harmonics of Langmuir waves increases with the density inhomogeneities. The harmonics are generated by the interaction of beam-generated Langmuir waves and their harmonics. In addition, we also find evidence for transverse harmonic electromagnetic wave emissions due to a coalescence of beam-generated and fundamental Langmuir waves with a vanishing wavevector. We investigate the effects of density inhomogeneities on the conversion process of the free energy of the electron beams to electrostatic and electromagnetic waves and the frequency shift of electron resonances caused by perpendicular gradients in the beam velocity distribution function. Our findings explain the observation of Langmuir waves and their harmonics in solar radio bursts and of the observed frequency shifts in these emissions.

The formation of collisionless shock fronts is an ubiquitous phenomenon in space plasma environments. In the solar wind shocks might accompany coronal mass ejections, while even more violent events, such as supernovae, produce shock fronts traveling at relativistic speeds. While the basic concepts of shock formation and particle acceleration in their vicinity are known, many details on a micro-physical scope are still under discussion. In recent years the hybrid kinetic simulation approach has allowed to study the dynamics and acceleration of protons and heavier ions in great detail. However, Particle-in-Cell codes allow to study the process including also electron dynamics and the radiation pressure. Additionally a further numerical method allows for crosschecking results. We therefore investigate shock formation and particle acceleration with a fully kinetic particle-in-cell code. Besides electrons and protons we also include helium and carbon ions in our simulations of a quasi-parallel shock. We are able to reproduce characteristic features of the energy spectra of the particles, such as the temperature ratios of the different ion species in the downstream which scale with the ratio of particle mass to charge. We also find that approximately 12-15% of the energy of the unperturbed upstream is transferred to the accelerated particles escaping the shock.

We analyzed properties of waves excited by mildly relativistic electron beams propagating along magnetic field with a ring-shape perpendicular momentum distribution in neutral and current-free solar coronal plasmas. These plasmas are subject to both the beam and the electron cyclotron maser (ECM) instabilities driven by the positive momentum gradient of the ring-beam electron distribution in the directions parallel and perpendicular to the ambient magnetic field, respectively. To explore the related kinetic processes self-consistently, 2.5-dimensional fully kinetic particle-in-cell (PIC) simulations were carried out. To quantify excited wave properties in different coronal conditions, we investigated the dependence of their energy and polarization on the ring-beam electron density and magnetic field. In general, electrostatic waves dominate the energetics of waves and nonlinear waves are ubiquitous. In weakly magnetized plasmas, where the electron cyclotron frequency $\omega_{ce}$ is lower than the electron plasma frequency $\omega_{pe}$, it is difficult to produce escaping electromagnetic waves with frequency $\omega > \omega_{pe}$ and small refractive index $|c k / \omega| < 1$ ($k$ and $c$ are the wavenumber and the light speed, respectively). Highly polarized and anisotropic escaping electromagnetic waves can, however, be effectively excited in strongly magnetized plasmas with $\omega_{ce}/\omega_{pe} \geq 1$. The anisotropy of the energy, circular polarization degree (CPD), and spectrogram of these escaping electromagnetic waves strongly depend on the number density ratio of the ring-beam electrons to the background electrons. In particular, their CPDs can vary from left-handed to right-handed with the decrease of the ring-beam density, which may explain some observed properties of solar radio bursts (e.g., radio spikes) from the solar corona.

Although turbulence has been conjectured to be important for magnetic reconnection, still very little is known about its role in collisionless plasmas. Previous attempts to quantify the effect of turbulence on reconnection usually prescribed Alfv\'enic or other low-frequency fluctuations or investigated collisionless kinetic effects in just two-dimensional configurations and antiparallel magnetic fields. In view of this, we analyzed the kinetic turbulence self-generated by three-dimensional guide-field reconnection through force-free current sheets in frequency and wavenumber spaces, utilizing 3D particle-in-cell code numerical simulations. Our investigations reveal reconnection rates and kinetic turbulence with features similar to those obtained by current in-situ spacecraft observations of MMS as well as in the laboratory reconnection experiments MRX, VTF and \textsc{Vineta}-II. In particular we found that the kinetic turbulence developing in the course of 3D guide-field reconnection exhibits a broadband power-law spectrum extending beyond the lower-hybrid frequency and up to the electron frequencies. In the frequency space the spectral index of the turbulence appeared to be close to -2.8 at the reconnection X-line. In the wavenumber space it also becomes -2.8 as soon as the normalized reconnection rate reaches 0.1. The broadband kinetic turbulence is mainly due to current-streaming and electron-flow-shear instabilities excited in the sufficiently thin current sheets of kinetic reconnection. The growth of the kinetic turbulence corresponds to high reconnection rates which exceed those of fast laminar, non-turbulent reconnection.

We report a newly found two-stage mechanism of electron acceleration near X-lines of 3D collisionless guide-field magnetic reconnection in the non-relativistic regime typical, e.g., for stellar coronae. We found that after electrons are first pre-accelerated during the linear growth of reconnection, they become additionally accelerated in the course of the nonlinear stage of 3D guide-field magnetic reconnection. This additional acceleration is due to the filamentation of electric and magnetic fields caused by streaming instabilities. In addition to enhanced parallel electric fields, the filamentation leads to additional curvature-driven electron acceleration in the guide-field direction. As a result, part of the the accelerated electron spectra becomes a power law with a spectral index of $\sim-1.6$ near the X-line. This second stage of acceleration due to nonlinear reconnection is relevant for the production of energetic electrons in, e.g., thin current sheets of stellar coronae.

Aims: We investigate the electron acceleration in convective electric fields of cascading magnetic reconnection in a flaring solar corona and show the resulting hard X-ray (HXR) radiation spectra caused by Bremsstrahlung for the coronal source. Methods: We perform test particle calculation of electron motions in the framework of a guiding center approximation. The electromagnetic fields and their derivatives along electron trajectories are obtained by linearly interpolating the results of high-resolution adaptive mesh refinement (AMR) MHD simulations of cascading magnetic reconnection. Hard X-ray (HXR) spectra are calculated using an optically thin Bremsstrahlung model. Results: Magnetic gradients and curvatures in cascading reconnection current sheet accelerate electrons: trapped in magnetic islands, precipitating to the chromosphere and ejected into the interplanetary space. The final location of an electron is determined by its initial position, pitch angle and velocity. These initial conditions also influence electron acceleration efficiency. Most of electrons have enhanced perpendicular energy. Trapped electrons are considered to cause the observed bright spots along coronal mass ejection CME-trailing current sheets as well as the flare loop-top HXR emissions.

The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores.

We have proposed a method to calculate the relative magnetic helicity in a finite volume as given the magnetic field in the former paper (Yang et al. {\it Solar Physics}, {\bf 283}, 369, 2013). This method requires that the magnetic flux to be balanced on all the side boundaries of the considered volume. In this paper, we propose a scheme to obtain the vector potentials at the boundaries to remove the above restriction. We also used a theoretical model (Low and Lou, {\it Astrophys. J.} {\bf 352}, 343, 1990) to test our scheme.

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.

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.

We present a newly developed approach to solar coronal magnetic field extrapolation from vector magnetograms, based on the Principle of Minimum Dissipation Rate (MDR). The MDR system was derived from a variational problem that is more suitable for an open and externally driven system, like the solar corona. The resulting magnetic field equation is more general than force-free. Its solution can be expressed as the superposition of two linear (constant-$\alpha$) force-free fields (LFFFs) with distinct $\alpha$ parameters, and one potential field. Thus the original extrapolation problem is decomposed into three LFFF extrapolations, utilizing boundary data. The full MDR-based approach requires two layers of vector magnetograph measurements on solar surface, while a slightly modified practical approach only requires one. We test both approaches against 3D MHD simulation data in a finite volume. Both yield quantitatively good results. The errors in the magnetic energy estimate are within a few percents. In particular, the main features of relatively strong perpendicular current density structures, representative of the non-force freeness of the solution, are well recovered.