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.

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.

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.

This paper presents a novel approach for simulating plasma instabilities in tokamak plasmas using the piecewise field-aligned finite element method in combination with the particle-in-cell method. Our method traditionally aligns the computational grid, but defines the basis functions in piecewise field-aligned coordinates to avoid grid deformation while naturally representing the field-aligned mode structures. This scheme is formulated and implemented numerically. It also applies to the unstructured triangular meshes in principle. We have conducted linear benchmark tests, which agree well with previous results and traditional schemes. Furthermore, multiple-$n$ simulations are also carried out as a proof of principle, demonstrating the efficiency of this scheme in nonlinear turbulence simulations within the framework of the finite element method.

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.

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.

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.

Dusty plasmas typically contain various species of dust particles, though most studies have focused on homogeneous systems. This paper investigates the propagation of dust acoustic waves in an inhomogeneous dusty plasma with an interface, analysing how plasma inhomogeneity influences wave behaviour. Using scattering and reductive perturbation methods, we show that both transmitted and reflected waves depend strongly on the mass ratio between regions. Dust acoustic waves cannot propagate through a dust lattice when the wavelength is smaller than the lattice constant. At a discontinuous interface, at least one transmitted solitary wave is generated, with its amplitude determined by the mass ratio, while at most one reflected solitary wave can exist. These results underscore the critical role of the mass ratio in wave propagation and suggest a method for estimating dust particle masses and properties by analysing the incident, transmitted and reflected waves.

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Stability of the ‘rigid’ ($m = 1$) ballooning mode in a mirror axisymmetric trap is studied for the case of oblique neutral beam injection (NBI), which creates an anisotropic population of fast sloshing ions. Since small-scale modes with azimuthal numbers $m>1$ in long thin (paraxial) mirror traps are easily stabilized by finite-Larmor-radius (FLR) effects, suppression of the rigid ballooning and flute modes would mean stabilization of all magnetohydrodynamic (MHD) modes, with the exception of the mirror and firehose disturbances, which are intensively studied in geophysics, but have not yet been identified in mirror traps. Large-scale ballooning mode can, in principle, be suppressed either by the lateral perfectly conducting wall, or by the end MHD anchors such as the cusp, by biased limiters or by a combination of these two methods. The effects of the wall shape, vacuum gap width between the plasma column and the lateral wall, angle of oblique NBI, radial profile of the plasma pressure and axial profile of the vacuum magnetic field are studied. It is confirmed that the lateral conducting wall still creates the upper stability zone, where the ratio $\beta$ of the plasma pressure to the pressure of vacuum magnetic field exceeds the second critical value $\beta_{\text{cr2}}$, $\beta >\beta _{\text {cr2}}$. However, in many cases the upper zone is clamped from above by mirror instability. When the lateral wall is combined with end MHD anchors, a lower stability zone $\beta <\beta _{\text {cr}1}$ appears, where $\beta$ is below the first critical value $\beta_{\text{cr1}}$. These two zones can overlap in the case of a sufficiently smooth radial pressure profile, and/or a sufficiently low mirror ratio and/or a sufficiently narrow vacuum gap between the plasma column and the lateral wall. However, even in this case, the range of permissible values of beta is limited from above by the threshold of mirror instability $\beta_{\text{mm}}$, so that $\beta <\beta _{\text {mm}}<1$, in contrast to the case of transversal NBI, when neutral beams are injected perpendicularly to the magnetic field.

We present a quantum algorithm based on repeated measurement to solve initial-value problems for nonlinear ordinary differential equations (ODEs), which may be generated from partial differential equations in plasma physics. We map a dynamical system to a Hamiltonian form, where the Hamiltonian matrix is a function of dynamical variables. To advance in time, we measure expectation values from the previous time step and evaluate the Hamiltonian function classically, which introduces stochasticity into the dynamics. We then perform standard quantum Hamiltonian simulation over a short time, using the evaluated constant Hamiltonian matrix. This approach requires evolving an ensemble of quantum states, which are consumed each step to measure the required observables. We apply this approach to the classic logistic and Lorenz systems, in both integrable and chaotic regimes. Our analysis shows that the solutions’ accuracy is influenced by both the stochastic sampling rate and the nature of the dynamical system.

Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalised by particle collisions, we formulate a detailed phenomenological theory for the influence of perpendicular flow shear on magnetised-plasma turbulence. Our theory introduces two distinct regimes, called the weak-shear and strong-shear regimes, each with its own set of scaling laws for the scale and amplitude of the fluctuations and for the level of turbulent heat transport. We discover that the ratio of the typical radial and poloidal wavenumbers of the fluctuations (i.e. their aspect ratio) at the outer scale plays a central role in determining the dependence of the turbulent transport on the imposed flow shear. Our theoretical predictions are found to be in excellent agreement with numerical simulations of two paradigmatic models of fusion-relevant plasma turbulence: (i) an electrostatic fluid model of slab electron-scale turbulence, and (ii) Cyclone-base-case gyrokinetic ion-scale turbulence. Additionally, our theory envisions a potential mechanism for the suppression of electron-scale turbulence by perpendicular ion-scale flows based on the role of the aforementioned aspect ratio of the electron-scale fluctuations.

