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.

We consider the tilting instability of a magnetically confined spheromak using three-dimensional magnetohydrodynamic and relativistic particle-in-cell calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter-alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially, the spheromak rotates – tilts – trying to lower its magnetic potential energy. As a result, a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on a few Alfvén time scales. In the case of a higher-order (second-order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfvén waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetar giant flares, fast development of tilting instabilities and no stabilization of the higher-order spheromaks make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.

We apply field–particle correlations – a technique that tracks the time-averaged velocity-space structure of the energy density transfer rate between electromagnetic fields and plasma particles – to data drawn from a hybrid Vlasov–Maxwell simulation of Alfvén-ion cyclotron turbulence. Energy transfer in this system is expected to include both Landau and cyclotron wave–particle resonances, unlike previous systems to which the field–particle correlation technique has been applied. In this simulation, the energy transfer rate mediated by the parallel electric field $E_{\Vert }$ comprises approximately 60 % of the total rate, with the remainder mediated by the perpendicular electric field $E_{\bot }$. The parallel electric field resonantly couples to protons, with the canonical bipolar velocity-space signature of Landau damping identified at many points throughout the simulation. The energy transfer mediated by $E_{\bot }$ preferentially couples to particles with $v_{tp}\lesssim v_{\bot }\lesssim 3v_{tp}$, where $v_{tp}$ is the proton thermal speed, in agreement with the expected formation of a cyclotron diffusion plateau. Our results demonstrate clearly that the field–particle correlation technique can distinguish distinct channels of energy transfer using single-point measurements, even at points in which multiple channels act simultaneously, and can be used to determine quantitatively the rates of particle energization in each channel.

We present the elements required to construct two devices used in an undergraduate plasma physics laboratory. The materials and construction costs of the sources, the vacuum systems and probe drives and electrical circuits are presented in detail in the text and the first appendix. We also provide the software for probe motion and data acquisition as well as the electrical schematics for key components. Experiments which have been performed are listed and two (resonance cones and whistler waves) are described in greater detail. The machines are flexible and original research is possible.

Numerical simulations of Alfvén modes driven by energetic particles are performed with the gyrokinetic (GK) global particle-in-cell code ORB5. A reversed shear equilibrium magnetic field is adopted. A simplified configuration with circular flux surfaces and large aspect ratio is considered. The nonlinear saturation of beta-induced Alfvén eigenmodes (BAE) is investigated. The roles of the wave–particle nonlinearity of the different species, i.e. thermal ions, electrons and energetic ions are described, in particular for their role in the saturation of the BAE and the generation of zonal flows. The nonlinear redistribution of the electron population is found to be important in increasing the BAE saturation level and the zonal flow amplitude.