Using transport models, the impacts of trapped electrons on zonal flows and turbulence in helical field configurations are studied. The effect of the trapped electrons on the characteristic quantities of the linear response for zonal flows is investigated for two different field configurations in the Large Helical Device. The turbulent potential fluctuation, zonal flow potential fluctuation and ion energy transport are quickly predicted by the reduced models for which the linear and nonlinear simulation results are used to determine dimensionless parameters related to turbulent saturation levels and typical zonal flow wavenumbers. The effects of zonal flows on the turbulent transport for the case of the kinetic electron response are much smaller than or comparable to those in an adiabatic electron condition for the two different field configurations. It is clarified that the effect of zonal flows on the turbulent transport due to the trapped electrons changes, depending on the field configurations.

We focus our attention on some relevant aspects of the beam–plasma instability in order to refine some features of the linear and nonlinear dynamics. After a re-analysis of the Poisson equation and of the assumption dealing with the background plasma in the form of a linear dielectric, we study the non-perturbative properties of the linear dispersion relation, showing the necessity for a better characterization of the mode growth rate in those flat regions of the distribution function where the Landau formula is no longer predictive. We then upgrade the original $N$-body approach in O'Neil et al. (Phys. Fluids, vol. 14, 1971, pp. 1204–1212), in order to include a return current in the background plasma. This correction term is responsible for smaller saturation levels and growth rates of the Langmuir modes, as result of the energy density transferred to the plasma via the return current. Finally, we include friction effects, as those due to the collective influence of all the plasma charges on the motion of the beam particles. The resulting force induces a progressive resonance detuning, because particles are losing energy and decreasing their velocity. This friction phenomenon gives rise to a deformation of the distribution function, associated with a significant growth of the less energetic particle population. The merit of this work is to show how a fine analysis of the beam–plasma instability outlines a number of subtleties about the linear, intermediate and late dynamics which can be of relevance when such a system is addressed as a paradigm to describe relevant nonlinear wave–particle phenomena (Chen & Zonca, Rev. Mod. Phys., vol. 88, 2016, 015008).

Magnetised plasma turbulence can have a multiscale character: instabilities driven by mean temperature gradients drive turbulence at the disparate scales of the ion and the electron gyroradii. Simulations of multiscale turbulence, using equations valid in the limit of infinite scale separation, reveal novel cross-scale interaction mechanisms in these plasmas. In the case that both long-wavelength (ion-gyroradius-scale) and short-wavelength (electron-gyroradius-scale) linear instabilities are driven far from marginal stability, we show that the short-wavelength instabilities are suppressed by interactions with long-wavelength turbulence. Two novel effects contributed to the suppression: parallel-to-the-field-line shearing by the long-wavelength ${{\boldsymbol {E}} \times \boldsymbol {B}}$ flows, and the modification of the background density gradient by the piece of the long-wavelength electron adiabatic response with parallel-to-the-field-line variation. In contrast, simulations of multiscale turbulence where instabilities at both scales are driven near marginal stability demonstrate that when the long-wavelength turbulence is sufficiently collisional and zonally dominated the effect of cross-scale interaction can be parameterised solely in terms of the local modifications to the mean density and temperature gradients. We discuss physical arguments that qualitatively explain how a change in equilibrium drive leads to the observed transition in the impact of the cross-scale interactions.

When mass falls on the polar regions of a neutron star in a binary X-ray source system, it tends to spread out over the entire surface. A long-standing question in research on this problem is: will the mass be anchored on the magnetic field lines and drag the field with it or is there a special mechanism that allows the mass to slip through the magnetic field lines, leading to much less distortion? As the amount of mass falling on the neutron star can actually be comparable with the neutron star mass, the question of which alternative holds is very important. We suggest an efficient mechanism that will allow the mass to slip through the lines. The mechanism is based on a strong ideal Schwarzschild (Structure and Evolution of the Stars. Princeton University Press, 1958) instability. As the instability itself is ideal, it cannot directly force the mass to slip though the lines. However, it can create a cascade of eddies whose scale extends down to a resistive scale, at the same time mixing the field lines up without breaking them. On this scale the mass can cross the lines. This instability is efficient enough that it can produce a mass flow in the plasma without growing to a large amplitude but saturates at a small one. The instability determines the mass per flux distribution of the accumulated material on different lines so that the equilibrium is marginal to the instability on every line. This makes the equilibrium unique. Thus, as the extra mass on the neutron star grows, the state of the outer shell proceeds through a sequence of unique critically unstable equilibria. In an appendix, an attempt is made to track the critical equilibria over long times.

