The bootstrap current in stellarators can be presented as a sum of a collisionless value given by the Shaing–Callen asymptotic formula and an off-set current, which non-trivially depends on plasma collisionality and radial electric field. Using NEO-2 modeling, analytical estimates and semi-analytical studies with the help of a propagator method, it is shown that the off-set current in the $1/\nu$ regime does not converge with decreasing collisionality $\nu _\ast$ but rather shows oscillations over $\log \nu _\ast$ with an amplitude of the order of the bootstrap current in an equivalent tokamak. The convergence to the Shaing–Callen limit appears in regimes with significant orbit precession, in particular, due to a finite radial electric field, where the off-set current decreases as $\nu _\ast ^{3/5}$. The off-set current strongly increases in case of nearly aligned magnetic field maxima on the field line where it diverges as $\nu _\ast ^{-1/2}$ in the $1/\nu$ regime and saturates due to the precession at a level exceeding the equivalent tokamak value by ${v_E^\ast }^{-1/2}$, where $v_E^\ast$ is the perpendicular Mach number. The latter off-set, however, can be minimized by further aligning the local magnetic field maxima and by fulfilling an extra integral condition of “equivalent ripples” for the magnetic field. A criterion for the accuracy of this alignment and of ripple equivalence is derived. In addition, the possibility of the bootstrap effect at the magnetic axis caused by the above off-set is also discussed.

Landau damping is one of the cornerstones of plasma physics. Based on the initial-value approach adopted by Landau in his original derivation of Landau damping, we examine the solutions of the linear Vlasov–Poisson system for different equilibrium distribution functions $f_0(v)$, going beyond the traditional focus on the root with largest imaginary part and investigating the full set of roots that the dispersion relation of the system generally admits. Specifically, we provide analytical insights into the number and the structure of the roots for entire and meromorphic functions $f_0(v)$, such as Maxwellian and $\kappa$ distributions, we discuss the potential issues related to the redefinition of $\partial{f}_0/\partial{v}$ as a complex variable function and we show how different sigmoids affect the root structure associated with non-meromorphic cut-off distribution functions. Finally, based on the comparison of the several root structures considered, we wonder if the multiple roots might hint at a deeper understanding of the Landau damping phenomenon.

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

In magnetized, stratified environments such as the Sun's corona and solar wind, Alfvénic fluctuations ‘reflect’ from background gradients, enabling nonlinear interactions that allow their energy to dissipate into heat. This process, termed ‘reflection-driven turbulence’, likely plays a key role in coronal heating and solar-wind acceleration, explaining a range of detailed observational correlations and constraints. Building on previous works focused on the inner heliosphere, here we study the basic physics of reflection-driven turbulence using reduced magnetohydrodynamics in an expanding box – the simplest model that can capture local turbulent plasma dynamics in the super-Alfvénic solar wind. Although idealized, our high-resolution simulations and simple theory reveal a rich phenomenology that is consistent with a diverse range of observations. Outwards-propagating fluctuations, which initially have high imbalance (high cross-helicity), decay nonlinearly to heat the plasma, becoming more balanced and magnetically dominated. Despite the high imbalance, the turbulence is strong because Elsässer collisions are suppressed by reflection, leading to ‘anomalous coherence’ between the two Elsässer fields. This coherence, together with linear effects, causes the growth of ‘anastrophy’ (squared magnetic potential) as the turbulence decays, forcing the energy to rush to larger scales and forming a ‘$1/f$-range’ energy spectrum in the process. Eventually, expansion overcomes the nonlinear and Alfvénic physics, forming isolated, magnetically dominated ‘Alfvén vortices’ with minimal nonlinear dissipation. These results can plausibly explain the observed radial and wind-speed dependence of turbulence imbalance (cross-helicity), residual energy, fluctuation amplitudes, plasma heating and fluctuation spectra, as well as making a variety of testable predictions for future observations.

Turbulence self-organization is studied in the flux-driven regime by means of the reduced model Tokam1D. Derived in the electrostatic and isothermal limit but keeping finite electron and ion temperatures, it features two instabilities that are suspected to dominate turbulent transport at the edge of L-mode tokamak plasmas: interchange (a reduced version of the resistive ballooning modes) and collisional drift waves, governed respectively by an effective gravity parameter $g$ and the adiabaticity parameter $C$. The usual properties of these two instabilities are recovered in the linear regime. The nonlinear study focuses on the self-organization of collisional drift-wave turbulence at $g=0$. It is found that the energy stored in zonal flows (ZFs) decreases smoothly at small $C$ due to the reduction of both electric and diamagnetic stresses. Conversely to gradient-driven simulations, no sharp collapse is observed due to the self-consistent evolution of the equilibrium density profile. The ZFs are found to structure into staircases at small and large $C$. These structures exhibit a rich variety of dynamics but are found to be robust to large perturbations. Their nucleation is found to be critically governed by the phase dynamics. Finally, the staircase structures are lost in the gradient-driven regime, when the system is prevented from storing turbulent energy into the equilibrium density (pressure) profile.

