Steady, helical perturbations known as ‘density snakes’ with poloidal and toroidal mode numbers $m=1$, $n=1$ have been studied in several tokamak experiments. These three-dimensional, helical states are interesting due to their stability and persistence, including their coexistence with the sawtooth cycle. Presented here are studies of density snakes in tokamak plasmas in the Madison Symmetric Torus (MST) device. They are diagnosed using an 11-chord interferometer, internal and edge magnetic coils and impurity ion spectroscopy. Compared with observations in other tokamak plasmas, snakes in MST form with relatively high resistivity and low edge safety factor, $ q(a) \geqslant 2.2$, which moves the $q=1$ resonant surface outward in radius and probably forms a large magnetic island. As a result, the density perturbation associated with the snake is larger, the structure occupies a broader span of minor radius and the snakes are somewhat less stable. The helical structure and distribution of snake events are characterized, including whether they are best described as ideal or resistive kink modes. Finally, an analysis of their perturbation or destruction during sawtooth crashes is given.

This work tackles a significant challenge in dynamo theory: the possibility of long-term amplification and maintenance of an axisymmetric magnetic field. We introduce a novel model that allows for non-trivial axially symmetric steady-state solutions for the magnetic field, particularly when the dynamo operates primarily within a ‘nearly spherical’ toroidal volume inside a fluid shell surrounding a solid core. In this model, Ohm’s law is generalised to include the dissipative force, arising from electron collisions, that tends to align the velocity of the shell with the rotational speed of the inner core and outer mantle. Our findings reveal that, in this context, Cowling’s theorem and the neutral point argument are modified, leading to magnetic energy growth for a suitable choice of toroidal flow. The global equilibrium magnetic field that emerges from our model exhibits a dipolar character. The central insight of the model developed here is that if an additional force is incorporated into Ohm’s law, symmetric dynamos become possible.

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.

This study investigates the impact of local thermal non-equilibrium on the stability analysis of partially ionized plasma within a porous medium. The plasma, heated from below, is enclosed by various combinations of bounding surfaces. Both nonlinear (via the energy method) and linear (utilizing the normal mode analysis method) analyses are performed. Eigenvalue problems for both analyses are formulated and solved using the Galerkin method. The study also explores the effects of compressibility, medium permeability and magnetic fields on system stability. The collisional frequency among plasma components and the thermal diffusivity ratio significantly influence energy decay. The results reveal that the Rayleigh–Darcy number is identical for both nonlinear and linear analyses, thus eliminating the possibility of a subcritical region and confirming global stability. The principle of exchange of stabilities is validated, indicating the absence of oscillatory convection modes. Medium permeability, heat-transfer coefficient and compressibility delay the onset of convection, demonstrating stabilizing effects. Conversely, the porosity-modified conductivity ratio hastens the convection process, indicating destabilizing effects. Rigid–rigid bounding surfaces are found to be thermally more stable for confining the partially ionized plasma. Additionally, the magnetic field exerts a stabilizing influence.

We augment the ‘quasisisymmetric stellarator repository’ (QUASR) to include vacuum field stellarators with quasihelical symmetry using a globalized optimization workflow. The database now has over 300 000 quasisaxisymmetry and quasihelically symmetric devices along with coil sets, optimized for a variety of aspect ratios, rotational transforms and discrete rotational symmetries. This paper outlines a couple of ways to explore and characterize the data set. We plot devices on a near-axis quasisymmetry (QS) landscape, revealing close correspondence to this predicted landscape. We also use principal component analysis (PCA) to reduce the dimensionality of the data so that it can easily be visualized in two or three dimensions. The PCA also gives a mechanism to compare the new devices here with previously published ones in the literature. We are able to characterize the structure of the data, observe clusters and visualize the progression of devices in these clusters. The topology of the data are governed by the interplay of the design constraints and valleys of the QS objective. These techniques reveal that the data has structure, and that typically one, two or three principal components are sufficient to characterize it. The latest version of QUASR is archived at https://zenodo.org/doi/10.5281/zenodo.10050655 and can be explored online at quasr.flatironinstitute.org.

Insight into plasma dynamics under usual pulsed laser deposition (PLD) conditions for NiO thin film growth is provided by implementing angle- and time-resolved Langmuir probe (LP) methods. The selective separation generated an acceleration region that separates ions based on nature and ionisation state. A maximum of the kinetic energy for most plasma components was found for 0.5–2 Pa Ar, while the time-resolved analysis revealed a multipeak evolution of the electron temperature, which widened and shifted with increasing pressure. Evidence of two temperature structures for NiO plasma is presented, and the estimation of the accelerating field generated between the two plasma structures reveals selective in acceleration in the first microsecond. The acceleration field has a maximum value for the O2 atmosphere at approximately 2 Pa, which shows the separation between drift-dominated kinetics and reaction-based dynamics. Further investigation in this 2 Pa region revealed the appearance of a perturbation consistent with the formation of a plasma fireball on the probe. The dynamics of these perturbations is affected by the nature of the gas having different incubation times.

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