The concept of the available energy of a collisionless plasma is discussed in the context of magnetic confinement. The available energy quantifies how much of the plasma energy can be converted into fluctuations (including nonlinear ones) and is thus a measure of plasma stability, which can be used to derive linear and nonlinear stability criteria without solving an eigenvalue problem. In a magnetically confined plasma, the available energy is determined by the density and temperature profiles as well as the magnetic geometry. It also depends on what constraints limit the possible forms of plasma motion, such as the conservation of adiabatic invariants and the requirement that the transport be ambipolar. A general method based on Lagrange multipliers is devised to incorporate such constraints in the calculation of the available energy, and several particular cases are discussed for which it can be calculated explicitly. In particular, it is shown that it is impossible to confine a plasma in a Maxwellian ground state relative to perturbations with frequencies exceeding the ion bounce frequency.

The linear MHD (magnetohydrodynamic) stability for high beta plasmas in the inward shifted Large Helical Device (LHD) configurations has been investigated for a wide range of magnetic Reynolds numbers $S$ using numerical simulations based on the kinetic MHD model with kinetic thermal ions where the beta is the ratio of the plasma pressure to the magnetic pressure. It is found that the dependence of the linear growth rate of the resistive ballooning modes on the $S$ number changes from $\unicode[STIX]{x1D6FE}\propto S^{-1/3}$ to $\unicode[STIX]{x1D6FE}\propto S^{-1}$ by the kinetic thermal ion effects so that the resistive ballooning modes are significantly suppressed as the $S$ number increases. For a high $S$ number comparable to experimental values, the most unstable modes are interchange modes. The kinetic thermal ion effects change the most unstable interchange mode from the ideal mode to the resistive mode. This transition of the interchange modes by kinetic thermal ion effects is consistent with the shift of the marginal stability boundary for the ideal interchange modes observed in the LHD experiments.

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight into the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect ratio of the ideal magnetohydrodynamics equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet–Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren–Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest-order solution and the W7-X magnetic field is performed.

In this paper, we discuss the results of a new particle pusher in realistic ultra-strong electromagnetic fields such as those encountered around rotating neutron stars. After presenting the results of this algorithm in simple fields and comparing them to expected exact analytical solutions, we present new simulations for a rotating magnetic dipole in vacuum for a millisecond pulsar by using the Deutsch solution. Particles are injected within the magnetosphere, neglecting radiation reaction, interaction among them and their feedback on the fields. Our simulations are therefore not yet fully self-consistent because the Maxwell equations are not solved according to the current produced by these particles. The code highlights the symmetrical behaviour of particles of opposite charge to mass ratio, $q/m$, with respect to the north and south hemispheres. The relativistic Lorentz factor $\gamma$ of the accelerated particles is proportional to this ratio $q/m$: protons reach up to $\gamma _p \simeq 10^{10.7}$, whereas electrons reach up to $\gamma _e \simeq 10^{14}$. Our simulations show that particles could either be captured by the neutron star, trapped around it or ejected far from it, well outside the light cylinder. Actually, for a given charge to mass ratio, particles follow similar trajectories. These particle orbits show some depleted directions, especially at high magnetic inclination with respect to the rotation axis for positive charges and at low inclination for negative charges because of symmetry. Other directions are preferred and loaded with a high density of particles, some directions concentrating the highest or lowest acceleration efficiencies.

