Laser-plasma instabilities (LPI) play a detrimental role in energy coupling to the target in inertial confinement fusion (ICF). The recent development of applied strong magnetic fields for use in ICF and laboratory astrophysics experiments has opened opportunities to investigate the role of external magnetic fields on LPIs. Recent numerical studies have shown that stimulated Raman scattering (SRS) can be mitigated by external magnetic fields in the kinetic regime of the instability and warrant systematic experimental studies to validate modelling. To this end, we design experiments at the OMEGA-EP laser facility to investigate the effect of an external perpendicular $B$-field of 5–30 T on the backscattered light spectrum from a gas-jet target. We present measurements from a proof-of-principle experiment, where the backward-SRS (BSRS) is in the kinetic regime, for which the magnetization is expected to produce the greatest effects on instability growth. New simulations motivated by this experiment are used to inform the proposal of an upgraded experimental design. Our simulation predictions show that the new design is suited to experimentally demonstrating BSRS mitigation by an external magnetic field, despite the magnetization effects on the hydrodynamics, as well as the inherent temporal and spatial variations in plasma conditions.

A collisionless shock is a self-organized structure where fields and particle distributions are mutually adjusted to ensure a stable mass, momentum and energy transfer from the upstream to the downstream region. This adjustment may involve rippling, reformation or whatever else is needed to maintain the shock. The fields inside the shock front are produced due to the motion of charged particles, which is in turn governed by the fields. The overshoot arises due to the deceleration of the ion flow by the increasing magnetic field, so that the drop of the dynamic pressure should be compensated by the increase of the magnetic pressure. The role of the overshoot is to regulate ion reflection, thus properly adjusting the downstream ion temperature and kinetic pressure and also speeding up the collisionless relaxation and reducing the anisotropy of the eventually gyrotropized distributions.

The transport phenomena of dust particles have been widely observed in fusion plasmas. In this article, we report the observations of dust fragmentations in the Experimental Advanced Superconducting Tokamak (EAST). A dust particle splits into two daughter particles and their motions are recorded before and after the breakup with a fast video camera. The trajectories of the daughter particles in the experiment are consistent with equation-of-motion simulations. The stability of a rotating charged particle in the plasma is briefly discussed.

In this article, we study nonlinear waves propagating along the background magnetic field in relativistic electron–positron plasmas. Using the reductive perturbation method, we derive a three-dimensional equation describing these waves. When the perturbations do not vary in the directions orthogonal to the background magnetic field this equation reduces to the vector modified Kortewed–de Vries equation. We present solutions of the obtained equation in the form of planar solitary waves and describe the results of study of their stability with respect to transverse perturbations. We also study numerically non-planar solitary waves.

The problem of the electron dynamics on a closed magnetic field line passing through a high-$Z$ plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the $H$ theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide insights into the behaviour of the disruption dynamics for representative scenarios. The simulations are conducted using the numerical framework Dream (Disruption Runaway Electron Analysis Model). We show that, irrespective of the quantities of the material deposition, multi-megaampere runaway currents will be produced in the deuterium–tritium phase of operations, even in the optimal scenarios. However, the severity of the outcome can be influenced by tailoring the radial profile of the injected material; in particular, if the injected neon is deposited at the edge region it leads to a significant reduction of both the final runaway current and the transported heat losses. The Bayesian approach allows us to map the parameter space efficiently, with more accuracy in favourable parameter regions, thereby providing us with information about the robustness of the optima.

An analytical theory of wave turbulence is developed for pure compressible magnetohydrodynamics in the small $\beta$ limit. In contrast to previous works where the multiple scale method was not mentioned and slow magneto-acoustic waves were included, we present here a theory for fast magneto-acoustic waves for which only an asymptotic closure is possible in three dimensions. We introduce the compressible Elsässer fields (canonical variables) and show their linear relationship with the mass density and the compressible velocity. The kinetic equations of wave turbulence for three-wave interactions are obtained and the detailed conservation is shown for the two invariants, energy and momentum (cross-helicity). An exact stationary solution (Kolmogorov-Zakharov spectrum) exists only for the energy. We find a $k^{-3/2}$ energy spectrum compatible with the Iroshnikov–Kraichnan (IK) phenomenological prediction; this leads to a mass density spectrum with the same scaling. Despite the presence of a relatively strong uniform magnetic field, this turbulence is characterized by an energy spectrum with a power index that is independent of the angular direction; its amplitude, however, shows an angular dependence. We prove the existence of the IK solution using the locality condition, show that the energy flux is positive and hence the cascade direct and find the Kolmogorov constant. This theory offers a plausible explanation for recent observations in the solar wind at small $\beta$ where isotropic spectra with a $-3/2$ power-law index are found and associated with fast magneto-acoustic waves. This theory may also be used to explain the IK spectrum often observed near the Sun. Besides, it provides a rigorous theoretical basis for the well-known phenomenological IK spectrum, which coincides with the Zakharov–Sagdeev spectrum for acoustic wave turbulence.

