An overview of the neutral beam injectors developed at the Budker Institute of Nuclear Physics in Novosibirsk during the last 10 years is presented. These neutral injectors are used for plasma diagnostics, heating and current drive in modern fusion devices with magnetic confinement. An arc or a radio-frequency (RF) discharge generates a plasma in the ion sources of the injectors, and a positive hydrogen or deuterium ion beam is extracted and accelerated by a multiaperture ion-optical system (IOS). The accelerated ion beam is converted into a neutral one in a gas target. The precision multiaperture IOS with spherically concave electrodes provides ballistic focusing of the neutral beam. The high-energy, high-power beam injector based on negative ions, which is currently under development, is described as well. It comprises a RF negative ion source and a wide-aperture electrostatic accelerator separated from the source by a low-energy beam transport line, thereby improving the injector reliability.

We present a stellarator configuration optimized for a large threshold (‘critical gradient’) for the onset of the ion temperature gradient (ITG) driven mode, which achieves the largest critical gradient we have seen in any stellarator. Above this threshold, gyrokinetic simulations show that the configuration has low turbulence levels over an experimentally relevant range of the drive strength. The applied optimization seeks to maximize the drift curvature, leading to enhanced local-shear stabilization of toroidal ITG modes, and the associated turbulence. These benefits are combined with excellent quasi-symmetry, yielding low neoclassical transport and vanishingly small alpha particle losses. Analysis of the resulting configuration suggests a trade-off between magnetohydrodynamic (MHD) and ITG stability, as the new configuration possesses a vacuum magnetic hill.

Linearized gravity around a rotating black hole or compact object introduces the concept of a gravitomagnetic field, which originates from the matter–current in the rotating object. Plasma in proximity to this object is subsequently subjected to motion guided by this gravitomagnetic term (where mass serves as the effective charge) in addition to the conventional magnetic field. Such an interplay of fields complicates the accessible plasma waves of the system and thus merits exploration to delineate this interplay, to identify any observable signatures of the gravitomagnetic field in weakly magnetized systems, and to motivate future numerical work in a fully relativistic setting where the effects may be stronger. In this work we analyse the dispersion of the upper and lower hybrid electrostatic waves in a plasma immersed in both magnetic and gravitomagnetic fields. In particular, we discuss the effective augmentation or cancellation of the two fields under the right conditions for the upper hybrid wave. In contrast, the lower hybrid wave experiences a frequency up-shift from the gravitomagnetic field regardless of whether it is parallel or antiparallel to the magnetic field for the studied field strengths.

Our earth is immersed in the near-earth space plasma environment, which plays a vital role in protecting our planet against the solar-wind impact and influencing space activities. It is significant to investigate the physical processes dominating the environment, for deepening our scientific understanding of it and improving the ability to forecast the space weather. As a crucial part of the National Major Scientific and Technological Infrastructure–Space Environment Simulation Research Infrastructure (SESRI) in Harbin, the Space Plasma Environment Research Facility (SPERF) builds a system to replicate the near-earth space plasma environment in the laboratory. The system aims to simulate the three-dimensional (3-D) structure and processes of the terrestrial magnetosphere for the first time in the world, providing a unique platform to reveal the physics of the 3-D asymmetric magnetic reconnection relevant to the earth's magnetopause, wave–particle interaction in the earth's radiation belt, particles’ dynamics during the geomagnetic storm, etc. The paper will present the engineering design and construction of the near-earth space plasma simulation system of the SPERF, with a focus on the critical technologies that have been resolved to achieve the scientific goals. Meanwhile, the possible physical issues that can be studied based on the apparatus are sketched briefly. The earth-based system is of great value in understanding the space plasma environment and supporting space exploration.

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

Collisions are an important dissipation mechanism in plasmas. When approximating collision operators numerically, it is important to preserve their mathematical structure in order to retain the laws of thermodynamics at the discrete level. This is particularly challenging when considering particle methods. A simple but commonly used collision operator is the Lenard–Bernstein operator, or its modified energy- and momentum-conserving counterpart. In this work, we present a macro-particle discretisation of this operator that is provably energy and momentum preserving.

