The decay of a turbulent magnetic field is slower with helicity than without. Furthermore, the magnetic correlation length grows faster for a helical than a non-helical field. Both helical and non-helical decay laws involve conserved quantities: the mean magnetic helicity density and the Hosking integral. Using direct numerical simulations in a triply periodic domain, we show quantitatively that in the fractionally helical case the mean magnetic energy density and correlation length are approximately given by the maximum of the values for the purely helical and purely non-helical cases. The time of switchover from one to the other decay law can be obtained on dimensional grounds and is approximately given by $I_{H}^{1/2}I_{M}^{-3/2}$, where $I_{H}$ is the Hosking integral and $I_{M}$ is the mean magnetic helicity density. An earlier approach based on the decay time is found to agree with our new result and suggests that the Hosking integral exceeds naive estimates by the square of the same resistivity-dependent factor by which also the turbulent decay time exceeds the Alfvén time. In the presence of an applied magnetic field, the mean magnetic helicity density is known to be not conserved, and we show that then also the Hosking integral is not conserved.

We analyse distributions of the spatial scales of coherent intermittent structures – current sheets – obtained from fully kinetic, two-dimensional simulations of relativistic turbulence in a collisionless pair plasma using unsupervised machine-learning data dissection. We find that the distribution functions of sheet length $\ell$ (longest scale of the analysed structure in the direction perpendicular to the dominant guide field) and curvature $r_c$ (radius of a circle fitted to the structures) can be well-approximated by power-law distributions, indicating self-similarity of the structures. The distribution for the sheet width $w$ (shortest scale of the structure) peaks at the kinetic scales and decays exponentially at larger values. The data shows little or no correlation between $w$ and $\ell$, as expected from theoretical considerations. The typical $r_c$ depends linearly on $\ell$, which indicates that the sheets all have a similar curvature relative to their sizes. We find a weak correlation between $r_c$ and $w$. Our results can be used to inform realistic magnetohydrodynamic subgrid models for plasma turbulence in high-energy astrophysics.

It is common wisdom that collisionless shocks become non-planar and non-stationary at sufficiently high Mach numbers. Whatever the shock structure, the upstream and downstream fluxes of the mass, momentum and energy should be equal. At low Mach numbers, these conservation laws are satisfied when the shock front is planar and stationary. When this becomes impossible, inhomogeneity and time dependence, presumably in the form of rippling, develop. In this study, we show that the shock structure changes as a kind of ‘phase transition’ when the Mach number is increased while other shock parameters are kept constant.

Standard approaches to neoclassical theory do not extend into regions of strong gradients in tokamaks such as the pedestal and internal transport barriers. Here, we calculate the modifications to neoclassical electron physics inside strong gradient regions of large aspect ratio tokamaks in the banana regime. We show that these modifications are due to the different ion flow and the strong poloidal variation of the potential. We also provide a physical interpretation of the mechanisms that drive poloidal asymmetries and hence a poloidal electric field. We apply our model to two specific example cases of pedestal profiles, calculating the neoclassical electron flux and the bootstrap current. We find that, depending on the ion flow, weak gradient neoclassical theory overestimates or underestimates the neoclassical electron transport and the bootstrap current in regions with strong gradients. We show that the determination of the mean parallel flow is more complex than in weak gradient neoclassical theory. For vanishing turbulence, we can determine the radial electric field for a given flow profile in the pedestal.

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.

We investigate a model of turbulent magnetic reconnection introduced by (Higashimori, Yokoi and Hoshino 2013 Phys. Rev. Lett. 110, 255001) and show that the classic two-dimensional, steady-state Sweet–Parker and Petschek reconnection solutions are supported. We present evidence that these are the only two steady-state reconnection solutions, and we determine the criterion for their selection. Sweet–Parker reconnection occurs when there is no growth in turbulent energy, whereas Petschek reconnection occurs when the current density in the reconnecting current sheet is able to surpass a critical value, allowing for the growth of turbulent energy that creates the diffusion region. Further, we show that the Petschek solutions are self-similar, depending on the value of the turbulent time scale, and produce a universal steady reconnection rate. The self-consistent development of Petschek reconnection through turbulence, within the model, is an example of fast and steady magnetic reconnection without an explicit need for the collisionless terms in an extended Ohm’s law.