A new method named B-TIS (Bourgeois & Davoine, J. Comput. Phys., vol. 413, 2020, 109426) has recently been proposed for suppressing the influence of numerical Cherenkov radiation that appears in particle-in-cell (PIC) simulation of laser wakefield acceleration (LWFA). However, while this method provides good results when applied to the already accelerated electrons, we show here that it cannot model correctly most of the plasma electron bulk interacting with the laser field. We thus investigate in this paper the origins of this limitation and propose an improved method for which this limitation is removed. This new method, named B-TIS3, can now be applied to a much broader variety of problems and improve the performance in comparison with the standard PIC algorithm. We show that, for an electron interacting directly with a laser pulse, this new technique offers greater accuracy in terms of momentum and motion than the conventional scheme used in many PIC codes. These improvements translate into more faithful energy spectrum and electric charge for the accelerated beam in simulations of vacuum laser acceleration (VLA) or LWFA involving direct laser acceleration (DLA) at low plasma density. This new method, easy to implement and not computationally demanding, should then prove useful to study in depth and help develop novel VLA, DLA and LWFA techniques.

In the present scheme of work, the Hermite–Gaussian (HG) laser beam dynamics has been investigated under the influence of an upward density ramp inside magnetized plasma, where both relativistic and ponderomotive nonlinearities are operative. One can achieve self-focusing of laser beam due to the change in the medium's dielectric function, which comes into operation due to the expulsion of plasma electrons from the high intensity to the low-intensity region by ponderomotive force and their motion at relativistic speeds. The dynamics of the laser beam and terahertz generation have been investigated by using the moment theory approach. It has been observed from the present analysis that the dynamics of the laser beam and the production of terahertz radiations strongly depends upon the HG laser beam and plasma parameters. In addition to this, the effect of density ramp and magnetic field has also been investigated on the efficiency of terahertz generation. It has been observed that higher-order modes of the HG laser beam play a dominant role in the production of terahertz radiations.

We study the parametric decay instability of parallel-propagating Alfvén waves in a low-beta plasma using one-dimensional fully kinetic simulations. We focus for the first time on the conversion of the energy stored in the initial Alfvén wave into particle internal energy, and on its partition between particle species. We show that compressible fluctuations generated by the decay of the pump wave into a secondary ion-acoustic mode and a reflected Alfvén wave contribute to the gain of internal energy via two distinct mechanisms. First, the ion-acoustic mode leads nonlinearly to proton trapping and proton phase-space mixing, in agreement with previous work based on hybrid simulations. Second, during the nonlinear stage, a compressible front of the fast type develops at the steepened edge of the backward Alfvén wave leading to a field-aligned proton beam propagating backwards at the Alfvén speed. We find that parametric decay heats preferentially protons, which gain approximately 50 % of the pump wave energy in the form of internal energy. However, we find that electrons are also energized and that they contribute to the total energy balance by gaining 10 % of the pump wave energy. By investigating energy partition and particle heating during parametric decay, our results contribute to the determination of the role of compressible and kinetic effects in wave-driven models of the solar wind.

A high-intensity laser pulse propagating through a gas target disturbs the uniform plasma distribution. Plasma density structures, created by high-order Bessel–Gauss beams for guiding the accelerating Gaussian beam and laser wakefield acceleration of electrons, are analysed using Wake-T and Fourier–Bessel particle-in-cell (FBPIC) simulation tools. The use of Bessel–Gauss doughnut beams increases the acceleration distance and energy of accelerated electrons up to 2.3 times at a 2 mm distance relative to the Gaussian beam of the same intensity.

Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.

The application scopes of two different reductive perturbation methods to derive the Korteweg–de Vries (KdV) equation and coupled KdV (CKdV) equation in two-temperature-ion dusty plasma are given by using the particle-in-cell (PIC) numerical method in the present paper. It suggests that the reductive perturbation method (RPM) is valid if the amplitude of the CKdV solitary wave is small enough. However, for the KdV solitary wave, RPM is valid not only if the amplitude of the KdV solitary wave is small enough, but also if the nonlinear coefficient of the KdV equation is not tending to zero.

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.

