In a magnetised plasma on scales well above ion kinetic scales, any constant-magnitude magnetic field, accompanied by parallel Alfvénic flows, forms a nonlinear solution in an isobaric, constant-density background. These structures, which are also known as spherically polarised Alfvén waves, are observed ubiquitously in the solar wind, presumably created by the growth of small-amplitude fluctuations as they propagate outwards in the corona. Here, we present a computational method to construct such solutions of arbitrary amplitude with general multidimensional structure, and explore some of their properties. The difficulty lies in computing a zero-divergence, constant-magnitude magnetic field, which leaves a single, quasi-free function to define the solution, while requiring strong constraints on any individual component of the field. Motivated by the physical process of wave growth in the solar wind, our method circumvents this issue by starting from low-amplitude Alfvénic fluctuations dominated by a strong mean field, then ‘growing’ magnetic perturbations into the large-amplitude regime. We present example solutions with non-trivial structure in one, two and three dimensions, demonstrating a clear tendency to form very sharp gradients or discontinuities, unless the solution is one-dimensional. As well as being useful as an input for other calculations, particularly the study of parametric decay, our results provide a natural explanation for the extremely sharp field discontinuities observed across magnetic field switchbacks in the low solar wind.

Weakly collisional plasmas are subject to nonlinear relaxation processes, which can operate at rates much faster than the particle collision frequencies. This causes the plasma to respond like a magnetised fluid despite having long particle mean free paths. In this Letter the effective collisional mechanisms are modelled in the plasma kinetic equation to produce density, pressure and magnetic-field responses to compare with spacecraft measurements of the solar wind compressive fluctuations at 1 AU. This enables a measurement of the effective mean free path of the solar wind protons, found to be ${\approx }4 \times 10^{5}$ km, which is approximately $10^{3}$ times shorter than the collisional mean free path. These measurements are shown to support the effective fluid behaviour of the solar wind at scales above the proton gyroradius and demonstrate that effective collision processes alter the thermodynamics and transport of weakly collisional plasmas.

We examine the turbulence driven by the ion and electron temperature gradients in selected magnetic configurations of the Wendelstein 7-X (W7-X) stellarator. The inherent flexibility in the configuration space of W7-X enables us to find candidate configurations manifesting low turbulent transport. We follow insights gained by stellarator optimization techniques, in order to identify key geometric features, which are directly related to the ion and electron heat fluxes produced by plasma turbulence. One such a feature is the flux expansion at locations where the curvature is particularly unfavourable. Starting from a configuration routinely used in the W7-X experiment, we end up with an optimized configuration. Based on this equilibrium, we select a configuration from W7-X configuration database with a similar feature as the optimized one. With the help of nonlinear gyrokinetic simulations, we show that the heat flux in this configuration is less stiff than in the initial configuration, both for ion temperature gradient and electron temperature gradient turbulence.

Dissipative processes cause collisionless plasmas in many systems to develop non-thermal particle distributions with broad power-law tails. The prevalence of power-law energy distributions in space/astrophysical observations and kinetic simulations of systems with a variety of acceleration and trapping (or escape) mechanisms poses a deep mystery. We consider the possibility that such distributions can be modelled from maximum-entropy principles, when accounting for generalizations beyond the Boltzmann–Gibbs entropy. Using a dimensional representation of entropy (related to the Renyi and Tsallis entropies), we derive generalized maximum-entropy distributions with a power-law tail determined by the characteristic energy scale at which irreversible dissipation occurs. By assuming that particles are typically energized by an amount comparable to the free energy (per particle) before equilibrating, we derive a formula for the power-law index as a function of plasma parameters for magnetic dissipation in systems with sufficiently complex topologies. The model reproduces several results from kinetic simulations of relativistic turbulence and magnetic reconnection.

We report a generation of energetic protons by the interaction of a high-energy electron driving beam with an underdense plasma slab. After an interaction period of approximately 4000 fs, a proton beam with maximum energy greater than 250 MeV can be achieved by applying a driving beam with energy 1.0 GeV to a 200 $\mathrm {\mu }$m plasma slab. Our two-dimensional particle-in-cell simulations also show that the proton acceleration process can be divided into two stages. In the first stage, a strong positive longitudinal electric field appears near the rear boundary of the plasma slab after the driving beam has passed through it. This acceleration process is similar to the target normal sheath acceleration scheme by the interaction between intense pulsed laser with overdense plasma targets. In the second stage, the accelerated protons experience a long-range acceleration process with a two-stream instability between the high-energy driving beam and the proton beam. Further analyses show that this accelerated proton beam is equipped with the property of good collimation and high energy. This scheme presents a new way for proton or ion acceleration on some special occasions.

