This paper discusses thermal equilibrium states of single-species plasmas, such as pure electron plasmas and pure positron plasmas, that are confined in a dipole trap. Thermal equilibrium states for such plasmas are routinely realized in the homogeneous magnetic field of Penning–Malmberg traps. We generalize the theory of these states to include inhomogeneous magnetic dipole fields. The approach to thermal equilibrium takes place in two stages with well separated time scales. On the collision time scale, thermal equilibrium is established along each magnetic field line. On the much longer transport time scale, heat conduction and viscosity bring the plasmas on different flux contours into thermal equilibrium, we call this a state of global thermal equilibrium. We present numerical results for local and global thermal equilibria. These results agree with the analytic predictions for charge collections that are large compared with the Debye length. There is, in principle, no limit to the confinement time of a single-species plasma in a global thermal equilibrium state. Experiments with hot electron–ion plasmas performed in the LDX and RT1 devices give us confidence that, in contrast to a Penning–Malmberg trap, a magnetic dipole field can also confine cold quasi-neutral electron–positron pair plasmas on the time scale of the phenomena of interest. Such pair plasmas are assumed to form in the magnetosphere of neutron stars but have so far not been realized in a laboratory. The creation of an electron–positron pair plasma is the main goal of the APEX collaboration.

In tokamaks and neoclassically optimised stellarators, like Wendelstein 7-X (W7-X) and the Helically Symmetric Experiment, turbulent transport is expected to be the dominant transport mechanism. Among the electrostatic instabilities that drive turbulence, the trapped-electron mode (TEM) has been shown both analytically and in simulations to be absent over large ranges of parameter space in quasi-isodynamic stellarator configurations with the maximum-$J$ property. It has been proposed that the reduction of the linear TEM growth rate in such configurations may lead to the passing-electron-driven universal instability, which is often subdominant to the TEM, becoming the fastest-growing instability over some range of parameter space. Here, we show through gyrokinetic simulations using the Gene code, that the universal instability is dominant in a variety of stellarator geometries over a range of parameter space typically occupied by the TEM, but most consequentially in devices which possess beneficial TEM stability properties like W7-X, which locally satisfies the maximum-$J$ property for deeply trapped particles in the regions of worst curvature. We find that the universal instability exists at long perpendicular wavelengths and, as a result, dominates the potential fluctuation amplitude in nonlinear simulations. In W7-X, universal modes are found to differ in parallel mode structure from trapped-particle modes, which may impact turbulence localisation in experiments.

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

We summarize recent attempts to unravel the role of plasma kinetic effects in radiation mediated shocks. Such shocks form in all strong stellar explosions and are responsible for the early electromagnetic emission released from these events. A key issue that has been overlooked in all previous works is the nature of the coupling between the charged leptons, that mediate the radiation force, and the ions, which are the dominant carriers of the shock energy. Our preliminary investigation indicates that in the case of relativistic shocks, as well as Newtonian shocks in multi-ion plasma, this coupling is driven by either, transverse magnetic fields of a sufficiently magnetized upstream medium, or plasma microturbulence if strong enough magnetic fields are absent. We discuss the implications for the shock breakout signal, as well as abundance evolution and kilonova emission in binary neutron star mergers.

We upgrade the three-dimensional linear stability CASTOR3D code by taking into account gyro-viscosity in the momentum equation, as well as pressure and Hall terms in Ohm's law. Making use of the drift approximation, the resulting simplified version of Ohm's law is implemented. The ion and electron diamagnetic drift effects described by these newly implemented terms are investigated for ideal and resistive, low- and medium-$n$ internal modes depending on the ratio of ion and electron pressures, and plasma resistivity. Stability studies are performed for two- and three-dimensional test equilibria. The results confirm previous findings, and provide new insights especially in the case of the three-dimensional equilibrium (e.g. reduction of the coupling of toroidal harmonics).

The oscillating relativistic electrons in the accelerating/focusing wakefields of plasma accelerators emit electromagnetic radiation known as betatron radiation (BR). The proton-driven plasma wakefield acceleration has been demonstrated in the Advanced Wakefield Experiment (AWAKE) at CERN; however, its accompanying radiation emission is less explored compared with those in the laser- and electron beam-driven plasma accelerators. In this paper, a detailed simulation study of BR in the AWAKE is presented. Considering the new set-up of the AWAKE Run 2 (2021–), the radiation emission from both the witness electron beam and the seeding electron beam is investigated using particle-in-cell simulations. The influence of radial size mismatch and misaligned off-axis injection on the witness beam dynamics, as well as the spectral features of the relevant BR are studied. These non-ideal electron injections are likely to occur in experiment. The proton self-modulation stage is also investigated with a close look at the seeding electron beam dynamics and its BR. As a footprint of the emitting particles, BR can provide valuable information about the beam dynamics. Some practical challenges to implement the betatron diagnostic in the AWAKE Run 2 experiment are also addressed.

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by dHybridR, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise a new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field–particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. Together, the field–particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation.

Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are computationally expensive. Here, we present a one-dimensional semi-analytical model for the radial displacement of ablated material after pellet injection, taking into account both the Alfvén and ohmic currents which shortcircuit the charge separation creating the drift. The model is suitable for rapid calculation of the radial drift displacement, and can be useful for e.g. modelling of disruption mitigation via pellet injection.

