We present a novel framework for finite element particle-in-cell methods based on the discretization of the underlying Hamiltonian structure of the Vlasov–Maxwell system. We derive a semi-discrete Poisson bracket, which retains the defining properties of a bracket, anti-symmetry and the Jacobi identity, as well as conservation of its Casimir invariants, implying that the semi-discrete system is still a Hamiltonian system. In order to obtain a fully discrete Poisson integrator, the semi-discrete bracket is used in conjunction with Hamiltonian splitting methods for integration in time. Techniques from finite element exterior calculus ensure conservation of the divergence of the magnetic field and Gauss’ law as well as stability of the field solver. The resulting methods are gauge invariant, feature exact charge conservation and show excellent long-time energy and momentum behaviour. Due to the generality of our framework, these conservation properties are guaranteed independently of a particular choice of the finite element basis, as long as the corresponding finite element spaces satisfy certain compatibility conditions.

We experimentally study the excitation and propagation of acoustic solitary waves in a one-dimensional dusty plasma (i.e. a Yukawa chain) with $n=65$ particles interacting through a screened Coulomb potential. The lattice constant $a=1.02\pm 0.02$ mm. Waves are launched by applying a 100 mW laser pulse to one end of the chain for laser pulse durations from 0.10 to 2.0 s. We observe damped solitary waves which propagate for distances ${\gtrsim}30a$ with an acoustic speed $c_{s}=11.5\pm 0.2~\text{mm}~\text{s}^{-1}$. The maximum velocity perturbation increases with laser pulse duration for durations ${\leqslant}0.5$ s and then saturates at ${\approx}15\,\%$. The wave speed is found to be independent of the maximum amplitude, indicating that the formation of nonlinear solitons is prevented by neutral-gas damping.

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.

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.

Five-dimensional gyrokinetic continuum simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using a full-$f$ discontinuous-Galerkin approach implemented in the Gkeyll code. While various simplifications have been used for now, such as long-wavelength approximations in the gyrokinetic Poisson equation and the Hamiltonian, these simulations include the basic elements of a fusion-device scrape-off layer: localised sources to model plasma outflow from the core, cross-field turbulent transport, parallel flow along magnetic field lines, and parallel losses at the limiter or divertor with sheath-model boundary conditions. The set of sheath-model boundary conditions used in the model allows currents to flow through the walls. In addition to details of the numerical approach, results from numerical simulations of turbulence in the Large Plasma Device, a linear device featuring straight magnetic field lines, are presented.

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.

Although, for current laser pulse energies, the weakly nonlinear regime of laser wakefield acceleration is known to be the optimal for reaching the highest possible electron energies, the capabilities of upcoming large laser systems will provide the possibility of running highly nonlinear regimes of laser pulse propagation in underdense or near-critical plasmas. Using an extended particle-in-cell (PIC) model that takes into account all the relevant physics, we show that such regimes can be implemented with external guiding for a relatively long distance of propagation and allow for the stable transformation of laser energy into other types of energy, including the kinetic energy of a large number of high energy electrons and their incoherent emission of photons. This is despite the fact that the high intensity of the laser pulse triggers a number of new mechanisms of energy depletion, which we investigate systematically.

Relativistic electron beam transport through a high-density, magnetized plasma is studied numerically and theoretically. An electron beam injected into a cold plasma excites Weibel and two-stream instabilities that heat the beam and saturate. In the absence of an applied magnetic field, the heated beam continues to propagate. However, when a magnetic field of particular strength is applied along the direction of beam propagation, a secondary instability of off-angle whistler modes is excited. These modes then couple nonlinearly creating a large amplitude parallel-propagating whistler that stops the beam. Here, we will show these phenomena in detail and explain the mechanism of whistler mediated beam stagnation.

The multiple colliding laser pulse concept formulated by Bulanov et al. (Phys. Rev. Lett., vol. 104, 2010b, 220404) is beneficial for achieving an extremely high amplitude of coherent electromagnetic field. Since the topology of electric and magnetic fields of multiple colliding laser pulses oscillating in time is far from trivial and the radiation friction effects are significant in the high field limit, the dynamics of charged particles interacting with the multiple colliding laser pulses demonstrates remarkable features corresponding to random walk trajectories, limit circles, attractors, regular patterns and Lévy flights. Under extremely high intensity conditions the nonlinear dissipation mechanism stabilizes the particle motion resulting in the charged particle trajectory being located within narrow regions and in the occurrence of a new class of regular patterns made by the particle ensembles.

High-frequency photons travelling in plasma exhibit a linear polarizability that can influence the dispersion of linear plasma waves. We present a detailed calculation of this effect for Langmuir waves as a characteristic example. Two alternative formulations are given. In the first formulation, we calculate the modified dispersion of Langmuir waves by solving the governing equations for the electron fluid, where the photon contribution enters as a ponderomotive force. In the second formulation, we provide a derivation based on the photon polarizability. Then, the calculation of ponderomotive forces is not needed, and the result is more general.

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.

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.

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.

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.

