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.

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.

The basic physics of the processes playing the most important role in divertor plasma detachment is reviewed. The models used in the two-dimensional edge plasma transport codes that are widely used to address different issues of the edge plasma physics and to simulate the experimental data, as well as the numerical schemes and convergence issues, are described. The processes leading to ultimate divertor plasma detachment, the transition to and the stability of the detached regime, as well as the impact of the magnetic configuration and divertor geometry on detachment, are considered. A consistent, integral physical picture of ultimate detachment of a tokamak divertor plasma is developed.

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.

We compare, in an extensive and systematic way, linear theory results obtained with the hybrid (ion kinetic and electron fluid), the gyrokinetic and the fully kinetic plasma models. We present a test case with parameters that are relevant for solar wind turbulence at small scales, which is a topic now recognized to need a kinetic treatment, to a certain extent. We comment on the comparison of low-frequency single modes (Alfvén/ion-cyclotron, ion-acoustic and fast modes) for a wide range of propagation angles, and on the overall spectral properties of the linear operators, for quasi-perpendicular propagation. The methodology and the results presented in this paper will be valuable when choosing which model should be used in regimes where the assumptions of each model are not trivially satisfied.

The results obtained by the plasma physics community for the validation and the prediction of turbulence and transport in magnetized plasmas come mainly from the use of very central processing unit (CPU)-consuming particle-in-cell or (gyro)kinetic codes which naturally include non-Maxwellian kinetic effects. To date, fluid codes are not considered to be relevant for the description of these kinetic effects. Here, after revisiting the limitations of the current fluid theory developed in the 19th century, we generalize the fluid theory including kinetic effects such as non-Maxwellian super-thermal tails with as few fluid equations as possible. The collisionless and collisional fluid closures from the nonlinear Landau Fokker–Planck collision operator are shown for an arbitrary collisionality. Indeed, the first fluid models associated with two examples of collisionless fluid closures are obtained by assuming an analytic non-Maxwellian distribution function (e.g. the INMDF (Izacard, O. 2016b Kinetic corrections from analytic non-Maxwellian distribution functions in magnetized plasmas. Phys. Plasmas 23, 082504) that stands for interpreted non-Maxwellian distribution function). One of the main differences with the literature is our analytic representation of the distribution function in the velocity phase space with as few hidden variables as possible thanks to the use of non-orthogonal basis sets. These new non-Maxwellian fluid equations could initiate the next generation of fluid codes including kinetic effects and can be expanded to other scientific disciplines such as astrophysics, condensed matter or hydrodynamics. As a validation test, we perform a numerical simulation based on a minimal reduced INMDF fluid model. The result of this test is the discovery of the origin of particle and heat diffusion. The diffusion is due to the competition between a growing INMDF on short time scales due to spatial gradients and the thermalization on longer time scales. The results shown here could provide the insights to break some of the unsolved puzzles of turbulence.

We investigate the question of how plasma currents circulate and close in the scrape-off-layer (SOL) of convection-limited tokamak plasmas. A simplified two-fluid model describes how currents must evacuate charge at the sheaths due to cross-field currents that are not divergence-free. These include turbulence-driven polarization currents and poloidally asymmetric equilibrium diamagnetic currents. The theory provides an estimate for the radial profile of the floating potential, which reveals a dipolar structure like the one observed experimentally. Simulations with a fluid turbulence code provide evidence for the predicted behaviour of currents and floating potential.

The long-standing challenge to describing charged particle dynamics in strong classical electromagnetic fields is how to incorporate classical radiation, classical radiation reaction and quantized photon emission into a consistent unified framework. The current, semiclassical methods to describe the dynamics of quantum particles in strong classical fields also provide the theoretical framework for fundamental questions in gravity and hadron–hadron collisions, including Hawking radiation, cosmological particle production and thermalization of particles created in heavy-ion collisions. However, as we show, these methods break down for highly relativistic particles propagating in strong fields. They must therefore be improved and adapted for the description of laser–plasma experiments that typically involve the acceleration of electrons. Theory developed from quantum electrodynamics, together with dedicated experimental efforts, offer the best controllable context to establish a robust, experimentally validated foundation for the fundamental theory of quantum effects in strong classical potentials.