We propose a theoretical scheme for the study of the nonlinear interaction of drift-wave-like turbulence and energetic particles in stellarators. The approach is based on gyrokinetics, and features a separation of time and scales, for electromagnetic fluctuations, inspired by linear ballooning theory. Two specific moments of the gyrokinetic equation constitute the main equations of the system, which requires a full kinetic nonlinear solution. This is found iteratively, expanding in the smallness of the bounce-average radial drift frequency, and nonlinear $\boldsymbol {E}\times \boldsymbol {B}$ drift frequency, compared with the inverse time scales of the resonantly interacting energetic particles. Our analysis is therefore valid for neoclassically optimised stellators. The resummation of all iterative and perturbative nonlinear kinetic solutions is discussed in terms of Feynman diagrams. Particular emphasis is put on the role of collisionlessly undamped large-scale structures in phase space, the kinetic equivalent of zonal flows, i.e. phase-space zonal structures, and on wave-like fluctuations generated by energetic particles.

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

The growing amount of data produced by simulations and observations of space physics processes encourages the use of methods rooted in machine learning for data analysis and physical discovery. We apply a clustering method based on self-organizing maps to fully kinetic simulations of plasmoid instability, with the aim of assessing their suitability as a reliable analysis tool for both simulated and observed data. We obtain clusters that map well, a posteriori, to our knowledge of the process; the clusters clearly identify the inflow region, the inner plasmoid region, the separatrices and regions associated with plasmoid merging. Self-organizing map-specific analysis tools, such as feature maps and the unified distance matrix, provide us with valuable insights into both the physics at work and specific spatial regions of interest. The method appears as a promising option for the analysis of data, both from simulations and from observations, and could also potentially be used to trigger the switch to different simulation models or resolution in coupled codes for space simulations.

We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the ion poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes.

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-$\beta$ collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma $\beta$ is increased. For mode amplitudes $|\delta B/B_0| \gtrsim 0.3$, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies $|\delta B/B_0| \gtrsim 2\beta ^{-1}$. The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-$\beta$ plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

A new perturbation and continuation method is presented for computing and analysing stellarator equilibria. The method is formally derived from a series expansion about the equilibrium condition $\boldsymbol {F} \equiv \boldsymbol {J}\times \boldsymbol {B} - \boldsymbol {\nabla } p = 0$, and an efficient algorithm for computing solutions to second- and third-order perturbations is developed. The method has been implemented in the DESC stellarator equilibrium code, using automatic differentiation to compute the required derivatives. Examples are shown demonstrating its use for computing complicated equilibria, perturbing a tokamak into a stellarator and performing parameter scans in pressure, rotational transform and boundary shape in a fraction of the time required for a full solution.

Global electromagnetic turbulence is simulated in stellarator geometry using the gyrokinetic particle-in-cell code EUTERPE. The evolution of the turbulent electromagnetic field and the plasma profiles is considered at different values of the plasma beta and for different magnetic configurations. It is found that turbulence is linearly driven at relatively high toroidal mode numbers. In the nonlinear regime, lower toroidal mode numbers, including zonal flows, are excited resulting in a quench of the linear instability drive. The turbulent heat flux is outward and leads to the nonlinear relaxation of the plasma temperature profile. The particle flux is inward for the parameters considered. The effect of the parallel perturbation of the magnetic field on the stellarator turbulence is addressed.

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

The impurity transport at Wendelstein 7-X during its most recent campaign is characterized and documented for a variety of different plasma scenarios. An overview of its dependence on several quantities is given, which allows identification of transport regimes and the major driver for impurity transport. Beyond this, a comparison with the impurity behavior in other fusion devices is now possible. In contrast to other stellarators, no density dependence of the impurity transport has been found. Additionally, the influence of the turbulence contribution to the overall transport is reflected in the dependence on various parameters, e.g. turbulent diffusion and density fluctuation amplitudes. With this database approach, one can now also apply scaling laws to make extrapolations about the impurity confinement in future plasma scenarios.

Good fast ion confinement in quasi-isodynamic stellarators still has to be demonstrated experimentally. In the absence of fusion alphas, auxiliary heating systems can be used to generate fast ions. For this purpose, W7-X has neutral beam injection and its newly installed ion cyclotron heating system. In practice, both systems can run simultaneously, the synergetic running of both is the focus of the current modelling development. With the upgrades to the SCENIC code described here, we find that, compared with pure neutral beam injection (NBI), the fast ion tail is significantly enhanced when radio frequency (RF) heating is enabled. But the effectiveness of heating and fast ion tail generation strongly depends on the plasma parameters and antenna frequency, which calls for further study, both numerical and experimental.

Starting from the equations of collisionless linear gyrokinetics for magnetised plasmas with an imposed inhomogeneous magnetic field, we present the first known analytical, closed-form solution for the resulting velocity-space integrals in the presence of resonances due to both parallel streaming and constant magnetic drifts. These integrals are written in terms of the well-known plasma dispersion function (Faddeeva & Terent'ev, Tables of Values of the Function $w(z)=\exp (-z^2)(1+2i/\sqrt {\pi }\int _0^z \exp (t^2) \,\mathrm {d} t)$ for Complex Argument, 1954. Gostekhizdat. English translation: Pergamon Press, 1961; Fried & Conte, The Plasma Dispersion Function, 1961. Academic Press), rendering the subsequent expressions simpler to treat analytically and more efficient to compute numerically. We demonstrate that our results converge to the well-known ones in the straight-magnetic-field and two-dimensional limits, and show good agreement with the numerical solver by Gürcan (J. Comput. Phys., vol. 269, 2014, p. 156). By way of example, we calculate the exact dispersion relation for a simple electrostatic, ion-temperature-gradient-driven instability, and compare it with approximate kinetic and fluid models.