These lecture notes were presented by Allan N. Kaufman in his graduate plasma theory course and a follow-on special topics course (Physics 242A, B, C and Physics 250 at the University of California Berkeley). The notes follow the order of the lectures. The equations and derivations are as Kaufman presented, but the text is a reconstruction of Kaufman’s discussion and commentary. The notes were transcribed by Bruce I. Cohen in 1971 and 1972, and word processed, edited and illustrations added by Cohen in 2017 and 2018. The series of lectures is divided into four major parts: (i) collisionless Vlasov plasmas (linear theory of waves and instabilities with and without an applied magnetic field, Vlasov–Poisson and Vlasov–Maxwell systems, Wentzel–Kramers–Brillouin–Jeffreys (WKBJ) eikonal theory of wave propagation); (ii) nonlinear Vlasov plasmas and miscellaneous topics (the plasma dispersion function, singular solutions of the Vlasov–Poisson system, pulse-response solutions for initial-value problems, Gardner’s stability theorem, gyroresonant effects, nonlinear waves, particle trapping in waves, quasilinear theory, nonlinear three-wave interactions); (iii) plasma collisional and discreteness phenomena (test-particle theory of dynamic friction and wave emission, classical resistivity, extension of test-particle theory to many-particle phenomena and the derivation of the Boltzmann and Lenard–Balescu equations, the Fokker–Planck collision operator, a general scattering theory, nonlinear Landau damping, radiation transport and Dupree’s theory of clumps); (iv) non-uniform plasmas (adiabatic invariance, guiding-centre drifts, hydromagnetic theory, introduction to drift-wave stability theory).

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modelling the turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non-Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the fluid model of Chew–Goldberger–Low (CGL) in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the fourth-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the ‘normal’ closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the ‘normal’ closure, and beyond the fourth-order moment, Landau fluid closures are required.

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(\unicode[STIX]{x1D701})$, and the associated plasma response function $R(\unicode[STIX]{x1D701})=1+\unicode[STIX]{x1D701}Z(\unicode[STIX]{x1D701})=-Z^{\prime }(\unicode[STIX]{x1D701})/2$. We then consider a one-dimensional (1-D) (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the fourth-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1-D closures at higher-order moments than the fourth order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3-D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the fourth-order moment level. In appendix A, we provide comprehensive tables with Padé approximants of $R(\unicode[STIX]{x1D701})$ up to the eighth-pole order, with many given in an analytic form.

Integrated modelling of electron runaway requires computationally expensive kinetic models that are self-consistently coupled to the evolution of the background plasma parameters. The computational expense can be reduced by using parameterized runaway generation rates rather than solving the full kinetic problem. However, currently available generation rates neglect several important effects; in particular, they are not valid in the presence of partially ionized impurities. In this work, we construct a multilayer neural network for the Dreicer runaway generation rate which is trained on data obtained from kinetic simulations performed for a wide range of plasma parameters and impurities. The neural network accurately reproduces the Dreicer runaway generation rate obtained by the kinetic solver. By implementing it in a fluid runaway-electron modelling tool, we show that the improved generation rates lead to significant differences in the self-consistent runaway dynamics as compared to the results using the previously available formulas for the runaway generation rate.

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle $\unicode[STIX]{x1D6FC}$ between the wall and the magnetic field $\boldsymbol{B}$ is oblique. Here, we consider the fusion-relevant case of a shallow-angle, $\unicode[STIX]{x1D6FC}\ll 1$, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for $\unicode[STIX]{x1D6FC}/\sqrt{\unicode[STIX]{x1D70F}+1}\gg \sqrt{m_{\text{e}}/m_{\text{i}}}$, where $m_{\text{e}}$ is the electron mass, $m_{\text{i}}$ is the ion mass, $\unicode[STIX]{x1D70F}=T_{\text{i}}/ZT_{\text{e}}$, $T_{\text{e}}$ is the electron temperature, $T_{\text{i}}$ is the ion temperature and $Z$ is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii $\unicode[STIX]{x1D70C}_{\text{s}}=\sqrt{m_{\text{i}}(ZT_{\text{e}}+T_{\text{i}})}/ZeB$, where e is the proton charge and $B=|\boldsymbol{B}|$ is the magnitude of the magnetic field. We study the dependence on $\unicode[STIX]{x1D70F}$ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by $\unicode[STIX]{x1D70F}$. The kinetic model is shown to be asymptotically equivalent to Chodura’s fluid model at small ion temperature, $\unicode[STIX]{x1D70F}\ll 1$, for $|\text{ln}\,\unicode[STIX]{x1D6FC}|>3|\text{ln}\,\unicode[STIX]{x1D70F}|\gg 1$. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, $\unicode[STIX]{x1D70F}\gg 1$, relevant because $T_{\text{i}}$ is measured to be a few times larger than $T_{\text{e}}$ near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrt{\unicode[STIX]{x1D6FC}}$ or $1/\sqrt{\unicode[STIX]{x1D70F}}$, depending on which is largest.

Boundary layers in space and astrophysical plasmas are the location of complex dynamics where different mechanisms coexist and compete, eventually leading to plasma mixing. In this work, we present fully kinetic particle-in-cell simulations of different boundary layers characterized by the following main ingredients: a velocity shear, a density gradient and a magnetic gradient localized at the same position. In particular, the presence of a density gradient drives the development of the lower-hybrid drift instability (LHDI), which competes with the Kelvin–Helmholtz instability (KHI) in the development of the boundary layer. Depending on the density gradient, the LHDI can even dominate the dynamics of the layer. Because these two instabilities grow on different spatial and temporal scales, when the LHDI develops faster than the KHI an inverse cascade is generated, at least in two dimensions. This inverse cascade, starting at the LHDI kinetic scales, generates structures at scale lengths at which the KHI would typically develop. When that is the case, those structures can suppress the KHI itself because they significantly affect the underlying velocity shear gradient. We conclude that, depending on the density gradient, the velocity jump and the width of the boundary layer, the LHDI in its nonlinear phase can become the primary instability for plasma mixing. These numerical simulations show that the LHDI is likely to be a dominant process at the magnetopause of Mercury. These results are expected to be of direct impact to the interpretation of the forthcoming BepiColombo observations.

A method is given to rapidly compute quasisymmetric stellarator magnetic fields for plasma confinement, without the need to call a three-dimensional magnetohydrodynamic equilibrium code inside an optimization iteration. The method is based on direct solution of the equations of magnetohydrodynamic equilibrium and quasisymmetry using Garren & Boozer’s expansion about the magnetic axis (Phys Fluids B, vol. 3, 1991, p. 2805), and it is several orders of magnitude faster than the conventional optimization approach. The work here extends the method of Landreman et al. (J. Plasma Phys., vol. 85, 2019, 905850103), which was limited to flux surfaces with elliptical cross-section, to higher order in the aspect-ratio expansion. As a result, configurations can be generated with strong shaping that achieve quasisymmetry to high accuracy. Using this construction, we give the first numerical demonstrations of Garren and Boozer’s ideal scaling of quasisymmetry breaking with the cube of the inverse aspect ratio. We also demonstrate a strongly non-axisymmetric configuration (vacuum rotational transform $\unicode[STIX]{x1D704}>0.4$) in which symmetry-breaking mode amplitudes throughout a finite volume are ${<}2\times 10^{-7}$, the smallest ever reported. To generate boundary shapes of finite-minor-radius configurations, a careful analysis is given of the effect of substituting a finite minor radius into the near-axis expansion. The approach here can provide analytic insight into the space of possible quasisymmetric stellarator configurations, and it can be used to generate good initial conditions for conventional stellarator optimization.

The stability of conducting Taylor–Couette flows under the presence of toroidal magnetic background fields is considered. For strong enough magnetic amplitudes such magnetohydrodynamic flows are unstable against non-axisymmetric perturbations which may also transport angular momentum. In accordance with the often used diffusion approximation, one expects the angular momentum transport to be vanishing for rigid rotation. In the sense of a non-diffusive $\unicode[STIX]{x1D6EC}$ effect, however, even for rigidly rotating $z$-pinches, an axisymmetric angular momentum flux appears which is directed outward (inward) for large (small) magnetic Mach numbers. The internal rotation in a magnetized rotating tank can thus never be uniform. Those particular rotation laws are used to estimate the value of the instability-induced eddy viscosity for which the non-diffusive $\unicode[STIX]{x1D6EC}$ effect and the diffusive shear-induced transport compensate each other. The results provide the Shakura & Sunyaev viscosity ansatz leading to numerical values linearly growing with the applied magnetic field.

The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren & Boozer (Phys. Fluids B, vol. 3 (10), 1991a, pp. 2805–2821). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman & Sengupta (J. Plasma Phys., vol. 84 (6), 2018, 905840616); Landreman et al. (J. Plasma Phys., vol. 85 (1), 2019, 905850103), enables their efficient numerical construction.

Beam-driven instabilities are considered in a pulsar plasma assuming that both the background plasma and the beam are relativistic Jüttner distributions. In the rest frame of the background, the only waves that can satisfy the resonance condition are in a tiny range of slightly subluminal phase speeds. The growth rate for the kinetic (or maser) version of the weak-beam instability is much smaller than has been estimated for a relativistically streaming Gaussian distribution, and the reasons for this are discussed. The growth rate for the reactive version of the weak-beam instability is treated in a conventional way. We compare the results with exact calculations, and find that the approximate solutions are not consistent with the exact results. We conclude that, for plausible parameters, there is no reactive version of the instability. The growth rate in the pulsar frame is smaller than that in the rest frame of the background plasma by a factor $2\unicode[STIX]{x1D6FE}_{\text{s}}$, where $\unicode[STIX]{x1D6FE}_{\text{s}}=10^{2}{-}10^{3}$ is the Lorentz factor of the bulk motion of the background plasma, placing a further constraint on effective wave growth. Based on these results, we argue that beam-driven wave growth probably plays no role in pulsar radio emission.

A gyrokinetic Coulomb collision operator is derived, which is particularly useful to describe the plasma dynamics at the periphery region of magnetic confinement fusion devices. The derived operator is able to describe collisions occurring in distribution functions arbitrarily far from equilibrium with variations on spatial scales at and below the particle Larmor radius. A multipole expansion of the Rosenbluth potentials is used in order to derive the dependence of the full Coulomb collision operator on the particle gyroangle. The full Coulomb collision operator is then expressed in gyrocentre phase-space coordinates, and a closed formula for its gyroaverage in terms of the moments of the gyrocentre distribution function in a form ready to be numerically implemented is provided. Furthermore, the collision operator is projected onto a Hermite–Laguerre velocity space polynomial basis and expansions in the small electron-to-ion mass ratio are provided.

This paper discusses a detailed computational analysis that illustrated the influences of the magnetic field and external potential on the performance of a high-temperature Hall-effect electric thruster. Uniform and non-uniform magnetic field configurations were examined. The Lorentz force in the $x$ direction, acting on the plasma, was shown to substantially enhance the flow velocity in the non-uniform magnetic field, which indicated that the non-uniform magnetic field was more suitable for Hall-effect electromagnetic acceleration. The static temperature increased with the external potential, especially near the region of cathode. This increment in gas temperature, together with the effect of the Lorentz force, results in the enhancement of the velocity at the front and back of the cathode. However, the Mach number and gas density decreased due to static temperature increases caused by the conversion of more electric power into internal energy. The thrust increased eventually with the increase of the average exit velocity.

We derive a dispersion relation for the damping of acoustic waves in equi-molar deuterium–tritium (DT) gas due to radiation coupling and electron thermal conduction and discuss its significance for inertial confinement fusion (ICF) targets with high-Z shells surrounding a central DT fuel region. As the shell implodes around DT fuel in such a target, shocks and waves are transmitted through the DT gas. If the shell is perturbed due to drive non-uniformity or manufacturing imperfection, these shocks and waves may be perturbed as well, and can potentially re-perturb the shell. This can complicate calculation of shell stability and implosion asymmetry and in general make the target less robust against implosion non-uniformity. Damping of perturbations in DT gas can alleviate these complications. Also, damping of low-order modes, which is primarily due to radiation coupling, can drive the DT gas to an isobaric and isothermal ‘equilibrium’ configuration during ignition. We find that for the range of common ignition temperatures in targets with high-Z shells, $2.5\lesssim T_{ig}\lesssim 3.5$  keV, damping of low-order modes is significant for areal densities ($\unicode[STIX]{x1D70C}r$) in the broad range of $0.6\lesssim \unicode[STIX]{x1D70C}r\lesssim 1.8~\text{g}~\text{cm}^{-2}$. This suggests it is advantageous to design these targets to achieve areal densities at ignition within this range. Furthermore, we derive a simple constraint between areal density and temperature, $\unicode[STIX]{x1D70C}r=0.34T_{o}$ where $T_{o}$ is in keV, such that DT gas undergoing equilibrium ignition is optimally robust against non-uniformity.

We show analytically that for $\unicode[STIX]{x1D704}$-profiles similar to the one of the Wendelstein 7-X stellarator, where $\unicode[STIX]{x1D704}$ is the rotational transform of the equilibrium magnetic field, a highly conducting toroidal plasma is unstable to kinetically mediated pressure-driven long-wavelength reconnecting modes, of the infernal type. The modes are destabilized either by the electron temperature gradient or by a small amount of current, depending on how far from unity the average value of $\unicode[STIX]{x1D704}$ is, which is assumed to be slowly varying. We argue that, for W7-X, a broad mode with toroidal and poloidal mode numbers $(n,m)=(1,1)$ can be destabilized due to the strong geometric side-band coupling of the resonant kinetic electron response at locations where $\unicode[STIX]{x1D704}$ is rational for harmonics that belong to the mode family of the $(n,m)=(1,1)$ mode itself. In many regimes, the growth rate is insensitive to the plasma density, thus it is likely to persist in high performance W7-X discharges. For a peaked electron temperature, with a maximum of $T_{e}=5~\text{keV}$, larger than the ion temperature, $T_{i}=2.5~\text{keV}$, and a density $n_{0}=10^{19}~\text{m}^{-3}$, instability is found in regimes which show plasma sawtooth activity, with growth rates of the order of tens of kiloHertz. Frequencies are either electron diamagnetic or of the ideal magnetohydrodynamic type, but sub-Alfvénic. The kinetic infernal mode is thus a good candidate for the explanation of sawtooth oscillations in present-day stellarators and poses a new challenge to the problem of stellarator reactor optimization.

The multipolar representation of the magnetic field has, for the lowest-order term, a magnetic dipole that dominates the far field. Thus the far-field representation of the magnetic field of the Earth, Sun and other celestial bodies is a dipole. Since these bodies consist of or are surrounded by plasma, which can support Alfvén waves, their propagation along dipole magnetic field lines is considered using a new coordinate system: dipolar coordinates. The present paper introduces multipolar coordinates, which are an example of conformal coordinates; conformal coordinates are orthogonal with equal scale factors, and can be extended from the plane to space, for instance as cylindrical or spherical dipolar coordinates. The application considered is to Alfvén waves propagating along a circle, that is a magnetic field line of a dipole, with transverse velocity and magnetic field perturbations; the various forms of the wave equation are linear second-order differential equations, with variable coefficients, specified by a background magnetic field, which is force free. The absence of a background magnetic force leads to a mean state of hydrostatic equilibrium, specified by the balance of gravity against the pressure gradient, for a perfect gas or incompressible liquid. The wave equation is simplified to a Gaussian hypergeometric type in the case of zero frequency, otherwise, for non-zero frequency, an extended Gaussian hypergeometric equation is obtained. The solution of the latter specifies the magnetic field perturbation spectrum, and also, via a polarisation relation, the velocity perturbation spectrum; both are plotted, over half a circle, for three values of the dimensionless frequency.

The linear theory stability of different collisionless plasma sheath structures, including the classic sheath, inverse sheath and space-charge limited (SCL) sheath, is investigated as a typical eigenvalue problem. The three background plasma sheaths formed between a Maxwellian plasma source and a dielectric wall with a fully self-consistent secondary electron emission condition are solved by recent developed 1D3V (one-dimensional space and three-dimensional velocities), steady-state, collisionless kinetic sheath model, within a wide range of Maxwellian plasma electron temperature $T_{e}$. Then, the eigenvalue equations of sheath plasma fluctuations through the three sheaths are numerically solved, and the corresponding damping and growth rates $\unicode[STIX]{x1D6FE}$ are found: (i) under the classic sheath structure (i.e. $T_{e}<T_{ec}$ (the first threshold)), there are three damping solutions (i.e. $\unicode[STIX]{x1D6FE}_{1}$, $\unicode[STIX]{x1D6FE}_{2}$ and $\unicode[STIX]{x1D6FE}_{3}$, $0>\unicode[STIX]{x1D6FE}_{1}>\unicode[STIX]{x1D6FE}_{2}>\unicode[STIX]{x1D6FE}_{3}$) for most cases, but there is only one growth-rate solution $\unicode[STIX]{x1D6FE}$ when $T_{e}\rightarrow T_{ec}$ due to the inhomogeneity of sheath being very weak; (ii) under the inverse sheath structure, which arises when $T_{e}>T_{ec}$, there are no background ions in the sheath so that the fluctuations are stable; (iii) under the SCL sheath conditions (i.e. $T_{e}\geqslant T_{e\text{SCL}}$, the second threshold), the obvious ion streaming through the sheath region again emerges and the three damping solutions are again found.

Simulations of one-dimensional Vlasov–Maxwell solutions with non-periodic boundary conditions are discussed. Results obtained using a recently developed filtered flux-balance simulation system are compared to those obtained using a filtered, Fourier–Fourier transformed system. Excellent agreement is confirmed except for the appearance of the Gibbs phenomenon on the discontinuous simulated solutions of the transformed system. Recovery of the flux-balance results from the Fourier transformed results using the inverse polynomial reconstruction method is demonstrated.

The propagation of dust-acoustic (DA) shock waves (SWs) is studied in a four-component dusty plasma system containing viscous dust fluids of opposite polarity, Schamel distributed ions and Boltzmann distributed electrons. The reductive perturbation method is employed to derive a modified Burgers equation which gives rise to the DA shock waves with stronger nonlinearity. The viscous force acting in the dust fluids is identified as a source of dissipation, and is responsible for the formation of the DA shock waves. The basic characteristics (viz., speed, amplitude, width) of the DA shock waves are found to be significantly modified by the combined effects of opposite polarity dust fluids and trapped ions. The applications of this investigation in different space plasma environments and laboratory devices are pinpointed.

The results of a basic electron heat transport experiment using multiple localized heat sources in close proximity and embedded in a large magnetized plasma are presented. The set-up consists of three biased probe-mounted crystal cathodes, arranged in a triangular spatial pattern, that inject low energy electrons along a strong magnetic field into a pre-existing, cold afterglow plasma, forming electron temperature filaments. When the three sources are activated and placed within a few collisionless electron skin depths of each other, a non-azimuthally symmetric wave pattern emerges due to interference of the drift-Alfvén modes that form on each filament’s temperature gradient. Enhanced cross-field transport from chaotic ($\boldsymbol{E}\times \boldsymbol{B}$, where $\boldsymbol{E}$ is the electric field and $\boldsymbol{B}$ the magnetic field) mixing rapidly relaxes the gradients in the inner triangular region of the filaments and leads to growth of a global nonlinear drift-Alfvén mode that is driven by the thermal gradient in the outer region of the triangle. Azimuthal flow shear arising from the emissive cathode sources modifies the linear eigenmode stability and convective pattern. A steady-current model with emissive sheath boundary predicts the plasma potential and shear flow contribution from the sources.