Comparisons are made of two different particle simulations: one for the study of plasma-based accelerators (Gueroult & Fisch, Phys. Plasmas, vol. 23, 2016, 032113) and the other for the study of shock formation in the interstellar medium (Yamauchi & Ohsawa, Phys. Plasmas, vol. 14, 2007, 053110). In the former, shock waves used for plasma density control create ion beams by reflection. In the latter, a fast and dense beam of exploding ions penetrates a surrounding plasma. In both simulations, magnetic bumps are generated from the motion of ion beams perpendicular to a magnetic field. Despite the apparent differences of their purposes, configurations and spatial scales, the two simulations show strong similarities in the generation processes and effects of the bumps, suggesting that these are not rare plasma phenomena. The bump created by the exploding ions develops into backward and forward magnetosonic pulses.

This is a tutorial-style selective review explaining basic concepts of forced magnetic reconnection. It is based on a celebrated model of forced reconnection suggested by J. B. Taylor. The standard magnetohydrodynamic (MHD) theory of this process has been pioneered by Hahm & Kulsrud (Phys. Fluids, vol. 28, 1985, p. 2412). Here we also discuss several more recent developments related to this problem. These include energetics of forced reconnection, its Hall-mediated regime, and nonlinear effects with the associated onset of the secondary tearing (plasmoid) instability.

Determining the physical mechanisms that extract energy from turbulent fluctuations in weakly collisional magnetized plasmas is necessary for a more complete characterization of the behaviour of a variety of space and astrophysical plasmas. Such a determination is complicated by the complex nature of the turbulence as well as observational constraints, chiefly that in situ measurements of such plasmas are typically only available at a single point in space. Recent work has shown that correlations between electric fields and particle velocity distributions constructed from single-point measurements produce a velocity-dependent signature of the collisionless damping mechanism. We extend this work by constructing field–particle correlations using data sets drawn from single points in strongly driven, turbulent, electromagnetic gyrokinetic simulations to demonstrate that this technique can identify the collisionless mechanisms operating in such systems. The velocity-space structure of the correlation between proton distributions and parallel electric fields agrees with expectations of resonant mechanisms transferring energy collisionlessly in turbulent systems. This work motivates the eventual application of field–particle correlations to spacecraft measurements in the solar wind, with the ultimate goal to determine the physical mechanisms that dissipate magnetized plasma turbulence.

While electron kinetic effects are well known to be of fundamental importance in several situations, the electron mean-flow inertia is often neglected when length scales below the electron skin depth become irrelevant. This has led to the formulation of different reduced models, where electron inertia terms are discarded while retaining some or all kinetic effects. Upon considering general full-orbit particle trajectories, this paper compares the dispersion relations emerging from such models in the case of the Weibel instability. As a result, the question of how length scales below the electron skin depth can be neglected in a kinetic treatment emerges as an unsolved problem, since all current theories suffer from drawbacks of different nature. Alternatively, we discuss fully kinetic theories that remove all these drawbacks by restricting to frequencies well below the plasma frequency of both ions and electrons. By giving up on the length scale restrictions appearing in previous works, these models are obtained by assuming quasi-neutrality in the full Vlasov–Maxwell system.

Measurements of local plasma parameters in dusty plasma are crucial for understanding the physics issues related to such systems. The Langmuir probe, a small electrode immersed in the plasma, provides such measurements. However, designing of a Langmuir probe system in a dusty plasma environment demands special consideration. First, the probe has to be miniaturized enough so that its perturbation on the ambient dust structure is minimal. At the same time, the probe dimensions must be such that a well-defined theory exists for interpretation of its characteristics. The associated instrumentation must also support the measurement of current collected by the probe with high signal to noise ratio. The most important consideration, of course, comes from the fact that the probes are prone to dust contamination, as the dust particles tend to stick to the probe surface and alter the current collecting area in unpredictable ways. This article describes the design and operation of a Langmuir probe system that resolves these challenging issues in dusty plasma. In doing so, first, different theories that are used to interpret the probe characteristics in collisionless as well as in collisional regimes are discussed, with special emphasis on application. The critical issues associated with the current–voltage characteristics of Langmuir probe obtained in different operating regimes are discussed. Then, an algorithm for processing these characteristics efficiently in presence of ion-neutral collisions in the probe sheath is presented.

Subsequent to the recent rigorous derivation of an energetically consistent gyrokinetic collision operator in the so-called Landau representation, this paper investigates the possibility of finding a differential formulation of the gyrokinetic Landau collision operator. It is observed that, while a differential formulation is possible in the gyrokinetic phase space, reduction of the resulting system of partial differential equations to five dimensions via gyroaveraging poses a challenge. Based on the present work, it is likely that the gyrocentre analogues of the Rosenbluth–MacDonald–Judd potential functions must be kept gyroangle dependent.

