The paper studies the influence of the background plasma density fluctuations on the dynamics of the Langmuir turbulence generated by electron beams, for parameters typical for solar type III beams and plasmas near 1 AU. A self-consistent Hamiltonian model based on the Zakharov and the Newton equations is used, which presents several advantages compared to the Vlasov approach. Beams generating Langmuir turbulence can be accelerated as a result of wave transformation effects or/and decay cascade processes; in both cases, the beam-driven Langmuir waves transfer part of their energy to waves of smaller wavenumbers, which can be reabsorbed later on by beam particles of higher velocities. As a consequence, beams can conserve a large part of their initial kinetic energy while propagating and radiating wave turbulence over long distances in inhomogeneous plasmas. Beam particles can also be accelerated in quasi-homogeneous plasmas due to the second cascade of wave decay, the wave transformation processes being very weak in this case. The net gains and losses of energy of a beam and the wave turbulence it radiates are calculated as a function of the average level of plasma density fluctuations and the beam parameters. The results obtained provide relevant information on the mechanism of energy reabsorption by beams radiating Langmuir turbulence in solar wind plasmas.

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two-dimensional ‘ABC’ structures) and zero-total-current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse, 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 plasma magnetization $\unicode[STIX]{x1D70E}\leqslant 10^{2}$ the spectrum power-law index is $p>2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\unicode[STIX]{x1D70E}\geqslant 10^{2}$, the spectra are hard, $p<2$, yet the maximal energy $\unicode[STIX]{x1D6FE}_{\text{max}}$ can still exceed the average magnetic energy per particle, ${\sim}\unicode[STIX]{x1D70E}$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

Accurate models of X-ray absorption and re-emission in partly stripped ions are necessary to calculate the structure of stars, the performance of hohlraums for inertial confinement fusion and many other systems in high-energy-density plasma physics. Despite theoretical progress, a persistent discrepancy exists with recent experiments at the Sandia Z facility studying iron in conditions characteristic of the solar radiative–convective transition region. The increased iron opacity measured at Z could help resolve a longstanding issue with the standard solar model, but requires a radical departure for opacity theory. To replicate the Z measurements, an opacity experiment has been designed for the National Facility (NIF). The design uses established techniques scaled to NIF. A laser-heated hohlraum will produce X-ray-heated uniform iron plasmas in local thermodynamic equilibrium (LTE) at temperatures ${\geqslant}150$ eV and electron densities ${\geqslant}7\times 10^{21}~\text{cm}^{-3}$. The iron will be probed using continuum X-rays emitted in a ${\sim}200$ ps, ${\sim}200~\unicode[STIX]{x03BC}\text{m}$ diameter source from a 2 mm diameter polystyrene (CH) capsule implosion. In this design, $2/3$ of the NIF beams deliver 500 kJ to the ${\sim}6$ mm diameter hohlraum, and the remaining $1/3$ directly drive the CH capsule with 200 kJ. Calculations indicate this capsule backlighter should outshine the iron sample, delivering a point-projection transmission opacity measurement to a time-integrated X-ray spectrometer viewing down the hohlraum axis. Preliminary experiments to develop the backlighter and hohlraum are underway, informing simulated measurements to guide the final design.

We propose a new semi-Lagrangian Vlasov–Poisson solver. It employs metric elements to follow locally the flow and its deformation, allowing one to find quickly and accurately the initial phase-space position $\boldsymbol{Q}(\boldsymbol{P})$ of any test particle $\boldsymbol{P}$ , by expanding at second order the geometry of the motion in the vicinity of the closest element. It is thus possible to reconstruct accurately the phase-space distribution function at any time $t$ and position $\boldsymbol{P}$ by proper interpolation of initial conditions, following Liouville theorem. When distortion of the elements of metric becomes too large, it is necessary to create new initial conditions along with isotropic elements and repeat the procedure again until next resampling. To speed up the process, interpolation of the phase-space distribution is performed at second order during the transport phase, while third-order splines are used at the moments of remapping. We also show how to compute accurately the region of influence of each element of metric with the proper percolation scheme. The algorithm is tested here in the framework of one-dimensional gravitational dynamics but is implemented in such a way that it can be extended easily to four- or six-dimensional phase space. It can also be trivially generalised to plasmas.

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}$.

Penetration of a magnetic field into plasma that is faster than resistive diffusion can be induced by the Hall electric field in a non-uniform plasma. This mechanism explained successfully the measured velocity of the magnetic field penetration into pulsed plasmas. Major related issues have not yet been resolved. Such is the theoretically predicted, but so far not verified experimentally, high magnetic energy dissipation, as well as the correlation between the directions of the density gradient and of the field penetration.