This study reports on a set of experiments designed to clarify the impact of the rotational transform on confinement quality at the TJ-II stellarator. For this purpose, the net plasma current is controlled using external coils, resulting in the modification of the rotational transform profile. Significant and systematic variations of the edge electron density gradients (up to $50\,\%{-}60\,\%$) and the plasma energy content ($20\,\%{-}30\,\%$) are achieved. The explanation of this behaviour relies on the placement of low-order rational surfaces in relation to the edge gradient region, which affect local turbulence fluctuation levels, facilitating the formation of zonal flows and concomitant transport barriers. This hypothesis is confirmed experimentally on the basis of a broad array of diagnostic measurements. Calculations based on a resistive magnetohydrodynamic turbulence model provide qualitative support for this hypothesis, clarifying the impact on confinement of specific rational surfaces and highlighting the complex nature of magnetically confined fusion plasmas.

Magnetic reconnection, a fundamental plasma process, is pivotal in understanding energy conversion and particle acceleration in astrophysical systems. While extensively studied in two-dimensional (2-D) configurations, the dynamics of reconnection in three-dimensional (3-D) systems remains under-explored. In this work, we extend the classical tearing mode instability to three dimensions by introducing a modulation along the otherwise uniform direction in a 2-D equilibrium, given by $g(y)$, mimicking a flux-tube-like configuration. We perform linear stability analysis (both analytically and numerically) and direct numerical simulations to investigate the effects of three-dimensionality. Remarkably, we find that a tearing-like instability arises in three dimensions as well, even without the presence of guide fields. Further, our findings reveal that the 3-D tearing instability exhibits reduced growth rates compared with two dimensions by a factor of $\int g(y)^{1/2} {\rm d}y\,/\int {\rm d}y$, with the dispersion relation maintaining similar scaling characteristics. We show that the modulation introduces spatially varying resistive layer properties, which influence the reconnection dynamics.

An ionized gas medium (plasma state) turns to a complex state of plasma or dusty plasma if micrometre- to submicrometre-sized solid dust particles are introduced in it. The dusty plasma medium exhibits fluid- as well as solid-like characteristics at different background plasma conditions. It supports various linear and nonlinear dynamical structures because of external perturbation and internal instabilities. The vortical or coherent structure in the dusty plasma medium is a kind of self-sustained dynamical structure that is formed either by instabilities or by external forcing. In this review article, the author discusses the past theoretical, experimental and computational investigations on vortical and coherent structures in unmagnetized and magnetized dusty plasmas. The possible mechanisms of the formation of vortices in a dust-grain medium are discussed in detail. The studies on the evolution of vortices and their correlation with turbulence are also reviewed.

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ${\sim }$14 times faster and $\sim$7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

Optimal-mode theory (Landreman et al. 2015 J. Plasma Phys. 81, 905810501) can be used to derive upper bounds on growth rates of local gyrokinetic instabilities (Helander & Plunk 2021 Phys. Rev. Lett. 127, 155001). These bounds follow from thermodynamic principles (specifically on the Helmholtz free energy) (Helander & Plunk Phys. Rev. Lett. 127, 2021, p. 155001), and thus apply to any instability and geometry, independently of many plasma parameters. In this work, we compare these upper bounds with the growth rates of linear gyrokinetic eigenmodes. Experimentally relevant scenarios of density-gradient- and ion-temperature-gradient-driven instabilities are considered. The difference between the upper bounds and the numerically computed growth rates is always positive, as it must be, but depends strongly on the instability in question and on the geometry of the magnetic field. The nature of this difference can be analysed by examining the contributions of optimal modes to gyrokinetic eigenmodes. This approach exploits the completeness and orthogonality properties of optimal modes.