We carry out three-dimensional and two-dimensional particle-in-cell simulations of the expansion of a magnetized plasma that initially uniformly fills a half-space and contains a semicylindrical region of heated electrons elongated along the surface of the plasma boundary. This geometry is related, for instance, to ablation of a plane target by a femtosecond laser beam under quasi-cylindrical focusing. We find that a decay of the inhomogeneous plasma–vacuum discontinuity is strongly affected by an external magnetic field parallel to its boundary. We observe various transient phenomena, including an anisotropic scattering of electrons and an accompanying Weibel instability, and reveal various spatial structures of the arising magnetic field and current, including multiple flying-apart filaments of a Z-pinch type and slowly evolving current sheets with different orientations. The magnitude of the self-generated magnetic field can be of the order of, or significantly exceed that of, the external one. Such phenomena are expected in the laser and cosmic plasmas, including the explosive processes in the planetary magnetospheres and stellar coronal arches.

Fluctuation measurements reveal the outward electromagnetic energy flux needed to drive the dynamo electromotive force supporting magnetic self-organization in a reversed-field pinch plasma. The radial Poynting flux due to tearing mode fluctuations is measured with an insertable probe during magnetic relaxation. This flux corresponds to transient power levels much larger than the input power and comparable to the global equilibrium magnetic energy transient loss rate. The probe measurements of this flux are roughly as predicted by a simple Poynting's theorem model upon substitution of equilibrium measurement data.

In the presence of a stationary gravitomagnetic field, a weakly rotating self-gravitating fluid relaxes into equilibrium flow configurations which admit both large- and small-scale structures. Unlike the traditional Beltrami equilibria in fluid, the equilibrium states in a rotating self-gravitating fluid in a gravitomagnetic field have structure similar to that of double curl Beltrami equilibria in plasma where the generalized vortical lines are aligned with the flow field. Different equilibrium flow configurations in the rotating fluid can be distinguished by the ratio between total energy and helicity. However, these fluid equilibria do not exhibit diamagnetic behaviours as observed in multi-species plasma equilibria.

We show that large-amplitude, non-planar, Alfvén-wave (AW) packets are exact nonlinear solutions of the relativistic magnetohydrodynamic equations when the total magnetic-field strength in the local fluid rest frame ($b$) is a constant. We derive analytic expressions relating the components of the fluctuating velocity and magnetic field. We also show that these constant-$b$ AWs propagate without distortion at the relativistic Alfvén velocity and never steepen into shocks. These findings and the observed abundance of large-amplitude, constant-$b$ AWs in the solar wind suggest that such waves may be present in relativistic outflows around compact astrophysical objects.

Magnetic reconnection is explored on the Terrestrial Reconnection Experiment (TREX) for asymmetric inflow conditions and in a configuration where the absolute rate of reconnection is set by an external drive. Magnetic pileup enhances the upstream magnetic field of the high-density inflow, leading to an increased upstream Alfvén speed and helping to lower the normalized reconnection rate to values expected from theoretical consideration. In addition, a shock interface between the far upstream supersonic plasma inflow and the region of magnetic flux pileup is observed, important to the overall force balance of the system, thereby demonstrating the role of shock formation for configurations including a supersonically driven inflow. Despite the specialized geometry where a strong reconnection drive is applied from only one side of the reconnection layer, previous numerical and theoretical results remain robust and are shown to accurately predict the normalized rate of reconnection for the range of system sizes considered. This experimental rate of reconnection is dependent on system size, reaching values as high as 0.8 at the smallest normalized system size applied.

In an early study of the properties and capabilities of the multiregion, relaxed magnetohydrodynamic model, Hole, Hudson & Dewar claim that they are able to construct a multiregion stepped pressure cylindrical equilibrium which does not require the existence of surface currents. We present a brief argument showing that this claim is incorrect, and clarify the meaning of their statement. Furthermore, even with the statement clarified, we demonstrate that it is not possible to find solutions to reproduce the equilibrium corresponding to the parameters given in the article. We invite the authors to provide a corrigendum with the correct values of the equilibrium they constructed.

Finding less complicated coils that have adequately low field errors is a crucial step in stellarator development. One coil metric that is of high importance is the maximum curvature of the coil centreline, or coil single filament. Conductors cannot be bent below some threshold minimum radius of curvature. High coil curvatures can cause strains to exceed acceptable levels, especially in superconducting coils. We investigate three ways to optimize coil curvature and find that applying penalty functions to the coil curvature solves for coils that have a constrained maximum curvature and low field error. Penalty functions are implemented in FOCUS and coil solutions optimized for an HSX-like ‘plasma boundary’ are presented.

