A two-field Hamiltonian gyrofluid model for kinetic Alfvén waves retaining ion finite Larmor radius corrections, parallel magnetic field fluctuations and electron inertia, is used to study turbulent cascades from the magnetohydrodynamic (MHD) to the sub-ion scales. Special attention is paid to the case of imbalance between waves propagating along or opposite to the ambient magnetic field. For weak turbulence in the absence of electron inertia, kinetic equations for the spectral density of the conserved quantities (total energy and generalized cross-helicity) are obtained. They provide a global description, matching between the regimes of reduced MHD at large scales and electron reduced MHD at small scales, previously considered in the literature. In the limit of ultra-local interactions, Leith-type nonlinear diffusion equations in the Fourier space are derived and heuristically extended to the strong turbulence regime by modifying the transfer time appropriately. Relations with existing phenomenological models for imbalanced MHD and balanced sub-ion turbulence are discussed. It turns out that in the presence of dispersive effects, the dynamics is sensitive on the way turbulence is maintained in a steady state. Furthermore, the total energy spectrum at sub-ion scales becomes steeper as the generalized cross-helicity flux is increased.

To explain the possible destabilization of a two-dimensional magnetic equilibrium such as the near-Earth magnetotail, we developed a kinetic model describing the resonant interaction of electromagnetic fluctuations and bouncing electrons trapped in the magnetic bottle. A small-$\unicode[STIX]{x1D6FD}$ approximation (i.e., the plasma pressure is lower than the magnetic pressure) is used in agreement with a small field line curvature. The linearized gyro-kinetic Vlasov equation is integrated along the unperturbed particle trajectories, including cyclotron and bounce motions. The dispersion relation for drift-Alfvèn waves is obtained through the plasma quasi-neutrality condition and Ampere’s law for the parallel current. It has been found that for a quasi-dipolar configuration ($\text{L}\sim 8$ corresponds to the set of the Earth’s magnetic field lines, crossing the Earths magnetic equator at 8 Earth radii), unstable electromagnetic modes may develop in the current sheet with a growth rate of the order of a few tenths of a second provided that the typical scale of density gradient slope responsible for the diamagnetic drift effects is over one Earth radius or less. This instability growth rate is large enough to destabilize the current sheet on time scales of 2–4 minutes as observed during substorm onset.

We have combined global magnetohydrodynamic (MHD) simulations of the solar wind and magnetosphere interaction with an implicit particle-in-cell simulation (PIC) and used this approach to model magnetic reconnection at both the dayside magnetopause and in the magnetotail plasma sheet. In this approach, we first model the magnetospheric configuration driven by the solar wind using the MHD simulation. At a time of interest (usually when a thin current sheet has formed in the MHD simulation), we load a large particle-in-cell simulation with plasma and fields based on the MHD state. We use the MHD results to set the boundary conditions on the PIC simulation. The coupling between the two models is one way – the PIC results do not change the MHD results. In these calculations, we use the UCLA global MHD code and the iPic3D implicit particle-in-cell code. In this paper we describe this technique in detail. As an example of this approach, we present PIC results on reconnection in the magnetotail during a magnetospheric substorm.

Stochastic heating refers to an increase in the average magnetic moment of charged particles interacting with electromagnetic fluctuations whose frequencies are smaller than the particles’ cyclotron frequencies. This type of heating arises when the amplitude of the gyroscale fluctuations exceeds a certain threshold, causing particle orbits in the plane perpendicular to the magnetic field to become stochastic rather than nearly periodic. We consider the stochastic heating of protons by Alfvén-wave (AW) and kinetic-Alfvén-wave (KAW) turbulence, which may make an important contribution to the heating of the solar wind. Using phenomenological arguments, we derive the stochastic-proton-heating rate in plasmas in which $\unicode[STIX]{x1D6FD}_{\text{p}}\sim 1$–30, where $\unicode[STIX]{x1D6FD}_{\text{p}}$ is the ratio of the proton pressure to the magnetic pressure. (We do not consider the $\unicode[STIX]{x1D6FD}_{\text{p}}\gtrsim 30$ regime, in which KAWs at the proton gyroscale become non-propagating.) We test our formula for the stochastic-heating rate by numerically tracking test-particle protons interacting with a spectrum of randomly phased AWs and KAWs. Previous studies have demonstrated that at $\unicode[STIX]{x1D6FD}_{\text{p}}\lesssim 1$, particles are energized primarily by time variations in the electrostatic potential and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the electrostatic potential. In contrast, at $\unicode[STIX]{x1D6FD}_{\text{p}}\gtrsim 1$, particles are energized primarily by the solenoidal component of the electric field and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the magnetic field.

