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.

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.

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.

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.

The first nonlinear gyrokinetic simulations obtained using a moment approach based on the Hermite–Laguerre decomposition of the distribution function are presented, implementing advanced models for the collision operator. Turbulence in a two-dimensional Z-pinch is considered within a flux-tube configuration. In the collisionless regime, our gyromoment approach shows very good agreement with nonlinear simulations carried out with the continuum gyrokinetic code GENE, even with fewer gyromoments than required for the convergence of the linear growth rate. By using advanced linear collision operators, the role of collisions in setting the level of turbulent transport is then analysed. The choice of collision operator model is shown to have a crucial impact when turbulence is quenched by the presence of zonal flows. The convergence properties of the gyromoment approach improve when collisions are included.

How much does the cross-section of a toroidal magnetic field configuration tell us about its magnetohydrodynamic (MHD) stability? It is generally believed that positive triangularity (typically leading to bean-shaped cross-sections with their indentation on the inboard side in stellarators) contributes positively to MHD stability. In this paper, we explore the basis of this statement within a near-axis description for axisymmetric and quasisymmetric magnetic configurations. In agreement with the existing literature, we show that positive triangularity stabilises vertically elongated tokamaks. In quasisymmetric stellarators, the toroidal asymmetry of flux surfaces modifies this relation. The behaviour of stellarator-symmetric, quasisymmetric stellarators can still be described in terms of the shape of one of their up–down symmetric cross-sections. However, we show that for a sample of quasisymmetric configurations, the positive-bean-shaped cross-sections do not contribute positively to stability. Unlike in the axisymmetric case, we also learn that finite $\beta$ can improve stability even without magnetic shear.

Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

A high-intensity laser pulse propagating through a gas target disturbs the uniform plasma distribution. Plasma density structures, created by high-order Bessel–Gauss beams for guiding the accelerating Gaussian beam and laser wakefield acceleration of electrons, are analysed using Wake-T and Fourier–Bessel particle-in-cell (FBPIC) simulation tools. The use of Bessel–Gauss doughnut beams increases the acceleration distance and energy of accelerated electrons up to 2.3 times at a 2 mm distance relative to the Gaussian beam of the same intensity.

We study the parametric decay instability of parallel-propagating Alfvén waves in a low-beta plasma using one-dimensional fully kinetic simulations. We focus for the first time on the conversion of the energy stored in the initial Alfvén wave into particle internal energy, and on its partition between particle species. We show that compressible fluctuations generated by the decay of the pump wave into a secondary ion-acoustic mode and a reflected Alfvén wave contribute to the gain of internal energy via two distinct mechanisms. First, the ion-acoustic mode leads nonlinearly to proton trapping and proton phase-space mixing, in agreement with previous work based on hybrid simulations. Second, during the nonlinear stage, a compressible front of the fast type develops at the steepened edge of the backward Alfvén wave leading to a field-aligned proton beam propagating backwards at the Alfvén speed. We find that parametric decay heats preferentially protons, which gain approximately 50 % of the pump wave energy in the form of internal energy. However, we find that electrons are also energized and that they contribute to the total energy balance by gaining 10 % of the pump wave energy. By investigating energy partition and particle heating during parametric decay, our results contribute to the determination of the role of compressible and kinetic effects in wave-driven models of the solar wind.