A magnetic field configuration of an island divertor for a quasi-axisymmetric stellarator (CFQS) is proposed. The configuration incorporates large islands surrounding the core confinement region. The interface between the core region and the peripheral region of the island divertor is a clear magnetic separatrix similar to a tokamak divertor. The structure of divertor magnetic field lines is very regular without stochasticity and the connection length is sufficiently long for good divertor performance. Such a divertor configuration is produced in the magnetic field configuration for the CFQS device, which is now under construction in China.

A fluid system is derived to describe electrostatic magnetized plasma turbulence at scales somewhat larger than the Larmor radius of a given species. It is related to the Hasegawa–Mima equation, but does not conserve enstrophy, and, as a result, exhibits a forward cascade of energy, to small scales. The inertial-range energy spectrum is argued to be shallower than a $-11/3$ power law, as compared to the $-5$ law of the Hasegawa–Mima enstrophy cascade. This property, confirmed here by direct numerical simulations of the fluid system, may help explain the fluctuation spectrum observed in gyrokinetic simulations of streamer-dominated electron-temperature-gradient driven turbulence (Plunk et al., Phys. Rev. Lett., vol. 122, 2019, 035002), and also possibly some cases of ion-temperature-gradient driven turbulence where zonal flows are suppressed (Plunk et al., Phys. Rev. Lett., vol. 118, 2017, 105002).

Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.

We develop a simple model for the kinematics of charged particles in regions of magnetic turbulence. We approximate the local magnetic field as smoothly varying in strength and direction, where adiabatic invariance prevails, or as presenting rapid changes in direction or ‘kinks’. Particles execute guiding centre gyromotion around a field line. However, in analogy to kinetic theory for collisional environments, when the particle undergoes a rapid change in direction by some angle $\unicode[STIX]{x1D703}$, it would instantaneously transition to Larmor motion around the new field line. This mimics Brownian motion wherein we replace collisions with other particles by rapid transitions or ‘collisions’ with other field lines. Using standard methods drawn from Brownian motion, we follow the evolution of the parallel and perpendicular components of the velocity, namely $v_{\Vert }$ and $v_{\bot }$, and rigorously show that kinetic energy isotropization necessarily emerges.

Mechanisms setting the density decay in the scrape-off layer (SOL) at the outer midplane of a tokamak plasma are disentangled using two-fluid numerical simulations in a double-null magnetic configuration and analytical estimates. Typical experimental observations are retrieved, in particular increasing intermittency of the turbulence going from the near to the far SOL, which is reflected in two different density decay lengths. The decay length of the near SOL is well described as the result of transport driven by a nonlinearly saturated ballooning instability, while in the far SOL, the density decay length is described using a model of intermittent transport mediated by blobs. The analytical estimates of the decay lengths agree well with the simulation results and typical experimental values and can therefore be used to guide tokamak design and operation.

Plasma terminating disruptions in tokamaks may result in relativistic runaway electron beams with potentially serious consequences for future devices with large plasma currents. In this paper, we investigate the effect of plasma elongation on the coupled dynamics of runaway generation and resistive diffusion of the electric field. We find that elongated plasmas are less likely to produce large runaway currents, partly due to the lower induced electric fields associated with larger plasmas, and partly due to direct shaping effects, which mainly lead to a reduction in the runaway avalanche gain.

In toroidal magnetic confinement devices, such as tokamaks and stellarators, neoclassical transport is usually an order of magnitude larger than its classical counterpart. However, when a high-collisionality species is present in a stellarator optimized for low Pfirsch–Schlüter current, its classical transport can be comparable to the neoclassical transport. In this letter, we compare neoclassical and classical fluxes and transport coefficients calculated for Wendelstein 7-X (W7-X) and Large Helical Device (LHD) cases. In W7-X, we find that the classical transport of a collisional impurity is comparable to the neoclassical transport for all radii, while it is negligible in the LHD cases, except in the vicinity of radii where the neoclassical transport changes sign. In the LHD case, electrostatic potential variations on the flux surface significantly enhance the neoclassical impurity transport, while the classical transport is largely insensitive to this effect in the cases studied.

We study the damping of collisionless Alfvénic turbulence in a strongly magnetised plasma by two mechanisms: stochastic heating (whose efficiency depends on the local turbulence amplitude $\unicode[STIX]{x1D6FF}z_{\unicode[STIX]{x1D706}}$) and linear Landau damping (whose efficiency is independent of $\unicode[STIX]{x1D6FF}z_{\unicode[STIX]{x1D706}}$), describing in detail how they affect and are affected by intermittency. The overall efficiency of linear Landau damping is not affected by intermittency in critically balanced turbulence, while stochastic heating is much more efficient in the presence of intermittent turbulence. Moreover, stochastic heating leads to a drop in the scale-dependent kurtosis over a narrow range of scales around the ion gyroscale.