We have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We find that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can eventually trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq en_{\text{CR}}v_{\text{sh}}$, where $en_{\text{CR}}$ is the seed charge density and $v_{\text{sh}}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than ${\sim}10~\text{yr}$, thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.

Starting from expressions in Connor et al. (Phys. Fluids, vol. 31, 1988, p. 577), we derive a one-dimensional tearing equation similar to the approximate equation obtained by Hegna & Callen (Phys. Plasmas, vol. 1, 1994, p. 2308) and Nishimura et al. (Phys. Plasmas, vol. 5, 1998, p. 4292), but for more realistic toroidal equilibria. The intention is to use this approximation to explore the role of steep profiles, bootstrap currents and strong shaping in the vicinity of a separatrix, on the stability of tearing modes which are resonant in the H-mode pedestal region of finite aspect ratio, shaped cross-section tokamaks, e.g. the Joint European Torus (JET). We discuss how this one-dimensional model for tearing modes, which assumes a single poloidal harmonic for the perturbed poloidal flux, compares with a model that includes poloidal coupling Fitzpatrick et al. (Nucl. Fusion, vol. 33, 1993, p. 1533).

The homogeneous dynamo effect is at the root of magnetic field generation in cosmic bodies, including planets, stars and galaxies. While the underlying theory had increasingly flourished since the middle of the 20th century, hydromagnetic dynamos were not realized in the laboratory until 1999. On 11 November 1999, this situation changed with the first observation of a kinematic dynamo in the Riga experiment. Since that time, a series of experimental campaigns has provided a wealth of data on the kinematic and the saturated regime. This paper is intended to give a comprehensive survey about these experiments, to summarize their main results and to compare them with numerical simulations.

Mean field electrodynamics (MFE) facilitates practical modelling of secular, large scale properties of astrophysical or laboratory systems with fluctuations. Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages. Often however, real systems do not exhibit such scale separation. This raises two questions: (I) What are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) How precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel. Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise. The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.

We discuss a mean-field theory of the generation of large-scale vorticity in a rotating density stratified developed turbulence with inhomogeneous kinetic helicity. We show that the large-scale non-uniform flow is produced due to either a combined action of a density stratified rotating turbulence and uniform kinetic helicity or a combined effect of a rotating incompressible turbulence and inhomogeneous kinetic helicity. These effects result in the formation of a large-scale shear, and in turn its interaction with the small-scale turbulence causes an excitation of the large-scale instability (known as a vorticity dynamo) due to a combined effect of the large-scale shear and Reynolds stress-induced generation of the mean vorticity. The latter is due to the effect of large-scale shear on the Reynolds stress. A fast rotation suppresses this large-scale instability.

In most numerical simulations of the Earth’s core the dynamo is located outside the tangent cylinder and, in a zero-order sense, takes the form of a classical $\unicode[STIX]{x1D6FC}^{2}$ dynamo. Such a dynamo usually requires a distribution of helicity, $h$, which is asymmetric about the equator and in the simulations it is observed that, outside the tangent cylinder, the helicity is predominantly negative in the north and positive in the south. If we are to extrapolate the results of these simulations to the planets, we must understand how this asymmetry in helicity is established and ask if the same mechanism is likely to operate in a planet. In some of the early numerical dynamos, which were too viscous by a factor of at least $10^{9}$, as measured by the Ekman number, the asymmetric helicity distribution was attributed to Ekman pumping. However, Ekman pumping plays much less of a role in more recent, and less viscous, numerical dynamos, and almost certainly plays no significant role in the core of a planet. So the question remains: what establishes the asymmetric helicity distribution in the simulations and is this mechanism likely to carry over to planetary cores? In this paper we review the evidence that planetary dynamos, and their numerical analogues, might be maintained by helical waves, especially inertial waves, excited in and around the equatorial regions. This cartoon arises from the observation that there tends to be a statistical bias in the buoyancy flux towards the equatorial regions, and so waves are preferentially excited there. Moreover, upward (downward) propagating inertial waves carry negative (positive) helicity, which leads naturally to a segregation in $h$.

Using a variational optimization method we find the smallest flow perturbations that can trigger kinematic dynamo action in Kolmogorov flow. In comparison to previous work, a second-order mean field dynamo model is used to track down the optimal dynamos in the high magnetic Reynolds number limit ($Rm$). The magnitude of minimal perturbation flows decays inversely proportional to the magnetic Reynolds number. We reveal the asymptotic high-$Rm$ structure of the optimal flow perturbation and the magnetic eigenmode. We identify the optimal dynamo as of $\unicode[STIX]{x1D6FC}{-}\unicode[STIX]{x1D6FA}$ type, with magnetic fluctuations that localize on a critical layer.

We apply a nonlinear mean-field dynamo model which includes a budget equation for the dynamics of Wolf numbers to predict solar activity. This dynamo model takes into account the algebraic and dynamic nonlinearities of the $\unicode[STIX]{x1D6FC}$ effect, where the equation for the dynamic nonlinearity is derived from the conservation law for the magnetic helicity. The budget equation for the evolution of the Wolf number is based on a formation mechanism of sunspots related to the negative effective magnetic pressure instability. This instability redistributes the magnetic flux produced by the mean-field dynamo. To predict solar activity on the time scale of one month we use a method based on a combination of the numerical solution of the nonlinear mean-field dynamo equations and the artificial neural network. A comparison of the results of the prediction of the solar activity with the observed Wolf numbers demonstrates a good agreement between the forecast and observations.