Particle transport, acceleration and energization are phenomena of major importance for both space and laboratory plasmas. Despite years of study, an accurate theoretical description of these effects is still lacking. Validating models with self-consistent, kinetic simulations represents today a new challenge for the description of weakly collisional, turbulent plasmas. We perform simulations of steady state turbulence in the 2.5-dimensional approximation (three-dimensional fields that depend only on two-dimensional spatial directions). The chosen plasma parameters allow to span different systems, going from the solar corona to the solar wind, from the Earth’s magnetosheath to confinement devices. To describe the ion diffusion we adapted the nonlinear guiding centre (NLGC) theory to the two-dimensional case. Finally, we investigated the local influence of coherent structures on particle energization and acceleration: current sheets play an important role if the ions’ Larmor radii are of the order of the current sheet’s size. This resonance-like process leads to the violation of the magnetic moment conservation, eventually enhancing the velocity-space diffusion.

We consider the one-dimensional equilibrium problem of a shear-flow boundary layer within an ‘extended-fluid model’ of a plasma that includes the Hall and the electron pressure terms in Ohm’s law, as well as dynamic equations for anisotropic pressure for each species and first-order finite-Larmor-radius (FLR) corrections to the ion dynamics. We provide a generalized version of the analytic expressions for the equilibrium configuration given in Cerri et al., (Phys. Plasmas, vol. 20 (11), 2013, 112112), highlighting their intrinsic asymmetry due to the relative orientation of the magnetic field $\boldsymbol{B}$, $\boldsymbol{b}=\boldsymbol{B}/|\boldsymbol{B}|$, and the fluid vorticity $\unicode[STIX]{x1D74E}=\unicode[STIX]{x1D735}\times \boldsymbol{u}$ (‘$\unicode[STIX]{x1D74E}\boldsymbol{b}$ asymmetry’). Finally, we show that FLR effects can modify the Chapman–Ferraro current layer at the flank magnetopause in a way that is consistent with the observed structure reported by Haaland et al., (J. Geophys. Res. (Space Phys.), vol. 119, 2014, pp. 9019–9037). In particular, we are able to qualitatively reproduce the following key features: (i) the dusk–dawn asymmetry of the current layer, (ii) a double-peak feature in the current profiles and (iii) adjacent current sheets having thicknesses of several ion Larmor radii and with different current directions.

An analytic theory is presented of magnetic structures in collisionless, high-$\unicode[STIX]{x1D6FD}$ plasmas. Using a reduced magnetohydrodynamic model, a stationary nonlinear solution is constructed in the form of a Kelvin–Stuart cat’s eyes chain of magnetic islands, on the spatial scale that exceeds the characteristic ion lengths. The solution is imbedded in a background current sheet and possesses both a significant plasma density perturbation trapped inside the magnetic islands and a compressional magnetic field component that is driven mostly by a current loop located at the separatrix of the islands. This theory may provide an explanation for the magnetic structures observed in the solar wind close to the Earth by the Cluster spacecraft.

Properties of the electron mirror instability and its competition with the usually dominant whistler (electron cyclotron) instability driven by the electron perpendicular temperature anisotropy are investigated at the linear level using a Vlasov linear solver and at the nonlinear level using a two-dimensional full particle code. The simulation results show that the linearly subdominant electron mirror instability may compete at the nonlinear level with the whistler instability and may even eventually become the dominant mode that generates robust non-propagating sub-ion-scale coherent structures in the form of magnetic peaks.

Observational evidence in space and astrophysical plasmas with a long collisional mean free path suggests that more massive charged particles may be preferentially heated. One possible mechanism for this is the turbulent cascade of energy from injection to dissipation scales, where the energy is converted to heat. Here we consider a simple system consisting of a magnetized plasma slab of electrons and a single ion species with a cross-field density gradient. We show that such a system is subject to an electron drift wave instability, known as the universal instability, which is stabilized only when the electron and ion thermal speeds are equal. For unequal thermal speeds, we find from quasilinear analysis and nonlinear simulations that the instability gives rise to turbulent energy exchange between ions and electrons that acts to equalize the thermal speeds. Consequently, this turbulent heating tends to equalize the component temperatures of pair plasmas and to heat ions to much higher temperatures than electrons for conventional mass-ratio plasmas.

