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.

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.

The ideal Chew–Goldberger–Low (CGL) plasma equations, including the double adiabatic conservation laws for the parallel ($p_\parallel$) and perpendicular pressure ($p_\perp$), are investigated using a Lagrangian variational principle. An Euler–Poincaré variational principle is developed and the non-canonical Poisson bracket is obtained, in which the non-canonical variables consist of the mass flux ${\boldsymbol {M}}$, the density $\rho$, the entropy variable $\sigma =\rho S$ and the magnetic induction ${\boldsymbol {B}}$. Conservation laws of the CGL plasma equations are derived via Noether's theorem. The Galilean group leads to conservation of energy, momentum, centre of mass and angular momentum. Cross-helicity conservation arises from a fluid relabelling symmetry, and is local or non-local depending on whether the gradient of $S$ is perpendicular to ${\boldsymbol {B}}$ or otherwise. The point Lie symmetries of the CGL system are shown to comprise the Galilean transformations and scalings.

This paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the $\smash {H}$-theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence (Phys. Fluids, vol. 16, 1973, p. 1102) is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

Dissipative processes cause collisionless plasmas in many systems to develop non-thermal particle distributions with broad power-law tails. The prevalence of power-law energy distributions in space/astrophysical observations and kinetic simulations of systems with a variety of acceleration and trapping (or escape) mechanisms poses a deep mystery. We consider the possibility that such distributions can be modelled from maximum-entropy principles, when accounting for generalizations beyond the Boltzmann–Gibbs entropy. Using a dimensional representation of entropy (related to the Renyi and Tsallis entropies), we derive generalized maximum-entropy distributions with a power-law tail determined by the characteristic energy scale at which irreversible dissipation occurs. By assuming that particles are typically energized by an amount comparable to the free energy (per particle) before equilibrating, we derive a formula for the power-law index as a function of plasma parameters for magnetic dissipation in systems with sufficiently complex topologies. The model reproduces several results from kinetic simulations of relativistic turbulence and magnetic reconnection.

Following indirect-drive experiments which demonstrated promising performance for low convergence ratios (below 17), previous direct-drive simulations identified a fusion-relevant regime which is expected to be robust to hydrodynamic instability growth. This paper expands these results with simulated implosions at lower energies of 100 and 270 kJ, and ‘hydrodynamic equivalent’ capsules which demonstrate comparable convergence ratio, implosion velocity and in-flight aspect ratio without the need for cryogenic cooling, which would allow the assumptions of one-dimensional-like performance to be tested on current facilities. A range of techniques to improve performance within this regime are then investigated, including the use of two-colour and deep ultraviolet laser pulses. Finally, further simulations demonstrate that the deposition of electron energy into the hotspot of a low convergence ratio implosion through auxiliary heating also leads to significant increases in yield. Results include break even for 1.1 MJ of total energy input (including an estimated 370 kJ of short-pulse laser energy to produce electron beams for the auxiliary heating), but are found to be highly dependent upon the efficiency with which electron beams can be created and transported to the hotspot to drive the heating mechanism.

An optimized compact stellarator with four simple coils is obtained from direct optimization via a coil shape. The new stellarator consists of two interlocking coils and two vertical field coils similar to those of the Columbia Non-neutral Torus (CNT) (Pedersen et al., Phys. Rev. Lett., vol. 88, 2002, pp. 205002). The optimized configuration has a global magnetic well and a low helical ripple level comparable to that of Wendelstein 7-X (W7-X) (Wolf et al., Nucl. Fusion, vol. 57, 2017, pp. 102020). The two interlocking coils have a smooth three-dimensional shape much simpler than those of advanced stellarators such as W7-X. This result opens up possibilities of future stellarator reactors with simplified coils.

We study the time-dependent formation and evolution of a current sheet (CS) in a magnetised, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field increases, and adiabatic invariance in the inflowing fluid elements produces a field-biased pressure anisotropy with excess perpendicular pressure. At large plasma beta, this anisotropy excites the mirror instability, which deforms the reconnecting field on ion-Larmor scales and dramatically reduces the effective thickness of the CS. Tearing modes whose wavelengths are comparable to that of the mirrors then become unstable, triggering reconnection on smaller scales and at earlier times than would have occurred if the thinning CS were to have retained its Harris profile. A novel method for identifying and tracking X-points is introduced, yielding X-point separations that are initially intermediate between the perpendicular and parallel mirror wavelengths in the upstream plasma. These mirror-stimulated tearing modes ultimately grow and merge to produce island widths comparable to the CS thickness, an outcome we verify across a range of CS formation timescales and initial CS widths. Our results may find their most immediate application in the tearing disruption of magnetic folds generated by turbulent dynamo in weakly collisional, high-beta, astrophysical plasmas.

