We study the longitudinal stability of beam–plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam–plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts.

The creation of well-thermalized, hot and dense plasmas is attractive for warm dense matter studies. We investigate collisionally induced energy absorption of an ultraintense and ultrashort laser pulse in a solid copper target using particle-in-cell simulations. We find that, upon irradiation by a $2\times 10^{20}~\text{W}\,\text{cm}^{-2}$ intensity, 60 fs duration, circularly polarized laser pulse, the electrons in the collisional simulation rapidly reach a well-thermalized distribution with ${\sim}3.5~\text{keV}$ temperature, while in the collisionless simulation the absorption is several orders of magnitude weaker. Circular polarization inhibits the generation of suprathermal electrons, while ensuring efficient bulk heating through inverse bremsstrahlung, a mechanism usually overlooked at relativistic laser intensity. An additional simulation, taking account of both collisional and field ionization, yields similar results: the bulk electrons are heated to ${\sim}2.5~\text{keV}$, but with a somewhat lower degree of thermalization than in the pre-set, fixed-ionization case. The collisional absorption mechanism is found to be robust against variations in the laser parameters. At fixed laser pulse energy, increasing the pulse duration rather than the intensity leads to a higher electron temperature.

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. It is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

Accelerated statistical computation of collisionless fusion alpha particle losses in stellarator configurations is presented based on direct guiding-centre orbit tracing. The approach relies on the combination of recently developed symplectic integrators in canonicalized magnetic flux coordinates and early classification into regular and chaotic orbit types. Only chaotic orbits have to be traced up to the end, as their behaviour is unpredictable. An implementation of this technique is provided in the code SIMPLE (symplectic integration methods for particle loss estimation, Albert et al., 2020b, doi:10.5281/zenodo.3666820). Reliable results were obtained for an ensemble of 1000 orbits in a quasi-isodynamic, a quasi-helical and a quasi-axisymmetric configuration. Overall, a computational speed up of approximately one order of magnitude is achieved compared to direct integration via adaptive Runge–Kutta methods. This reduces run times to the range of typical magnetic equilibrium computations and makes direct alpha particle loss computation adequate for use within a stellarator optimization loop.

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

Density functional theory (DFT) is a powerful theoretical tool widely used in such diverse fields as computational condensed-matter physics, atomic physics and quantum chemistry. DFT establishes that a system of $N$ interacting electrons can be described uniquely by its single-particle density $n(\boldsymbol{r})$, instead of the $N$-body wave function, yielding an enormous gain in terms of computational speed and memory storage space. Here, we use time-dependent DFT to show that a classical collisionless plasma can always, in principle, be described by a set of fluid equations for the single-particle density and current. The results of DFT guarantee that an exact closure relation, fully reproducing the Vlasov dynamics, necessarily exists, although it may be complicated (non-local in space and time, for instance) and difficult to obtain in practice. This goes against the common wisdom in plasma physics that the Vlasov and fluid descriptions are mutually incompatible, with the latter inevitably missing some ‘purely kinetic’ effects.

The evolution of a linearly polarized, long-wavelength Alfvén wave propagating in a collisionless magnetized plasma with a sheared parallel-directed velocity flow is here studied by means of two-dimensional hybrid Vlasov–Maxwell (HVM) simulations. The unperturbed sheared flow has been represented by an exact solution of the HVM set of equations of (Malara et al., Phys. Rev. E, vol. 97, 2018, 053212), thus avoiding spurious oscillations that would arise from the non-stationary initial state and inevitably affect the dynamics of the system. We have considered the evolution of both a small and a moderate amplitude Alfvén wave, in order to separate linear wave–shear flow couplings from kinetic effects, the latter being more relevant for larger wave amplitudes. The phase mixing generated by the shear flow modifies the initial perturbation, leading to the formation of small-scale transverse fluctuations at scales comparable with the proton inertial length/Larmor radius. By analysing both the polarization and group velocity of perturbations in the shear regions, we identify them as kinetic Alfvén waves (KAWs). In the moderate amplitude run, kinetic effects distort the proton distribution function in the shear region. This leads to the formation of a proton beam, at the Alfvén speed and parallel to the magnetic field. Such a feature, due to the parallel electric field associated with KAWs, positively compares with solar wind observations of suprathermal ion populations, suggesting that it may be related to the presence of ion-scale KAW-like fluctuations.