The nature of the cross-scale connections between the inertial-range turbulent energy cascade and the small-scale kinetic processes in collisionless plasmas is explored through the analysis of two-dimensional hybrid Vlasov–Maxwell numerical simulation (HVM), with $\unicode[STIX]{x1D6FC}$ particles, and through a proxy of the turbulent energy transfer rate, namely the local energy transfer (LET) rate. Correlations between pairs of variables, including those related to kinetic processes and to deviation from Maxwellian distributions, are first evidenced. Then, the general properties and the statistical scaling laws of the LET are described, confirming its reliability for the description of the turbulent cascade and revealing its textured topology. Finally, the connection between such proxy and the diagnostic variables is explored using conditional averaging, showing that several quantities are enhanced in the presence of large positive energy flux, and reduced near sites of negative flux. These observations can help in determining which processes are involved in the dissipation of energy at small scales, as for example the ion-cyclotron or mirror instabilities typically associated with perpendicular anisotropy of temperature.

The structure of a whistler precursor in a quasi-perpendicular shock is studied within two-fluid approach in one-dimensional case. The complete set of equations is reduced to the KdV equation, if no dissipation is included. With a phenomenological resistive dissipation the structure is described with the KdV–Burgers equation. The shock profile is intrinsically time dependent. For sufficiently strong dissipation, temporal evolution of a steepening profile results in generation of a stationary decaying whistler ahead of the shock front. With the decrease of the dissipation parameter, whistler wave trains begin to detach and propagate toward the upstream and the ramp is weakly time dependent. In the weakly dissipative regime the shock front experiences reformation.

Reconnection outflows have been under intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from an investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.

In this paper, weak-turbulence theory is used to investigate the nonlinear evolution of the parametric instability in three-dimensional low-$\unicode[STIX]{x1D6FD}$ plasmas at wavelengths much greater than the ion inertial length under the assumption that slow magnetosonic waves are strongly damped. It is shown analytically that the parametric instability leads to an inverse cascade of Alfvén wave quanta, and several exact solutions to the wave kinetic equations are presented. The main results of the paper concern the parametric decay of Alfvén waves that initially satisfy $e^{+}\gg e^{-}$, where $e^{+}$ and $e^{-}$ are the frequency ($f$) spectra of Alfvén waves propagating in opposite directions along the magnetic field lines. If $e^{+}$ initially has a peak frequency $f_{0}$ (at which $fe^{+}$ is maximized) and an ‘infrared’ scaling $f^{p}$ at smaller $f$ with $-1<p<1$, then $e^{+}$ acquires an $f^{-1}$ scaling throughout a range of frequencies that spreads out in both directions from $f_{0}$. At the same time, $e^{-}$ acquires an $f^{-2}$ scaling within this same frequency range. If the plasma parameters and infrared $e^{+}$ spectrum are chosen to match conditions in the fast solar wind at a heliocentric distance of 0.3 astronomical units (AU), then the nonlinear evolution of the parametric instability leads to an $e^{+}$ spectrum that matches fast-wind measurements from the Helios spacecraft at 0.3 AU, including the observed $f^{-1}$ scaling at $f\gtrsim 3\times 10^{-4}~\text{Hz}$. The results of this paper suggest that the $f^{-1}$ spectrum seen by Helios in the fast solar wind at $f\gtrsim 3\times 10^{-4}~\text{Hz}$ is produced in situ by parametric decay and that the $f^{-1}$ range of $e^{+}$ extends over an increasingly narrow range of frequencies as $r$ decreases below 0.3 AU. This prediction will be tested by measurements from the Parker Solar Probe.

Understanding the removal of energy from turbulent fluctuations in a magnetized plasma and the consequent energization of the constituent plasma particles is a major goal of heliophysics and astrophysics. Previous work has shown that nonlinear interactions among counterpropagating Alfvén waves – or Alfvén wave collisions – are the fundamental building block of astrophysical plasma turbulence and naturally generate current sheets in the strongly nonlinear limit. A nonlinear gyrokinetic simulation of a strong Alfvén wave collision is used to examine the damping of the electromagnetic fluctuations and the associated energization of particles that occurs in self-consistently generated current sheets. A simple model explains the flow of energy due to the collisionless damping and the associated particle energization, as well as the subsequent thermalization of the particle energy by collisions. The net particle energization by the parallel electric field is shown to be spatially localized, and the nonlinear evolution is essential in enabling spatial non-uniformity. Using the recently developed field–particle correlation technique, we show that particles resonant with the Alfvén waves in the simulation dominate the energy transfer, demonstrating conclusively that Landau damping plays a key role in the spatially localized damping of the electromagnetic fluctuations and consequent energization of the particles in this strongly nonlinear simulation.

The Kelvin–Helmholtz instability, proposed a long time ago for its role in and impact on the transport properties at magnetospheric flanks, has been widely investigated in the Earth’s magnetosphere context. This review covers more than fifty years of theoretical and numerical efforts in investigating the evolution of Kelvin–Helmholtz vortices and how the rich nonlinear dynamics they drive allow solar wind plasma bubbles to enter into the magnetosphere. Special care is devoted to pointing out the main advantages and weak points of the different plasma models that can be adopted for describing the collisionless magnetospheric medium and in underlying the important role of the three-dimensional geometry of the system.