Mean-field dynamo equations are addressed with the aid of the path integral method. The evolution of magnetic field is treated as a three-dimensional Wiener random process, and the mean magnetic-field equations are obtained with the Wiener integrals taken over all the trajectories of the fluid particles. The form of the equations is just the same as the conventional mean-field equations, but here the equations are derived with the velocity field realisation affected by the force exerted by the magnetic field. In this sense, we derive nonlinear dynamo equations.

Stellar radiative zones are typically assumed to be motionless in standard models of stellar structure but there is sound theoretical and observational evidence that this cannot be the case. We investigate by direct numerical simulations a three-dimensional and time-dependent model of stellar radiation zones consisting of an electrically conductive and stably stratified anelastic fluid confined to a rotating spherical shell and driven by a baroclinic torque. As the baroclinic driving is gradually increased a sequence of transitions from an axisymmetric and equatorially symmetric time-independent flow to flows with a strong poloidal component and lesser symmetry are found. It is shown that all flow regimes characterised by significant non-axisymmetric components are capable of generating a self-sustained magnetic field. As the value of the Prandtl number is decreased and the value of the Ekman number is decreased, flows become strongly time dependent with progressively complex spatial structure and dynamos can be generated at lower values of the magnetic Prandtl number.

Hamiltonian extended magnetohydrodynamics (XMHD) is restricted to respect helical symmetry by reducing the Poisson bracket for the three-dimensional dynamics to a helically symmetric one, as an extension of the previous study for translationally symmetric XMHD (Kaltsas et al., Phys. Plasmas, vol. 24, 2017, 092504). Four families of Casimir invariants are obtained directly from the symmetric Poisson bracket and they are used to construct Energy–Casimir variational principles for deriving generalized XMHD equilibrium equations with arbitrary macroscopic flows. The system is then cast into the form of Grad–Shafranov–Bernoulli equilibrium equations. The axisymmetric and the translationally symmetric formulations can be retrieved as geometric reductions of the helically symmetric one. As special cases, the derivation of the corresponding equilibrium equations for incompressible plasmas is discussed and the helically symmetric equilibrium equations for the Hall MHD system are obtained upon neglecting electron inertia. An example of an incompressible double-Beltrami equilibrium is presented in connection with a magnetic configuration having non-planar helical magnetic axis.

This paper addresses one aspect of the problem of the suppression of tearing mode magnetic islands by electron cyclotron current drive (ECCD) injection, formulating the problem as the converse of a forced reconnection problem. New physical conditions are discussed which should be considered in the technical approach towards a robust control strategy. Limits on the ECCD deposition are determined to avoid driving the system into regimes where secondary instabilities develop. Numerical simulations confirming the theory are also presented.

Kinetic ballooning modes in magnetically confined toroidal plasmas are investigated putting emphasis on specific stellarator features. In particular, we propose a Mercier criterion which is purposely designed to allow for direct comparison with local flux-tube gyrokinetics simulations. We investigate the influence on the marginal frequency of the mode of a magnetic curvature which is inhomogeneous on the magnetic flux surface due to the fieldline-label dependence. This is a typical (surface) global effect present in non-axisymmetry. Finally, we propose an artificial equilibrium model that explicitly retains the fieldline-label dependence in the magnetic drift, and analyse the stability of the system by introducing a representation of the perturbations similar to the flux-bundle model of Sugama et al. (Plasma Fusion Res., vol. 7, 2012, 2403094). The coupling of flux bundles is shown to have a stabilising effect on the most unstable local flux-tube mode.

First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in a pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as ‘ABC fields’ that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density and very consistent profiles of the scalar product $\boldsymbol{E}\boldsymbol{\cdot }\boldsymbol{B}$ of electric field $\boldsymbol{E}$ and magnetic field $\boldsymbol{B}$. We decompose the profiles of $\boldsymbol{E}\boldsymbol{\cdot }\boldsymbol{B}$ with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.

Transport of test particles in two-dimensional weak turbulence with waves propagating along the poloidal direction is studied using a reduced model. Finite Larmor radius (FLR) effects are included by gyroaveraging over one particle orbit. For low wave amplitudes the motion is mostly regular with particles trapped in the potential wells. As the amplitude increases the trajectories become chaotic and the Larmor radius modifies the orbits. For a thermal distribution of Finite Larmor radii the particle distribution function (PDF) is Gaussian for small $\unicode[STIX]{x1D70C}_{th}$ (thermal gyroradius) but becomes non-Gaussian for large $\unicode[STIX]{x1D70C}_{th}$. However, the time scaling of transport is diffusive, as characterized by a linear dependence of the variance of the PDF with time. An explanation for this behaviour is presented that provides an expression for an effective diffusion coefficient and reproduces the numerical results for large wave amplitudes which implies generalized chaos. When a shear flow is added in the direction of wave propagation, a modified model is obtained that produces free-streaming particle trajectories in addition to trapped ones; these contribute to ballistic transport for low wave amplitude but produce super-ballistic transport in the chaotic regime. As in the previous case, the PDF is Gaussian for low $\unicode[STIX]{x1D70C}_{th}$ becoming non-Gaussian as it increases. The perpendicular transport presents the same behaviour as in the case with no flow but the diffusion is faster in the presence of the flow.