We present a theoretical analysis of electron pitch-angle scattering by ion-acoustic electrostatic fluctuations present in the Earth's bow shock and, presumably, collisionless shocks in general. We numerically simulate electron interaction with a single wave packet to demonstrate the scattering through phase bunching and phase trapping and quantify electron pitch-angle scattering in dependence on the wave amplitude and wave normal angle to the local magnetic field. The iterative mapping technique is used to model pitch-angle scattering of electrons by a large number of wave packets, which have been reported in the Earth's bow shock. Assuming that successive electron scatterings are not correlated, we revealed that the long-term dynamics of electrons is diffusive. The diffusion coefficient depends on the ratio $\varPhi _0/W$ between the wave packet amplitude and electron energy, $D\propto (\varPhi _0/W)^{\nu }$. A quasi-linear scaling ($\nu \approx 2$) is observed for sufficiently small wave amplitudes, $\varPhi _0\lesssim 10^{-3}W$, while the diffusion is nonlinear ($1<\nu <2$) above this threshold. We show that pitch-angle diffusion of ${\lesssim }1$ keV electrons in the Earth's bow shock can be nonlinear. The corresponding diffusion coefficient scales with the intensity $E_{w}^{2}$ of the electrostatic fluctuations in a nonlinear fashion, $D\propto E_{w}^{\nu }$ with $\nu <2$, while its expected values in the Earth's bow shock are $D\sim 0.1\unicode{x2013}100$ $(T_{e}/W)^{\nu -1/2}\,{\rm rad}^{2}\,{\rm s}^{-1}$. We speculate that in the Earth's quasi-perpendicular bow shock the stochastic shock drift acceleration mechanism with pitch-angle scattering provided by the electrostatic fluctuations can contribute to the acceleration of thermal electrons up to approximately 1 keV. The potential effects of a finite perpendicular coherence scale of the wave packets on the efficiency of electron scattering are discussed.

The main emphasis of the Laser Ion Generation, Handling and Transport (LIGHT) beamline at GSI Helmholtzzentrum für Schwerionenforschung GmbH are phase-space manipulations of laser-generated ion beams. In recent years, the LIGHT collaboration has successfully generated and focused intense proton bunches with an energy of 8 MeV and a temporal duration shorter than 1 ns (FWHM). An interesting area of application that exploits the short ion bunch properties of LIGHT is the study of ion-stopping power in plasmas, a key process in inertial confinement fusion for understanding energy deposition in dense plasmas. The most challenging regime is found when the projectile velocity closely approaches the thermal plasma electron velocity ($v_{i}\approx v_{e,\text {th}}$), for which existing theories show high discrepancies. Since conclusive experimental data are scarce in this regime, we plan to conduct experiments on laser-generated plasma probed with ions generated with LIGHT at a higher temporal resolution than previously achievable. The high temporal resolution is important because the parameters of laser-generated plasmas are changing on the nanosecond time scale. To meet this goal, our recent studies have dealt with ions of lower kinetic energies. In 2021, laser accelerated carbon ions were transported with two solenoids and focused temporally with LIGHT's radio frequency cavity. A bunch length of 1.2 ns (FWHM) at an energy of 0.6 MeV u$^{-1}$ was achieved. In 2022, protons with an energy of 0.6 MeV were transported and temporally compressed to a bunch length of 0.8 ns. The proton beam was used to measure the energy loss in a cold foil. Both the ion and proton beams will also be employed for energy loss measurements in a plasma target.

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high-$\beta$ collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high-$\beta$ means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a ‘global-to-local’ approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighbourhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform and quality of quasi-axisymmetry. The database of stellarators, which comprises approximately 200 000 coil sets, is available online and is called QUASR, for ‘quasi-symmetric stellarator repository’.

This paper establishes a mean-field equation set and an energy theorem to provide a theoretical basis in view of the development of self-consistent, physics-based turbulent transport models for mean-field transport codes. A rigorous averaging procedure identifies the exact form of the perpendicular turbulent fluxes which are modelled by ad hoc diffusive terms in mean-field transport codes, next to other closure terms which are not commonly considered. Earlier work suggested that the turbulent $E\times B$ particle and heat fluxes, which are thus identified to be important closure terms, can be modelled to reasonable accuracy using the kinetic energy in the $E\times B$ velocity fluctuations ($k_{E}$). The related enstrophy led to further modelling improvements in an initial study, although further analysis is required. To support this modelling approach, transport equations are derived analytically for both quantities. In particular, an energy theorem is established in which the various source and sink terms of $k_{E}$ are shown to couple to mean-field and turbulent parallel kinetic energy, kinetic energy in the other perpendicular velocity components, the thermal energy and the magnetic energy. This provides expressions for the interchange, drift-wave and Reynolds stress terms amongst others. Note that most terms in these energy equations are in turn closure terms. It is suggested to evaluate these terms using reference data from detailed turbulence code simulations in future work.