The concept of the available energy of a collisionless plasma is discussed in the context of magnetic confinement. The available energy quantifies how much of the plasma energy can be converted into fluctuations (including nonlinear ones) and is thus a measure of plasma stability, which can be used to derive linear and nonlinear stability criteria without solving an eigenvalue problem. In a magnetically confined plasma, the available energy is determined by the density and temperature profiles as well as the magnetic geometry. It also depends on what constraints limit the possible forms of plasma motion, such as the conservation of adiabatic invariants and the requirement that the transport be ambipolar. A general method based on Lagrange multipliers is devised to incorporate such constraints in the calculation of the available energy, and several particular cases are discussed for which it can be calculated explicitly. In particular, it is shown that it is impossible to confine a plasma in a Maxwellian ground state relative to perturbations with frequencies exceeding the ion bounce frequency.

The classical collision cross-sections of a proton with an argon atom as well as the thermal transport coefficients and rate constant of the colliding $\text{H}^{+}-\text{Ar}$ system are evaluated at the kinetic temperature $T\in [100,10\,000]~(\text{K})$ by means of the asymptotically correct analytical potential constructed for the ground $X^{1}\unicode[STIX]{x1D6F4}^{+}$ state of the ArH+ cation from the highly accurate ab initio data available in the entire range of internuclear distances (Terashkevich et al., J. Quant. Spectrosc. Radiat. Transfer, vol. 234, 2019, pp. 139–146). The results can be useful to estimate thermodynamic, transport and kinetic properties of the Ar/H2 plasma in a wide temperature range.

The dependence of confinement on input power for a tokamak plasma with regions having a stiff temperature profile is explored. The resilience of the confinement of the core energy to increasing power loss by core radiation from impurities in such situations, as it is anticipated will be required in a demonstration fusion reactor (DEMO) design, is examined.

The problem of plasma expansion into a vacuum is revisited with the addition of a finite boundary condition; an electrically insulated surface. As plasma expands towards a charge-accumulating surface, the leading electron cloud charges the surface negatively, which in turn repels electrons and attracts ions. This plasma–surface interaction is shown to result in a feedback process which accelerates the plasma expansion. In addition, we examine the decrease in (negative) surface potential and associated near-surface electron density. To investigate this plasma coupling with an electrically floating surface, we develop an analytic model including four neighbouring plasma regions: (i) undisturbed plasma, (ii) quasi-neutral self-similar expansion, (iii) ion front boundary layer and (iv) electron cloud. A key innovation in our approach is a self-contained analytic approximation of the ion front boundary layer, providing a spatially continuous electric field model for the early phase of bounded plasma expansion.

A gyrokinetic model is presented that can properly describe large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and the arbitrarily large deviations from thermal equilibrium that are present in the plasma periphery of tokamak devices. The formulation of the gyrokinetic model is based on a second-order accurate description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behaviour of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function. The collisional effects are included by a nonlinear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a Hermite–Laguerre velocity-space polynomial basis. Then, the gyrokinetic equation is projected onto the same basis yielding the spatial and temporal evolution of the Hermite–Laguerre expansion coefficients. A closed set of fluid equations for the lowest-order coefficients is presented. The Hermite–Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwell equations derived from a variational principle where the velocity integrals are explicitly evaluated.

The effects of the electron energy distribution function (EEDF) on the structure of a dusty plasma sheath are investigated. Here, it is assumed that the electrons obey a Druyvesteyn-type distribution with a parameter $x$ controlling the shape of the EEDF. The Druyvesteyn-like distribution tends to a Maxwellian distribution as $x$ varies from 2 to 1. Using the orbital motion limited theory, the incident electron current on the dust is evaluated for a given $x$. The results of numerical simulations are compared with those of a Maxwellian distribution. It was found that the sheath dynamics depends strongly on the magnitude of $x$. The sheath thickness increases monotonically with increasing $x$. However, the absolute dust charge decreases and, as a result, the accelerating ion drag force is weakened and thus the dust number density is enhanced. For a plasma with a Druyvesteyn-like distribution, the Bohm speed is a function of $x$ and increases with increasing $x$.