When analysing stellarator configurations, it is common to perform an asymptotic expansion about the magnetic axis. This so-called near-axis expansion is convenient for the same reason asymptotic expansions often are, namely, it reduces the dimension of the problem. This leads to convenient and quickly computed expressions of physical quantities, such as quasisymmetry and stability criteria, which can be used to gain further insight. However, it has been repeatedly found that the expansion diverges at high orders in the distance from axis, limiting the physics the expansion can describe. In this paper, we show that the near-axis expansion diverges in vacuum due to ill-posedness and that it can be regularised to improve its convergence. Then, using realistic stellarator coil sets, we demonstrate numerical convergence of the vacuum magnetic field and flux surfaces to the true values as the order increases. We numerically find that the regularisation improves the solutions of the near-axis expansion under perturbation, and we demonstrate that the radius of convergence of the vacuum near-axis expansion is correlated with the distance from the axis to the coils.

We investigate the global stability properties of an electron–positron pair plasma in the linear regime. The plasma is confined by the magnetic field of an infinitely long wire. This configuration is the large-aspect-ratio limit of the levitated dipole experiment of the APEX collaboration. The stability is governed by the diocotron mode and the interchange mode. The diocotron mode dominates in the case of a cold, non-neutral plasma. For specific density profiles we find analytic solutions. We derive a necessary condition for instability and find unstable solutions if the plasma forms a thin shell around the wire. Solutions for arbitrary density profiles with finite temperature are obtained numerically. We find that finite-temperature effects stabilise the diocotron mode. The interchange mode, on the other hand, dominates if the plasma is neutral and has a finite temperature. This mode becomes unstable for a steep-enough density gradient, that is aligned with the gradient of the magnetic field strength and is stabilised by the equilibrium $E\times B$ drift of a non-neutral plasma.

High-beta magnetised plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel temperature gradients can excite whistler waves also in collisional plasmas. Here, we develop a Wigner–Moyal theory for the collisional whistler instability starting from Braginskii-like fluid equations in a slab geometry. This formalism is necessary because, for a large region in parameter space, the fastest-growing whistler waves have wavelengths comparable to the background temperature gradients. We find additional damping terms in the expression for the instability growth rate involving inhomogeneous Nernst advection and resistivity. They (i) enable whistler waves to re-arrange the electron temperature profile via growth, propagation and subsequent dissipation, and (ii) allow non-constant temperature profiles to exist stably. For high-beta plasmas, the marginally stable solutions take the form of a temperature staircase along the magnetic field lines. The electron heat flux can also be suppressed by the Ettingshausen effect when the whistler intensity profile is sufficiently peaked and oriented opposite the background temperature gradient. This mechanism allows cold fronts without magnetic draping, might reduce parallel heat losses in inertial fusion experiments and generally demonstrates that whistler waves can regulate transport even in the collisional limit.

For small-shear helical-axis stellarators, linear ideal-magnetohydrodynamic (MHD) stability calculations and full-torus, nonlinear, electromagnetic gyrokinetic (GK) simulations (the latter with this unprecedented combination of objectives in stellarator GKs) in their linear phase are shown to yield well agreeing spatio-temporal structures of unstable, globally extended perturbations. Likewise, good agreement is found for their dependence on the plasma pressure and the vacuum-field magnetic well in plasma equilibria with identical gradient lengths of the temperature and density profiles. In the nonlinear phase, these perturbations with MHD signatures entail deformations of the magnetic surfaces, growing magnetic islands which rotate in the electron diamagnetic direction and, eventually, lead to ergodisation of a larger part of the magnetic surfaces.

This paper presents a theoretical model of plasma equilibrium in the diamagnetic confinement mode in an axisymmetric mirror device with neutral beam injection. The hot ionic component is described within the framework of the kinetic theory, since the Larmor radius of the injected ions appears to be comparable to or even larger than the characteristic scale of the magnetic field inhomogeneity. The electron drag of the hot ions is taken into account, while the angular scattering of the hot ions due to Coulomb collisions is neglected. The background warm plasma, on the contrary, is considered to be in local thermal equilibrium, i.e. has a Maxwellian distribution function and is described in terms of magnetohydrodynamics. The density of the hot ions is assumed to be negligible compared with that of the warm plasma. Both the conventional gas-dynamic loss and the non-adiabatic loss specific to the diamagnetic confinement mode are taken into account. In this work, we do not consider the effects of the warm plasma rotation as well as the inhomogeneity of the electrostatic potential. A self-consistent theoretical model of the plasma equilibrium is constructed. In the case of the cylindrical bubble, this model is reduced to a simpler one. The numerical solutions in the limit of a thin transition layer of the diamagnetic bubble are found. Examples of the equilibria corresponding to the gas-dynamic multiple-mirror trap device are considered.

