The strong, sharp flow structures that are seen frequently in tokamak cores, and large amplitude spontaneous global toroidal rotation are both surprising in light of current theories where toroidal flow evolution is dominantly diffusive. Mechanisms for spontaneously generating strong poloidal shear flows have been extensively investigated, but these processes were thought not to apply to toroidal flows. We confirm, however, that there is a regime with near-zero toroidal momentum diffusivity, where toroidal flow structures are spontaneously generated, as shown in earlier global gyrokinetic simulations. This allows strong rotation with negligible applied external torque, and the regime where this occurs is also favourable for strong turbulence stabilisation by rotation. The transition to low momentum diffusivity occurs for tight aspect ratio, low safety factor and at low levels of heat flux, in agreement with a simple zero-dimensional model for the flow–turbulence interaction.

We study the damping of collisionless Alfvénic turbulence in a strongly magnetised plasma by two mechanisms: stochastic heating (whose efficiency depends on the local turbulence amplitude $\unicode[STIX]{x1D6FF}z_{\unicode[STIX]{x1D706}}$) and linear Landau damping (whose efficiency is independent of $\unicode[STIX]{x1D6FF}z_{\unicode[STIX]{x1D706}}$), describing in detail how they affect and are affected by intermittency. The overall efficiency of linear Landau damping is not affected by intermittency in critically balanced turbulence, while stochastic heating is much more efficient in the presence of intermittent turbulence. Moreover, stochastic heating leads to a drop in the scale-dependent kurtosis over a narrow range of scales around the ion gyroscale.

A two-field Hamiltonian gyrofluid model for kinetic Alfvén waves retaining ion finite Larmor radius corrections, parallel magnetic field fluctuations and electron inertia, is used to study turbulent cascades from the magnetohydrodynamic (MHD) to the sub-ion scales. Special attention is paid to the case of imbalance between waves propagating along or opposite to the ambient magnetic field. For weak turbulence in the absence of electron inertia, kinetic equations for the spectral density of the conserved quantities (total energy and generalized cross-helicity) are obtained. They provide a global description, matching between the regimes of reduced MHD at large scales and electron reduced MHD at small scales, previously considered in the literature. In the limit of ultra-local interactions, Leith-type nonlinear diffusion equations in the Fourier space are derived and heuristically extended to the strong turbulence regime by modifying the transfer time appropriately. Relations with existing phenomenological models for imbalanced MHD and balanced sub-ion turbulence are discussed. It turns out that in the presence of dispersive effects, the dynamics is sensitive on the way turbulence is maintained in a steady state. Furthermore, the total energy spectrum at sub-ion scales becomes steeper as the generalized cross-helicity flux is increased.

The model of the global gyrokinetic particle-in-cell code ORB5 has been extended for the study of pair plasmas. This has been done by including the physics of the Debye shielding, by including the electron polarization density and by retaining the effects of the electron finite Larmor radius. This model is verified against previous numerical results for the cyclone base case tokamak scenario in deuterium plasmas, and for local pair plasma simulations. The linear dynamics of temperature-gradient driven instabilities and geodesic acoustic modes is investigated. Mass dependencies for different Debye lengths are studied.

It is shown that in low-beta, weakly collisional plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae and some laboratory plasmas, Alfvénic fluctuations produce no ion heating within the gyrokinetic approximation, i.e. as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion-Larmor frequency (even though their spatial scales can be above or below the ion Larmor scale). Thus, all low-frequency ion heating in such plasmas is due to compressive fluctuations (‘slow modes’): density perturbations and non-Maxwellian perturbations of the ion distribution function. Because these fluctuations energetically decouple from the Alfvénic ones already in the inertial range, the above conclusion means that the energy partition between ions and electrons in low-beta plasmas is decided at the outer scale, where turbulence is launched, and can be determined from magnetohydrodynamic (MHD) models of the relevant astrophysical systems. Any additional ion heating must come from non-gyrokinetic mechanisms such as cyclotron heating or the stochastic heating owing to distortions of ions’ Larmor orbits. An exception to these conclusions occurs in the Hall limit, i.e. when the ratio of the ion to electron temperatures is as low as the ion beta (equivalently, the electron beta is order unity). In this regime, slow modes couple to Alfvénic ones well above the Larmor scale (viz., at the ion inertial or ion sound scale), so the Alfvénic and compressive cascades join and then separate again into two cascades of fluctuations that linearly resemble kinetic Alfvén and ion-cyclotron waves, with the former heating electrons and the latter ions. The two cascades are shown to decouple, scalings for them are derived and it is argued physically that the two species will be heated by them at approximately equal rates.