In this work, the first nonlinear particle-in-cell simulations carried out in a stellarator with the global gyrokinetic code EUTERPE using adiabatic electrons and realistic plasma parameters are reported. Several studies are conducted with the aim of enabling reliable nonlinear simulations in stellarators with this code. First, EUTERPE is benchmarked against ORB5 in both linear and nonlinear settings in a tokamak configuration. Next, the use of noise control and stabilization tools, a Krook-type collision operator, markers’ weight smoothing and heating sources is investigated. It is studied in detail how these tools influence the linear growth rate of instabilities in both tokamak and stellarator geometries, and their influence on the linear zonal flow evolution in a stellarator. Then, it is studied how these tools allow improvement of the quality of the results in a set of nonlinear simulations of electrostatic turbulence in a stellarator configuration. Finally, these tools are applied to a W7-X magnetic configuration using experimental plasma parameters.

The classical collision cross-sections of a proton with an argon atom as well as the thermal transport coefficients and rate constant of the colliding $\text{H}^{+}-\text{Ar}$ system are evaluated at the kinetic temperature $T\in [100,10\,000]~(\text{K})$ by means of the asymptotically correct analytical potential constructed for the ground $X^{1}\unicode[STIX]{x1D6F4}^{+}$ state of the ArH+ cation from the highly accurate ab initio data available in the entire range of internuclear distances (Terashkevich et al., J. Quant. Spectrosc. Radiat. Transfer, vol. 234, 2019, pp. 139–146). The results can be useful to estimate thermodynamic, transport and kinetic properties of the Ar/H2 plasma in a wide temperature range.

This paper investigates the motion of three-dimensional ideal magnetohydrodynamics with incompressible flows. The governing equation is performed at steady state, with the magnetic field parallel to the plasma flow. The equations of stationary equilibrium are derived and described mathematically in Cartesian space. Two approaches for derivation of general three-dimensional solutions for Alfvénic and non-Alfvénic flows at constant and variable fluid densities are constructed. The general vector and scalar potentials of the velocity field are used to derive general formulas of general three-dimensional solutions for Alfvénic and non-Alfvénic flows. To verify the general results we have obtained, some examples are presented. An application that may be of interest for coronal loops and solar prominences is presented.

