Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul (Phys. Rev. Lett., vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$. An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$. Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. (Phys. Plasma, vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$.

The main emphasis of the Laser Ion Generation, Handling and Transport (LIGHT) beamline at GSI Helmholtzzentrum für Schwerionenforschung GmbH are phase-space manipulations of laser-generated ion beams. In recent years, the LIGHT collaboration has successfully generated and focused intense proton bunches with an energy of 8 MeV and a temporal duration shorter than 1 ns (FWHM). An interesting area of application that exploits the short ion bunch properties of LIGHT is the study of ion-stopping power in plasmas, a key process in inertial confinement fusion for understanding energy deposition in dense plasmas. The most challenging regime is found when the projectile velocity closely approaches the thermal plasma electron velocity ($v_{i}\approx v_{e,\text {th}}$), for which existing theories show high discrepancies. Since conclusive experimental data are scarce in this regime, we plan to conduct experiments on laser-generated plasma probed with ions generated with LIGHT at a higher temporal resolution than previously achievable. The high temporal resolution is important because the parameters of laser-generated plasmas are changing on the nanosecond time scale. To meet this goal, our recent studies have dealt with ions of lower kinetic energies. In 2021, laser accelerated carbon ions were transported with two solenoids and focused temporally with LIGHT's radio frequency cavity. A bunch length of 1.2 ns (FWHM) at an energy of 0.6 MeV u$^{-1}$ was achieved. In 2022, protons with an energy of 0.6 MeV were transported and temporally compressed to a bunch length of 0.8 ns. The proton beam was used to measure the energy loss in a cold foil. Both the ion and proton beams will also be employed for energy loss measurements in a plasma target.

We present a stellarator configuration optimized for a large threshold (‘critical gradient’) for the onset of the ion temperature gradient (ITG) driven mode, which achieves the largest critical gradient we have seen in any stellarator. Above this threshold, gyrokinetic simulations show that the configuration has low turbulence levels over an experimentally relevant range of the drive strength. The applied optimization seeks to maximize the drift curvature, leading to enhanced local-shear stabilization of toroidal ITG modes, and the associated turbulence. These benefits are combined with excellent quasi-symmetry, yielding low neoclassical transport and vanishingly small alpha particle losses. Analysis of the resulting configuration suggests a trade-off between magnetohydrodynamic (MHD) and ITG stability, as the new configuration possesses a vacuum magnetic hill.

Laser-driven ion acceleration is well established using solid targets mainly in the target normal sheath acceleration regime. To follow the increasing repetition rate available on high-intensity lasers, the use of high-density gas targets has been explored in the past decade. When interacting with targets reaching densities close to the critical one, the laser pulse can trigger different acceleration mechanisms such as Collisionless Shock Acceleration (CSA) or hole boring. Particle-in-cell simulations using ideal target profiles show that CSA can accelerate a collimated, narrow energy spread and few hundreds of megaelectronvolts ion beam on the laser axis. Nevertheless, in real experiments, the laser will not only interact with an overcritical, thin plasma slab with sharp density gradients, but also with lower density regions surrounding the core of the gas jet, extending to several hundreds of micrometres. The interaction of the laser with these lower density wings will lead to nonlinear effects that will reduce the available energy to drive the shock in the high-density region of the target. Optically tailoring this target could mitigate that issue. Recent experiments conducted on different laser facilities aimed at testing several tailoring configurations. We first tested a scheme with a copropagating picosecond prepulse to create a lower density plasma channel to facilitate the propagation of the main pulse, while the second one was a transverse tailoring driven by nanosecond laser pulses to generate blast waves and form a high-density plasma slab. The main results will be presented here and the methods compared.

In the limit of sufficiently fast rotation, rotating mirror traps are known to be stable against the loss-cone modes associated with conventional (non-rotating) mirrors. This paper calculates how quickly a mirror configuration must rotate in order for several of these modes to be stabilized (in particular, the high-frequency convective loss cone, drift cyclotron loss cone and Dory–Guest–Harris modes). Commonalities in the stabilization conditions for these modes then motivate a modified formulation of the Gardner free energy and diffusively accessible free energy to be used for systems in which the important modes have wavevectors that are orthogonal or nearly orthogonal to the magnetic field, as well as a modification to include the effects of a loss region in phase space.

This paper establishes a mean-field equation set and an energy theorem to provide a theoretical basis in view of the development of self-consistent, physics-based turbulent transport models for mean-field transport codes. A rigorous averaging procedure identifies the exact form of the perpendicular turbulent fluxes which are modelled by ad hoc diffusive terms in mean-field transport codes, next to other closure terms which are not commonly considered. Earlier work suggested that the turbulent $E\times B$ particle and heat fluxes, which are thus identified to be important closure terms, can be modelled to reasonable accuracy using the kinetic energy in the $E\times B$ velocity fluctuations ($k_{E}$). The related enstrophy led to further modelling improvements in an initial study, although further analysis is required. To support this modelling approach, transport equations are derived analytically for both quantities. In particular, an energy theorem is established in which the various source and sink terms of $k_{E}$ are shown to couple to mean-field and turbulent parallel kinetic energy, kinetic energy in the other perpendicular velocity components, the thermal energy and the magnetic energy. This provides expressions for the interchange, drift-wave and Reynolds stress terms amongst others. Note that most terms in these energy equations are in turn closure terms. It is suggested to evaluate these terms using reference data from detailed turbulence code simulations in future work.

Generalising the Elsässer variables, we introduce the $Q$-variables. These are more flexible than the Elsässer variables, because they also allow us to track waves with phase speeds different than the Alfvén speed. We rewrite the magnetohydrodynamics (MHD) equations with these $Q$-variables. We consider also the linearised version of the resulting MHD equations in a uniform plasma, and recover the classical Alfvén waves, but also separate the fast and slow magnetosonic waves into upward- and downward-propagating waves. Moreover, we show that the $Q$-variables may also track the upward- and downward-propagating surface Alfvén waves in a non-uniform plasma, displaying the power of our generalisation. In the end, we lay the mathematical framework for driving solar wind models with a multitude of wave drivers.

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high-$\beta$ collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high-$\beta$ means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

When a particle crosses a region of space where the curvature radius of the magnetic field line shrinks below its gyroradius $r_{g}$, it experiences a non-adiabatic (magnetic moment violating) change in pitch angle. The present paper carries that observation into magnetohydrodynamic (MHD) turbulence to examine the influence of intermittent, sharp bends of the magnetic field lines on particle transport. On the basis of dedicated measurements in a simulation of incompressible turbulence, it is argued that regions of sufficiently large curvature exist in sufficient numbers on all scales to promote scattering. The parallel mean free path predicted by the power-law statistics of the curvature strength scales in proportion to $r_{g}^{0.3}\,\ell _{c}^{0.7}$ ($\ell _{c}$ is the coherence scale of the turbulence), which is of direct interest for cosmic-ray phenomenology. Particle tracking in that numerical simulation confirms that the magnetic moment diffuses through localized, violent interactions, in agreement with the above picture. Correspondingly, the overall transport process is non-Brownian up to length scales $\gtrsim \ell _{c}$.

The OMFIT STEP (Meneghini et al., Nucl. Fusion, vol. 10, 2020, p. 1088) workflow has been used to develop inductive and steady-state H-mode core plasma scenario use cases for a $B_0 = 8 \, {\rm T}$, $R_0 = 4 \, {\rm m}$ machine to help guide and inform future higher-fidelity studies of core transport and confinement in compact tokamak reactors. Both use cases are designed to produce 200 MW or more of net electric power in an up-down symmetric plasma with minor radius $a = 1.4 \, {\rm m}$, elongation $\kappa = 2.0$, triangularity $\delta = 0.5$ and effective charge $Z_{{\rm eff}} \simeq 2$. Additional considerations based on the need for compatibility of the core with reactor-relevant power exhaust solutions and external actuators were used to guide and constrain the use case development. An extensive characterization of core transport in both scenarios is presented, the most important feature of which is the extreme sensitivity of the results to the quantitative stiffness level of the transport model used as well as the predicted critical gradients. This sensitivity is shown to arise from different levels of transport stiffness exhibited by the models, combined with the gyroBohm-normalized fluxes of the predictions being an order of magnitude larger than other H-mode plasmas. Additionally, it is shown that although heating in both plasmas is predominantly to the electrons and collisionality is low, the plasmas remain sufficiently well coupled for the ions to carry a significant fraction of the thermal transport. As neoclassical transport is negligible in these conditions, this situation inherently requires long-wavelength ion gyroradius-scale turbulence to be the dominant transport mechanism in both plasmas. These results are combined with other basic considerations to propose a simple heuristic model of transport in reactor-relevant plasmas, along with simple metrics to quantify coupling and core transport properties across burning and non-burning plasmas.

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

We summarize recent attempts to unravel the role of plasma kinetic effects in radiation mediated shocks. Such shocks form in all strong stellar explosions and are responsible for the early electromagnetic emission released from these events. A key issue that has been overlooked in all previous works is the nature of the coupling between the charged leptons, that mediate the radiation force, and the ions, which are the dominant carriers of the shock energy. Our preliminary investigation indicates that in the case of relativistic shocks, as well as Newtonian shocks in multi-ion plasma, this coupling is driven by either, transverse magnetic fields of a sufficiently magnetized upstream medium, or plasma microturbulence if strong enough magnetic fields are absent. We discuss the implications for the shock breakout signal, as well as abundance evolution and kilonova emission in binary neutron star mergers.

The oscillating relativistic electrons in the accelerating/focusing wakefields of plasma accelerators emit electromagnetic radiation known as betatron radiation (BR). The proton-driven plasma wakefield acceleration has been demonstrated in the Advanced Wakefield Experiment (AWAKE) at CERN; however, its accompanying radiation emission is less explored compared with those in the laser- and electron beam-driven plasma accelerators. In this paper, a detailed simulation study of BR in the AWAKE is presented. Considering the new set-up of the AWAKE Run 2 (2021–), the radiation emission from both the witness electron beam and the seeding electron beam is investigated using particle-in-cell simulations. The influence of radial size mismatch and misaligned off-axis injection on the witness beam dynamics, as well as the spectral features of the relevant BR are studied. These non-ideal electron injections are likely to occur in experiment. The proton self-modulation stage is also investigated with a close look at the seeding electron beam dynamics and its BR. As a footprint of the emitting particles, BR can provide valuable information about the beam dynamics. Some practical challenges to implement the betatron diagnostic in the AWAKE Run 2 experiment are also addressed.

Global electromagnetic turbulence is simulated in stellarator geometry using the gyrokinetic particle-in-cell code EUTERPE. The evolution of the turbulent electromagnetic field and the plasma profiles is considered at different values of the plasma beta and for different magnetic configurations. It is found that turbulence is linearly driven at relatively high toroidal mode numbers. In the nonlinear regime, lower toroidal mode numbers, including zonal flows, are excited resulting in a quench of the linear instability drive. The turbulent heat flux is outward and leads to the nonlinear relaxation of the plasma temperature profile. The particle flux is inward for the parameters considered. The effect of the parallel perturbation of the magnetic field on the stellarator turbulence is addressed.

How much does the cross-section of a toroidal magnetic field configuration tell us about its magnetohydrodynamic (MHD) stability? It is generally believed that positive triangularity (typically leading to bean-shaped cross-sections with their indentation on the inboard side in stellarators) contributes positively to MHD stability. In this paper, we explore the basis of this statement within a near-axis description for axisymmetric and quasisymmetric magnetic configurations. In agreement with the existing literature, we show that positive triangularity stabilises vertically elongated tokamaks. In quasisymmetric stellarators, the toroidal asymmetry of flux surfaces modifies this relation. The behaviour of stellarator-symmetric, quasisymmetric stellarators can still be described in terms of the shape of one of their up–down symmetric cross-sections. However, we show that for a sample of quasisymmetric configurations, the positive-bean-shaped cross-sections do not contribute positively to stability. Unlike in the axisymmetric case, we also learn that finite $\beta$ can improve stability even without magnetic shear.

We present a general method for optimizing the configuration of an experimental diagnostic to minimize uncertainty and bias in inferred quantities from experimental data. The method relies on Bayesian inference to sample the posterior using a physical model of the experiment and instrument. The mean squared error (MSE) of posterior samples relative to true values obtained from a high fidelity model (HFM) across multiple configurations is used as the optimization metric. The method is demonstrated on a common problem in dense plasma research, the use of radiation detectors to estimate physical properties of the plasma. We optimize a set of filtered photoconducting diamond detectors to minimize the MSE in the inferred X-ray spectrum, from which we can derive quantities like the electron temperature. In the optimization we self-consistently account for uncertainties in the instrument response with appropriate prior probabilities. We also develop a penalty term, acting as a soft constraint on the optimization, to produce results that avoid negative instrumental effects. We show results of the optimization and compare with two other reference instrument configurations to demonstrate the improvement. The MSE with respect to the total inferred X-ray spectrum is reduced by more than an order of magnitude using our optimized configuration compared with the two reference cases. We also extract multiple other quantities from the inference and compare with the HFM, showing an overall improvement in multiple inferred quantities like the electron temperature, the peak in the X-ray spectrum and the total radiated energy.

The paper presents experimental results from the SMOLA device on the testing of the helical mirror confinement hypothesis. Helical mirror confinement is the technique of an active control of axial plasma losses from a confinement zone by multiple magnetic mirrors that move along the axis in the reference frame of the plasma that experiences $\boldsymbol{E} \times \boldsymbol{B}$ rotation due to an applied radial electric field. Theory predicts that a helical mirror will provide an axial force that modifies the plasma flow and, simultaneously, density pinching to the axis. The force direction depends on the plasma rotation direction. Experimental data on the axial plasma losses at different direction of the magnetic mirror movement are presented. If the trapped ions move in the direction opposite to the direction of the axial losses, then the particle flux reduces in the broad range of the plasma density. The confinement improves with the increase of the fraction of the trapped particles (effective mirror ratio was up to $R_{{\rm eff}}=5.8\pm 1.4$). If the trapped ions move in the same direction as the axial losses, then the flux depends on density. At intermediate densities, the integral flux through the transport section rises compared to the plasma flowing through the straight magnetic field. The effective mirror ratio is lower and does not significantly depend on the fraction of the trapped particles (effective mirror ratio at intermediate density was $R_{{\rm eff}}=3.3\pm 0.8$).

Plasma supports collective modes and particle–wave interactions that lead to complex behaviour in, for example, inertial fusion energy applications. While plasma can sometimes be modelled as a charged fluid, a kinetic description is often crucial for studying nonlinear effects in the higher-dimensional momentum–position phase space that describes the full complexity of the plasma dynamics. We create a differentiable solver for the three-dimensional partial-differential equation describing the plasma kinetics and introduce a domain-specific objective function. Using this framework, we perform gradient-based optimization of neural networks that provide forcing function parameters to the differentiable solver given a set of initial conditions. We apply this to an inertial-fusion-relevant configuration and find that the optimization process exploits a novel physical effect.

The effect of sheared E × B flow on the blob dynamics in the scrape-off layer (SOL) of HL-2A tokamak has been studied during the plasma current ramp-up in ohmically heated deuterium plasmas by the combination of poloidal and radial Langmuir probe arrays. The experimental results indicate that the SOL sheared E × B flow is substantially enhanced as the plasma current exceeds a certain value and the strong sheared E × B flow has the ability to slow the blob radial motion via stretching its poloidal correlation length. The locally accumulated blobs are suggested to be responsible for the increase of plasma density just outside the Last Closed Flux Surface (LCFS) observed in this experiment. The results presented here reveal the significant role played by the strong sheared E × B flow on the blob dynamics, which provides a potential method to control the SOL width by modifying the sheared E × B flow in future tokamak plasmas.

The energy functional of ideal magnetohydrodynamics (MHD) is extended to the regime of finite resistivity and transformed into an intuitive form, which separates pressure-gradient, current-density drive and stabilising phenomena, analogously to the ideal version of the functional derived by Greene and Johnson (Plasma Phys., vol. 10, issue 8, 1968, pp. 729–745). The resulting resistive energy functional is not in the form of a solvable eigenvalue problem but is intended to analyse the energetic composition of eigenfunctions which were calculated in advance by solving the MHD equations. The new resistive energy functional provides information on the different energetic drives and stabilising phenomena of an eigenfunction. The new energy functional has been implemented in the CASTOR3D code in general 3D curvilinear coordinates and successfully validated across three different coordinate systems for a numerical test case.