The maser instability associated with the loss-cone distribution has been widely invoked to explain the radio bursts observed in the astrophysical plasma environment, such as aurora and corona. In the laboratory plasma of a tokamak, events reminiscent of these radio bursts have also been frequently observed as an electron cyclotron emission (ECE) burst in the microwave range ($\mathrm{\sim }2{f_{\textrm{ce}}}$ near the last closed flux surface) during transient magnetohydrodynamic events. These bursts have a short duration of ~10 μs and display a radiation spectrum corresponding to a radiation temperature ${T_{e,\textrm{rad}}}$ of over $30\ \textrm{keV}$ while the edge thermal electron temperature ${T_e}$ is only in the range of $1\ \textrm{keV}$. Suprathermal electrons can be generated through magnetic reconnection, and a loss-cone distribution can be generated through open stochastic field lines in the magnetic mirror of the near-edge region of a tokamak plasma. Radiation modelling shows that a sharp distribution gradient $\partial f/\partial {v_ \bot } > 0$ at the loss-cone boundary can cause a negative absorption of ECE radiation through the maser instability. The negative absorption then amplifies the radiation so that the microwave intensity is significantly stronger than the thermal value. The significant ${T_{e,\textrm{rad}}}$ from the simulations suggests the potential role of the loss-cone maser instability in generating the ECE burst in a tokamak.

A fusion neutron source (FNS) based on the gas-dynamic trap (GDT, Budker Institute, Novosibirsk) is considered for confinement of two-species plasma heated by neutral beam injection in a regime where the fast ion distribution function is far from Maxwellian. Kinetic instabilities are expected to develop in this regime, and in this paper we investigate the ion-cyclotron instability evolving in moderate densities of pure hydrogen and mixed deuterium–hydrogen target plasmas. The properties of the studied unstable mode, such as its azimuthal wavenumbers, propagation direction and its being affected by changes in the bulk plasma density and composition, allow us to identify it as the drift cyclotron loss cone (DCLC) instability. This mode scatters fast ions and thereby leads to drops in diamagnetic flux signals and increases longitudinal energy and particle losses, with the average energy of the lost ions estimated to be far above the temperature of warm Maxwellian ions. Our interpretation is that the unstable wave grows due to interaction with the fast ions located near the loss cone in the velocity space and scatters them. Applying the method of suppressing the DCLC instability by filling the loss cone with warm plasma, we have determined the values of plasma density and deuterium percentage that allow us to suppress the DCLC instability in the GDT. These findings justify using mixed bulk plasmas in fusion neutron source operation.

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.

The response tensor is derived for a relativistically streaming, strongly magnetized, one-dimensional Jüttner distribution of electrons and positrons, referred to as a pulsar plasma. This is used to produce a general treatment of wave dispersion in a pulsar plasma. Specifically, relativistic streaming, the spread in Lorentz factors in a pulsar rest frame and cyclotron resonances are taken into account. Approximations to the response tensor are derived by making approximations to relativistic plasma dispersion functions appearing in the general form of the response tensor. The cold-plasma limit, the highly relativistic limit and limits related to cyclotron resonances are considered. The theory developed in this paper has applications to generalized Faraday rotation in pulsars and magnetars.

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.

Free normal-flow (NF) conditions at the plasma boundary are shown to be essential for three-dimensional magnetohydrodynamic (MHD) simulations to agree with linear stability theory. A comparative verification study is presented between two different formulations of the boundary conditions (BCs) for velocity perturbations: (i) fully consistent free NF implementation and (ii) fixed NF formulation, neglecting flow perturbations at the numerical boundary. Numerical results are compared with consolidated figures of merit from the linear theory of external kink modes. We consider two classes of initial equilibria presenting different numerical challenges: a uniform current channel surrounded by pure vacuum and a shaped Wesson-like equilibrium, with smooth (polynomial) radial dependency. Only the fully consistent free NF formulation is invariably accurate in modelling the plasma interface at the numerical boundary, without the need of enforcing a pseudovacuum region at the edge of the simulation domain, as in most analogous past studies. This study employs the cylindrical code SPECYL (Cappello & Biskamp, Nucl. Fusion, vol. 36, no. 5, 1996, p. 571) that solves a full-MHD model without pressure gradients, whose fully consistent resistive wall module with free NF BCs was recently successfully verified against the independent code PIXIE3D (Spinicci et al., AIP Adv., vol. 13, no. 9, 2023, p. 095111).

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.

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.

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the $\delta f$ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

This paper investigates the influence of compressibility on the internal wave modes and the density gradient-driven Rayleigh–Taylor (R–T) instability in magnetized strongly coupled dusty plasmas. The dusty magnetohydrodynamic model is formulated for compressible fluids, accounting for the effects of weakly coupled electrons/ions and strongly coupled dust particles under the influence of the gravitational field. The effect of the magnetic field on the dust dynamics has been incorporated through the magnetic force on the electrons and ions in quasineutral dusty plasmas. A dispersion relation of the R–T instability has been derived which has been modified because of the compressibility effect, dust acoustic wave speed and viscoelastic coefficients. The shear Alfvén and compressional viscoelastic wave modes become coupled in the dispersion characteristics. The modified R–T instability criterion is derived in terms of the Alfvén speed, viscoelastic effects and dust grain parameters. The graphical illustrations show that the growth rate of R–T instability has been suppressed due to compressibility, viscoelastic coefficients and dust acoustic speed. The results are useful to discuss the development of R–T instability in magnetized astrophysical dusty plasmas.

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.

