The Arbitrary Linear Plasma Solver (ALPS) is a parallelised numerical code that solves the dispersion relation in a hot (even relativistic) magnetised plasma with an arbitrary number of particle species with arbitrary gyrotropic equilibrium distribution functions for any direction of wave propagation with respect to the background field. ALPS reads the background momentum distributions as tables of values on a $(p_{\bot },p_{\Vert })$ grid, where $p_{\bot }$ and $p_{\Vert }$ are the momentum coordinates in the directions perpendicular and parallel to the background magnetic field, respectively. We present the mathematical and numerical approach used by ALPS and introduce our algorithms for the handling of poles and the analytic continuation for the Landau contour integral. We then show test calculations of dispersion relations for a selection of stable and unstable configurations in Maxwellian, bi-Maxwellian, $\unicode[STIX]{x1D705}$-distributed and Jüttner-distributed plasmas. These tests demonstrate that ALPS derives reliable plasma dispersion relations. ALPS will make it possible to determine the properties of waves and instabilities in the non-equilibrium plasmas that are frequently found in space, laboratory experiments and numerical simulations.

Azimuthal surface waves are eigenmodes of cylindrical plasma–dielectric–metal structures both in the presence of and without an axial static magnetic field. They are actively studied due to possible applications in plasma electronics, nanotechnologies and biomedical diagnostics. Higher radial modes are known to propagate at higher frequencies and shorter wavelengths compared to those of the zeroth mode, a feature which is of interest for practical applications. To gain the advantage of the excitation of higher radial modes of azimuthal surface waves one has first to know their dispersion properties. This paper generalizes the results of earlier papers by including a static axial magnetic field and considering the higher radial modes. The presence of the constant axial magnetic field removes the degeneracy in the wave spectrum with respect to the sign of the azimuthal wavenumber.

First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in a pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as ‘ABC fields’ that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density and very consistent profiles of the scalar product $\boldsymbol{E}\boldsymbol{\cdot }\boldsymbol{B}$ of electric field $\boldsymbol{E}$ and magnetic field $\boldsymbol{B}$. We decompose the profiles of $\boldsymbol{E}\boldsymbol{\cdot }\boldsymbol{B}$ with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.

In this paper we have studied the effect of the density and temperature of negative ions on the nonlinear dust-acoustic wave propagation in a Lorentzian dusty plasma. We have considered both adiabatic and non-adiabatic dust charge variation. The presence of both low and high populations of negative ions are considered. Separate models have been developed because the two populations give rise to opposite polarity of grain charges. In both models electrons are assumed to follow a kappa velocity distribution while the positive and negative ions satisfy a Maxwellian velocity distribution. Adiabatic dust charge variation shows the propagation of a dust-acoustic soliton in cases of both a high and low population of negative ions whose amplitude depends on the negative ion temperature and negative ion density. On the other hand, non-adiabatic dust charge variation generates a stable oscillatory dust-acoustic shock when the negative ion population is low. An unstable potential has been predicted from this analysis when the negative ion population is high and the dust charge variation is non-adiabatic.

A drift-kinetic calculation in an axisymmetric tokamak, with super-diamagnetic flows, is presented to elucidate the relation between the radial electric field, $E_{r}$, zonal flows and the rapid precession of the trapped particle (TP) population. It has been shown earlier (Rosenbluth & Hinton, Phys. Rev. Lett., vol. 80(4), 1998, p. 724, hereafter RH) that an initial radial electric field results in geodesic acoustic mode oscillations which subsequently Landau damp, resulting in a much smaller final residual electric field, and accompanying parallel zonal flows. We observe an apparent paradox: the final angular momentum in the RH parallel zonal flow is much smaller than the angular momentum expected from the well-known rapid precession of the trapped particle population in the RH residual electric field. We reconcile this paradox by illuminating the presence of a population of reverse circulating particle flows that, dominantly, are equal and opposite to the rapid TP precession. Mathematically, the calculation is facilitated by transforming to an energy coordinate shifted from conventional by an amount proportional to $E_{r}$. We also discuss the well-known RH coefficient in the context of effective mass and show how the TP precession and the opposite circulating flows contribute to this mass. Finally, we show that in the long wavelength limit, the RH flows and RH coefficient arise as a natural consequence of conservation of toroidal angular momentum and the second adiabatic invariant.