The energy budget of a collisionless plasma subject to electrostatic fluctuations is considered, and the excess of thermal energy over the minimum accessible to it under various constraints that limit the possible forms of plasma motion is calculated. This excess measures how much thermal energy is ‘available’ for conversion into plasma instabilities, and therefore constitutes a nonlinear measure of plasma stability. A distribution function with zero available energy defines a ‘ground state’ in the sense that its energy cannot decrease by any linear or nonlinear plasma motion. In a Vlasov plasma with small density and temperature fluctuations, the available energy is proportional to the mean square of these quantities, and exceeds the corresponding energy in ideal or resistive magnetohydrodynamics. If the first or second adiabatic invariant is conserved, ground states generally have inhomogeneous density and temperature. Magnetically confined plasmas are usually not in any ground state, but certain types of stellarator plasmas are so with respect to fluctuations that conserve both these adiabatic invariants, making the plasma linearly and nonlinearly stable to such fluctuations. Similar stability properties can also be enjoyed by plasmas confined by a dipole magnetic field.

The paper investigates the dynamics of magnetized non-neutral (electron) plasmas subjected to external electric field perturbations. A two-dimensional (2-D) particle-in-cell code is effectively exploited to model this system with a special attention to the role that non-axisymmetric, multipolar radio frequency (RF) drives applied to the cylindrical (circular) boundary play on the insurgence of azimuthal instabilities and the subsequent formation of coherent structures preventing the relaxation to a fully developed turbulent state, when the RF fields are chosen in the frequency range of the low-order fluid modes themselves. The isomorphism of such system with a 2-D inviscid incompressible fluid offers an insight into the details of forced 2-D fluid turbulence. The choice of different initial density (i.e. fluid vorticity) distributions allows for a selection of conditions where different levels of turbulence and intermittency are expected and a range of final states is achieved. Integral and spectral quantities of interest are computed along the flow using a multiresolution analysis based on a wavelet decomposition of both enstrophy and energy 2-D maps. The analysis of a variety of cases shows that the qualitative features of turbulent relaxation are similar in conditions of both free and forced evolution; at the same time, fine details of the flow beyond the self-similarity turbulence properties are highlighted in particular in the formation of structures and their timing, where the influence of the initial conditions and the effect of the external forcing can be distinguished.

A numerical investigation is carried out to understand the equilibrium $\unicode[STIX]{x1D6FD}$-limit in a classical stellarator. The stepped-pressure equilibrium code (Hudson et al., Phys. Plasmas, vol. 19 (11), 2012) is used in order to assess whether or not magnetic islands and stochastic field-lines can emerge at high $\unicode[STIX]{x1D6FD}$. Two modes of operation are considered: a zero-net-current stellarator and a fixed-iota stellarator. Despite the fact that relaxation is allowed (Taylor, Rev. Mod. Phys., vol. 58 (3), 1986, pp. 741–763), the former is shown to maintain good flux surfaces up to the equilibrium $\unicode[STIX]{x1D6FD}$-limit predicted by ideal-magnetohydrodynamics (MHD), above which a separatrix forms. The latter, which has no ideal equilibrium $\unicode[STIX]{x1D6FD}$-limit, is shown to develop regions of magnetic islands and chaos at sufficiently high $\unicode[STIX]{x1D6FD}$, thereby providing a ‘non-ideal $\unicode[STIX]{x1D6FD}$-limit’. Perhaps surprisingly, however, the value of $\unicode[STIX]{x1D6FD}$ at which the Shafranov shift of the axis reaches a fraction of the minor radius follows in all cases the scaling laws predicted by ideal-MHD. We compare our results to the High-Beta-Stellarator theory of Freidberg (Ideal MHD, 2014, Cambridge University Press) and derive a new prediction for the non-ideal equilibrium $\unicode[STIX]{x1D6FD}$-limit above which chaos emerges.

A long-standing problem in magnetic reconnection is to explain why it tends to proceed at or below a normalized rate of 0.1. This article gives a review of observational and numerical evidence for this rate and discusses recent theoretical work addressing this problem. Some remaining open questions are summarized.

The semi-Lagrangian discontinuous Galerkin method, coupled with a splitting approach in time, has recently been introduced for the Vlasov–Poisson equation. Since these methods are conservative, local in space and able to limit numerical diffusion, they are considered a promising alternative to more traditional semi-Lagrangian schemes. In this paper we study the conservation of important physical invariants and the long-time behaviour of the semi-Lagrangian discontinuous Galerkin method. To that end we conduct a theoretical analysis and perform a number of numerical simulations. In particular, we find that the entropy is non-decreasing for the discontinuous Galerkin scheme, while unphysical oscillations in the entropy are observed for the traditional method based on cubic spline interpolation.

