Upper bounds on the growth of instabilities in gyrokinetic systems have recently been derived by considering the optimal perturbations that maximise the growth of a chosen energy norm. This technique has previously been applied to two-species gyrokinetic systems with fully kinetic ions and electrons. However, in tokamaks and stellarators, the expectation from linear instability analyses is that the most important kinetic electron contribution to ion-scale modes often comes from the trapped electrons, which bounce faster than the time scale upon which instabilities evolve. As a result, a fully kinetic electron response is not required to describe unstable modes in many cases. Here, we apply the optimal mode analysis to a reduced two-species system consisting of fully gyrokinetic ions and bounce-averaged electrons, with the aim of finding a tighter bound on ion-scale instabilities in toroidal geometry. This analysis yields bounds that are greatly reduced in comparison with the earlier two-species result. Moreover, if the energy norm is properly chosen, wave–particle resonance effects can be captured, reproducing the stabilisation of density-gradient-driven instabilities in maximum-$J$ devices. The optimal mode analysis also reveals that the maximum-$J$ property has an additional stabilising effect on ion-temperature-gradient-driven instabilities, even in the absence of an electron free energy source. This effect is explained in terms of the concept of mode inertia, making it distinct from other mechanisms.

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.

Ion cyclotron resonance heating is a versatile heating method that has been demonstrated to be able to efficiently couple power directly to the ions via the fast magnetosonic wave. However, at temperatures relevant for reactor grade devices such as DEMO, electron damping becomes increasingly important. To reduce electron damping, it is possible to use an antenna with a power spectrum dominated by low parallel wavenumbers. Moreover, using an antenna with a unidirectional spectrum, such as a travelling wave array antenna, the parallel wavenumber can be downshifted by mounting the antenna in an elevated position relative to the equatorial plane. This downshift can potentially enhance ion heating as well as fast wave current drive efficiency. Thus, such a system could benefit ion heating during the ramp-up phase and be used for current drive during flat-top operation. To test this principle, both ion heating and current drive have been simulated in a DEMO-like plasma for a few different mounting positions of the antenna using the FEMIC code. We find that moving the antenna off the equatorial plane makes ion heating more efficient for all considered plasma temperatures at the expense of on-axis heating. Moreover, although current drive efficiency is enhanced, electron damping is reduced for lower mode numbers, thus reducing the driven current in this part of the spectrum.

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.

Using two counter-propagating ultra-intense laser interactions with a solid target, we conducted a study on the generation of electron-positron pairs via the multi-photon Breit–Wheeler (BW) process and trident process. These processes were simulated using the particle-in-cell (PIC) code EPOCH. Our proposed scheme involves irradiating two targets with two counter-propagating lasers. High-energy photons are produced when hot electrons collide with the reflected laser pulse at the target's front, leading to electron and positron pair production. In the single-target scenario, electron bunches are extracted from the target by the p-polarized laser electromagnetic field and accelerated by the laser ponderomotive force before colliding with the counter-propagating laser. However, using two targets enhances pair creation compared with the single-target set-up. We observed that in two-target configurations, the increased number of high-energy gamma-rays contributes to higher-energy electron–positron generation. Additionally, the generation of hot electrons is also more pronounced in this scheme. Consequently, the laser demonstrates higher efficiency in generating gamma photons and positrons in the dual-target set-up, which is beneficial for investigating high-energy pair production and gamma-ray emission. The generated positrons exhibit a density of the order of $10^{27}\,\text {m}^{-3}$ and can be accelerated to energies of 1.5 GeV. The involvement of hot electrons in the target is crucial for generating high-energy photons and positrons. The maximum pair yield reaches $8 \times 10^9$ for the BW process and $10^8$ for the trident process. Notably, the total laser energy conversion efficiencies to electrons, $\gamma$-rays and positrons show improvement in the dual-target configuration. Specifically, the laser energy absorbed by positrons increases from 11.62 % in Case A to 13.12 % in Case B. These enhancements in conversion efficiency and electron/positron density have significant practical implications in experimental set-ups. In both the BW and trident processes, the two-target set-up dominates, highlighting its effectiveness. We also compared the strengths of both approaches, suggesting that these simple models of implementing two targets can be used in experiments as well.

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 present the first experimental observations of the dust acoustic wave where the wave was observed to propagate in the directions of gravity and magnetic field when these directions were not aligned. The experiments were conducted in the Magnetized Dusty Plasma eXperiment facility using a novel electrode system that allows for the angle between gravity and the magnetic field to be varied in a controlled way. This letter reports on measurements in an rf glow discharge argon plasma environment where the angle between direction of gravity and the magnetic field is 45$^{\circ }$. When there was no applied magnetic field, the wave was observed to propagate in the direction of gravity. However, as the magnetic field increased and the ions transitioned from flowing in the direction of gravity to the direction of the magnetic field, a second wave emerged and two distinct waves were observed to simultaneously propagate, one in the direction of gravity and one in the direction of the magnetic field. As the magnetic field was further increased, the wave that propagated in the direction of gravity was suppressed and the wave was only observed to propagate in the direction of the applied magnetic field. We also observe that the speed and the kinetic temperature of the dust for the mode that propagated in the direction of gravity decreased with increasing magnetic field while the speed and the kinetic temperature of the dust for the mode that propagated in the direction of the magnetic field increased with increasing magnetic field. These measurements suggest that an ion-dust streaming instability is at least partially responsible for the high temperatures that have previously been observed in dusty plasmas when the dust acoustic wave is present.

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.

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.

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.

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 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 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.

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.

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.

Progress in understanding multi-scale collisionless plasma phenomena requires employing tools which balance computational efficiency and physics fidelity. Collisionless fluid models are able to resolve spatio-temporal scales that are unfeasible with fully kinetic models. However, constructing such models requires truncating the infinite hierarchy of moment equations and supplying an appropriate closure to approximate the unresolved physics. Data-driven methods have recently begun to see increased application to this end, enabling a systematic approach to constructing closures. Here, we use sparse regression to search for heat flux closures for one-dimensional electrostatic plasma phenomena. We examine OSIRIS particle-in-cell simulation data of Landau-damped Langmuir waves and two-stream instabilities. Sparse regression consistently identifies six terms as physically relevant, together regularly accounting for more than 95 % of the variation in the heat flux. We further quantify the relative importance of these terms under various circumstances and examine their dependence on parameters such as thermal speed and growth/damping rate. The results are discussed in the context of previously known collisionless closures and linear collisionless theory.