I analyse and numerically evaluate the radiation field generated by an experimentally realized embodiment of an electric polarization current whose rotating distribution pattern moves with linear speeds exceeding the speed of light in vacuum. I find that the flux density of the resulting emission (i) has a dominant value and is linearly polarized within a sharply delineated radiation beam whose orientation and polar width are determined by the range of values of the linear speeds of the rotating source distribution, and (ii) decays with the distance $d$ from the source as $d^{-\unicode[STIX]{x1D6FC}}$ in which the value of $\unicode[STIX]{x1D6FC}$ lies between $1$ and $2$ (instead of being equal to $2$ as in a conventional radiation) across the beam. In that the rate at which boundaries of the retarded distribution of such a source change with time depends on its duration monotonically, this is an intrinsically transient emission process: temporal rate of change of the energy density of the radiation generated by it has a time-averaged value that is negative (instead of being zero as in a conventional radiation) at points where the envelopes of the wave fronts emanating from the constituent volume elements of the source distribution are cusped. The difference in the fluxes of power across any two spheres centred on the source is in this case balanced by the change with time of the energy contained inside the shell bounded by those spheres. These results are relevant not only to long-range transmitters in communications technology but also to astrophysical objects containing rapidly rotating neutron stars (such as pulsars) and to the interpretation of the energetics of the multi-wavelength emissions from sources that lie at cosmological distances (such as radio and gamma-ray bursts). The analysis presented in this paper is self-contained and supersedes my earlier works on this problem.

Wave dispersion in a pulsar plasma (a one-dimensional, strongly magnetized, pair plasma streaming highly relativistically with a large spread in Lorentz factors in its rest frame) is discussed, motivated by interest in beam-driven wave turbulence and the pulsar radio emission mechanism. In the rest frame of the pulsar plasma there are three wave modes in the low-frequency, non-gyrotropic approximation. For parallel propagation (wave angle $\unicode[STIX]{x1D703}=0$) these are referred to as the X, A and L modes, with the X and A modes having dispersion relation $|z|=z_{\text{A}}\approx 1-1/2\unicode[STIX]{x1D6FD}_{\text{A}}^{2}$, where $z=\unicode[STIX]{x1D714}/k_{\Vert }c$ is the phase speed and $\unicode[STIX]{x1D6FD}_{\text{A}}c$ is the Alfvén speed. The L mode dispersion relation is determined by a relativistic plasma dispersion function, $z^{2}W(z)$, which is negative for $|z|<z_{0}$ and has a sharp maximum at $|z|=z_{\text{m}}$, with $1-z_{\text{m}}<1-z_{0}\ll 1$. We give numerical estimates for the maximum of $z^{2}W(z)$ and for $z_{\text{m}}$ and $z_{0}$ for a one-dimensional Jüttner distribution. The L and A modes reconnect, for $z_{\text{A}}>z_{0}$, to form the O and Alfvén modes for oblique propagation ($\unicode[STIX]{x1D703}\neq 0$). For $z_{\text{A}}<z_{0}$ the Alfvén and O mode curves reconnect forming a new mode that exists only for $\tan ^{2}\unicode[STIX]{x1D703}\gtrsim z_{0}^{2}-z_{\text{A}}^{2}$. The L mode is the nearest counterpart to Langmuir waves in a non-relativistic plasma, but we argue that there are no ‘Langmuir-like’ waves in a pulsar plasma, identifying three features of the L mode (dispersion relation, ratio of electric to total energy and group speed) that are not Langmuir like. A beam-driven instability requires a beam speed equal to the phase speed of the wave. This resonance condition can be satisfied for the O mode, but only for an implausibly energetic beam and only for a tiny range of angles for the O mode around $\unicode[STIX]{x1D703}\approx 0$. The resonance is also possible for the Alfvén mode but only near a turnover frequency that has no counterpart for Alfvén waves in a non-relativistic plasma.