The dynamics of a turbulent plasma not only manifests the transport of energy from large to small scales, but also can lead to a tangling of the magnetic field that threads through the plasma. The resulting magnetic field line wander can have a large impact on a number of other important processes, such as the propagation of energetic particles through the turbulent plasma. Here we explore the saturation of the turbulent cascade, the development of stochasticity due to turbulent tangling of the magnetic field lines and the separation of field lines through the turbulent dynamics using nonlinear gyrokinetic simulations of weakly collisional plasma turbulence, relevant to many turbulent space and astrophysical plasma environments. We determine the characteristic time $t_{2}$ for the saturation of the turbulent perpendicular magnetic energy spectrum. We find that the turbulent magnetic field becomes completely stochastic at time $t\lesssim t_{2}$ for strong turbulence, and at $t\gtrsim t_{2}$ for weak turbulence. However, when the nonlinearity parameter of the turbulence, a dimensionless measure of the amplitude of the turbulence, reaches a threshold value (within the regime of weak turbulence) the magnetic field stochasticity does not fully develop, at least within the evolution time interval $t_{2}<t\leqslant 13t_{2}$. Finally, we quantify the mean square displacement of magnetic field lines in the turbulent magnetic field with a functional form $\langle (\unicode[STIX]{x1D6FF}r)^{2}\rangle =A(z/L_{\Vert })^{p}$ ($L_{\Vert }$ is the correlation length parallel to the magnetic background field $\boldsymbol{B}_{\mathbf{0}}$, $z$ is the distance along $\boldsymbol{B}_{\mathbf{0}}$ direction), providing functional forms of the amplitude coefficient $A$ and power-law exponent $p$ as a function of the nonlinearity parameter.

Reduced fluid models for collisionless plasmas including electron inertia and finite Larmor radius corrections are derived for scales ranging from the ion to the electron gyroradii. Based either on pressure balance or on the incompressibility of the electron fluid, they respectively capture kinetic Alfvén waves (KAWs) or whistler waves (WWs), and can provide suitable tools for reconnection and turbulence studies. Both isothermal regimes and Landau fluid closures permitting anisotropic pressure fluctuations are considered. For small values of the electron beta parameter $\unicode[STIX]{x1D6FD}_{e}$ , a perturbative computation of the gyroviscous force valid at scales comparable to the electron inertial length is performed at order $O(\unicode[STIX]{x1D6FD}_{e})$ , which requires second-order contributions in a scale expansion. Comparisons with kinetic theory are performed in the linear regime. The spectrum of transverse magnetic fluctuations for strong and weak turbulence energy cascades is also phenomenologically predicted for both types of waves. In the case of moderate ion to electron temperature ratio, a new regime of KAW turbulence at scales smaller than the electron inertial length is obtained, where the magnetic energy spectrum decays like $k_{\bot }^{-13/3}$ , thus faster than the $k_{\bot }^{-11/3}$ spectrum of WW turbulence.

The properties of a non-relativistic magnetised low-beta electron–positron plasma in slab geometry are investigated. The two species are taken to be drift kinetic while we retain Larmor radius effects in quasi-neutrality, and inertia in Ohm’s law. A linear analysis shows that, for small magnetic perturbations, Alfvénic perturbations travel at the electron Alfvén speed, which is based on the electron mass. We discuss the role of the displacement current when Larmor-scale and Debye-scale effects are both retained. We predict the existence of a kinetic electron Alfvén wave which connects to the K-modes of Mishchenko et al. (J. Plasma Phys., 2017 (submitted)) in the electrostatic limit. It is found that linear drift waves are not supported by the system if the two species have the same temperature. Tearing modes can be driven unstable by equilibrium current density gradients. Also in this case, the characteristic time is based on the electron Alfvén speed. Nonlinear hybrid fluid-kinetic equations are also derived. It is shown that each species is described, to leading order, by the kinetic reduced electron heating model (KREHM) kinetic equation of Zocco & Schekochihin (Phys. Plasmas, vol. 18, 2011, 102309). The model is extended to retain first-order Larmor radius effects. It supports collisionless dispersive waves, which can greatly impact nonlinear magnetic reconnection. Diamagnetic effects enter the nonlinear equations via the first-order magnetic compressibility. A minimal nonlinear model for two-dimensional low-frequency isothermal pair plasmas is derived.