Starting from the equations of collisionless linear gyrokinetics for magnetised plasmas with an imposed inhomogeneous magnetic field, we present the first known analytical, closed-form solution for the resulting velocity-space integrals in the presence of resonances due to both parallel streaming and constant magnetic drifts. These integrals are written in terms of the well-known plasma dispersion function (Faddeeva & Terent'ev, Tables of Values of the Function $w(z)=\exp (-z^2)(1+2i/\sqrt {\pi }\int _0^z \exp (t^2) \,\mathrm {d} t)$ for Complex Argument, 1954. Gostekhizdat. English translation: Pergamon Press, 1961; Fried & Conte, The Plasma Dispersion Function, 1961. Academic Press), rendering the subsequent expressions simpler to treat analytically and more efficient to compute numerically. We demonstrate that our results converge to the well-known ones in the straight-magnetic-field and two-dimensional limits, and show good agreement with the numerical solver by Gürcan (J. Comput. Phys., vol. 269, 2014, p. 156). By way of example, we calculate the exact dispersion relation for a simple electrostatic, ion-temperature-gradient-driven instability, and compare it with approximate kinetic and fluid models.

The first nonlinear gyrokinetic simulations obtained using a moment approach based on the Hermite–Laguerre decomposition of the distribution function are presented, implementing advanced models for the collision operator. Turbulence in a two-dimensional Z-pinch is considered within a flux-tube configuration. In the collisionless regime, our gyromoment approach shows very good agreement with nonlinear simulations carried out with the continuum gyrokinetic code GENE, even with fewer gyromoments than required for the convergence of the linear growth rate. By using advanced linear collision operators, the role of collisions in setting the level of turbulent transport is then analysed. The choice of collision operator model is shown to have a crucial impact when turbulence is quenched by the presence of zonal flows. The convergence properties of the gyromoment approach improve when collisions are included.

We present a numerical method to solve the Vlasov–Maxwell equations for spin-1/2 particles, in a semiclassical approximation where the orbital motion is treated classically while the spin variable is fully quantum. Unlike the spinless case, the phase-space distribution function is a $2\times 2$ matrix, which can also be represented, in the Pauli basis, as one scalar function $f_0$ and one three-component vector function $\boldsymbol f$. The relationship between this ‘vectorial’ representation and the fully scalar representation on an extended phase space first proposed by Brodin et al. (Phys. Rev. Lett., vol. 101, 2008, p. 245002) is analysed in detail. By means of suitable approximations and symmetries, the vectorial spin-Vlasov–Maxwell model can be reduced to two-dimensions in the phase space, which is amenable to numerical solutions using a high-order grid-based Eulerian method. The vectorial model enjoys a Poisson structure that paves the way to accurate Hamiltonian split-time integrators. As an example, we study the stimulated Raman scattering of an electromagnetic wave interacting with an underdense plasma, and compare the results with those obtained earlier with the scalar spin-Vlasov–Maxwell model and a particle-in-cell code.

The explanation of the mechanisms of anomalous resistance and diffusion of plasma, as well as the generation of fast particles, is given from a unified position on the basis of experimental results obtained in plasma in a magnetically isolated diode of an electron accelerator with a microsecond pulse duration. It is shown that oscillations of the potential and current in plasma with current occur by charge separation and the formation of electric domains with a strong field. Electric domains are generated by means of runaway electrons resulting from the redistribution of the current density and its channelling in the current-carrying plasma. An increase in plasma resistance occurs due to a decrease in the number of electrons in the current-conducting state, since some of the plasma electrons pass into layers of excess negative charge of electric domains. The speed of movement of electrons in the composition of the electric domain is much lower than the velocity of electrons in the conducting state. The nucleation of domains is accompanied by microwave emission. Quasi-neutrals in whole electrical domains are injected into the region of the anode wall. The destruction of domains leads to the appearance of plasma channels between the anode and cathode plasma. The transverse electromagnetic waves generated during the nucleation of domains capture charged particles and inject them in a direction that is perpendicular to the insulating longitudinal magnetic field. Plasma turbulence occurs due to charge separation and the generation of electrical domains that create channels between the plasma and the chamber wall, as well as generating fluxes of fast particles.

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.