The regulation of electron heat transport in high-$\beta$, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as $\beta _e^{-1}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of $\beta _e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected $\beta _e^{-1}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known $90^{\circ }$ scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

In hydrodynamic (HD) turbulence, an exact decomposition of the energy flux across scales has been derived that identifies the contributions associated with vortex stretching and strain self-amplification (Johnson, Phys. Rev. Lett., vol. 124, 2020 104501; J. Fluid Mech., vol. 922, 2021, A3) to the energy flux across scales. Here, we extend this methodology to general coupled advection–diffusion equations and, in particular, to homogeneous magnetohydrodynamic (MHD) turbulence. We show that several MHD subfluxes are related to each other by kinematic constraints akin to the Betchov relation in HD. Applied to data from direct numerical simulations, this decomposition allows for an identification of physical processes and for the quantification of their respective contributions to the energy cascade, as well as a quantitative assessment of their multi-scale nature through a further decomposition into single- and multi-scale terms. We find that vortex stretching is strongly depleted in MHD compared with HD, and the kinetic energy is transferred from large to small scales almost exclusively by the generation of regions of small-scale intense strain induced by the Lorentz force. In regions of large strain, current sheets are stretched by large-scale straining motion into regions of magnetic shear. This magnetic shear in turn drives extensional flows at smaller scales. Magnetic energy is transferred from large to small scales predominantly by the aforementioned current-sheet thinning in regions of high strain. The contributions from current-filament stretching – the analogue to vortex stretching – and from bending of magnetic field-lines into current filaments by vortical motion are both almost negligible, although the latter induces strong backscatter of magnetic energy. Consequences of these results for subgrid-scale turbulence modelling are discussed.