A numerical study of the homo-interactions between two falling droplets and between two rising bubbles in a strongly coupled dusty plasma medium is presented in this article. The strongly coupled dusty plasma is considered as a viscoelastic fluid using the generalized hydrodynamic fluid model formalism. Two factors that affect homo-interactions are taken into account: the initial spacing and the coupling strength of the medium. Three different spacings between two droplets are simulated: widely, medium and closely. In each case, the coupling strength has been given as mild–strong and strong. It is shown that the overall dynamic is governed by the competition between the acceleration of two droplets/bubbles due to gravity and the interaction due to the closeness of the droplets/bubbles. Particularly in viscoelastic fluids, apart from the initial separation, shear waves originating from rotating vortices are responsible for the closeness of two droplets or bubbles. Several two-dimensional simulations have been carried out. This work is a continuation of the work done in Parts 1 (Dharodi & Das, J. Plasma Phys., vol. 87, issue 2, 2021, 905870216) and 2 (Dharodi, J. Plasma Phys., vol. 87, issue 4, 2021, 905870402).

The available energy of a plasma is defined as the maximum amount by which the plasma energy can be lowered by volume-preserving rearrangements in phase space, so-called Gardner restacking. A general expression is derived for the available energy of a nearly homogeneous plasma and is shown to be closely related to the Helmholtz free energy, which it can never exceed. A number of explicit examples are given.

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$, and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$, the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma (viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$, the electron inertial scale $d_e$ and the plasma $\beta$. This work is of general consequence in emphasising the fact that high-$\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$. In this situation, the determination of transport coefficients via the standard CE approach is not valid.

This paper is devoted to experimental studies of plasma confinement with high relative pressure ($\beta$) in the Gas Dynamic Trap (BINP, Novosibirsk). In previous high-$\beta$ confinement studies a maximum local $\beta = 0.6$ was achieved in the fast-ion turning point, contributed to by a beam-driven population of fast ions with an anisotropic distribution function. In this study the axial magnetic field profile was modified to bring the turning points closer to one another, which effectively increased the energy density of plasma and pushed the $\beta$ value higher. Experiments were performed for two non-standard magnetic configurations, where the axial fast-ion confinement region length was reduced by 1.5 and 2 times compared with the standard configuration. The average values of $\langle \beta _{\perp } \rangle$ over the plasma central cross-section were found to be 0.1 and 0.18, respectively, for the two configurations, with the latter value significantly exceeding the $\langle \beta _{\perp } \rangle =0.08$ of the standard configuration, in which the previous record was set. Moreover, halving the fast ion confinement region almost doubled the D–D fusion proton flux from the trap centre compared with the standard configuration. The electron temperature in both new magnetic configurations was only slightly smaller than in the standard configuration. In addition, an effect of Alfvén ion–cyclotron instability (AICI) development on the pressure in the turning points is discussed. Presumably, with some decrease in magnetic field an evolving AICI does not result in considerable pressure axial redistribution, so the pressure maximum is in the turning points’ vicinity despite the instability.

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and $P_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number $M_A \approx 1.5$), supersonic (Sonic Mach number $M_S \approx 4-5$) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio $L/\delta \approx 100$, Lundquist number $S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = $\tau _{\text {cool}}^{-1}/\tau _{H}^{-1}\approx 100$); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

Large-amplitude current-driven plasma instabilities, which can transition to the Buneman instability, were observed in one-dimensional simulations to generate high-energy back-streaming ions. We investigate the saturation of multi-dimensional plasma instabilities and its effects on energetic ion formation. Such ions directly impact spacecraft thruster lifetimes and are associated with magnetic reconnection and cosmic ray inception. An Eulerian Vlasov–Poisson solver employing the grid-based direct kinetic method is used to study the growth and saturation of 2D2V collisionless, electrostatic current-driven instabilities spanning two dimensions each in the configuration (D) and velocity (V) spaces supporting ion and electron phase-space transport. Four stages characterise the electric potential evolution in such instabilities: linear modal growth, harmonic growth, accelerated growth via quasi-linear mechanisms alongside nonlinear fill-in and saturated turbulence. Its transition and isotropisation process bears considerable similarities to the development of hydrodynamic turbulence. While a tendency to isotropy is observed in the plasma waves, followed by electron and then ion phase spaces after several ion-acoustic periods, the formation of energetic back-streaming ions is more limited in the 2D2V than in the 1D1V simulations. Plasma waves formed by two-dimensional electrostatic kinetic instabilities can propagate in the direction perpendicular to the net electron drift. Thus, large-amplitude multi-dimensional waves generate high-energy transverse-streaming ions and eventually limit energetic backward-streaming ions along the longitudinal direction. The multi-dimensional study sheds light on interactions between longitudinal and transverse electrostatic plasma instabilities, as well as fundamental characteristics of the inception and sustenance of unmagnetised plasma turbulence.

We solve the linearized gyrokinetic equation, quasineutrality condition and Ampere's law to obtain the dispersion relation of microtearing modes (MTMs) in collisionless low-$\beta$ toroidal plasmas. Consistent with past studies, we find that MTMs are driven unstable by the electron temperature gradient and that this instability drive is mediated by magnetic drifts. The dispersion relation that we derive can be evaluated numerically very quickly and may prove useful for devising strategies to mitigate MTM instability in fusion devices.