We employ the $\lambda ^{3}$ regime where a near-single-cycle laser pulse is tightly focused, thus providing the highest possible intensity for the minimal energy at a certain laser power. The quantum electrodynamics processes in the course of the interaction of an ultra-intense laser with a solid target are studied via three-dimensional particle-in-cell simulations, revealing the generation of copious $\gamma$-photons and electron–positron pairs. A parametric study of the laser polarisation, target thickness and electron number density shows that a radially polarised laser provides the optimal regime for $\gamma$-photon generation. By varying the laser power in the range of 1 to 300 PW we find the scaling of the laser to $\gamma$-photon energy conversion efficiency. The laser-generated $\gamma$-photon interaction with a high-$Z$ target is further studied using Monte Carlo simulations revealing further electron–positron pair generation and radioactive nuclide creation.

The paper presents experimental results from the SMOLA device that is the first facility with a helical mirror section of the magnetic system. This device was built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea that is the technique of an active control of axial losses from a confinement zone. Theory predicts that, with rotating plasma, a helical mirror will provide suppression of the axial plasma flow and, simultaneously, density pinching to the axis. Experiments demonstrated the increase in plasma density in the entrance trap by a factor of 1.6 in the helical configuration. The integral axial flux from the transport section drops severalfold. The effective mirror ratio of the helical section was $R_{eff} > 10$. Particle flux returning by the helical mirror section towards the confinement zone was observed. At high corrugation ratios, the axial flux direction is different at the magnetic axis and in the periphery of the plasma in the helical section. All axial fluxes scale linearly with the plasma density, even if the ion mean free path is comparable to the total length of the helical section. Good agreement of the experimental results with theoretical predictions is found.

We show that large-amplitude, non-planar, Alfvén-wave (AW) packets are exact nonlinear solutions of the relativistic magnetohydrodynamic equations when the total magnetic-field strength in the local fluid rest frame ($b$) is a constant. We derive analytic expressions relating the components of the fluctuating velocity and magnetic field. We also show that these constant-$b$ AWs propagate without distortion at the relativistic Alfvén velocity and never steepen into shocks. These findings and the observed abundance of large-amplitude, constant-$b$ AWs in the solar wind suggest that such waves may be present in relativistic outflows around compact astrophysical objects.

The GENE-3D code, the global stellarator version of the established GENE framework, has been extended to an electromagnetic gyrokinetic code. This paper outlines the basic structure of the algorithm, highlighting the treatment of the electromagnetic terms. The numerical implementation is verified against the radially global GENE code in linear and nonlinear tokamak simulations, recovering excellent agreement between both codes. As a first application to stellarator plasmas, linear and nonlinear global simulations with kinetic electrons of ion temperature gradient (ITG) turbulence in Wendelstein 7-X were performed, showing a decrease of ITG activity through the introduction of electromagnetic effects via a finite plasma-$\beta$. The upgrade makes it possible to study a large variety of new physical scenarios, including kinetic electron and electromagnetic effects, reducing the gap between gyrokinetic models and physically realistic systems.

Radio frequency (RF) sheaths occur under a wide variety of conditions when RF waves, material surfaces and plasma coexist. RF sheaths are of special importance in describing the interaction of ion cyclotron range of frequency (ICRF) waves with the boundary plasma in tokamaks, stellarators and other magnetic confinement devices. In this article the basic physics of RF sheaths is discussed in the context of magnetic fusion research. Techniques for modelling RF sheaths, their interaction with RF wave fields and the resulting consequences are highlighted. The article is intended as a guide for the early-career ICRF researcher, but it may equally well serve to provide an overview of basic RF sheath concepts and modelling directions for any interested fusion scientist.

Electron energization during merging of magnetized plasmas is studied using the OMEGA and OMEGA EP laser facilities by colliding two plasma plumes, each containing a Biermann-battery self-generated magnetic field. Two neighbouring plasma plumes are produced by intense laser beams, and the anti-parallel Biermann fields merge and reconnect in the process of the plumes’ expansion and collision. To isolate the merging as an acceleration source, the electron energy spectra obtained from two-plume collision shots are compared with the spectra from single-plume shots. Single-plume shots exhibit an energized electron tail with energies up to ${\sim }250\ \textrm {keV}$. The electrons in merging experiments are additionally accelerated by ${\sim }50\text {--}100$ keV compared to single-plume shots.

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the $\boldsymbol {E\times B}$ rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the $\boldsymbol {E\times B}$ velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