Even if the magnetic field in a stellarator is integrable, phase-space integrability for energetic particle guiding-center trajectories is not guaranteed. Both trapped and passing particle trajectories can experience convective losses, caused by wide phase-space island formation, and diffusive losses, caused by phase-space island overlap. By locating trajectories that are closed in the angle coordinate but not necessarily closed in the radial coordinate, we can quantify the magnitude of the perturbation that results in island formation. We characterize island width and island overlap in quasihelical (QH) and quasiaxisymmetric (QA) equilibria with finite plasma pressure $\beta$ for both trapped and passing energetic particles. For trapped particles in QH, low-shear toroidal precession frequency profiles near zero result in wide island formation. While QA transit frequencies do not cross through the zero resonance, we observe that island overlap is more likely since higher shear results in the crossing of more low-order resonances.

Taking advantage of the invariance of the generalised canonical momentum associated with a translational symmetry along a given direction, we describe the dynamics of a plasma by solving an ensemble of $N$ relativistic reduced Vlasov equations coupled in a self-consistent way with the Maxwell equations. This approach, hereafter referred to as the multi-stream model, allows for a drastic reduction in the computational time compared with the full kinetic Vlasov–Maxwell approach. It is also well adapted to a parallel environment. In addition, we extend the model to a two-dimensional geometry in the configuration space, which makes it possible to treat the interaction between several instabilities of beam–plasma and Weibel type, with a relatively small number of streams. The model provides an exact description of current densities perpendicular to a cyclic coordinate, which are responsible, at both fundamental and microscopic levels, for key features of energy transfer, plasma heating and magnetic reconnection processes in collisionless plasmas.

A charged particle in a suitably strong magnetic field spirals along the field lines while slowly drifting transversely. This note provides a brief derivation of an effective Lagrangian formulation for the guiding-centre approximation that captures this dynamics without resolving the gyro motion. It also explains how the effective Lagrangian may, for special magnetic fields, admit a ‘quasi-symmetry’ which can give rise to a conserved quantity helpful for plasma confinement in fields lacking a geometric isometry. The aim of this note is to offer a pedagogical introduction and some perspectives on this well-established subject.

The intense applied lower hybrid electric fields used to drive current in tokamaks can result in the formation of velocity space island structure. When this happens the lower hybrid current drive efficiency can be calculated for a monochromatic wave and is shown to be below the quasilinear level.

The change in direction of the wavevector and group velocity experienced by a wave refracted at the interface of an anisotropic medium in uniform linear motion are determined analytically. These transmission conditions, which are shown to be consistent with the generalised Snell’s law written in the laboratory frame, are then used to examine the effect of motion on waves incident on a magnetised plasma. For an incident wave in the plane perpendicular to the magnetic field the motion is observed to lead to non-negligible deviation of the low-frequency X-mode, as well as to non-symmetrical total reflection angles. These effects are shown to be further complicated when the magnetic field is in the plane formed by the incident wavevector and the medium’s velocity, as the anisotropy now competes with the motion-induced drag. Although obtained in simplified configurations, these results suggest that accounting for motion when modelling plasma wave trajectories could be important under certain conditions, calling for a more detailed quantification of the effect of motion in actual diagnostics and plasma control schemes.

Kinetic theory of particles near resonances is a current topic of discussion in plasma physics and astrophysics. We extend this discussion to the kinetic theory of the interaction between alpha particles (energetic particles predicted to exist in large quantities in next-generation fusion experiments) and a neoclassical tearing mode (NTM), a resistively-driven perturbation which sometimes exists in a tokamak. We develop a quasilinear treatment of the interaction between alpha particles and an NTM, showing why an NTM can be a source of significant passing alpha particle transport in tokamaks. The limitations on quasilinear theory constrain our theory's applicability to small amplitude NTMs, highlighting the importance of nonlinear studies.

In this work, we present a detailed assessment of fusion-born alpha-particle confinement, their wall loads and stability of Alfvén eigenmodes driven by these energetic particles in the Infinity Two Fusion Pilot Plant baseline plasma design, a four-field-period quasi-isodynamic stellarator to operate in deuterium–tritium fusion conditions. Using the Monte Carlo codes, SIMPLE, ASCOT5 and KORC-T, we study the collisionless and collisional dynamics of guiding-centre and full-orbit alpha-particles in the core plasma. We find that core energy losses to the wall are less than 4 %. Our simulations shows that peak power loads on the wall of this configuration are approximately 2.5 MW m-$^2$ and are spatially localised, toroidally and poloidaly, in the vicinity of x-points of the magnetic island chain $n/m = 4/5$ outside the plasma volume. Also, an exploratory analysis using various simplified walls shows that shaping and distance of the wall from the plasma volume can help reduce peak power loads. Our stability assessment of Alfvén eigenmodes using the STELLGAP and FAR3d codes shows the absence of unstable modes driven by alpha-particles in Infinity Two due to the relatively low alpha-particle beta at the envisioned 800 MW operating scenario.