Between the base of the solar corona at $r=r_\textrm {b}$ and the Alfvén critical point at $r=r_\textrm {A}$, where $r$ is heliocentric distance, the solar-wind density decreases by a factor $ \mathop > \limits_\sim 10^5$, but the plasma temperature varies by a factor of only a few. In this paper, I show that such quasi-isothermal evolution out to $r=r_\textrm {A}$ is a generic property of outflows powered by reflection-driven Alfvén-wave (AW) turbulence, in which outward-propagating AWs partially reflect, and counter-propagating AWs interact to produce a cascade of fluctuation energy to small scales, which leads to turbulent heating. Approximating the sub-Alfvénic region as isothermal, I first present a brief, simplified calculation showing that in a solar or stellar wind powered by AW turbulence with minimal conductive losses, $\dot {M} \simeq P_\textrm {AW}(r_\textrm {b})/v_\textrm {esc}^2$, $U_{\infty } \simeq v_\textrm {esc}$, and $T\simeq m_\textrm {p} v_\textrm {esc}^2/[8 k_\textrm {B} \ln (v_\textrm {esc}/\delta v_\textrm {b})]$, where $\dot {M}$ is the mass outflow rate, $U_{\infty }$ is the asymptotic wind speed, $T$ is the coronal temperature, $v_\textrm {esc}$ is the escape velocity of the Sun, $\delta v_\textrm {b}$ is the fluctuating velocity at $r_\textrm {b}$, $P_\textrm {AW}$ is the power carried by outward-propagating AWs, $k_\textrm {B}$ is the Boltzmann constant, and $m_\textrm {p}$ is the proton mass. I then develop a more detailed model of the transition region, corona, and solar wind that accounts for the heat flux $q_\textrm {b}$ from the coronal base into the transition region and momentum deposition by AWs. I solve analytically for $q_\textrm {b}$ by balancing conductive heating against internal-energy losses from radiation, $p\,\textrm {d} V$ work, and advection within the transition region. The density at $r_\textrm {b}$ is determined by balancing turbulent heating and radiative cooling at $r_\textrm {b}$. I solve the equations of the model analytically in two different parameter regimes. In one of these regimes, the leading-order analytic solution reproduces the results of the aforementioned simplified calculation of $\dot {M}$, $U_\infty$, and $T$. Analytic and numerical solutions to the model equations match a number of observations.

Tokamak equilibria have been derived that are analytic solutions to the Grad–Shafranov equation. This paper, Part 1, describes a wide range of such equilibria including smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to include edge pedestals and toroidal flow.

We use three-dimensional (3-D) fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfvén waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic flux ropes as part of the turbulent cascade. The large magnetic structures present in the simulation domain retain the initial anisotropy while the small-scale structures produced by the turbulent cascade are less anisotropic. To quantify the occurrence of reconnection in our simulation domain, we develop a new set of indicators based on intensity thresholds to identify reconnection events in which both ions and electrons are heated and accelerated in 3-D particle-in-cell simulations. According to the application of these indicators, we identify the occurrence of reconnection events in the simulation domain and analyse one of these events in detail. The event is related to the reconnection of two flux ropes, and the associated ion and electron exhausts exhibit a complex 3-D structure. We study the profiles of plasma and magnetic-field fluctuations recorded along artificial-spacecraft trajectories passing near and through the reconnection region. Our results suggest the presence of particle heating and acceleration related to small-scale reconnection events within magnetic flux ropes produced by the anisotropic Alfvénic turbulent cascade in the solar wind. These events are related to current structures of the order of a few ion inertial lengths in size.

The formation of a substantial postdisruption runaway electron current in ASDEX Upgrade material injection experiments is determined by avalanche multiplication of a small seed population of runaway electrons. For the investigation of these scenarios, the runaway electron description of the coupled 1.5-D transport solvers ASTRA-STRAHL is amended by a fluid model describing electron runaway caused by the hot-tail mechanism. Applied in simulations of combined background plasma evolution, material injection and runaway electron generation in ASDEX Upgrade discharge #33108, both the Dreicer and hot-tail mechanism for electron runaway produce only ${\sim }$ 3 kA of runaway current. In colder plasmas with core electron temperatures $T_\textrm {e,c}$ below 9 keV, the postdisruption runaway current is predicted to be insensitive to the initial temperature, in agreement with experimental observations. Yet in hotter plasmas with $T_\textrm {e,c}$ above 10 keV, hot-tail runaway can be increased by up to an order of magnitude, contributing considerably to the total postdisruption runaway current. In ASDEX Upgrade high-temperature runaway experiments, however, no runaway current is observed at the end of the disruption, despite favourable conditions for both primary and secondary runaway.

