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.

We have derived a necessary condition for the achievement of breakeven in axisymmetric plasmas with zero toroidal field and confined by a dipole magnetic field (B. Lehnert, Nature, vol. 181, 1958, p. 4605; A. Hasegawa, Comments Plasma Phys. Control. Fusion, vol. 11, no. 3, 1987). Excellent MHD stability, high values of $\beta$ (up to $26$ %) and good confinement properties awaken the interest of private investors after years of neglect due to lack of public funding and competing alternative lines of research like the tokamak. Starting from a requirement of self-consistency between the balances of momentum and energy in a dipole-confined, two-species plasma and assuming a Maxwellian distribution function for ions and electrons, we derive a necessary condition for breakeven. This condition is more stringent than the Lawson criterion because of the lack of a stabilizing toroidal field. For a given current flowing across the toroidal coil internal to the plasma, the crucial factor at stake is the ratio between the radius of the main toroidal coil and the radius of the vacuum chamber.

This study presents observations of coherent modes (CMs) in a spherical tokamak using a microwave interferometer near the midplane. The CMs within the 30–60 kHz frequency range were observed during electron cyclotron resonance heating only, and the frequency of the CMs increased proportionally with the square root of the electron temperature near $R = 0.7m$. Generally, these modes displayed bursting and chirping signatures with strong density rise and fall. Their appearance indicated an increase in the intensity of hard x rays, suggesting a deterioration in energetic electron confinement. Furthermore, the effect of CMs on the intensity of energetic electron-driven whistler waves was observed. They decreased when CMs were present and gradually increased with the decrease in CM intensity. The CMs may influence the intensity of whistler waves by affecting the energetic electron confinement.

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.

The paper considers the electron transport in toroidal systems taking into account relativistic effects for electrons. The treatment is based on the relativistic drift-kinetic equation with the thermodynamic equilibrium given by the relativistic Maxwell–Jüttner distribution function. The definition of relativistic fluxes is given in a classic-like form using the same set of thermodynamic forces as in the classical (non-relativistic) approach. Such a formulation allows us to apply the currently used non-relativistic solvers for calculation of relativistic mono-energetic transport coefficients. As an example, the procedure for calculating electron fluxes is proposed, in which relativistic effects are taken into account using the DKES code. The model can be easily implemented in various transport codes, developed for the non-relativistic limit, making them accurate also for hot plasmas with non-negligible relativistic effects.

This study investigates the dynamics of various particles within the plasma sheath, focusing on the influence of secondary emissions from charged dust particles. The research concentrates on backscattered electron emission (BEE), inelastic reflection emission (IRE) and true-secondary electron emission (TEE) as key contributors to the behaviour of dust particles within the plasma sheath. Employing the semi-empirical model of Furman and Pivi (F-P model), the study defines the total emission of secondary electrons (EES), comprising these three types. The analysis aims to enhance our understanding of the complex interplay between secondary emission phenomena and the dynamics of charged particles within the plasma sheath, contributing valuable insights to the field. Furthermore, a comparative study has been conducted between the results obtained from the emission of secondary electrons according to the Sternglass theory and the emission of secondary electrons obtained using the F-P model. It is observed that the secondary electron emission (SEE) from the dust, based on the F-P model, demonstrates more pronounced effects on the sheath characteristics, particularly when considering lower values of the superextensive electron parameter ‘$q$’.

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.

Quasilinear treatments are widely used for tokamaks to evaluate radio frequency (rf) heating and current drive. Even though the core of a tokamak plasma is weakly collisional, the solution of the linearized kinetic equation is evaluated using unperturbed collisionless trajectories while often treating successive poloidal circuits of the passing (and trapped) particles as uncorrelated or nearly so. In addition, the most important effect of tokamak geometry, the mirror force, is usually mistreated or ignored when obtaining the solution. These concerning aspects of rf treatments are clarified by considering lower hybrid heating and current drive to illustrate that the electrons in resonance with the applied rf are enclosed by narrow collisional boundary layers, and that tokamak geometry makes it necessary to retain poloidal variation when solving a weakly collisional linearized kinetic equation. Other aspects such as collisional boundary layers at the trapped–passing boundary, cyclotron resonances, and the limitations of quasilinear theory are also considered. The new insights lead to a fundamentally different formulation and interpretation of the solution of the linearized Fokker–Planck equation used for rf quasilinear theory in a tokamak, while retaining many of the features that have contributed to its successful application to rf heating and current drive.

Electron heating and acceleration in collisionless shocks is a long-standing problem. Rapid isotropization of heated electrons cannot be explained solely by the cross-shock potential. In addition, the macroscopic cross-shock potential prevents efficient reflection and injection into the diffusive acceleration regime. Recent observations have shown that small-scale electric fields are present in the shock front, together with the large-scale cross-shock potential. These small-scale fields have been found also in the upstream and downstream regions. Electron heating in shocks is produced by the combined action of the large- and small-scale fields. The large-scale potential determines the energy transferred to the electrons. The small-scale electrostatic fields scatter electrons. Here we study the scattering of electrons on the typical waveforms, namely solitary bipolar spikes and wavepackets. The main effect is the generation of backstreaming electrons with large pitch angles. It is found that wavepackets are more efficient in electron reflection in the interaction of electrons both with a single spike and with multiple spikes.

We derive an advanced surrogate model for predicting turbulent transport at the edge of tokamaks driven by electron temperature gradient (ETG) modes. Our derivation is based on a recently developed sensitivity-driven sparse grid interpolation approach for uncertainty quantification and sensitivity analysis at scale, which informs the set of parameters that define the surrogate model as a scaling law. Our model reveals that ETG-driven electron heat flux is influenced by the safety factor $q$, electron beta $\beta _e$ and normalized electron Debye length $\lambda _D$, in addition to well-established parameters such as the electron temperature and density gradients. To assess the trustworthiness of our model's predictions beyond training, we compute prediction intervals using bootstrapping. The surrogate model's predictive power is tested across a wide range of parameter values, including within-distribution testing parameters (to verify our model) as well as out-of-bounds and out-of-distribution testing (to validate the proposed model). Overall, validation efforts show that our model competes well with, or can even outperform, existing scaling laws in predicting ETG-driven transport.

The formation and evolution of unconfined counter-helicity spheromaks merging have been experimentally investigated by using a magnetized coaxial plasma gun. By comparing the time-dependent photodiode signals and plasma radiation images of counter-helicity spheromaks merging and plasma jets merging, it is found that the field-reversed configuration (FRC) plasma formed by counter-helicity spheromaks merging has a distinct contour and a long maintenance time. For plasma jets merging, the resulting plasma has no discernible contours and a shorter lifetime. In addition, it is inferred from these data that stagnation heating and magnetic reconnection events occur during the counter-helicity spheromaks merging, causing a rapid rise in plasma pressure at the merging midplane and sharp kinks in the field lines near the merger region. By changing different operating parameters and observing the impact on the merger characteristics, it is suggested that the qualitative dynamics of the FRC plasma depends on the balance between the plasma pressure and the magnetic pressure. The high discharge voltage breaks the equilibrium in the merged body, while the large gas-puffed mass just weakens the compression effect of the merged body. These results give us an intuitive understanding of the counter-helicity spheromak merger process and its dependence on discharge parameters, and also provide a distinct perspective for the optimal design of FRC.

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.