The Weibel instability, driven by a temperature anisotropy, is investigated within different kinetic descriptions based on the semi-Lagrangian full kinetic and relativistic Vlasov–Maxwell model, on the multi-stream approach, which is based on a Hamiltonian reduction technique, and finally, with the full pressure tensor fluid-type description. Dispersion relations of the Weibel instability are derived using the three different models. A qualitatively different regime is observed in Vlasov numerical experiments depending on the excitation of a longitudinal plasma electric field driven initially by the combined action of the stream symmetry breaking and weak relativistic effects, in contrast with the existing theories of the Weibel instability based on their purely transverse characters. The multi-stream model offers an alternate way to simulate easily the coupling with the longitudinal electric field and particularly the nonlinear regime of saturation, making numerical experiments more tractable, when only a few moments of the distribution are considered. Thus a numerical comparison between the reduced Hamiltonian model (the multi-stream model) and full kinetic (relativistic) Vlasov simulations has been investigated in that regime. Although nonlinear simulations of the fluid model, including the dynamics of the pressure tensor, have not been carried out here, the model is strongly relevant even in the three-dimensional case.

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.

Magnetic fields play an important role in astrophysical accretion discs and in the interstellar and intergalactic medium. They drive jets, suppress fragmentation in star-forming clouds and can have a significant impact on the accretion rate of stars. However, the exact amplification mechanisms of cosmic magnetic fields remain relatively poorly understood. Here, I start by reviewing recent advances in the numerical and theoretical modelling of the turbulent dynamo, which may explain the origin of galactic and intergalactic magnetic fields. While dynamo action was previously investigated in great detail for incompressible plasmas, I here place particular emphasis on highly compressible astrophysical plasmas, which are characterised by strong density fluctuations and shocks, such as the interstellar medium. I find that dynamo action works not only in subsonic plasmas, but also in highly supersonic, compressible plasmas, as well as for low and high magnetic Prandtl numbers. I further present new numerical simulations from which I determine the growth of the turbulent (un-ordered) magnetic field component ($B_{turb}$) in the presence of weak and strong guide fields ($B_{0}$). I vary $B_{0}$ over five orders of magnitude and find that the dependence of $B_{turb}$ on $B_{0}$ is relatively weak, and can be explained with a simple theoretical model in which the turbulence provides the energy to amplify $B_{turb}$. Finally, I discuss some important implications of magnetic fields for the structure of accretion discs, the launching of jets and the star-formation rate of interstellar clouds.

The existence of a potential energy functional in the zero-Larmor-radius collisionless plasma theory of Kruskal & Oberman (Phys. Fluids, vol. 1, 1958 p. 275), Rosenbluth & Rostoker (Phys. Fluids, vol. 2, 1959, p. 23) allows us to derive easily sufficient conditions for linear stability. However, this kinetic magnetohydrodynamics (KMHD) theory does not have a self-adjointness property, making it difficult to derive necessary conditions. In particular, the standard methods to prove that an instability follows if some trial perturbation makes the incremental potential energy negative, which rely on the self-adjointness of the force operator or on the existence of a complete basis of normal modes, are not applicable to KMHD. This paper investigates KMHD linear stability criteria based on the time evolution of initial-value solutions, without recourse to the classic bounds or comparison theorems of Kruskal–Oberman and Rosenbluth–Rostoker for the KMHD potential energy. The adopted approach does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Most importantly, the investigation of a necessary condition for stability does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes. It is thereby shown that stability in isothermal ideal-MHD is a sufficient condition for stability in KMHD and that, with a proviso on the long-time behaviour of oscillations about stable equilibria, stability in the double-adiabatic fluid theory, including the variation of the parallel fluid displacement, would be a necessary condition for stability in KMHD.

The self-organisation of turbulence into regular zonal flows can be fruitfully investigated with quasi-linear methods and statistical descriptions. A wave-kinetic equation that assumes asymptotically large-scale zonal flows leads to ultraviolet divergence. From an exact description of quasi-linear dynamics emerges two better geometrical optics approximations. These involve not only the mean flow shear but also the second and third derivative of the mean flow. One approximation takes the form of a new wave-kinetic equation, but is only valid when the zonal flow is quasi-static and wave action is conserved.

Dust acoustic (or dust density) waves have been observed in many laboratory dusty plasmas. These low-frequency waves involve the dynamics of highly charged and massive dust grains, and can be excited by the flow of ions relative to dust. In this paper, we consider the nonlinear development of the dust acoustic instability, excited by thermal ion flow, in a collisional plasma containing dust with high kinetic temperature (warm dust). It is shown that under certain conditions there may be a long-wavelength secondary instability in the nonlinear stage as dust gets heated by the waves. The characteristics of the nonlinear development are considered as a function of the relative charge density of the dust. Application to possible experimental parameters is discussed.