A growing body of evidence suggests that the solar wind is powered to a large extent by an Alfvén-wave (AW) energy flux. AWs energize the solar wind via two mechanisms: heating and work. We use high-resolution direct numerical simulations of reflection-driven AW turbulence (RDAWT) in a fast-solar-wind stream emanating from a coronal hole to investigate both mechanisms. In particular, we compute the fraction of the AW power at the coronal base ($P_\textrm {AWb}$) that is transferred to solar-wind particles via heating between the coronal base and heliocentric distance $r$, which we denote by $\chi _{H}(r)$, and the fraction that is transferred via work, which we denote by $\chi _{W}(r)$. We find that $\chi _{W}(r_{A})$ ranges from 0.15 to 0.3, where $r_{A}$ is the Alfvén critical point. This value is small compared with one because the Alfvén speed $v_{A}$ exceeds the outflow velocity $U$ at $r < r_{A}$, so the AWs race through the plasma without doing much work. At $r>r_{A}$, where $v_{A} < U$, the AWs are in an approximate sense ‘stuck to the plasma’, which helps them do pressure work as the plasma expands. However, much of the AW power has dissipated by the time the AWs reach $r=r_{A}$, so the total rate at which AWs do work on the plasma at $r>r_{A}$ is a modest fraction of $P_\textrm {AWb}$. We find that heating is more effective than work at $r < r_{A}$, with $\chi _{H}(r_{A})$ ranging from 0.5 to 0.7. The reason that $\chi _{H} \geq 0.5$ in our simulations is that an appreciable fraction of the local AW power dissipates within each Alfvén-speed scale height in RDAWT, and there are a few Alfvén-speed scale heights between the coronal base and $r_{A}$. A given amount of heating produces more magnetic moment in regions of weaker magnetic field. Thus, paradoxically, the average proton magnetic moment increases robustly with increasing $r$ at $r>r_{A}$, even though the total rate at which AW energy is transferred to particles at $r>r_{A}$ is a small fraction of $P_\textrm {AWb}$.

Recent discoveries have demonstrated that matter can be distinguished on the basis of topological considerations, giving rise to the concept of topological phase. Introduced originally in condensed matter physics, the physics of topological phase can also be fruitfully applied to plasmas. Here, the theory of topological phase is introduced, including a discussion of Berry phase, Berry connection, Berry curvature and Chern number. One of the clear physical manifestations of topological phase is the bulk-boundary correspondence, the existence of localized unidirectional modes at the interface between topologically distinct phases. These concepts are illustrated through examples, including the simple magnetized cold plasma. An outlook is provided for future theoretical developments and possible applications.

This work aims to analyse the possibility of directly driven imploding spherical targets in order to create a source of energetic particles (neutrons, protons, alphas, tritium and 3He ions) for the Laser MégaJoule facility. D3He gas-filled spherical SiO2 glass pellets, irradiated by an absorbed laser intensity of 1014 W cm−2 or 1015 W cm−2 have been considered. Depending on the absorbed laser intensity and the amount of the ablated glass layer two distinct regimes have been identified: a massive pusher and an exploding pusher. Both regimes are analysed in terms of hydrodynamics and fast particle spectra. Energetic particle time-resolved spectra are calculated and used to infer ionic temperatures and total areal densities. A parametric study has been performed by varying the shell thickness and target inner radius for both laser absorbed intensities.

Magnetic wave perturbations are observed in the solar wind and in the vicinity of Earth's bow shock. For such environments, recent work on magnetic pumping with electrons trapped in the magnetic perturbations has demonstrated the possibility of efficient energization of superthermal electrons. Here we also analyse the energization of such energetic electrons for which the transit time through the system is short compared with time scales associated with the magnetic field evolution. In particular, considering an idealized magnetic configuration we show how trapping/detrapping of energetic magnetized electrons can cause effective parallel velocity ($v_{\parallel }$-) diffusion. This parallel diffusion, combined with naturally occurring mechanisms known to cause pitch angle scattering, such as whistler waves, produces enhanced heating rates for magnetic pumping. We find that at low pitch angle scattering rates, the combined mechanism enhances the heating beyond the predictions of the recent theory for magnetic pumping with trapped electrons.