We provide a fundamentally new perspective on subcritical turbulence in plasmas, based on coherent structures, which are obtained and characterised via direct numerical solution. The domains where these coherent states exist appear to be closely connected to the those where related turbulent states can exist, so there may be a deep connection between the stability of these coherent structures and the domain where sustained turbulence is possible. In contrast to previous descriptions of turbulence in terms of a stochastic collection of linear waves, we present a fundamentally nonlinear representation based on more general classes of translating oscillatory nonlinear solutions. In turbulent tokamak plasmas, the transport can often be completely suppressed by introducing a background shear flow, whose amplitude is an important control parameter. As this parameter is decreased below a critical value, radially localised structures appear, becoming larger and more complex, in both gyrokinetic simulations and a simpler fluid model of the plasma. For the fluid model, we directly solve for a particular class of nonlinear solutions, relative periodic orbits, and determine their stability, thus explaining why these isolated structures appear in initial-value simulations. The increase of complexity as the flow shear is reduced is explained by a series of Hopf bifurcations of these nonlinear solutions, which we quantify via stability analysis. In gyrokinetic simulations, we are able to indirectly determine the underlying relative periodic orbits by imposing symmetry conditions on the simulations.

The presence of multi-component protons with their distinct features is confirmed by various space missions in the Earth’s outer magnetosphere regions. Isotropic cold protons and anisotropic hot protons significantly influence/modify the dispersion behaviour of various modes and instabilities and regulate the magnetospheric dynamics effectively. Our present study pays attention to the left-hand-polarised proton cyclotron mode, which gets unstable in the large proton temperature anisotropy condition, i.e. $T_{\perp p}\gt T_{\parallel p}$. Such favourable thermal conditions for protons are extensively observed during the compression of the solar wind against the Earth’s magnetic field. To reveal the wave dynamics in more detail, i.e. time-scale variations in the cold and hot proton temperatures and resulting wave-energy density, we further allow the time evolution of our model bi-Maxwellian distribution function in response to the proton cyclotron instability. Based on velocity-moment techniques, we formulated a set of equations comprising an instantaneous dispersion relation, dynamical perpendicular and parallel temperature relations and a wave-energy density equation. For the graphical illustrations of our mathematical results, we choose initial conditions that are relevant to magnetospheric space environments and reported in various experimental studies. Our exact numerical analysis shows the notable impact of hot proton temperature anisotropy and relative density on the real frequency, growth rate, evolution of initial distributions and wave-energy density of the proton cyclotron instability. Such detailed outcomes will be quite helpful for global/local magnetospheric experimental and simulation studies.

We construct a mean-field model that describes the nonlinear dynamics of a spin-polarized electron gas interacting with fixed, positively charged ions possessing a magnetic moment that evolves in time. The mobile electrons are modelled by a four-component distribution function in the two-dimensional phase space $(x,v)$, obeying a Vlasov–Poisson set of equations. The ions are modelled by a Landau–Lifshitz equation for their spin density, which contains ion–ion and electron–ion magnetic exchange terms. We perform a linear response study of the coupled Vlasov–Poisson–Landau–Lifshitz (VPLL) equations for the case of a Maxwell–Boltzmann equilibrium, focussing in particular on the spin dispersion relation. Conditions of stability or instability for the spin modes are identified, which depend essentially on the electron spin polarization rate $\eta$ and the electron–ion magnetic coupling constant $K$. We also develop an Eulerian grid-based computational code for the fully nonlinear VPLL equations, based on the geometric Hamiltonian method first developed by Crouseilles et al. (J. Plasma Phys., vol. 89, no. 2, 2023, p. 905890215). This technique allows us to achieve great accuracy for the conserved quantities, such as the modulus of the ion spin vector and the total energy. Numerical tests in the linear regime are in accordance with the estimations of the linear response theory. For two-stream equilibria, we study the interplay of instabilities occurring in both the charge and the spin sectors. The set of parameters used in the simulations, with densities close to those of solids (${\approx }10^{29}\ {\rm m}^{-3}$) and temperatures of the order of 10 eV, may be relevant to the warm dense matter regime appearing in some inertial fusion experiments.