Aiming at a fusion reactor, two issues must be solved for the lower hybrid current drive (LHCD), namely good lower hybrid wave (LHW)–plasma coupling and effective current drive at high density. For this goal, efforts have been made to improve LHW–plasma coupling and current drive capability at high density in experimental advanced superconducting tokamak (EAST). LHW–plasma coupling is improved by means of local gas puffing and gas puffing from the electron side is taken as a routine way for EAST to operate with LHCD. Studies of high density experiments suggest that low recycling and high lower hybrid (LH) frequency are preferred for LHCD experiments at high density, consistent with previous results in other machines. With the combination of 2.45 GHz and 4.6 GHz LH waves, a repeatable high confinement mode plasma with maximum density up to $4.5\times 10^{19}~\text{m}^{-3}$ was obtained by LHCD in EAST. In addition, in the first stage of LHCD cyclic operation, an alternative candidate for more economical fusion reactors has been demonstrated in EAST and further work will be continued.

Two signatures of quantum effects on radiation reaction in the collision of a ${\sim}$GeV electron beam with a high intensity (${>}3\times 10^{20}~\text{W}~\text{cm}^{-2}$) laser pulse have been considered. We show that the decrease in the average energy of the electron beam may be used to measure the Gaunt factor $g$ for synchrotron emission. We derive an equation for the evolution of the variance in the energy of the electron beam in the quantum regime, i.e. quantum efficiency parameter $\unicode[STIX]{x1D702}\not \ll 1$. We show that the evolution of the variance may be used as a direct measure of the quantum stochasticity of the radiation reaction and determine the parameter regime where this is observable. For example, stochastic emission results in a 25 % increase in the standard deviation of the energy spectrum of a GeV electron beam, 1 fs after it collides with a laser pulse of intensity $10^{21}~\text{W}~\text{cm}^{-2}$. This effect should therefore be measurable using current high-intensity laser systems.

In traditional electron/ion laboratory plasmas, the system size $L$ is much larger than both the plasma skin depth $l_{s}$ and the Debye length $\unicode[STIX]{x1D706}_{D}$. In current and planned efforts to create electron/positron plasmas in the laboratory, this is not necessarily the case. A low-temperature, low-density system may have $\unicode[STIX]{x1D706}_{D}<L<l_{s}$; a high-density, thermally relativistic system may have $l_{s}<L<\unicode[STIX]{x1D706}_{D}$. Here we consider the question of what plasma physics phenomena are accessible (and/or diagnostically exploitable) in these different regimes and how this depends on magnetization. While particularly relevant to ongoing pair plasma creation experiments, the transition from single-particle behaviour to collective, ‘plasma’ effects – and how the criterion for that threshold is different for different phenomena – is an important but often neglected topic in electron/ion systems as well.

A brief critique is presented of some different classes of magnetohydrodynamic equilibrium solutions based on their continuity properties and whether the magnetic field is integrable or not. A generalized energy functional is introduced that is comprised of alternating ideal regions, with nested flux surfaces with an irrational rotational transform, and Taylor-relaxed regions, possibly with magnetic islands and chaos. The equilibrium states have globally continuous magnetic fields, and may be constructed for arbitrary three-dimensional plasma boundaries and appropriately prescribed pressure and rotational-transform profiles.

In this paper, we report a new type of instability of electron holes (EHs) interacting with passing ions. The nonlinear interaction of EHs and ions is investigated using a new theory of hole kinematics. It is shown that the oscillation in the velocity of the EH parallel to the magnetic field direction becomes unstable when the hole velocity in the ion frame is slower than a few times the cold ion sound speed. This instability leads to the emission of ion-acoustic waves from the solitary hole and decay in its magnitude. The instability mechanism can drive significant perturbations in the ion density. The instability threshold, oscillation frequency and instability growth rate derived from the theory yield quantitative agreement with the observations from a novel high-fidelity hole-tracking particle-in-cell code.

The time evolution of slow sound waves in a homogeneous, collisionless and quasineutral plasma, in particular their Landau damping, is investigated using the kinetic-magnetohydrodynamics formulation of Ramos (J. Plasma Phys. vol. 81, 2015 p. 905810325; vol. 82, 2016 p. 905820607). In this approach, the electric field is eliminated from a closed, hybrid fluid-kinetic system that ensures automatically the fulfilment of the charge neutrality condition. Considering the time dependence of a spatial-Fourier-mode linear perturbation with wavevector parallel to the equilibrium magnetic field, this can be cast as a second-order self-adjoint problem with a continuum spectrum of real and positive squared frequencies. Therefore, a conventional resolution of the identity with a continuum basis of singular normal modes is guaranteed, which simplifies significantly a Van Kampen-like treatment of the Landau damping. The explicit form of such singular normal modes is obtained, along with their orthogonality relations. These are used to derive the damped time evolution of the fluid moments of solutions of initial-value problems, for the most general kinds of initial conditions. The non-zero parallel electric field is not used explicitly in this analysis, but it is calculated from any given solution after the later has been obtained.