It is shown that an electron, positron and ion plasma can be self-organized to a double Beltrami state – the superposition of two force-free states. The scale parameters which determine the nature of the self-organized structures are found to depend on the number densities of the plasma species. The loss of equilibrium in slowly evolving double Beltrmi states is investigated. The effects of density ratios, helicities, positron flows and energy on equilibrium are investigated. It is found that the double Beltrami state transforms to a single Beltrami state at the termination of equilibrium. It is also shown that much of the magnetic energy converts to the flow kinetic energy through catastrophic transformation.

A modified Vlasov equation is obtained by developing a covariant statistical mechanics for a system of electrons without considering the effects of the ions and including the Landau–Lifshitz equation of motion. General dispersion relations for the transverse and longitudinal modes for any temperature are expressed. The results are similar to those found by Hakim & Mangeney (Phys. Fluids, vol. 14, 1971, pp. 2751–2781) for both the modified Vlasov equation and the dispersion relations. However, for the longitudinal mode, unlike the development of Hakim and Mangeney, correct expansions are done in order to give a numerical approach to obtain the longitudinal relativistic dispersion relations for any value of the wavenumber. Accordingly, new loop solutions, with turning points, crossing the super-luminous region and the super-thermal region are found. Although the expressions for the Landau damping and the damping due to the radiation reaction force coincide with the Hakim and Mangeney results for some particular cases, in general they are different. A Landau anti-damping appears in the second branch of the loop in a small region between the cutoff point and the intersection with the super-thermal line. The analysis of this effect leads us to a kind of wave pulse. We will call them bipolar waves. The treatment contains the relativistic interactions between all the electrons in the system with retarded effects. This explain the differences with Zhang’s recent work (Phys. Plasmas, vol. 20, 2013, 092112–092132). It is shown that for low densities, the cutoff of the wave is due to the dispersion relations and not due to the radiation reaction force damping. While for both high densities and temperatures, the damping due to the radiation reaction force is important.

An analytical study of the small amplitude electron acoustic double layers in a magnetized plasma consisting of superthermal electrons and ions along with cold fluid electrons is discussed. The dispersion relation allows electron acoustic waves with the frequency within electron and ion gyro-frequency in the modelled plasma. In the process of study of the nonlinear structures, the Sagdeev pseudo-potential method for small amplitude regions is employed. The existence domains for the double layers are investigated in terms of the Mach numbers of the structures and the temperature ratios of the species for different ratios of their concentration. The effects of the compositional parameters on the nature and size of the double layers are also explored and it is observed that the plasma can support both compressive and rarefactive double layers depending on the values of those parameters and the Mach numbers.

The use of radio frequency (RF) waves in fusion plasmas for heating, for non-inductive current generation, for profile control and for diagnostics has been well established. The RF waves, excited by antenna structures placed near the wall of a fusion device, have to propagate through density fluctuations at the plasma edge. These fluctuations can modify the properties of the RF waves that propagate towards the core of the plasma. A full-wave electromagnetic computational code ScaRF based on the finite difference frequency domain (FDFD) method has been developed to study the effect of density turbulence on RF waves. The anisotropic plasma permittivity used in the scattering studies is that for a magnetized, cold plasma. The code is used to study the propagation of an RF plane wave through a modulated, spatially periodic density interface. Such an interface could arise in the edge region due to magnetohydrodynamic instability or drift waves. The frequency of the plane wave is taken to be in the range of the electron cyclotron frequency. The scattering analysis is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D. The effect of different density contrasts across the interface and of different spatial modulations are discussed. While ScaRF is used to study a periodic density fluctuation, the code is general enough to include different varieties of density fluctuations in the edge region – such as blobs and filaments, and spatially random fluctuations.