The neutral beam (NB) fast ion confinement in the Large Helical Device (LHD) is studied for several full field ($B_{t}\sim 2.75~\text{T}$) magnetic configurations by a combination of neutron measurement and simulations. To investigate the NB fast ion confinement, we have performed a series of short-pulse NB injection experiments. The experiment results are analysed by the integrated code TASK3D-a. From this investigation, the effective particle diffusion coefficients of the tangential and perpendicular NBs are approximately $D^{\text{eff}}\sim 0.1~\text{m}^{2}~\text{s}^{-1}$ and $D^{\text{eff}}\sim 1~\text{m}^{2}~\text{s}^{-1}$ in the standard configuration. It is clarified that the NB fast ion confinement improves when the plasmas are shifted inward. Moreover, it is also found that the simulation, which considers the deuteron dilution effect due to the presence of impurity ions, can describe a neutron emission rate consistent with the measurement.

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $\delta n/n_0 \gtrsim 2/\beta$ (where $\beta \gg {1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $\delta n/n_0 \gtrsim 2/\beta$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

Kinetic plasma simulations are nowadays commonly used to study a wealth of nonlinear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in ultra-strong electromagnetic fields produced by intense laser beams for the former or by rotating compact objects such as neutron stars and black holes for the latter. In these ultra-strong electromagnetic fields, the gyro-period is several orders of magnitude smaller than the time scale on which we desire to investigate the plasma evolution. Some approximations are required such as, for instance, artificially decreasing the electromagnetic field strength, which is certainly not satisfactory. The main flaw of this downscaling is that it cannot reproduce particle acceleration to ultra-relativistic speeds with a Lorentz factor above $\gamma \approx 10^3$–$10^4$. In this paper, we design a new algorithm able to catch particle motion and acceleration to a Lorentz factor of up to $10^{15}$ or even higher by using Lorentz boosts to special frames where the electric and magnetic field are parallel. Assuming that these fields are locally uniform in space and constant in time, we solve analytically the equation of motion in a tiny region smaller than the length scale of the spatial and temporal gradient of the field. This analytical integration of the orbit severely reduces the constraint on the time step, allowing us to use large time steps, avoiding resolving the ultra-high gyro-frequency. We performed simulations in ultra-strong spatially and time-dependent electromagnetic fields, showing that our particle pusher is able to follow accurately the exact analytical solution for very long times. This property is crucial to properly capture for instance lepton electrodynamics in electromagnetic waves produced by fast rotating neutron stars. We conclude with a simple implementation of our new pusher into a one-dimensional relativistic electromagnetic particle-in-cell code, testing it against plasma oscillations, two-stream instabilities and strongly magnetized relativistic shocks.

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$-pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al., Phys. Plasmas, vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al., J. Plasma Phys., vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

In this semi-analytic study we develop a mathematical model for determining the electromagnetic field due to a current-driven antenna immersed in a cold, magnetized plasma, valid for frequencies below the electron plasma frequency. At each point in the plasma, it is shown that the vacuum electric field of the antenna couples to the plasma conductivity tensor and acts as an infinitesimal source term to drive plasma currents – the total field is then found from the aggregate sum of these point sources, expressed as an integral across the vacuum field. A general solution is provided for both azimuthally symmetric cylindrical coordinates as well as a fully generalized Cartesian solution. As an example of how this general solution may be applied, we solve for the field due to an electric dipole antenna of length $\ell$, aligned along the background field, at frequencies below the ion cyclotron frequency. It is found that the near field decays exponentially with increasing $k_{\bot }z$, whereas the far field exhibits wave-like behaviour. The radiation zone exhibits propagation cones emanating from either end of the dipole, with a propagation angle that is consistent with past analytic studies of inertial Alfvén waves. The mathematical model presented here may be advantageous over other numerical methods, as it allows the user to solve parts of the problems analytically, thereby cutting down significantly on computation time, as well as offering physical insight into the system that may not be evident with other numerical solvers.

The dependence of confinement on input power for a tokamak plasma with regions having a stiff temperature profile is explored. The resilience of the confinement of the core energy to increasing power loss by core radiation from impurities in such situations, as it is anticipated will be required in a demonstration fusion reactor (DEMO) design, is examined.

The formation and control of strongly nonlinear standing plasma waves (SPWs) from a trivial equilibrium by a chirped frequency drive are discussed. If the drive amplitude exceeds a threshold, after passage through the linear resonance in this system, the excited wave preserves the phase locking with the drive, yielding a controlled growth of the wave amplitude. We illustrate these autoresonant waves via Vlasov–Poisson simulations, showing the formation of sharply peaked excitations with local electron density maxima significantly exceeding the unperturbed plasma density. The Whitham averaged variational approach applied to a simplified water bag model yields the weakly nonlinear evolution of the autoresonant SPWs and the autoresonance threshold. If the chirped driving frequency approaches some constant level, the driven SPW saturates at a target amplitude, avoiding the kinetic wave breaking.

The deviation from thermodynamic equilibrium of the ion velocity distribution functions (VDFs), as measured by the Magnetospheric Multiscale (MMS) mission in the Earth’s turbulent magnetosheath, is quantitatively investigated. Making use of the unprecedented high-resolution MMS ion data, and together with Vlasov–Maxwell simulations, this analysis aims at investigating the relationship between deviation from Maxwellian equilibrium and typical plasma parameters. Correlations of the non-Maxwellian features with plasma quantities such as electric fields, ion temperature, current density and ion vorticity are found to be similar in magnetosheath data and numerical experiments, with a poor correlation between distortions of ion VDFs and current density, evidence that questions the occurrence of VDF departure from Maxwellian at the current density peaks. Moreover, strong correlation has been observed with the magnitude of the electric field in the turbulent magnetosheath, while a certain degree of correlation has been found in the numerical simulations and during a magnetopause crossing by MMS. This work could help shed light on the influence of electrostatic waves on the distortion of the ion VDFs in space turbulent plasmas.