Observations, mainly of outbursts in dwarf novae, imply that the anomalous viscosity in highly ionized accretion discs is magnetic in origin and requires that the plasma ${\beta \sim 1}$. Until now, most simulations of the magnetic dynamo in accretion discs have used a local approximation (known as the shearing box). While these simulations demonstrate the possibility of a self-sustaining dynamo, the magnetic activity generated in these models saturates at $\beta \gg 1$. This long-standing discrepancy has previously been attributed to the local approximation itself. There have been recent attempts at simulating magnetic activity in global accretion discs with parameters relevant to the dwarf novae. These too find values of $\beta \gg 1$. We speculate that the tension between these models and the observations may be caused by numerical magnetic diffusivity. As a pedagogical example, we present exact time-dependent solutions for the evolution of weak magnetic fields in an incompressible fluid subject to linear shear and magnetic diffusivity. We find that the maximum factor by which the initial magnetic energy can be increased depends on the magnetic Reynolds number as ${\mathcal {R}}_{m}^{2/3}$. We estimate that current global numerical simulations of dwarf nova discs have numerical magnetic Reynolds numbers around six orders of magnitude less than the physical value found in dwarf nova discs of ${\mathcal {R}}_{m} \sim 10^{10}$. We suggest that, given the current limitations on computing power, expecting to be able to compute realistic dynamo action in observable accretion discs using numerical MHD is, for the time being, a step too far.

Beam–target reactions in fusion plasmas play an important role in both magnetic confinement fusion and inertial confinement fusion in the condition of low-density plasmas with high-velocity interactions. The traditional method for calculating beam–target reaction rate neglects the transport process of incident particles in inhomogeneous plasmas, leading to errors providing that the temperature and density in the transport path of incident particles vary obviously. An improved method considering the transport process is proposed in this paper to eliminate the deficiencies. Then the method is employed in high-speed plasma collision studies. When the initial plasma density and temperature are set to $0.5\,{\rm g}\,{\rm cm}^{-3}$ and 100 eV, it is found that the beam–target reaction rate calculated by the traditional method is almost identical to that by our method if the collision velocity is less than 600 km s$^{-1}$. However, the traditional method is not suitable for study as the collision velocity gets higher, inducing obvious differences, which can reach 70 % at 1000 km s$^{-1}$. The improved method will make large corrections to evaluate the importance of the non-negligible beam–target reaction for inertial confinement fusion schemes with large implosion velocity such as double-cone ignition and impact ignition, in which the high-speed plasmas collide with each other to realize plasma ignition.

In the limit of sufficiently fast rotation, rotating mirror traps are known to be stable against the loss-cone modes associated with conventional (non-rotating) mirrors. This paper calculates how quickly a mirror configuration must rotate in order for several of these modes to be stabilized (in particular, the high-frequency convective loss cone, drift cyclotron loss cone and Dory–Guest–Harris modes). Commonalities in the stabilization conditions for these modes then motivate a modified formulation of the Gardner free energy and diffusively accessible free energy to be used for systems in which the important modes have wavevectors that are orthogonal or nearly orthogonal to the magnetic field, as well as a modification to include the effects of a loss region in phase space.

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

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.

Overtaking collisions of large-amplitude solitons are investigated via fluid simulations for a plasma consisting of cold ions and Boltzmann-distributed electrons. To achieve this, a new fluid simulation code is presented. In addition, a novel approach to construct soliton initial conditions is developed. Using these ideas, initial conditions are combined that allows the simulation of overtaking collisions. It is shown that, in the small-amplitude regime, simulation results agree well with the two-soliton solution obtained from reductive perturbation theory. Interestingly, in the large amplitude regime, both the slow and fast solitons re-emerge after the collision with no significant change, showing that the collisions remain elastic. A comparison between reductive perturbation analysis and the simulations show that the only significant effect of higher order nonlinearities on overtaking collisions is a reduction in the magnitude of the phase shifts of both solitons.

Electron-cyclotron resonance heating (ECRH) is the main heating mechanism in the Wendelstein 7-X (W7-X) stellarator. Although second-harmonic ECRH (X2) has been used routinely for plasma startup, startup at third harmonic (X3) is known to be much more difficult. In this work, we investigate the energy gain of particles during nonlinear wave–particle interaction for conditions relevant to second- and third-harmonic startups in W7-X. We take into account both the beam and the ambient magnetic field inhomogeneities. The latter is shown to significantly increase the mean energy gain resulting from a single wave–particle resonant interaction. In W7-X-like conditions, the improvement in maximum gained energy is up to 4 times the analogous uniform magnetic field case. However, this improvement is not enough to ensure X3 startup. The optimal magnetic field inhomogeneity length scale for average energy gain and start up in W7-X-like conditions is found to be in the range of $1$ to $3\ {\rm km}^{-1}$. A possibility of using multiple beams with neighbouring resonances is also considered. A considerable enhancement of the energy gain is demonstrated.

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

The gyrokinetic theory of the residual flow, in the electrostatic limit, is revisited, with optimized stellarators in mind. We consider general initial conditions for the problem, and identify cases that lead to a non-zonal residual electrostatic potential, i.e. one having a significant component that varies within a flux surface. We investigate the behaviour of the ‘intermediate residual’ in stellarators, a measure of the flow that remains after geodesic acoustic modes have damped away, but before the action of the slower damping that is caused by unconfined particle orbits. The case of a quasi-isodynamic stellarator is identified as having a particularly large such residual, owing to the small orbit width achieved by optimization.