Magnetic field ripple is inherent in tokamaks since the toroidal magnetic field is generated by a finite number of toroidal field coils. The field ripple results in departures from axisymmetry that cause radial transport losses of particles and heat. These ripple losses are a serious concern for alphas near their birth speed $v_{0}$ since alpha heating of the background plasma is required to make fusion reactors into economical power plants. Ripple in tokamaks gives rise to at least two alpha transport regimes of concern. As the slowing down time $\unicode[STIX]{x1D70F}_{s}$ is much larger than the time for an alpha just born to make a toroidal transit, a regime referred to as the $1/\unicode[STIX]{x1D708}\propto \unicode[STIX]{x1D70F}_{s}$ regime can be encountered, with $\unicode[STIX]{x1D708}$ the appropriate alpha collision frequency. In this regime the radial transport losses increase as $v_{0}\unicode[STIX]{x1D70F}_{s}/R$, with $R$ the major radius of the tokamak. The deleterious effect of ripple transport is mitigated by electric and magnetic drifts within the flux surface. When drift tangent to the flux surface becomes significant another ripple regime, referred to as the $\sqrt{\unicode[STIX]{x1D708}}$ regime, is encountered where a collisional boundary layer due to the drift plays a key role. We evaluate the alpha transport in both regimes, taking account of the alphas having a slowing down rather than a Maxwellian distribution function and their being collisionally scattered by a collision operator appropriate for alphas. Alpha ripple transport is found to be in the $\sqrt{\unicode[STIX]{x1D708}}$ regime where it will be a serious issue for typical tokamak reactors as it will be well above the axisymmetric neoclassical level and can be large enough to deplete the alpha slowing down distribution function unless toroidal rotation is strong.

Three-dimensional direct numerical simulations are used to study the energy cascade rate in isothermal compressible magnetohydrodynamic turbulence. Our analysis is guided by a two-point exact law derived recently for this problem in which flux, source, hybrid and mixed terms are present. The relative importance of each term is studied for different initial subsonic Mach numbers $M_{S}$ and different magnetic guide fields $\boldsymbol{B}_{0}$ . The dominant contribution to the energy cascade rate comes from the compressible flux, which depends weakly on the magnetic guide field $\boldsymbol{B}_{0}$ , unlike the other terms whose moduli increase significantly with $M_{S}$ and $\boldsymbol{B}_{0}$ . In particular, for strong $\boldsymbol{B}_{0}$ the source and hybrid terms are dominant at small scales with almost the same amplitude but with a different sign. A statistical analysis undertaken with an isotropic decomposition based on the SO(3) rotation group is shown to generate spurious results in the presence of $\boldsymbol{B}_{0}$ , when compared with an axisymmetric decomposition better suited to the geometry of the problem. Our numerical results are compared with previous analyses made with in situ measurements in the solar wind and the terrestrial magnetosheath.

Transport of test particles in two-dimensional weak turbulence with waves propagating along the poloidal direction is studied using a reduced model. Finite Larmor radius (FLR) effects are included by gyroaveraging over one particle orbit. For low wave amplitudes the motion is mostly regular with particles trapped in the potential wells. As the amplitude increases the trajectories become chaotic and the Larmor radius modifies the orbits. For a thermal distribution of Finite Larmor radii the particle distribution function (PDF) is Gaussian for small $\unicode[STIX]{x1D70C}_{th}$ (thermal gyroradius) but becomes non-Gaussian for large $\unicode[STIX]{x1D70C}_{th}$. However, the time scaling of transport is diffusive, as characterized by a linear dependence of the variance of the PDF with time. An explanation for this behaviour is presented that provides an expression for an effective diffusion coefficient and reproduces the numerical results for large wave amplitudes which implies generalized chaos. When a shear flow is added in the direction of wave propagation, a modified model is obtained that produces free-streaming particle trajectories in addition to trapped ones; these contribute to ballistic transport for low wave amplitude but produce super-ballistic transport in the chaotic regime. As in the previous case, the PDF is Gaussian for low $\unicode[STIX]{x1D70C}_{th}$ becoming non-Gaussian as it increases. The perpendicular transport presents the same behaviour as in the case with no flow but the diffusion is faster in the presence of the flow.