We apply a continuation method to recently optimized stellarator equilibria with excellent quasi-axisymmetry to generate new equilibria with a wide range of rotational transform profiles. Using these equilibria, we investigate how the rotational transform affects fast-particle confinement, the maximum coil–plasma distance, the maximum growth rate in linear gyrokinetic ion-temperature gradient simulations and the ion heat flux in corresponding nonlinear simulations. We find values of two-term quasi-symmetry error comparable to or lower than those of the similar Landreman–Paul (Phys. Rev. Lett., vol. 128, 2022, 035001) configuration for values of the mean rotational transform $\bar {\iota }$ between $0.12$ and $0.75$. The fast-particle confinement improves with $\bar {\iota }$ until $\bar {\iota } = 0.73$, at which point the degradation in quasi-symmetry outweighs the benefits of further increasing $\bar {\iota }$. The required coil–plasma distance only varies by about ${\pm }10\,\%$ for the configurations under consideration, and is between $2.8$ and $3.3\ \mathrm {m}$ when the configuration is scaled up to reactor size (minor radius $a=1.7\ \mathrm {m}$ and volume-averaged magnetic field strength of $5.86\ \mathrm {T}$). The maximum growth rate from linear gyrokinetic simulations increases with $\bar {\iota }$, but also shifts towards higher $k_y$ values. The maximum linear growth rate is sensitive to the choice of flux tube at rational $\bar {\iota }$, but this can be compensated for by taking the maximum over several flux tubes. The corresponding ion heat fluxes from nonlinear simulations display a non-monotonic relation to $\bar {\iota }$. Sufficiently large positive shear is destabilizing. This is reflected in both linear growth rates and nonlinear heat fluxes.

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, $\beta \gg 1$, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-$\beta$, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (${\sim }5\,\%$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

Where and under what conditions the transfer of energy between electromagnetic fields and particles takes place in the solar wind remains an open question. We investigate the conditions that promote the growth of kinetic instabilities predicted by linear theory to infer how turbulence and temperature-anisotropy-driven instabilities are interrelated. Using a large dataset from Solar Orbiter, we introduce the radial rate of strain, a novel measure computed from single-spacecraft data, which we interpret as a proxy for the double-adiabatic strain rate. The solar wind exhibits high absolute values of the radial rate of strain at locations with large temperature anisotropy. We measure the kurtosis and skewness of the radial rate of strain from the statistical moments to show that it is non-Gaussian for unstable intervals and increasingly intermittent at smaller scales with a power-law scaling. We conclude that the velocity field fluctuations in the solar wind contribute to the presence of temperature anisotropy sufficient to create potentially unstable conditions.

Prompted by the relevant problem of temperature inversion (i.e. gradient of density anti-correlated to the gradient of temperature) in astrophysics, we introduce a novel method to model a gravitationally confined multi-component collisionless plasma in contact with a fluctuating thermal boundary. We focus on systems with anti-correlated (inverted) density and temperature profiles, with applications to solar physics. The dynamics of the plasma is analytically described via the coupling of an appropriated coarse-grained distribution function and a temporally coarse-grained Vlasov dynamics. We derive a stationary solution of the system and predict the inverted density and temperature profiles of the two species for scenarios relevant for the corona. We validate our method by comparing the analytical results with kinetic numerical simulations of the plasma dynamics in the context of the two-species Hamiltonian mean-field model. Finally, we apply our theoretical framework to the problem of the temperature inversion in the solar corona, obtaining density and temperature profiles in remarkably good agreement with the observations.

A specific set of dimensionless plasma and turbulence parameters is introduced to characterize the nature of turbulence and its dissipation in weakly collisional space and astrophysical plasmas. Key considerations are discussed for the development of predictive models of the turbulent plasma heating that characterize the partitioning of dissipated turbulent energy between the ion and electron species and between the perpendicular and parallel degrees of freedom for each species. Identifying the kinetic physical mechanisms that govern the damping of the turbulent fluctuations is a critical first step in constructing such turbulent heating models. A set of ten general plasma and turbulence parameters are defined, and reasonable approximations along with the exploitation of existing scaling theories for magnetohydrodynamic turbulence are used to reduce this general set of ten parameters to just three parameters in the isotropic temperature case: the ion plasma beta, the ion-to-electron temperature ratio and the isotropic driving wavenumber. A critical step forward in this study is to identify the dependence of all of the proposed kinetic mechanisms for turbulent damping in terms of the same set of fundamental plasma and turbulence parameters. Analytical estimations of the scaling of each damping mechanism on these fundamental parameters are presented. The power of this approach is illustrated in the development of the first phase diagram for the turbulent damping mechanisms as a function of the ion plasma beta and isotropic driving wavenumber for unity ion-to-electron temperature ratio, showing the regions of this two-dimensional parameter space in which ion Landau and transit-time damping, electron Landau and transit-time damping, ion cyclotron damping, ion stochastic heating, collisionless magnetic reconnection and kinetic ‘viscous’ heating play a role in the damping of the turbulent fluctuations.

Hot electrons produced by parametric instabilities can negatively affect the thermonuclear gain of a direct-drive inertial confinement fusion capsule. A Monte Carlo transport is coupled to the hydrodynamics MULTI-IFE code in order to study the hot electron transport in a high-gain target implosion. The thermonuclear energy produced by the directly driven implosion of a reference target drops from 27 MJ without hot electrons to 13.5 kJ when the hot electrons are taken into account. A clear relationship is established between the level of hot electrons produced and the degradation of the thermonuclear yield produced by the target. It is shown that adding a relatively short laser spike (shock ignition) to the main pulse restores the thermonuclear energy released by the implosion. Different shock-ignition windows are presented depending on the duration of the spike. It appears that the longer the duration of the spike, the lower the power required for shock ignition. Finally, a 300 ps-duration spike requires only 50 TW to restore a thermonuclear energy output of 20 MJ.

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a quasi-axisymmetric stellarator for favourable physics properties while preventing strong shaping of the plasma boundary, which can be difficult to create with external current sources.

Microinstabilities drive turbulent fluctuations in inhomogeneous, magnetised plasmas. In the context of magnetic confinement fusion devices, this leads to an enhanced transport of particles, momentum and energy, thereby degrading confinement. In this work, we describe an application of the adjoint method to efficiently determine variations of gyrokinetic linear growth rates on a general set of external parameters in the local $\delta f$-gyrokinetic model. We then offer numerical verification of this approach. When coupled with gradient-based techniques, this methodology can facilitate the optimisation process for the microstability of the confined plasmas across a high-dimensional parameter space. We present a numerical demonstration wherein the ion-temperature-gradient instability growth rate in a tokamak plasma is minimised with respect to flux surface shaping parameters. The adjoint method approach demonstrates a significant computational speed-up compared with a finite-difference gradient calculation.

This study documents several correlations observed during the first run of the plasma wakefield acceleration experiment E300 conducted at FACET-II, using a single drive electron bunch. The established correlations include those between the measured maximum energy loss of the drive electron beam and the integrated betatron X-ray signal, the calculated total beam energy deposited in the plasma and the integrated X-ray signal, among three visible light emission measuring cameras and between the visible plasma light and X-ray signal. The integrated X-ray signal correlates almost linearly with both the maximum energy loss of the drive beam and the energy deposited into the plasma, demonstrating its usability as a measure of energy transfer from the drive beam to the plasma. Visible plasma light is found to be a useful indicator of the presence of a wake at three locations that overall are two metres apart. Despite the complex dynamics and vastly different time scales, the X-ray radiation from the drive bunch and visible light emission from the plasma may prove to be effective non-invasive diagnostics for monitoring the energy transfer from the beam to the plasma in future high-repetition-rate experiments.

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul (Phys. Rev. Lett., vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$. An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$. Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. (Phys. Plasma, vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$.

The main emphasis of the Laser Ion Generation, Handling and Transport (LIGHT) beamline at GSI Helmholtzzentrum für Schwerionenforschung GmbH are phase-space manipulations of laser-generated ion beams. In recent years, the LIGHT collaboration has successfully generated and focused intense proton bunches with an energy of 8 MeV and a temporal duration shorter than 1 ns (FWHM). An interesting area of application that exploits the short ion bunch properties of LIGHT is the study of ion-stopping power in plasmas, a key process in inertial confinement fusion for understanding energy deposition in dense plasmas. The most challenging regime is found when the projectile velocity closely approaches the thermal plasma electron velocity ($v_{i}\approx v_{e,\text {th}}$), for which existing theories show high discrepancies. Since conclusive experimental data are scarce in this regime, we plan to conduct experiments on laser-generated plasma probed with ions generated with LIGHT at a higher temporal resolution than previously achievable. The high temporal resolution is important because the parameters of laser-generated plasmas are changing on the nanosecond time scale. To meet this goal, our recent studies have dealt with ions of lower kinetic energies. In 2021, laser accelerated carbon ions were transported with two solenoids and focused temporally with LIGHT's radio frequency cavity. A bunch length of 1.2 ns (FWHM) at an energy of 0.6 MeV u$^{-1}$ was achieved. In 2022, protons with an energy of 0.6 MeV were transported and temporally compressed to a bunch length of 0.8 ns. The proton beam was used to measure the energy loss in a cold foil. Both the ion and proton beams will also be employed for energy loss measurements in a plasma target.

We present a stellarator configuration optimized for a large threshold (‘critical gradient’) for the onset of the ion temperature gradient (ITG) driven mode, which achieves the largest critical gradient we have seen in any stellarator. Above this threshold, gyrokinetic simulations show that the configuration has low turbulence levels over an experimentally relevant range of the drive strength. The applied optimization seeks to maximize the drift curvature, leading to enhanced local-shear stabilization of toroidal ITG modes, and the associated turbulence. These benefits are combined with excellent quasi-symmetry, yielding low neoclassical transport and vanishingly small alpha particle losses. Analysis of the resulting configuration suggests a trade-off between magnetohydrodynamic (MHD) and ITG stability, as the new configuration possesses a vacuum magnetic hill.

Laser-driven ion acceleration is well established using solid targets mainly in the target normal sheath acceleration regime. To follow the increasing repetition rate available on high-intensity lasers, the use of high-density gas targets has been explored in the past decade. When interacting with targets reaching densities close to the critical one, the laser pulse can trigger different acceleration mechanisms such as Collisionless Shock Acceleration (CSA) or hole boring. Particle-in-cell simulations using ideal target profiles show that CSA can accelerate a collimated, narrow energy spread and few hundreds of megaelectronvolts ion beam on the laser axis. Nevertheless, in real experiments, the laser will not only interact with an overcritical, thin plasma slab with sharp density gradients, but also with lower density regions surrounding the core of the gas jet, extending to several hundreds of micrometres. The interaction of the laser with these lower density wings will lead to nonlinear effects that will reduce the available energy to drive the shock in the high-density region of the target. Optically tailoring this target could mitigate that issue. Recent experiments conducted on different laser facilities aimed at testing several tailoring configurations. We first tested a scheme with a copropagating picosecond prepulse to create a lower density plasma channel to facilitate the propagation of the main pulse, while the second one was a transverse tailoring driven by nanosecond laser pulses to generate blast waves and form a high-density plasma slab. The main results will be presented here and the methods compared.