Between the base of the solar corona at $r=r_\textrm {b}$ and the Alfvén critical point at $r=r_\textrm {A}$, where $r$ is heliocentric distance, the solar-wind density decreases by a factor $ \mathop > \limits_\sim 10^5$, but the plasma temperature varies by a factor of only a few. In this paper, I show that such quasi-isothermal evolution out to $r=r_\textrm {A}$ is a generic property of outflows powered by reflection-driven Alfvén-wave (AW) turbulence, in which outward-propagating AWs partially reflect, and counter-propagating AWs interact to produce a cascade of fluctuation energy to small scales, which leads to turbulent heating. Approximating the sub-Alfvénic region as isothermal, I first present a brief, simplified calculation showing that in a solar or stellar wind powered by AW turbulence with minimal conductive losses, $\dot {M} \simeq P_\textrm {AW}(r_\textrm {b})/v_\textrm {esc}^2$, $U_{\infty } \simeq v_\textrm {esc}$, and $T\simeq m_\textrm {p} v_\textrm {esc}^2/[8 k_\textrm {B} \ln (v_\textrm {esc}/\delta v_\textrm {b})]$, where $\dot {M}$ is the mass outflow rate, $U_{\infty }$ is the asymptotic wind speed, $T$ is the coronal temperature, $v_\textrm {esc}$ is the escape velocity of the Sun, $\delta v_\textrm {b}$ is the fluctuating velocity at $r_\textrm {b}$, $P_\textrm {AW}$ is the power carried by outward-propagating AWs, $k_\textrm {B}$ is the Boltzmann constant, and $m_\textrm {p}$ is the proton mass. I then develop a more detailed model of the transition region, corona, and solar wind that accounts for the heat flux $q_\textrm {b}$ from the coronal base into the transition region and momentum deposition by AWs. I solve analytically for $q_\textrm {b}$ by balancing conductive heating against internal-energy losses from radiation, $p\,\textrm {d} V$ work, and advection within the transition region. The density at $r_\textrm {b}$ is determined by balancing turbulent heating and radiative cooling at $r_\textrm {b}$. I solve the equations of the model analytically in two different parameter regimes. In one of these regimes, the leading-order analytic solution reproduces the results of the aforementioned simplified calculation of $\dot {M}$, $U_\infty$, and $T$. Analytic and numerical solutions to the model equations match a number of observations.

Tokamak equilibria have been derived that are analytic solutions to the Grad–Shafranov equation. This paper, Part 1, describes a wide range of such equilibria including smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to include edge pedestals and toroidal flow.

We use three-dimensional (3-D) fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfvén waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic flux ropes as part of the turbulent cascade. The large magnetic structures present in the simulation domain retain the initial anisotropy while the small-scale structures produced by the turbulent cascade are less anisotropic. To quantify the occurrence of reconnection in our simulation domain, we develop a new set of indicators based on intensity thresholds to identify reconnection events in which both ions and electrons are heated and accelerated in 3-D particle-in-cell simulations. According to the application of these indicators, we identify the occurrence of reconnection events in the simulation domain and analyse one of these events in detail. The event is related to the reconnection of two flux ropes, and the associated ion and electron exhausts exhibit a complex 3-D structure. We study the profiles of plasma and magnetic-field fluctuations recorded along artificial-spacecraft trajectories passing near and through the reconnection region. Our results suggest the presence of particle heating and acceleration related to small-scale reconnection events within magnetic flux ropes produced by the anisotropic Alfvénic turbulent cascade in the solar wind. These events are related to current structures of the order of a few ion inertial lengths in size.

The formation of a substantial postdisruption runaway electron current in ASDEX Upgrade material injection experiments is determined by avalanche multiplication of a small seed population of runaway electrons. For the investigation of these scenarios, the runaway electron description of the coupled 1.5-D transport solvers ASTRA-STRAHL is amended by a fluid model describing electron runaway caused by the hot-tail mechanism. Applied in simulations of combined background plasma evolution, material injection and runaway electron generation in ASDEX Upgrade discharge #33108, both the Dreicer and hot-tail mechanism for electron runaway produce only ${\sim }$ 3 kA of runaway current. In colder plasmas with core electron temperatures $T_\textrm {e,c}$ below 9 keV, the postdisruption runaway current is predicted to be insensitive to the initial temperature, in agreement with experimental observations. Yet in hotter plasmas with $T_\textrm {e,c}$ above 10 keV, hot-tail runaway can be increased by up to an order of magnitude, contributing considerably to the total postdisruption runaway current. In ASDEX Upgrade high-temperature runaway experiments, however, no runaway current is observed at the end of the disruption, despite favourable conditions for both primary and secondary runaway.