To properly describe heating in weakly collisional turbulent plasmas such as the solar wind, interparticle collisions should be taken into account. Collisions can convert ordered energy into heat by means of irreversible relaxation towards the thermal equilibrium. Recently, Pezzi et al. (Phys. Rev. Lett., vol. 116, 2016a, 145001) showed that the plasma collisionality is enhanced by the presence of fine structures in velocity space. Here, the analysis is extended by directly comparing the effects of the fully nonlinear Landau operator and a linearized Landau operator. By focusing on the relaxation towards the equilibrium of an out of equilibrium distribution function in a homogeneous force-free plasma, here it is pointed out that it is significant to retain nonlinearities in the collisional operator to quantify the importance of collisional effects. Although the presence of several characteristic times associated with the dissipation of different phase space structures is recovered in both the cases of the nonlinear and the linearized operators, the influence of these times is different in the two cases. In the linearized operator case, the recovered characteristic times are systematically larger than in the fully nonlinear operator case, this suggesting that fine velocity structures are dissipated more slowly if nonlinearities are neglected in the collisional operator.

Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma models is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning approximately two decades in wavenumber – from the magnetohydrodynamics inertial range to scales well below the ion gyroradius – with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell simulations of freely decaying Alfvénic turbulence, whereas the other consists of Eulerian hybrid Vlasov–Maxwell simulations of turbulence continuously driven with partially compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_{\bot }\unicode[STIX]{x1D70C}_{i}\gtrsim 1$ are observed, where $k_{\bot }$ is the fluctuations’ wavenumber perpendicular to the background magnetic field and $\unicode[STIX]{x1D70C}_{i}$ is the ion Larmor radius. The system indeed self-consistently ‘reprocesses’ the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter ($\unicode[STIX]{x1D6FD}$ is the ratio between the thermal and the magnetic pressures). Small-scale turbulence has been found to be mainly populated by kinetic Alfvén wave (KAW) fluctuations for $\unicode[STIX]{x1D6FD}\geqslant 1$, whereas KAW fluctuations are only sub-dominant for low-$\unicode[STIX]{x1D6FD}$.

Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For a quasi-radial ambient magnetic field the expansion drives $T_{p\Vert }>T_{p\bot }$ and the system becomes eventually unstable with respect to the dominant parallel fire hose instability. This instability is generally unable to counteract the induced anisotropization and the system typically becomes unstable with respect to the oblique fire hose instability later on. The oblique instability efficiently reduces the anisotropy and the system rapidly stabilizes, while a significant part of the generated electromagnetic fluctuations is damped to protons. As long as the magnetic field is in the quasi-radial direction, this evolution repeats itself and the electromagnetic fluctuations accumulate. For a sufficiently oblique magnetic field the expansion drives $T_{p\bot }>T_{p\Vert }$ and brings the system to the stable region with respect to the fire hose instabilities.

Solar flares are powerful radiations occurring in the Sun’s atmosphere. They are powered by magnetic reconnection, a phenomenon that can convert magnetic energy into other forms of energy such as heat and kinetic energy, and which is believed to be ubiquitous in the universe. With the ever increasing spatial and temporal resolutions of solar observations, as well as numerical simulations benefiting from increasing computer power, we can now probe into the nature and the characteristics of magnetic reconnection in three dimensions to better understand the phenomenon’s consequences during eruptive flares in our star’s atmosphere. We review in the following the efforts made on different fronts to approach the problem of magnetic reconnection. In particular, we will see how understanding the magnetic topology in three dimensions helps in locating the most probable regions for reconnection to occur, how the current layer evolves in three dimensions and how reconnection leads to the formation of flux ropes, plasmoids and flaring loops.

This paper summarises some of the recent progress that has been made in understanding astrophysical plasma turbulence in the solar wind, from in situ spacecraft observations. At large scales, where the turbulence is predominantly Alfvénic, measurements of critical balance, residual energy and three-dimensional structure are discussed, along with comparison to recent models of strong Alfvénic turbulence. At these scales, a few per cent of the energy is also in compressive fluctuations, and their nature, anisotropy and relation to the Alfvénic component is described. In the small-scale kinetic range, below the ion gyroscale, the turbulence becomes predominantly kinetic Alfvén in nature, and measurements of the spectra, anisotropy and intermittency of this turbulence are discussed with respect to recent cascade models. One of the major remaining questions is how the turbulent energy is dissipated, and some recent work on this question, in addition to future space missions which will help to answer it, are briefly discussed.