Collisionless shocks follow the Rankine–Hugoniot jump conditions to a good approximation. However, for a shock propagating parallel to a magnetic field, magnetohydrodynamics states that the shock properties are independent of the field strength, whereas recent particle-in-cell simulations reveal a significant departure from magnetohydrodynamics behaviour for such shocks in the collisionless regime. This departure is found to be caused by a field-driven anisotropy in the downstream pressure, but the functional dependence of this anisotropy on the field strength is yet to be determined. Here, we present a non-relativistic model of the plasma evolution through the shock front, allowing for a derivation of the downstream anisotropy in terms of the field strength. Our scenario assumes double adiabatic evolution of a pair plasma through the shock front. As a result, the perpendicular temperature is conserved. If the resulting downstream is firehose stable, then the plasma remains in this state. If unstable, it migrates towards the firehose stability threshold. In both cases, the conservation equations, together with the relevant hypothesis made on the temperature, allows a full determination of the downstream anisotropy in terms of the field strength.

The confinement of heat in the core of a magnetic fusion reactor is optimised using a multidimensional optimisation algorithm. For the first time in such a study, the loss of heat due to turbulence is modelled at every stage using first-principles nonlinear simulations which accurately capture the turbulent cascade and large-scale zonal flows. The simulations utilise a novel approach, with gyrofluid treatment of the small-scale drift waves and gyrokinetic treatment of the large-scale zonal flows. A simple near-circular equilibrium with standard parameters is chosen as the initial condition. The figure of merit, fusion power per unit volume, is calculated, and then two control parameters, the elongation and triangularity of the outer flux surface, are varied, with the algorithm seeking to optimise the chosen figure of merit. A twofold increase in the plasma power per unit volume is achieved by moving to higher elongation and strongly negative triangularity.

We compute the real and imaginary parts of the electric permittivities and magnetic permeabilities of relativistic electrons from quantum electrodynamics at finite temperatures and densities, for weak fields, neglecting electron–electron interactions. For non-zero temperatures, electromagnetic responses are reduced to one-dimensional integrals computed numerically. For zero temperature, we find analytic expressions for both their real/dispersive and imaginary/absorptive parts. As an application of our results, we obtain the dispersion relation for longitudinal electric plasmons. Present calculations support our recent claim that, at low frequencies and long wavelengths, the system will exhibit simultaneously negative electric and magnetic responses.

The current filamentation instability (CFI) is capable of generating strong magnetic fields relevant to the explanation of radiation processes in astrophysical objects and leads to the onset of particle acceleration in collisionless shocks. Probing such extreme scenarios in the laboratory is still an open challenge. In this work, we investigate the possibility of using neutral $e^{-}~e^{+}$ beams to explore the CFI with realistic parameters, by performing two-dimensional particle-in-cell simulations. We show that CFI can occur unless the rate at which the beam expands due to finite beam emittance is larger than the CFI growth rate and as long as the role of the competing electrostatic two-stream instability (TSI) is negligible. We also show that the longitudinal energy spread, typical of plasma-based accelerated electron–positron fireball beams, plays a minor role in the growth of CFI in these scenarios.

It is known that the filtered, doubly Fourier transformed, Eulerian Vlasov–Poisson simulation method does not exhibit any effects of recurrence. It has not been understood why this is so. It is shown that there is only one boundary in this simulation method where user input is allowed and that the appearance or not of recurrence depends on the character of the boundary data that are imposed. Examples are given that show that the choice of boundary data that leads to recurrence-free simulations does not degrade the accuracy of Landau damping simulations.

Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we generalize the Fokker–Planck operator in a fully ionized plasma by accounting for the effect of screening. We detail the derivation of this generalized operator, and calculate the effective ion length scales, needed in the components of the collision operator, for a number of ion species commonly appearing in fusion experiments. We show that for high electric fields, the secondary runaway growth rate can be substantially larger than in a fully ionized plasma with the same effective charge, although the growth rate is significantly reduced at near-critical electric fields. Furthermore, by comparison with the Boltzmann collision operator, we show that the Fokker–Planck formalism is accurate even for large impurity content.

The structure of a whistler precursor in a quasi-perpendicular shock is studied within two-fluid approach in one-dimensional case. The complete set of equations is reduced to the KdV equation, if no dissipation is included. With a phenomenological resistive dissipation the structure is described with the KdV–Burgers equation. The shock profile is intrinsically time dependent. For sufficiently strong dissipation, temporal evolution of a steepening profile results in generation of a stationary decaying whistler ahead of the shock front. With the decrease of the dissipation parameter, whistler wave trains begin to detach and propagate toward the upstream and the ramp is weakly time dependent. In the weakly dissipative regime the shock front experiences reformation.

