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.

Eulerian and Lagrangian tools are used to detect coherent structures in the velocity and magnetic fields of a mean--field dynamo, produced by direct numerical simulations of the three--dimensional compressible magnetohydrodynamic equations with an isotropic helical forcing and moderate Reynolds number. Two distinct stages of the dynamo are studied, the kinematic stage, where a seed magnetic field undergoes exponential growth, and the saturated regime. It is shown that the Lagrangian analysis detects structures with greater detail, besides providing information on the chaotic mixing properties of the flow and the magnetic fields. The traditional way of detecting Lagrangian coherent structures using finite--time Lyapunov exponents is compared with a recently developed method called function M. The latter is shown to produce clearer pictures which readily permit the identification of hyperbolic regions in the magnetic field, where chaotic transport/dispersion of magnetic field lines is highly enhanced.