Plasma supports collective modes and particle–wave interactions that lead to complex behaviour in, for example, inertial fusion energy applications. While plasma can sometimes be modelled as a charged fluid, a kinetic description is often crucial for studying nonlinear effects in the higher-dimensional momentum–position phase space that describes the full complexity of the plasma dynamics. We create a differentiable solver for the three-dimensional partial-differential equation describing the plasma kinetics and introduce a domain-specific objective function. Using this framework, we perform gradient-based optimization of neural networks that provide forcing function parameters to the differentiable solver given a set of initial conditions. We apply this to an inertial-fusion-relevant configuration and find that the optimization process exploits a novel physical effect.

We study the onset of electron heating in intense laser–solid interactions and its impact on the spectral quality of radiation pressure accelerated ions in both hole boring and light sail regimes. Two- and three-dimensional particle-in-cell (PIC) simulations are performed over a wide range of laser and target parameters and reveal how the pulse duration, profile, polarization and target surface stability control the electron heating, the dominant ion acceleration mechanisms and the ion spectra. We find that the onset of strong electron heating is associated with the growth of the Rayleigh–Taylor-like instability at the front surface and must be controlled to produce high-quality ion beams, even when circularly polarized lasers are employed. We define a threshold condition for the maximum duration of the laser pulse that allows mitigation of electron heating and radiation pressure acceleration of narrow energy spread ion beams. The model is validated by three-dimensional PIC simulations, and the few experimental studies that reported low energy spread radiation pressure accelerated ion beams appear to meet the derived criteria. The understanding provided by our work will be important in guiding future experimental developments, for example for the ultrashort laser pulses becoming available at state-of-the-art laser facilities, for which we predict that proton beams with $\sim$150–250 MeV, $\sim$30% energy spread, and a total laser-to-proton conversion efficiency of $\sim$20% can be produced.

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al., Phys. Plasmas, vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$, catalysed by $\boldsymbol {B}_{n-n'}$. It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

A new formalism for the derivation of the cnoidal wave solution is presented by introducing a new set of initial conditions and a subacoustic moving frame. The resulting set of solutions is also constructed in a way that ensures the conservation of number density. The solutions are illustrated in a variety of graphical representations, and the effect of the amplitude's magnitude on the cnoidal wave form is presented. Interestingly, it is shown that the wavelength of the solutions decreases with amplitude. In addition, it is shown that the small-amplitude cnoidal wave solutions converge to linear waves in the small-amplitude limit.

The Saffman helicity invariant of Hosking & Schekochihin (Phys. Rev. X, vol. 11, issue 4, 2021, 041005), which we here call the Hosking integral, has emerged as an important quantity that may govern the decay properties of magnetically dominated non-helical turbulence. Using a range of different computational methods, we confirm that this quantity is indeed gauge invariant and nearly perfectly conserved in the limit of large Lundquist numbers. For direct numerical simulations with ordinary viscosity and magnetic diffusivity operators, we find that the solution develops in a nearly self-similar fashion. In a diagram quantifying the instantaneous decay coefficients of magnetic energy and integral scale, we find that the solution evolves along a line that is indeed suggestive of the governing role of the Hosking integral. The solution settles near a line in this diagram that is expected for a self-similar evolution of the magnetic energy spectrum. The solution will settle in a slightly different position when the magnetic diffusivity decreases with time, which would be compatible with the decay being governed by the reconnection time scale rather than the Alfvén time.

The energy functional of ideal magnetohydrodynamics (MHD) is extended to the regime of finite resistivity and transformed into an intuitive form, which separates pressure-gradient, current-density drive and stabilising phenomena, analogously to the ideal version of the functional derived by Greene and Johnson (Plasma Phys., vol. 10, issue 8, 1968, pp. 729–745). The resulting resistive energy functional is not in the form of a solvable eigenvalue problem but is intended to analyse the energetic composition of eigenfunctions which were calculated in advance by solving the MHD equations. The new resistive energy functional provides information on the different energetic drives and stabilising phenomena of an eigenfunction. The new energy functional has been implemented in the CASTOR3D code in general 3D curvilinear coordinates and successfully validated across three different coordinate systems for a numerical test case.

Vertical displacement normal modes in shaped tokamak plasmas are studied analytically, based on the reduced ideal-magnetohydrodynamic model. With the help of quadratic forms, and using the appropriate eigenfunction for vertical displacements with toroidal mode number $n=0$ and dominant elliptical-angle mode number $m=1$, a dispersion relation is derived, including the effects of ideal or resistive walls through a single parameter, $D_w(\gamma )$, which is, in general, a function of the complex eigenfrequency $\gamma = -{\rm i}\omega$. For the resistive-wall case, the dispersion relation is cubic in $\gamma$. One root corresponds to the well-known, non-rotating resistive-wall vertical mode, growing on the resistive-wall time scale. The other two roots are weakly damped by wall resistivity, but oscillate with a frequency below the poloidal Alfvén frequency, which makes them immune to continuum damping, but subject to possible instability due to resonant interaction with fast ions.

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

The Vlasov–Fokker–Planck (VFP) equation, a variant of the Boltzmann equation, is frequently used to model the dynamics of laboratory and astrophysical plasmas. A common approach in solving the VFP equation in numerical and analytic studies is to expand the momentum part of the distribution function $f$ (i.e. the particle density in phase-space) in terms of known functions. The Cartesian tensor and spherical harmonic expansion have been widely used, leading to the question of how to convert between the different coefficients of the two expansions. This problem is also familiar in multipole expansions of an electrostatic (or gravitational) potential. The coefficients of the Cartesian tensor expansion of the potential are called (Cartesian) multipole moments and the ones of the spherical harmonic expansion are called spherical multipole moments. In this paper, we investigate the relation between the two kinds of multipole moments and provide a general formalism to convert between them. We subsequently apply this formalism to the coefficients of the expansions of the distribution function $f$. A free, open-source command-line tool which implements this formalism is provided.

A wave turbulence theory is developed for inertial electron magnetohydrodynamics (IEMHD) in the presence of a relatively strong and uniform external magnetic field $\boldsymbol {B_0} = B_0 \hat {\boldsymbol {e}}_\|$. This regime is relevant for scales smaller than the electron inertial length $d_e$. We derive the kinetic equations that describe the three-wave interactions between inertial whistler or kinetic Alfvén waves. We show that for both invariants, energy and momentum, the transfer is anisotropic (axisymmetric) with a direct cascade mainly in the direction perpendicular ($\perp$) to $\boldsymbol {B_0}$. The exact stationary solutions (Kolmogorov–Zakharov spectra) are obtained for which we prove the locality. We also found the Kolmogorov constant $C_K \simeq 8.474$. In the simplest case, the study reveals an energy spectrum in $k_\perp ^{-5/2} k_\|^{-1/2}$ (with k the wavenumber) and a momentum spectrum enslaved to the energy dynamics in $k_\perp ^{-3/2} k_\|^{-1/2}$. These solutions correspond to a magnetic energy spectrum ${\sim }k_\perp ^{-9/2}$, which is steeper than the EMHD prediction made for scales larger than $d_e$. We conclude with a discussion on the application of the theory to space plasmas.

New orbit-averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity-space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS (Velasco et al., J. Comput. Phys., vol. 418, 2020, 109512; Velasco et al., Nucl. Fusion, vol. 61, 2021, 116013). The equations are used to study the $1/\nu$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/\nu$ regime transitions directly into the $\nu$ regime, without passing through a $\sqrt {\nu }$ regime. An explicit formula for the neoclassical fluxes in the $\nu$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrt {\nu }$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by Calvo et al. (Plasma Phys. Control. Fusion, vol. 59, 2017, 055014). In these stellarators, the $\sqrt {\nu }$ regime exists in the collisionality interval $(a/R) \ln (R/a) \ll \nu _{ii} R a/\rho _i v_{ti} \ll R/a$. In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrt {\nu }$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $\delta$: $\delta ^{2} |\ln \delta | \ll \nu _{ii} R a/\rho _i v_{ti} \ll 1$. Here, $\nu _{ii}$ is the ion–ion collision frequency, $\rho _i$ and $v_{ti}$ are the ion gyroradius and thermal speed, and $a$ and $R$ are the minor and major radii.

We extend our previous work on the two-dimensional (2-D) Dimits transition in ion-scale turbulence (Ivanov et al., J. Plasma Phys., vol. 86, 2020, 855860502) to include variations along the magnetic field. We consider a three-field fluid model for the perturbations of electrostatic potential, ion temperature, and ion parallel flow in a constant-magnetic-curvature geometry without magnetic shear. It is derived in the cold-ion, long-wavelength asymptotic limit of the gyrokinetic theory. Just as in the 2-D model, a low-transport (Dimits) regime exists and is found to be dominated by a quasistatic staircase-like arrangement of strong zonal flows and zonal temperature. This zonal staircase is formed and maintained by a negative turbulent viscosity for the zonal flows. Unlike the 2-D model, the three-dimensional (3-D) one does not suffer from an unphysical blow up beyond the Dimits threshold where the staircase becomes nonlinearly unstable. Instead, a well-defined finite-amplitude saturated state is established. This qualitative difference between the 2-D and 3-D models is due to the appearance of small-scale ‘parasitic’ modes that exist only if we allow perturbations to vary along the magnetic field lines. These modes extract energy from the large-scale perturbations and provide an effective enhancement of large-scale thermal diffusion, thus aiding the energy transfer from large injection scales to small dissipative ones. We show that in our model, the parasitic modes always favour a zonal-flow-dominated state. In fact, a Dimits state with a zonal staircase is achieved regardless of the strength of the linear drive, provided the system is sufficiently extended along the magnetic field and sufficient parallel resolution is provided.

The goal of this work is to generate large statistically representative data sets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student $t$ process regression. We apply Student $t$ process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via colouring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics and classic machine learning clustering algorithms.

We develop the formalism of the first-order near-axis expansion of the magnetohydrodynamic equilibrium equations described by Garren & Boozer (Phys. Fluids B, vol. 3, issue 10, 1991, pp. 2805–2821) and Plunk et al. (J. Plasma Phys., vol. 85, issue 6, 2019; J. Plasma Phys., vol. 87, issue 6, 2021) for the case of a quasi-isodynamic, $N$-field-period, stellarator-symmetric, single-well magnetic field equilibrium. The importance of the magnetic axis shape is investigated, and we conclude that control of the curvature and torsion is crucial to obtain omnigenous configurations with finite aspect ratio and low effective ripple, especially for a higher number of field periods. For this reason a method is derived to construct classes of axis shapes with favourable curvature and torsion. Solutions are presented, including a three-field-period configuration constructed at an aspect ratio of $A=20$, with a maximum elongation of $e=3.2$ and an effective ripple under $1\,\%$, which demonstrates that high elongation is not a necessary feature of quasi-isodynamic stellarators.

A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configurations defined by a set of adjustable magnetic field parameters. This optimization, completed in 30 s on a single CPU core using the SIMSOPT code, yields a solution with excellent confinement, as measured by the conventional figure of merit for neoclassical transport, effective ripple, at a modest aspect ratio of eight. The optimization parameters that led to this configuration are described, its confinement properties are assessed and a set of magnetic field coils is found. The resulting transport at low collisionality is much smaller than that of W7-X, and the device needs significantly fewer coils because of the reduced number of field periods.

In a magnetised plasma on scales well above ion kinetic scales, any constant-magnitude magnetic field, accompanied by parallel Alfvénic flows, forms a nonlinear solution in an isobaric, constant-density background. These structures, which are also known as spherically polarised Alfvén waves, are observed ubiquitously in the solar wind, presumably created by the growth of small-amplitude fluctuations as they propagate outwards in the corona. Here, we present a computational method to construct such solutions of arbitrary amplitude with general multidimensional structure, and explore some of their properties. The difficulty lies in computing a zero-divergence, constant-magnitude magnetic field, which leaves a single, quasi-free function to define the solution, while requiring strong constraints on any individual component of the field. Motivated by the physical process of wave growth in the solar wind, our method circumvents this issue by starting from low-amplitude Alfvénic fluctuations dominated by a strong mean field, then ‘growing’ magnetic perturbations into the large-amplitude regime. We present example solutions with non-trivial structure in one, two and three dimensions, demonstrating a clear tendency to form very sharp gradients or discontinuities, unless the solution is one-dimensional. As well as being useful as an input for other calculations, particularly the study of parametric decay, our results provide a natural explanation for the extremely sharp field discontinuities observed across magnetic field switchbacks in the low solar wind.

We examine the turbulence driven by the ion and electron temperature gradients in selected magnetic configurations of the Wendelstein 7-X (W7-X) stellarator. The inherent flexibility in the configuration space of W7-X enables us to find candidate configurations manifesting low turbulent transport. We follow insights gained by stellarator optimization techniques, in order to identify key geometric features, which are directly related to the ion and electron heat fluxes produced by plasma turbulence. One such a feature is the flux expansion at locations where the curvature is particularly unfavourable. Starting from a configuration routinely used in the W7-X experiment, we end up with an optimized configuration. Based on this equilibrium, we select a configuration from W7-X configuration database with a similar feature as the optimized one. With the help of nonlinear gyrokinetic simulations, we show that the heat flux in this configuration is less stiff than in the initial configuration, both for ion temperature gradient and electron temperature gradient turbulence.

Weakly collisional plasmas are subject to nonlinear relaxation processes, which can operate at rates much faster than the particle collision frequencies. This causes the plasma to respond like a magnetised fluid despite having long particle mean free paths. In this Letter the effective collisional mechanisms are modelled in the plasma kinetic equation to produce density, pressure and magnetic-field responses to compare with spacecraft measurements of the solar wind compressive fluctuations at 1 AU. This enables a measurement of the effective mean free path of the solar wind protons, found to be ${\approx }4 \times 10^{5}$ km, which is approximately $10^{3}$ times shorter than the collisional mean free path. These measurements are shown to support the effective fluid behaviour of the solar wind at scales above the proton gyroradius and demonstrate that effective collision processes alter the thermodynamics and transport of weakly collisional plasmas.

Ion heating in collisionless shocks is non-adiabatic and efficient. The amount of heating and the downstream distributions depend on the shock parameters and on the incident ion distribution. The number of reflected ions and their distribution depend on the detailed shape of the tail of the distribution. In supercritical shocks the reflected ion contribution is significant. Kappa distributed ions are heated more strongly and have a larger fraction of reflected ions than Maxwellian distributed ions with the same upstream temperature and the same shock parameters. For kappa distributions the phase space dips are shallower.

The stability properties of kinetic ballooning modes are investigated for a number of magnetic configurations of the stellarator Wendelstein 7-X. In particular, we consider the effects of the vacuum rotational transform, ${\raise.1pt-\kern-6pt\iota}$, and the mirror ratio. The analysis sheds light on the interplay between global magnetohydrodynamic configuration properties and local gyrokinetic stability, and is instrumental in the design of high-$\beta$ (the ratio of kinetic to magnetic pressure) operation scenarios. In particular, it is demonstrated that some Wendelstein 7-X magnetic configurations have a relatively low kinetic ballooning mode threshold.