The work described in this article was carried out to investigate how permanent magnets (PM) affect the plasma confinement and ion beam properties in an inductively coupled plasma which expands from a helicon source. The cylindrical plasma device Njord has a 13 cm long and 20 cm wide stainless steel port connecting the source chamber and the diffusion chamber. The source chamber has an axial magnetic field produced by two coils, with magnetic field lines expanding into the diffusion chamber. Simulations have shown that the field lines leaving the edge of the source hit the port wall, causing a loss of electrons in this section. In the experiments performed in this work, PMs were added around the port walls near the exit of a plasma source and the effect was investigated experimentally by means of a retarding field energy analyser probe. The plasma potential, ion density and ion beam parameters were estimated, and the results with and without the PMs were compared. The results showed that the plasma density in the centre can in some cases be doubled, and the density at the edges of the plasma increased significantly with PMs in place. Although the plasma potential was slightly affected, and the beam velocity dropped by ${\sim}$10 %, the ion beam flux increased by a factor of 1.5.

Wave dispersion in a pulsar plasma is discussed emphasizing the relevance of different inertial frames, notably the plasma rest frame ${\mathcal{K}}$ and the pulsar frame ${\mathcal{K}}^{\prime }$ in which the plasma is streaming with speed $\unicode[STIX]{x1D6FD}_{\text{s}}$. The effect of a Lorentz transformation on both subluminal, $|z|<1$, and superluminal, $|z|>1$, waves is discussed. It is argued that the preferred choice for a relativistically streaming distribution should be a Lorentz-transformed Jüttner distribution; such a distribution is compared with other choices including a relativistically streaming Gaussian distribution. A Lorentz transformation of the dielectric tensor is written down, and used to derive an explicit relation between the relativistic plasma dispersion functions in ${\mathcal{K}}$ and ${\mathcal{K}}^{\prime }$. It is shown that the dispersion equation can be written in an invariant form, implying a one-to-one correspondence between wave modes in any two inertial frames. Although there are only three modes in the plasma rest frame, it is possible for backward-propagating or negative-frequency solutions in ${\mathcal{K}}$ to transform into additional forward-propagating, positive-frequency solutions in ${\mathcal{K}}^{\prime }$ that may be regarded as additional modes.

Thousands of gyrokinetic papers have been published since the introduction of the gyrokinetic change of variables. The intent here is not to review the field, but rather the goal is to present a tutorial providing insight into why gyrokinetics is an appropriate description of turbulent transport. The focus is on turbulent transport in axisymmetric tokamaks, but many of the ideas and techniques are applicable to stellarators and other magnetic fusion devices with nested surfaces of constant pressure. Besides the origins of gyrokinetics and recent insights, gyrokinetic orderings and gyrokinetic variables are summarized. Then a compact, but careful, derivation of the simplest electrostatic gyrokinetic equation for tokamaks is presented, along with a brief mention of gyrokinetics for stellarators. The advantages of assuming scale separation between the finer spatial scales and faster time variation of the turbulence and the global behaviour and slow temporal evolution of the background are stressed. Moreover, the procedure for removing the adiabatic or Maxwell–Boltzmann response is emphasized. Scale separation allows the near Maxwellian behaviour of the background and the rapid variation of the turbulence to be described by separate, local gyrokinetic equations. The turbulent fluctuations are found by solving the nonlinear gyrokinetic equation for the non-adiabatic portion of the fluctuating distribution function. Also, generalizations of the gyrokinetic equation so that it retains electromagnetic and flow modifications are considered in detail along with some symmetry properties. In addition, ambipolarity, particle transport and heat transport are addressed, along with a discussion of the complications associated with describing momentum transport and profile evolution.