A quasi-linear reduced transport model is developed from a database of high-$\beta$ electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

Isolated, undamped geodesic-acoustic-mode (GAM) packets have been demonstrated to obey a (focusing) nonlinear Schrödinger equation (NLSE) (E. Poli, Phys. Plasmas, 2021). This equation predicts susceptibility of GAM packets to the modulational instability (MI). The necessary conditions for this instability are analysed analytically and numerically using the NLSE model. The predictions of the NLSE are compared with gyrokinetic simulations performed with the global particle-in-cell code ORB5, where GAM packets are created from initial perturbations of the axisymmetric radial electric field $E_r$. An instability of the GAM packets with respect to modulations is observed both in cases in which an initial perturbation is imposed and when the instability develops spontaneously. However, significant differences in the dynamics of the small scales are discerned between the NLSE and gyrokinetic simulations. These discrepancies are mainly due to the radial dependence of the strength of the nonlinear term, which we do not retain in the solution of the NLSE, and to the damping of higher radial spectral components $k_r$. The damping of the high-$k_r$ components, which develop as a consequence of the nonlinearity, can be understood in terms of Landau damping. The influence of the ion Larmor radius $\rho _i$ as well as the perturbation wavevector $k_\text {pert}$ on this effect is studied. For the parameters considered here the aforementioned damping mechanism hinders the MI process significantly from developing to its full extent and is strong enough to stabilize some of the (according to the undamped NLSE model) unstable wavevectors.

The magnetohydrodynamic (MHD) equilibrium and stability properties of the Infinity Two fusion pilot plant baseline plasma physics design are presented. The configuration is a four-field period, aspect ratio $A = 10$ quasi-isodynamic stellarator optimised for excellent confinement at elevated density and high magnetic field $B = 9\,T$. Magnetic surfaces exist in the plasma core in vacuum and retain good equilibrium surface integrity from vacuum to an operational $\beta = 1.6 \,\%$, the ratio of the volume average of the plasma and magnetic pressures, corresponding to $800\ \textrm{MW}$ deuterium–tritium fusion operation. Neoclassical calculations show that a self-consistent bootstrap current of the order of ${\sim} 1\ \textrm{kA}$ slightly increases the rotational transform profile by less than 0.001. The configuration has a magnetic well across its entire radius. From vacuum to the operating point, the configuration exhibits good ballooning stability characteristics, exhibits good Mercier stability across most of its minor radius and it is stable against global low-n MHD instabilities up to $\beta = 3.2\,\%$.

In this paper, we calculate the available energy, an upper bound on the thermal energy released that may in turn drive turbulence, due to an ion-temperature-gradient instability driven by curvature in the presence of adiabatic electrons. This is done by choosing an appropriate set of invariants that neglect parallel dynamics, whilst keeping the density profile fixed. Conditions for vanishing available energy are derived and are found to be qualitatively similar to conditions for stability derived from gyrokinetic theory, including strong stabilisation if the ratio of the temperature and density gradient, $\mathrm{d} \ln T / \mathrm{d} \ln n =\eta$, falls in the range $0 \leqslant \eta \leqslant 2/3$. To assess the utility of the available energy, a database consisting of $6 \times 10^4$ local gyrokinetic simulations in randomly sampled tokamak geometries is constructed. Using this database and a similar one sampling stellarators (Landreman et al. 2025 J. Plasma Phys. vol. 91, E120), the available energy is shown to exhibit correlation with the ion energy flux as long as the parallel dynamics is unimportant. Overall it is found that available energy is good at predicting the energy flux variability due to the gradients in density and temperature, but performs worse when it comes to predicting its variability arising from geometry.

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.

A theoretical investigation on the space–time evolution of low-frequency dust acoustic waves (DAWs) in opposite polarity dusty plasmas reveals that they undergo phase mixing for arbitrary initial amplitudes, causing them to suffer a gradual loss in coherency. Both positively and negatively charged dynamical dust grains have been considered to coexist in the plasma, in addition to Maxwell–Boltzmann distributed hot electrons and ions. A perturbative analysis of the governing fluid-Maxwell equations leads us to conclude that the competing dynamics of the opposite polarity dust grains is what causes the DAWs to phase mix. An estimate for the phase-mixing time has also been obtained, which has been found to be profoundly influenced by the values of the various plasma parameters, such as the equilibrium densities of the plasma species, the masses of the opposite polarity dust grains and the electron and ion temperatures. The investigation has also been extended to include phase mixing of DAWs in electron-depleted dusty plasmas. The findings of this study are expected to have relevance in various astrophysical and laboratory plasma environments.

First numerical results from the newly developed pseudo-spectral code P-FLARE (Penalised FLux-driven Algorithm for REduced models) are presented. This flux-driven turbulence/transport code uses a pseudo-spectral formulation with the penalisation method to impose radial boundary conditions. Its concise, flexible structure allows implementing various quasi-two-dimensional reduced fluid models in flux-driven formulation. Here, results from simulations of the modified Hasegawa–Wakatani system are discussed, where particle transport and zonal flow formation, together with profile relaxation, are studied. It is shown that coupled spreading/profile relaxation that one obtains for this system is consistent with a simple one-dimensional model of coupled spreading/transport equations. Then, the effect of a particle source is investigated, which results in the observation of sandpile-like critical behaviour. The model displays profile stiffness for certain parameters, with very different input fluxes resulting in very similar mean density gradients. This is due to different zonal flow levels around the critical value for the control parameter (i.e. the ratio of the adiabaticity parameter to the mean gradient) and the existence for this system of a hysteresis loop for the transition from two-dimensional turbulence to a zonal flow dominated state.