Landau’s collisionless result for a weakly damped plasma wave is precisely recovered in a weakly collisional, steady state plasma by treating the physics of the narrow collisional boundary layer associated with the resonant electrons. To recover Landau’s results, the collision frequency must be large enough that islands are unable to form and/or the wave amplitude must be small enough to allow linearization. However, the Landau treatment fails once the collision frequency becomes too weak and/or the wave amplitude too large. Remarkably, Landau’s weakly damped plasma wave results require collisions and are shown to be inappropriate in the collisionless limit for a nonlinear, finite amplitude, steady state wave!

The propagation and absorption of the slow waves in the plasma of the Joint European Torus (JET) tokamak have been investigated by ray tracing. The study aims to obtain a qualitative notion of the penetration into the plasma and absorption of the slow wave excited by the A2 ITER-like antenna. The slow waves are radiated by antennas in the ion cyclotron resonance frequency inverted minority heating or mode conversion heating regimes. It has been discovered that the rays propagate in the toroidal direction over a significant distance, up to 6 × 103 cm, from the antenna. Spreading in the peripheral plasma, mainly between the separatrix and the wall, they slowly shift in the poloidal direction and can reach the divertor region. The change in equilibrium of the JET tokamak has a strong influence on both the propagation and absorption of slow waves. Absorption of the slow waves is caused by ion–electron collisions and Landau damping. In the minority heating regimes, the slow waves are strongly damped in the cyclotron resonance of minority ions even at very low minority density.

Using the ONEDFEL code we perform free electron laser simulations in the astrophysically important guide-field dominated regime. For wigglers’ (Alfvén waves) wavelengths of tens of kilometres and beam Lorentz factor ${\sim }10^3$, the resulting coherently emitted waves are in the centimetre range. Our simulations show a growth of the wave intensity over fourteen orders of magnitude, over the astrophysically relevant scale of approximately a few kilometres. The signal grows from noise (unseeded). The resulting spectrum shows fine spectral substructures, reminiscent of those observed in fast radio bursts.

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

One of the critical challenges in future high-current tokamaks is the avoidance of runaway electrons during disruptions. Here, we investigate disruptions mitigated with combined deuterium and noble gas injection in SPARC. We use multi-objective Bayesian optimisation of the densities of the injected material, taking into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The simulations are conducted using the numerical framework Dream (disruption runaway electron analysis model). We show that during deuterium operation, runaway generation can be avoided with material injection, even when we account for runaway electron generation from deuterium–deuterium induced Compton scattering. However, when including the latter, the region in the injected-material-density space corresponding to successful mitigation is reduced. During deuterium–tritium operation, acceptable levels of runaway current and transported heat losses are only obtainable at the highest levels of achievable injected deuterium densities. Furthermore, disruption mitigation is found to be more favourable when combining deuterium with neon, compared with deuterium combined with helium or argon.

The intensity of the turbulence in tokamaks and stellarators depends on its ability to excite and sustain zonal flows. Insight into this physics may be gained by studying the ‘residual’, i.e. the late-time linear response of the system to an initial perturbation. We investigate this zonal-flow residual in the limit of a small magnetic mirror ratio, where we find that the typical quadratic approximation to RH (Rosenbluth & Hinton, 1998 Phys. Rev. Lett. vol. 80, issue 4, pp. 724–727) breaks down. Barely passing particles are in this limit central in determining the resulting level of the residual, which we estimate analytically. The role played by the population with large orbit width provides valuable physical insight into the response of the residual beyond this limit. Applying this result to tokamak, quasi-symmetric and quasi-isodynamic equilibria, using a near-axis approximation, we identify the effect to be more relevant (although small) in the core of quasi-axisymmetric fields, where the residual is smallest. The analysis in the paper also clarifies the relationship between the residual and the geodesic acoustic mode, whose typical theoretical set-ups are similar.