We review multiwavelength properties of pulsar wind nebulae created by supersonically moving pulsars and the effects of pulsar motion on the pulsar wind nebulae morphologies and the ambient medium. Supersonic pulsar wind nebulae are characterized by bow-shaped shocks around the pulsar and/or cometary tails filled with the shocked pulsar wind. In the past several years significant advances in supersonic pulsar wind nebula studies have been made in deep observations with the Chandra and XMM-Newton X-ray observatories and the Hubble Space Telescope. In particular, these observations have revealed very diverse supersonic pulsar wind nebula morphologies in the pulsar vicinity, different spectral behaviours of long pulsar tails, the presence of puzzling outflows misaligned with the pulsar velocity and far-UV bow shocks. Here we review the current observational status focusing on recent developments and their implications.

In this article, self-focusing of an intense circularly polarized laser pulse in the presence of an external oblique magnetic field in hot magnetized plasma, using Maxwell’s equations and the relativistic fluid momentum equation, is studied. An envelope equation governing the spot size of the laser beam for both of left- and right-hand polarizations has been derived and the effects of the plasma temperature and oblique magnetic field on the electron density distribution of hot plasma with respect to variation of the normalized laser spot size has been investigated. Numerical results depict that in right-hand polarization, self-focusing of the laser pulse along the propagation direction in hot magnetized plasma becomes better and more compressed with increasing $\unicode[STIX]{x1D703}$. Inversely, in left-hand polarization, increase of $\unicode[STIX]{x1D703}$ in an oblique magnetic field leads to enhancement of the spot size and reduction self-focusing. Besides, in the plasma density profile, self-focusing of the laser pulse improves in comparison with no oblique magnetic field. Also it is shown that plasma temperature has a key role in the laser spot size, normalized laser output power and the variation of plasma density.

The extreme properties of the gamma-ray flares in the Crab nebula present a clear challenge to our ideas on the nature of particle acceleration in relativistic astrophysical plasma. It seems highly unlikely that standard mechanisms of stochastic type are at work here and hence the attention of theorists has switched to linear acceleration in magnetic reconnection events. In this series of papers, we attempt to develop a theory of explosive magnetic reconnection in highly magnetized relativistic plasma which can explain the extreme parameters of the Crab flares. In the first paper, we focus on the properties of the X-point collapse. Using analytical and numerical methods (fluid and particle-in-cell simulations) we extend Syrovatsky’s classical model of such collapse to the relativistic regime. We find that the collapse can lead to the reconnection rate approaching the speed of light on macroscopic scales. During the collapse, the plasma particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\unicode[STIX]{x1D70E}$. For sufficiently high magnetizations and vanishing guide field, the non-thermal particle spectrum consists of two components: a low-energy population with soft spectrum that dominates the number census; and a high-energy population with hard spectrum that possesses all the properties needed to explain the Crab flares.

The paper is devoted to the shapes of the solar wind ion flux fluctuation spectrum at the transition between the inertial and the kinetic range using in situ high-resolution measurements of the Russian mission Spektr-R. We analyse the variability of the transition region and select five typical types of spectral shapes: (i) spectra with two slopes and one break, (ii) spectra characterized by a nonlinear steepening in the kinetic range, (iii) spectra with flattening in the vicinity of the break, (iv) spectra with a bump in the vicinity of the break and (v) spectra without any steepening in the kinetic range. The most popular is the well-known type (i) observed in approximately half of the cases. The second most popular type of spectra is type (iii) occurring in approximately one third of the cases. The other three types are observed less often: type (ii) – in approximately 6 %; type (iv) in 3 % and type (v) in 6 % of cases. An analysis of typical plasma conditions for different types of spectra revealed that the last two type of spectra (iv) and (v) are generally observed in a very slow solar wind with a low proton density, (i) and (iii) are observed in the solar wind with rather typical conditions and (ii) is usually observed in high-speed streams. The effect of nonlinear steepening of the spectra in the kinetic range increases with the solar wind speed. We present also the analysis of statistical properties of the observed events and compare them with the predictions of several statistical turbulence models. We show that intermittency is always observed in the solar wind flow despite the presence of one or another shape of spectra. The log-Poisson model with a dominant contribution of filament-like structures shows the best parameterization of the experimentally observed scaling.

One-dimensional self-consistent kinetic models may describe some states of intense charged particle beams. In the collisionless approximation, which is appropriate for the short current pulse duration, the kinetic distribution function may be built as a function of the motion integrals. Two situations are considered corresponding to the wide charged particle flux with the sharp front and to the sheet continuous beam freely propagating in space. In both cases, one-dimensional Vlasov equation solutions are shown to exist, which are based on algebraic functions of the motion invariants.