In this work, we highlight an issue that may reduce the accuracy of many local nonlinear gyrokinetic simulations – turbulent self-interaction through the parallel boundary condition. Given a sufficiently long parallel correlation length, individual turbulent eddies can span the full domain and ‘bite their own tails’, thereby altering their statistical properties. Such self-interaction is only modelled accurately when the simulation domain corresponds to a full flux surface, otherwise it is artificially strong. For Cyclone Base Case parameters and typical domain sizes, we find that this mechanism modifies the heat flux by approximately 40 % and it can be even more important. The effect is largest when using kinetic electrons, low magnetic shear and strong turbulence drive (i.e. steep background gradients). It is found that parallel self-interaction can be eliminated by increasing the parallel length and/or the binormal width of the simulation domain until convergence is achieved.

Linear gyrokinetic simulations of magnetically confined electron–positron plasmas are performed for the first time in the geometry and parameter regimes likely to be relevant for upcoming laboratory experiments. In such plasmas, the density will be sufficiently small as to render the plasma effectively collisionless. The magnetic field will be very large, meaning that the Debye length will exceed the gyroradius by a few orders of magnitude. We show the results of linear simulations in flux tubes close to the current carrying ring and also in the bulk of the plasma, demonstrating the existence of entropy modes and interchange modes in pair plasmas. We study linear stability and show that in the relevant configurations, almost complete linear stability is attainable in large swathes of parameter space.

Understanding particle drifts in a non-symmetric magnetic field is of primary interest in designing optimized stellarators in order to minimize the neoclassical radial loss of particles. Quasisymmetry and omnigeneity, two distinct properties proposed to ensure radial localization of collisionless trapped particles in stellarators, have been explored almost exclusively for magnetic fields with nested flux surfaces. In this work, we examine radial particle confinement when all field lines are closed. We then study charged particle dynamics in the special case of a non-symmetric vacuum magnetic field with closed field lines obtained recently by Weitzner & Sengupta (Phys. Plasmas, vol. 27, 2020, 022509). These magnetic fields can be used to construct magnetohydrodynamic equilibria for low pressure. Expanding in the amplitude of the non-symmetric fields, we explicitly evaluate the omnigeneity and quasisymmetry constraints. We show that the magnetic field is omnigeneous in the sense that the drift surfaces coincide with the pressure surfaces. However, it is not quasisymmetric according to the standard definitions.

We find a new family of solutions for force-free magnetic structures in cylindrical geometry. These solutions have radial power-law dependence and are periodic but non-harmonic in the azimuthal direction; they generalize the vacuum $z$-independent potential fields to current-carrying configurations.

Magnetic field quantization is an important issue for degenerate environments such as neutron stars, radio pulsars and magnetars etc., due to the fact that these stars have a magnetic field higher than the quantum critical field strength of the order of $4.4\times 10^{13}~\text{G}$, accordingly, the cyclotron energy may be equal to or even much more than the Fermi energy of degenerate particles. We shall formulate here the exotic physics of strongly magnetized neutron stars, known as pulsars, specifically focusing on the outcomes of the quantized magnetic pressure. In this scenario, while following the modified quantum hydrodynamic model, we shall investigate both linear and nonlinear fast magnetosonic waves in a strongly magnetized, weakly ionized degenerate plasma consisting of neutrons and an electron–ion plasma in the atmosphere of a pulsar. Here, linear analysis depicts that sufficiently long, fast magnetosonic waves may exist in a weakly dispersive pulsar having finite phase speed at cutoff. To investigate one-dimensional nonlinear fast magnetosonic waves, a neutron density expression as a function of both the electron magnetic and neutron degenerate pressures, is derived with the aid of Riemann’s wave solution. Consequently, a modified Korteweg–de Vries equation is derived, having a rarefractive solitary wave solution. It is found that the basic properties such as amplitude, width and phase speed of the fast magnetoacoustic waves are significantly altered by the electron magnetic and the neutron degenerate pressures. The results of this theoretical investigation may be useful for understanding the formation and features of the solitary structures in astrophysical compact objects such as pulsars, magnetars and white dwarfs etc.