We develop a phenomenological model of strong imbalanced magnetohydrodynamic (MHD) turbulence that accounts for intermittency and the reflection of Alfvén waves by spatial variations in the Alfvén speed. Our model predicts the slopes of the inertial-range Elsasser power spectra, the scaling exponents of the higher-order Elsasser structure functions and the way in which the parallel (to the magnetic field) length scale of the fluctuations varies with the perpendicular length scale. These predictions agree reasonably well with measurements of solar-wind turbulence from the Parker Solar Probe (PSP). In contrast to previous models of intermittency in balanced MHD turbulence, we find that intermittency in reflection-driven MHD turbulence increases the parallel wavenumbers of the energetically dominant fluctuations at small perpendicular length scales. This, like the PSP measurements with which our model agrees, suggests that turbulence in the solar wind and solar corona may lead to more ion cyclotron heating than previously realized.

The near-axis description of optimised stellarators, at second order in the expansion, provides important information about the field, both of physical and practical importance for stellarator optimisation. It, however, remains relatively underdeveloped for an important class of such stellarators, called quasi-isodynamic (QI). In this paper we develop the theoretical and numerical framework, applying the second-order omnigeneity conditions derived in Rodríguez & Plunk (2023), to make explicit construction of equilibrium solutions. We find that the case of QI stellarators calls for the careful treatment of continuity, smoothness and periodicity of the various functions involved, especially for so-called half-helicity fields, which feature prominently in existing QI designs. The numerical implementation of necessary elements is described, and several examples are constructed and quantitatively verified in detail. This work establishes a basis for further systematic exploration of the space of QI stellarators, and the development of both theoretical and practical tools to facilitate effective optimisation of QI stellarators.

The phenomenon of focusing of microwave beams in a plasma near a turning-point caustic is discussed by exploiting the analytical solution to the Gaussian beam-tracing equations in the two-dimensional (2-D) linear-layer problem. The location of maximum beam focusing and the beam width at that location are studied in terms of the beam initial conditions. This focusing must be taken into account to interpret Doppler backscattering (DBS) measurements. We find that the filter function that characterises the scattering intensity contribution along the beam path through the plasma is inversely proportional to the beam width, predicting enhanced scattering from the beam focusing region. We show that the DBS signal enhancement for decreasing incident angles between the beam path and the density gradient is due to beam focusing and not due to forward scattering, as was originally proposed by (Gusakov et al., (Plasma Phys. Contr. Fusion, vol. 56, 2014, p. 0250092014, 2017); Plasma Phys. Rep. vol. 43(6), 2017, pp. 605–613). The analytic beam model is used to predict the measurement of the $k_y$ density-fluctuation wavenumber power spectrum via DBS, showing that, in an NSTX-inspired example, the spectral exponent of the turbulent, intermediate-to-high $k_y$ density-fluctuation spectrum might be quantitatively measurable via DBS, but not the spectral peak corresponding to the driving scale of the turbulent cascade.

We present a study on the solitons in strongly coupled Yukawa fluids using a simple fluid model (SFE), supplemented by an appropriate equation of state for the medium. The formulation covers a broad range of coupling ($\Gamma$) and screening ($\kappa$) parameters, showing an agreement with the nonlinear quasilocalized charged approximation and generalized hydrodynamic models in the weak screening regime of the solitons in Yukawa media. The results also show a quantitative agreement with the experimentally measured values of the width and Mach number with the normalized amplitude. It has also been observed that the amplitude and width of the soliton in the weak screening limit increase with $\Gamma$ up to $\Gamma \sim 10$, beyond which they remain independent of $\Gamma$ values. Molecular dynamics simulations also confirm that the localization begins to emerge beyond $\Gamma \sim 10$, showing no significant effects on the characteristics of the solitons in Yukawa media. Therefore, the SFE model is capable of predicting the impact of the onset of the localization on the solitons in Yukawa media. Additionally, the amplitude of the soliton increases while its width decreases with $\kappa$ values. The SFE model also explores the possibility of forming refractive soliton structures, whose intensity increases with $\kappa$ values and decreases with $\Gamma$ values.

We analyse the collisionless tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing shear flow, and is completely stabilised as the shear flow becomes Alfvénic. We also show that, in the presence of strong flow shear, the tearing mode growth rate decreases with increasing background ion-to-electron temperature ratio, the opposite behaviour to the tearing mode without flow shear. We find that even a relatively small flow shear is enough to dramatically alter the scaling behaviour of the mode, because the growth rate is small compared with the shear flow across the ion scales (but large compared with shear flow across the electron scales). Our results may explain the relative absence of reconnection events in the near-Sun Alfvénic solar wind observed recently by NASA’s Parker Solar Probe.