Accelerated particles are ubiquitous in the Cosmos and play a fundamental role in many processes governing the evolution of the Universe at all scales, from the sub-AU scale relevant for the formation and evolution of stars and planets to the Mpc scale involved in Galaxy assembly. We reveal the presence of energetic particles in many classes of astrophysical sources thanks to their production of non-thermal radiation, and we detect them directly at the Earth as cosmic rays. In the last two decades both direct and indirect observations have provided us a wealth of new, high-quality data about cosmic rays and their interactions both in sources and during propagation, in the Galaxy and in the Solar System. Some of the new data have confirmed existing theories about particle acceleration and propagation and their interplay with the environment in which they occur. Some others have brought about interesting surprises, whose interpretation is not straightforward within the standard framework and may require a change of paradigm in terms of our ideas about the origin of cosmic rays of different species or in different energy ranges. In this article, we focus on cosmic rays of galactic origin, namely with energies below a few petaelectronvolts, where a steepening is observed in the spectrum of energetic particles detected at the Earth. We review the recent observational findings and the current status of the theory about the origin and propagation of galactic cosmic rays.

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$, where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$, $e\phi \sim T_{e}$, where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$, the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$. Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$, that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$, and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$, the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the $\boldsymbol {\nabla }\boldsymbol {B}$ and $\boldsymbol {\nabla }\boldsymbol {\nabla }\boldsymbol {B}$ tensors are computed, which can be used for direct derivative-based optimization of electromagnetic coil shapes to achieve the desired magnetic configuration. Moreover, if the norm of these tensors is large compared with the field strength for a given magnetic field, the field must have a short length scale, suggesting it may be hard to produce with coils that are suitably far away. Second, we evaluate the minor radius at which the flux surface shapes would become singular, providing a lower bound on the achievable aspect ratio. This bound is also shown to be related to an equilibrium beta limit. Finally, for configurations that are constructed to achieve a desired magnetic field strength to first order in the expansion, we compute the error field that arises due to second-order terms.

Protoplanetary discs are made of gas and dust orbiting a young star. They are also the birth place of planetary systems, which motivates a large amount of observational and theoretical research. In these lecture notes, I present a review of the magnetic mechanisms applied to the outer regions ($R\gtrsim 1\ \mathrm {AU}$) of these discs, which are the planet-formation regions. In contrast to usual astrophysical plasmas, the gas in these regions is noticeably cold ($T < 300\ \mathrm {K}$) and dense, which implies a very low ionisation fraction close to the disc midplane. In these notes, I deliberately ignore the innermost $(R\sim 0.1\ \mathrm {AU})$ region, which is influenced by the star–disc interaction and various radiative effects. I start by presenting a short overview of the observational evidence for the dynamics of these objects. I then introduce the methods and approximations used to model these plasmas, including non-ideal magnetohydrodynamics, and the uncertainties associated with this approach. In this framework, I explain how the global dynamics of these discs is modelled, and I present a stability analysis of this plasma in the local approximation, introducing the non-ideal magneto-rotational instability. Following this mostly analytical part, I discuss numerical models that have been used to describe the saturation mechanisms of this instability, and the formation of large-scale structures by various saturation mechanisms. Finally, I show that local numerical models are insufficient because magnetised winds are also emitted from the surface of these objects. After a short introduction on wind physics, I present global models of protoplanetary discs, including both a large-scale wind and the non-ideal dynamics of the disc.

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. The generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$, controllable by the coil currents, are observed in the simulations.

While it is well known that every nearly periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely defined scalar. A less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged-particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed.

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 rapidly emerging technology of high-temperature superconductors (HTS) opens new opportunities for the development of non-planar non-insulated HTS magnets. This type of HTS magnet offers attractive features via its simplicity and robustness, and is well suited for modest size steady-state applications such as a mid-scale stellarator. In non-planar coil applications the HTS tape may be subject to severe hard-way bending strain ($\epsilon _{\textrm {bend}}$), torsional strains ($\epsilon _{\textrm {tor}}$) and magnetic field components transverse to the HTS tape plane ($B_{\perp }$), all of which can limit the magnet operating space. A novel method of winding angle optimization is here presented to overcome these limitations for fixed input non-planar coil filamentary geometry. Essentially, this method: (i) calculates the peak $\epsilon _{\textrm {bend}}$ and $B_{\perp }$ for arbitrary winding angle along an input coil filamentary trajectory, (ii) defines a cost function including both and then (iii) uses tensioned splines to define a winding angle that reduces $\epsilon _{\textrm {tor}}$ and optimizes the $\epsilon _{\textrm {bend}}$ and $B_{\perp }$ cost function. As strain limits are present even without $B_{\perp }$, this optimization is able to provide an assessment of the minimum buildable size of an arbitrary non-planar non-insulating HTS coil. This optimization finds that for standard 4 mm wide HTS tapes the minimum size coils of the existing HSX, NCSX and W7-X stellarator geometries are around 0.3–0.5 m in mean coil radius. Identifying the minimum size provides a path to specify a mid-scale stellarator capable of achieving high-field or high-temperature operation with minimal HTS tape length. For coils larger than this size, strain optimization allows use of wider (higher current capacity) HTS tapes or alternatively permitting a finite (yet tolerable) strain allows reduction of $B_{\perp }$. Reduced $B_{\perp }$ enables a reduction of the HTS tape length required to achieve a given design magnetic field or equivalently an increase in the achievable magnetic field for fixed HTS tape length. The distinct considerations for optimizing a stellarator coilset to further ease compatibility with non-insulated HTS magnets are also discussed, highlighting relaxed curvature limits and the introduction of limits to the allowable torsion.