A theoretical framework has been established to investigate the modulational instability of electromagnetic waves in magnetized electron–positron plasmas. The framework is capable of analyzing electromagnetic waves of any intensity and plasmas at any temperature. A fully relativistic hydrodynamic model, incorporating relativistic velocities and thermal effects, is used to describe the relativistic dynamics of particles in plasmas. Under the weakly magnetized approximation, a modified nonlinear Schrödinger equation, governing the dynamics of the envelope of electromagnetic waves in plasmas, is obtained. The growth rate of the modulational instability is then given both theoretically and numerically. By analyzing the dependence of the growth rate on some key physical parameters, the coupled interplay of relativistic effects, ponderomotive forces, thermal effects and magnetic fields on electromagnetic waves can be clarified. The findings demonstrate that specific combinations of physical parameters can significantly enhance modulational instability, providing a theoretical basis for controlling the propagation of electromagnetic waves in plasmas. This framework has broad applicability to most current laser–plasma experiments and high-energy radiation phenomena from stellar surfaces.

The heating effect of electromagnetic waves in ion cyclotron range of frequencies (ICRFs) in magnetic confinement fusion device is different in different plasma conditions. In order to evaluate the ICRF heating effect in different plasma conditions, we conducted a series of experiments and corresponding TRANSP simulations on the EAST tokamak. Both simulation and experimental results show that the effect of ICRF heating is poor at low core electron density. The decrease in electron density changes the left-handed electric field near the resonant layer, resulting in a significant decrease in the power absorbed by the hydrogen fundamental resonance. However, quite a few experiments must be performed in plasma conditions with low electron density. It is necessary to study how to make ICRF heating best in low electron density plasma. Through a series of simulation scans of the parallel refractive index (n//) of the ICRF antenna, it is concluded that the change of the ICRF antenna n// will lead to the change of the left-handed electric field, which will change the fundamental absorption of ICRF power by the hydrogen minority ions. Fully considering the coupling of ion cyclotron wave at the tokamak boundary and the absorption in the plasma core, optimizing the ICRF antenna structure and selecting appropriate parameters such as parallel refractive index, minority ion concentration, resonance layer position, plasma current and core electron temperature can ensure better heating effect in the ICRF heating experiments in the future EAST upgrade. These results have important implications for the enhancement of the auxiliary heating effect of EAST and other tokamaks.

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 apply a continuation method to recently optimized stellarator equilibria with excellent quasi-axisymmetry to generate new equilibria with a wide range of rotational transform profiles. Using these equilibria, we investigate how the rotational transform affects fast-particle confinement, the maximum coil–plasma distance, the maximum growth rate in linear gyrokinetic ion-temperature gradient simulations and the ion heat flux in corresponding nonlinear simulations. We find values of two-term quasi-symmetry error comparable to or lower than those of the similar Landreman–Paul (Phys. Rev. Lett., vol. 128, 2022, 035001) configuration for values of the mean rotational transform $\bar {\iota }$ between $0.12$ and $0.75$. The fast-particle confinement improves with $\bar {\iota }$ until $\bar {\iota } = 0.73$, at which point the degradation in quasi-symmetry outweighs the benefits of further increasing $\bar {\iota }$. The required coil–plasma distance only varies by about ${\pm }10\,\%$ for the configurations under consideration, and is between $2.8$ and $3.3\ \mathrm {m}$ when the configuration is scaled up to reactor size (minor radius $a=1.7\ \mathrm {m}$ and volume-averaged magnetic field strength of $5.86\ \mathrm {T}$). The maximum growth rate from linear gyrokinetic simulations increases with $\bar {\iota }$, but also shifts towards higher $k_y$ values. The maximum linear growth rate is sensitive to the choice of flux tube at rational $\bar {\iota }$, but this can be compensated for by taking the maximum over several flux tubes. The corresponding ion heat fluxes from nonlinear simulations display a non-monotonic relation to $\bar {\iota }$. Sufficiently large positive shear is destabilizing. This is reflected in both linear growth rates and nonlinear heat fluxes.

Electron energisation by magnetic reconnection has historically been studied in the Lagrangian guiding-centre framework. Insights from such studies include that Fermi acceleration in magnetic islands can accelerate electrons to high energies. An alternative Eulerian fluid formulation of electron energisation was recently used to study electron energisation during magnetic reconnection in the absence of magnetic islands. Here, we use particle-in-cell simulations to compare the Eulerian and Lagrangian models of electron energisation in a set-up where reconnection leads to magnetic island formation. We find the largest energisation at the edges of magnetic islands. There, energisation related to the diamagnetic drift dominates in the Eulerian model, while the Fermi related term dominates in the Lagrangian model. The models predict significantly different energisation rates locally. A better agreement is found after integrating over the simulation domain. We show that strong magnetic curvature can break the magnetic moment conservation assumed by the Lagrangian model, leading to erroneous results. The Eulerian fluid model is a complete fluid description and accurately models bulk energisation. However, local measurements of its constituent energisation terms need not reflect locations where plasma is heated or accelerated. The Lagrangian guiding centre model can accurately describe the energisation of particles, but it cannot describe the evolution of the fluid energy. We conclude that while both models can be valid, they describe two fundamentally different quantities, and care should be taken when choosing which model to use.

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.

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-$\varPsi$ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-$\varPsi$ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

Magnetic reconnection leads to the formation of island-shaped magnetic structure(s). Due to disagreement between theoretical evaluations of the characteristic reconnection time and observations, it is commonly accepted that the collisionality (or resistivity) is too low to explain magnetic reconnection phenomena in fusion plasmas. Thus, magnetic reconnection still raises many open questions. The work presented here aims to improve the fundamental knowledge about ‘the life of a magnetic island’. Here, in the light of the many works of the last 70 years, a new paradigm for understanding magnetic reconnection in fusion plasmas is proposed. The life of a magnetic island (whatever its scale) follows three phases: the origin, the growth and the saturation. The possible physical mechanisms at play in these three phases will be investigated. First, for the island origin, typical time scales in link with magnetic reconnection will be evaluated for three tokamaks of different sizes (TCV, WEST and JET) to verify if magnetic reconnection is such an unexplained phenomenon in fusion plasmas. Second, for the island drive, the richness of possible mechanisms leading to ‘rapid’ magnetic island growth in fusion devices will be presented for small and large scales. Third comes the island saturation step. Results on the prediction of a large island width at saturation are presented and discussed.

This paper presents a Hammir tandem mirror confinement performance analysis based on Realta Fusion’s first-of-a-kind model for axisymmetric magnetic mirror fusion performance. This model uses an integrated end plug simulation model including, heating, equilibrium and transport combined with a new formulation of the plasma operation contours (POPCONs) technique for the tandem mirror central cell. Using this model in concert with machine learning optimization techniques, it is shown that an end plug utilizing high temperature superconducting magnets and modern neutral beams enables a classical tandem mirror pilot plant producing a fusion gain Q > 5. The approach here represents an important advance in tandem mirror design. The high-fidelity end plug model enables calculations of heating and transport in the highly non-Maxwellian end plug to be made more accurately. The detailed end plug modelling performed in this work has highlighted the importance of classical radial transport and neutral beam absorption efficiency on end plug viability. The central cell POPCON technique allows consideration of a wide range of parameters in the relatively simple near-Maxwellian central cell, facilitating the selection of more optimal central cell plasmas. These advances make it possible to find more conservative classical tandem mirror fusion pilot plant operating points with lower temperatures, neutral beam energies and end plug performance requirements than designs in the literature. Despite being more conservative, it is shown that these operating points have sufficient confinement performance to serve as the basis of a viable fusion pilot plant provided that they can be stabilized against magnetohydrodynamic and trapped particle modes.

Eddy currents play a significant role in the evolution of tokamak plasmas and must therefore be correctly taken into account in time-dependent simulations. In this paper, a computational method for solving the evolution of tokamak plasma considering eddy currents utilising VMEC (Hirshman & Whitson, Phys. Fluids, vol. 26, 1983, pp. 3553–3568), a commonly used static magnetohydrodynamic equilibrium solver, is proposed. This method is convenient since it does not modify the equilibrium solver internally and achieves convergence calculation through external processing. By allowing the components of the magnetic field to be treated separately, this method provides convergence for cases with displacements in arbitrary directions, which has been difficult to achieve with the previous methods.

Cosmic-ray transport in turbulent astrophysical environments remains a multifaceted problem and, despite decades of study, the impact of complex magnetic field geometry – evident in simulations and observations – has only recently received more focussed attention. To understand how ensemble-averaged transport behaviour emerges from the intricate interactions between cosmic rays and structured magnetic turbulence, we run test-particle experiments in snapshots of a strongly turbulent magnetohydrodynamics simulation. We characterise particle–turbulence interactions via the gyro radii of particles and their experienced field-line curvatures, which reveals two distinct transport modes: magnetised motion, where particles are tightly bound to strong coherent flux tubes and undergo large-scale mirroring; and unmagnetised motion, characterised by chaotic scattering through weak and highly tangled regions of the magnetic field. We formulate an effective stochastic process for each mode: compound subdiffusion with long mean free paths for magnetised motion, and a Langevin process with short mean free paths for unmagnetised motion. A combined stochastic walker that alternates between these two modes accurately reproduces the mean squared displacements observed in the test-particle data. Our results emphasise the critical role of coherent magnetic structures in comprehensively understanding cosmic-ray transport and lay a foundation for developing a theory of geometry-mediated transport.

The regulation of electron heat transport in high-$\beta$, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as $\beta _e^{-1}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of $\beta _e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected $\beta _e^{-1}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known $90^{\circ }$ scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.