The current filamentation instability (CFI) is capable of generating strong magnetic fields relevant to the explanation of radiation processes in astrophysical objects and leads to the onset of particle acceleration in collisionless shocks. Probing such extreme scenarios in the laboratory is still an open challenge. In this work, we investigate the possibility of using neutral $e^{-}~e^{+}$ beams to explore the CFI with realistic parameters, by performing two-dimensional particle-in-cell simulations. We show that CFI can occur unless the rate at which the beam expands due to finite beam emittance is larger than the CFI growth rate and as long as the role of the competing electrostatic two-stream instability (TSI) is negligible. We also show that the longitudinal energy spread, typical of plasma-based accelerated electron–positron fireball beams, plays a minor role in the growth of CFI in these scenarios.

A new configuration for an antenna, which increases the absorption power and plasma density, is proposed for helicon plasma sources. The influence of the electromagnetic wave pattern symmetry on the plasma density and absorption power in a helicon plasma source with a common antenna (Nagoya) is analysed by using the standard COMSOL Multiphysics 5.3 software. In contrast to the theoretical model prediction, the electromagnetic wave does not represent a symmetric pattern for the common Nagoya antenna. In this work, a new configuration for an antenna is proposed which refines the asymmetries of the wave pattern in helicon plasma sources. The plasma parameters such as plasma density and absorption rate for a common Nagoya antenna and our proposed antenna under the same conditions are studied using simulations. In addition, the plasma density of seven operational helicon plasma source devices, having a common Nagoya antenna, is compared with the simulation results of our proposed antenna and the common Nagoya antenna. The simulation results show that the density of the plasma, which is produced by using our proposed antenna, is approximately twice in comparison to the plasma density produced by using the common Nagoya antenna. In fact, the simulation results indicate that the electric and magnetic fields symmetry of the helicon wave plays a vital role in increasing wave–particle coupling. As a result, wave–particle energy exchange and the plasma density of helicon plasma sources will be increased.

In this paper, we determine the growth rate of the exponential radiation amplification in the ion-channel laser, where a relativistic electron beam wiggles in a focusing ion channel that can be created in a wakefield accelerator. For the first time the radiation diffraction, which can limit the amplification, is taken into account. The electron beam quality requirements to obtain this amplification are also presented. It is shown that both the beam energy and wiggler parameter spreads should be limited. Two-dimensional and three-dimensional particle-in-cell simulations of the self-consistent ion-channel laser confirm our theoretical predictions.

A heat flux in a high-$\unicode[STIX]{x1D6FD}$ plasma with low collisionality triggers the whistler instability. Quasilinear theory predicts saturation of the instability in a marginal state characterized by a heat flux that is fully controlled by electron scattering off magnetic perturbations. This marginal heat flux does not depend on the temperature gradient and scales as $1/\unicode[STIX]{x1D6FD}$. We confirm this theoretical prediction by performing numerical particle-in-cell simulations of the instability. We further calculate the saturation level of magnetic perturbations and the electron scattering rate as functions of $\unicode[STIX]{x1D6FD}$ and the temperature gradient to identify the saturation mechanism as quasilinear. Suppression of the heat flux is caused by oblique whistlers with magnetic-energy density distributed over a wide range of propagation angles. This result can be applied to high-$\unicode[STIX]{x1D6FD}$ astrophysical plasmas, such as the intracluster medium, where thermal conduction at sharp temperature gradients along magnetic-field lines can be significantly suppressed. We provide a convenient expression for the amount of suppression of the heat flux relative to the classical Spitzer value as a function of the temperature gradient and $\unicode[STIX]{x1D6FD}$. For a turbulent plasma, the additional independent suppression by the mirror instability is capable of producing large total suppression factors (several tens in galaxy clusters) in regions with strong temperature gradients.

Observational evidence in space and astrophysical plasmas with a long collisional mean free path suggests that more massive charged particles may be preferentially heated. One possible mechanism for this is the turbulent cascade of energy from injection to dissipation scales, where the energy is converted to heat. Here we consider a simple system consisting of a magnetized plasma slab of electrons and a single ion species with a cross-field density gradient. We show that such a system is subject to an electron drift wave instability, known as the universal instability, which is stabilized only when the electron and ion thermal speeds are equal. For unequal thermal speeds, we find from quasilinear analysis and nonlinear simulations that the instability gives rise to turbulent energy exchange between ions and electrons that acts to equalize the thermal speeds. Consequently, this turbulent heating tends to equalize the component temperatures of pair plasmas and to heat ions to much higher temperatures than electrons for conventional mass-ratio plasmas.