A new configuration for an antenna, which increases the absorption power and plasma density, is proposed for helicon plasma sources. The influence of the electromagnetic wave pattern symmetry on the plasma density and absorption power in a helicon plasma source with a common antenna (Nagoya) is analysed by using the standard COMSOL Multiphysics 5.3 software. In contrast to the theoretical model prediction, the electromagnetic wave does not represent a symmetric pattern for the common Nagoya antenna. In this work, a new configuration for an antenna is proposed which refines the asymmetries of the wave pattern in helicon plasma sources. The plasma parameters such as plasma density and absorption rate for a common Nagoya antenna and our proposed antenna under the same conditions are studied using simulations. In addition, the plasma density of seven operational helicon plasma source devices, having a common Nagoya antenna, is compared with the simulation results of our proposed antenna and the common Nagoya antenna. The simulation results show that the density of the plasma, which is produced by using our proposed antenna, is approximately twice in comparison to the plasma density produced by using the common Nagoya antenna. In fact, the simulation results indicate that the electric and magnetic fields symmetry of the helicon wave plays a vital role in increasing wave–particle coupling. As a result, wave–particle energy exchange and the plasma density of helicon plasma sources will be increased.

An electron acoustic super solitary wave has been derived using the Sagdeev pseudopotential technique for a four component magnetized plasma consisting of the beam and bulk fluid electrons and two ions with Maxwell Boltzmann distributions. This is the first theoretical report of a super solitary wave in a magnetized plasma which has no direct association with the singularity of the pseudopotential. It shows a narrow and spiky subwell near the low potential which causes the lateral inversion of the wiggle for the bipolar electric field vis-á-vis the unmagnetized plasma. An analytical formalism was developed to identify these novel kinds of super solitary waves and their transition processes have been characterized. It was observed that the super solitary wave is directly influenced by the singularity of the pseudopotential lying in the vicinity of the solution. The first ever prediction of the electron acoustic super solitary wave raises the possibility of its application to the interpretation of the satellite observations of the electrostatic field data.

Gyrokinetic stability of plasmas in different magnetic geometries is studied numerically using the GENE code. We examine the stability of plasmas, varying the mass ratio between the positive and negative charge carriers, from conventional hydrogen plasmas through to electron–positron plasmas. Stability is studied for prescribed temperature and density gradients in different magnetic geometries: (i) An axisymmetric, circular flux surface, low $\unicode[STIX]{x1D6FD}$ (tokamak) configuration. (ii) A non-axisymmetric quasi-isodynamic (optimised stellarator) configuration using the geometry of the stellarator Wendelstein 7-X. We also present the analytic theory of trapped particle modes in electron–positron plasmas. We found similar behaviour of the growth rate and real frequency compared to previous studies on the tokamak case. We are able to identify two distinct regimes dominated by modes propagating in the electron diamagnetic direction and modes propagating in the ion/positron diamagnetic direction, depending on the mass ratio. In both the tokamak and the stellarator case we observe that the real frequency tends to zero as the mass ratio approaches unity and are able to explain this using gyrokinetic theory.

We have explored the thermodynamics of compressed magnetized plasmas in laboratory experiments and we call these studies ‘magnetothermodynamics’. The experiments are carried out in the Swarthmore Spheromak eXperiment device. In this device, a magnetized plasma source is located at one end and at the other end, a closed conducting can is installed. We generate parcels of magnetized plasma and observe their compression against the end wall of the conducting cylinder. The plasma parameters such as plasma density, temperature and magnetic field are measured during compression using HeNe laser interferometry, ion Doppler spectroscopy and a linear ${\dot{B}}$ probe array, respectively. To identify the instances of ion heating during compression, a PV diagram is constructed using measured density, temperature and a proxy for the volume of the magnetized plasma. Different equations of state are analysed to evaluate the adiabatic nature of the compressed plasma. A three-dimensional resistive magnetohydrodynamic code (NIMROD) is employed to simulate the twisted Taylor states and shows stagnation against the end wall of the closed conducting can. The simulation results are consistent to what we observe in our experiments.

In this paper, we determine the growth rate of the exponential radiation amplification in the ion-channel laser, where a relativistic electron beam wiggles in a focusing ion channel that can be created in a wakefield accelerator. For the first time the radiation diffraction, which can limit the amplification, is taken into account. The electron beam quality requirements to obtain this amplification are also presented. It is shown that both the beam energy and wiggler parameter spreads should be limited. Two-dimensional and three-dimensional particle-in-cell simulations of the self-consistent ion-channel laser confirm our theoretical predictions.