A Hamiltonian and action principle formalism for deriving three-dimensional gyroviscous magnetohydrodynamic models is presented. The uniqueness of the approach in constructing the gyroviscous tensor from first principles and its ability to explain the origin of the gyromap and the gyroviscous terms are highlighted. The procedure allows for the specification of free functions, which can be used to generate a wide range of gyroviscous models. Through the process of reduction, the noncanonical Hamiltonian bracket is obtained and briefly analysed.

The purpose of this paper is to bridge kinetic plasma descriptions and low-frequency single-fluid models. More specifically, the asymptotics leading to magnetohydrodynamic regimes starting from the Vlasov–Maxwell system are investigated. The analogy with the derivation, from the Vlasov–Poisson system, of a fluid representation for ions coupled to the Boltzmann relation for electrons is also outlined. The aim is to identify asymptotic parameters explaining the transitions from one microscopic description to a macroscopic low-frequency model. These investigations provide groundwork for the derivation of multi-scale numerical methods, model coupling or physics-based preconditioning.

The isotope dependence of plasma transport will have a significant impact on the performance of future D-T experiments in JET and ITER and eventually on the fusion gain and economics of future reactors. In preparation for future D-T operation on JET, dedicated experiments and comprehensive transport analyses were performed in H, D and H-D mixed plasmas. The analysis of the data has demonstrated an unexpectedly strong and favourable dependence of the global confinement of energy, momentum and particles in ELMy H-mode plasmas on the atomic mass of the main ion species, the energy confinement time scaling as ${\tau _E}\sim {A^{0.5}}$ (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045; JET Team, Nucl. Fusion, vol. 39, 1999, pp. 1227–1244), i.e. opposite to the expectations based only on local gyro-Bohm (GB) scaling, ${\tau _E}\sim {A^{ - 0.5}}$, and stronger than in the commonly used H-mode scaling for the energy confinement (Saibene et al., Nucl. Fusion, vol. 39, 1999, 1133; ITER Physics Basis, Nucl. Fusion, vol. 39, 1999, 2175). The scaling of momentum transport and particle confinement with isotope mass is very similar to that of energy transport. Nonlinear local GENE gyrokinetic analysis shows that the observed anti-GB heat flux is accounted for if collisions, E × B shear and plasma dilution with low-Z impurities (9Be) are included in the analysis (E and B are, respectively the electric and magnetic fields). For L-mode plasmas a weaker positive isotope scaling ${\tau _E}\sim {A^{0.14}}$ has been found in JET (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045), similar to ITER97-L scaling (Kaye et al., Nucl. Fusion, vol. 37, 1997, 1303). Flux-driven quasi-linear gyrofluid calculations using JETTO-TGLF in L-mode show that local GB scaling is not followed when stiff transport (as is generally the case for ion temperature gradient modes) is combined with an imposed boundary condition taken from the experiment, in this case predicting no isotope dependence. A dimensionless identity plasma pair in hydrogen and deuterium L-mode plasmas has demonstrated scale invariance, confirming that core transport physics is governed, as expected, by the 4 dimensionless parameters ρ*, ν*, β, q (normalised ion Larmor radius, collisionality, plasma pressure and safety factor) consistently with global quasi-linear gyrokinetic TGLF calculations (Maggi et al., Nucl. Fusion, vol. 59, 2019, 076028). We compare findings in JET with those in different devices and discuss the possible reasons for the different isotope scalings reported from different devices. The diversity of observations suggests that the differences may result not only from differences affecting the core, e.g. heating schemes, but are to a large part due to differences in device-specific edge and wall conditions, pointing to the importance of better understanding and controlling pedestal and edge processes.

It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria are discussed.