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.

Following indirect-drive experiments which demonstrated promising performance for low convergence ratios (below 17), previous direct-drive simulations identified a fusion-relevant regime which is expected to be robust to hydrodynamic instability growth. This paper expands these results with simulated implosions at lower energies of 100 and 270 kJ, and ‘hydrodynamic equivalent’ capsules which demonstrate comparable convergence ratio, implosion velocity and in-flight aspect ratio without the need for cryogenic cooling, which would allow the assumptions of one-dimensional-like performance to be tested on current facilities. A range of techniques to improve performance within this regime are then investigated, including the use of two-colour and deep ultraviolet laser pulses. Finally, further simulations demonstrate that the deposition of electron energy into the hotspot of a low convergence ratio implosion through auxiliary heating also leads to significant increases in yield. Results include break even for 1.1 MJ of total energy input (including an estimated 370 kJ of short-pulse laser energy to produce electron beams for the auxiliary heating), but are found to be highly dependent upon the efficiency with which electron beams can be created and transported to the hotspot to drive the heating mechanism.

An optimized compact stellarator with four simple coils is obtained from direct optimization via a coil shape. The new stellarator consists of two interlocking coils and two vertical field coils similar to those of the Columbia Non-neutral Torus (CNT) (Pedersen et al., Phys. Rev. Lett., vol. 88, 2002, pp. 205002). The optimized configuration has a global magnetic well and a low helical ripple level comparable to that of Wendelstein 7-X (W7-X) (Wolf et al., Nucl. Fusion, vol. 57, 2017, pp. 102020). The two interlocking coils have a smooth three-dimensional shape much simpler than those of advanced stellarators such as W7-X. This result opens up possibilities of future stellarator reactors with simplified coils.

We study the time-dependent formation and evolution of a current sheet (CS) in a magnetised, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field increases, and adiabatic invariance in the inflowing fluid elements produces a field-biased pressure anisotropy with excess perpendicular pressure. At large plasma beta, this anisotropy excites the mirror instability, which deforms the reconnecting field on ion-Larmor scales and dramatically reduces the effective thickness of the CS. Tearing modes whose wavelengths are comparable to that of the mirrors then become unstable, triggering reconnection on smaller scales and at earlier times than would have occurred if the thinning CS were to have retained its Harris profile. A novel method for identifying and tracking X-points is introduced, yielding X-point separations that are initially intermediate between the perpendicular and parallel mirror wavelengths in the upstream plasma. These mirror-stimulated tearing modes ultimately grow and merge to produce island widths comparable to the CS thickness, an outcome we verify across a range of CS formation timescales and initial CS widths. Our results may find their most immediate application in the tearing disruption of magnetic folds generated by turbulent dynamo in weakly collisional, high-beta, astrophysical plasmas.

We employ the $\lambda ^{3}$ regime where a near-single-cycle laser pulse is tightly focused, thus providing the highest possible intensity for the minimal energy at a certain laser power. The quantum electrodynamics processes in the course of the interaction of an ultra-intense laser with a solid target are studied via three-dimensional particle-in-cell simulations, revealing the generation of copious $\gamma$-photons and electron–positron pairs. A parametric study of the laser polarisation, target thickness and electron number density shows that a radially polarised laser provides the optimal regime for $\gamma$-photon generation. By varying the laser power in the range of 1 to 300 PW we find the scaling of the laser to $\gamma$-photon energy conversion efficiency. The laser-generated $\gamma$-photon interaction with a high-$Z$ target is further studied using Monte Carlo simulations revealing further electron–positron pair generation and radioactive nuclide creation.

The paper presents experimental results from the SMOLA device that is the first facility with a helical mirror section of the magnetic system. This device was built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea that is the technique of an active control of axial losses from a confinement zone. Theory predicts that, with rotating plasma, a helical mirror will provide suppression of the axial plasma flow and, simultaneously, density pinching to the axis. Experiments demonstrated the increase in plasma density in the entrance trap by a factor of 1.6 in the helical configuration. The integral axial flux from the transport section drops severalfold. The effective mirror ratio of the helical section was $R_{eff} > 10$. Particle flux returning by the helical mirror section towards the confinement zone was observed. At high corrugation ratios, the axial flux direction is different at the magnetic axis and in the periphery of the plasma in the helical section. All axial fluxes scale linearly with the plasma density, even if the ion mean free path is comparable to the total length of the helical section. Good agreement of the experimental results with theoretical predictions is found.

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.

The GENE-3D code, the global stellarator version of the established GENE framework, has been extended to an electromagnetic gyrokinetic code. This paper outlines the basic structure of the algorithm, highlighting the treatment of the electromagnetic terms. The numerical implementation is verified against the radially global GENE code in linear and nonlinear tokamak simulations, recovering excellent agreement between both codes. As a first application to stellarator plasmas, linear and nonlinear global simulations with kinetic electrons of ion temperature gradient (ITG) turbulence in Wendelstein 7-X were performed, showing a decrease of ITG activity through the introduction of electromagnetic effects via a finite plasma-$\beta$. The upgrade makes it possible to study a large variety of new physical scenarios, including kinetic electron and electromagnetic effects, reducing the gap between gyrokinetic models and physically realistic systems.

Radio frequency (RF) sheaths occur under a wide variety of conditions when RF waves, material surfaces and plasma coexist. RF sheaths are of special importance in describing the interaction of ion cyclotron range of frequency (ICRF) waves with the boundary plasma in tokamaks, stellarators and other magnetic confinement devices. In this article the basic physics of RF sheaths is discussed in the context of magnetic fusion research. Techniques for modelling RF sheaths, their interaction with RF wave fields and the resulting consequences are highlighted. The article is intended as a guide for the early-career ICRF researcher, but it may equally well serve to provide an overview of basic RF sheath concepts and modelling directions for any interested fusion scientist.

Electron energization during merging of magnetized plasmas is studied using the OMEGA and OMEGA EP laser facilities by colliding two plasma plumes, each containing a Biermann-battery self-generated magnetic field. Two neighbouring plasma plumes are produced by intense laser beams, and the anti-parallel Biermann fields merge and reconnect in the process of the plumes’ expansion and collision. To isolate the merging as an acceleration source, the electron energy spectra obtained from two-plume collision shots are compared with the spectra from single-plume shots. Single-plume shots exhibit an energized electron tail with energies up to ${\sim }250\ \textrm {keV}$. The electrons in merging experiments are additionally accelerated by ${\sim }50\text {--}100$ keV compared to single-plume shots.

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the $\boldsymbol {E\times B}$ rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the $\boldsymbol {E\times B}$ velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

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.