Adam Stahl

Synchrotron radiation from a runaway electron distribution in tokamaks

Analysis of the synchrotron radiation emitted by runaway electrons in a fusion plasma presents a direct route to knowledge of the particle momenta and pitch-angles of the runaway electron population, through the strong dependence of the synchrotron spectrum on these parameters. Information about the runaway density and its spatial distribution, as well as the time evolution of the above quantities can also be deduced. Previously, synchrotron spectra have been interpreted under the assumption that all runaways have the same energy and pitch-angle [1]. In practice, however, runaway electrons have a wide range of energies and pitch-angles, and that influences the synchrotron radiation spectrum. In this paper we present the synchrotron radiation spectra for typical avalanching runaway electron distributions. We compare the characteristics of the spectra obtained for a distribution of electrons to the one based on the emission of electrons having the same energy and pitch-angle. We examine the eects on the spectrum of including or neglecting magnetic field curvature in the calculation and analyse the sensitivity of the resulting spectrum to perturbations to the runaway distribution. We also discuss the implications for the deduced runaway parameters.

Numerical calculation of the runaway electron distribution function and associated synchrotron emission

Synchrotron emission from runaway electrons may be used to diagnose plasma conditions during a tokamak disruption, but solving this inverse problem requires rapid simulation of the electron distribution function and associated synchrotron emission as a function of plasma parameters. Here we detail a framework for this forward calculation, beginning with an efficient numerical method for solving the Fokker-Planck equation in the presence of an electric field of arbitrary strength. The approach is continuum (Eulerian), and we employ a relativistic collision operator, valid for arbitrary energies. Both primary and secondary runaway electron generation are included. For cases in which primary generation dominates, a time-independent formulation of the problem is described, requiring only the solution of a single sparse linear system. In the limit of dominant secondary generation, we present the first numerical verification of an analytic model for the distribution function. The numerical electron distribution function in the presence of both primary and secondary generation is then used for calculating the synchrotron emission spectrum of the runaways. It is found that the average synchrotron spectra emitted from realistic distribution functions are not well approximated by the emission of a single electron at the maximum energy.

Eddington capture sphere around luminous stars

Test particles infalling from infinity onto a compact spherical star with a mildly super-Eddington luminosity at its surface are typically trapped on the “Eddington capture sphere” and do not reach the surface of the star. The presence of a sphere on which radiation pressure balances gravity for static particles was first discovered some twenty five years ago. Subsequently, it was shown to be a capture sphere for particles in radial motion, and more recently also for particles in non-radial motion, in which the Poynting-Robertson radiation drag efficiently removes the orbital angular momentum of the particles, reducing it to zero. Here we develop this idea further, showing that “levitation” on the Eddington sphere (above the stellar surface) is a state of stable equilibrium, and discuss its implications for Hoyle-Lyttleton accretion onto a luminous star. When the Eddington sphere is present, the cross-section of a compact star for actual accretion is typically less than the geometrical cross-section πR2, direct infall onto the stellar surface only being possible for relativistic particles, with the required minimum particle velocity at infinity typically about half the speed of light. We further show that particles on typical trajectories in the vicinity of the stellar surface will also be trapped on the Eddington capture sphere.

Oscillations of the Eddington capture sphere

We present a toy model of mildly super-Eddington, optically thin accretion onto a compact star in the Schwarzschild metric, which predicts periodic variations of luminosity when matter is supplied to the system at a constant accretion rate. These are related to the periodic appearance and disappearance of the Eddington capture sphere. In the model the frequency is found to vary inversely with the luminosity. If the input accretion rate varies (strictly) periodically, the luminosity variation is quasi-periodic, and the quality factor is inversely proportional to the relative amplitude of mass accretion fluctuations, with its largest value Q ≈ 1/(10 |δṀ/Ṁ|) attained in oscillations at about 1 to 2 kHz frequencies for a 2  M⊙ star.

Effect of Partially Screened Nuclei on Fast-electron Dynamics

We analyze the dynamics of fast electrons in plasmas containing partially ionized impurity atoms, where the screening effect of bound electrons must be included. We derive analytical expressions for the deflection and slowing-down frequencies, and show that they are increased significantly compared to the results obtained with complete screening, already at subrelativistic electron energies. Furthermore, we show that the modifications to the deflection and slowing down frequencies are of equal importance in describing the runaway current evolution. Our results greatly affect fast-electron dynamics and have important implications, e.g., for the efficacy of mitigation strategies for runaway electrons in tokamak devices, and energy loss during relativistic breakdown in atmospheric discharges.

Numerical characterization of bump formation in the runaway electron tail

Runaway electrons are generated in a magnetized plasma when the parallel electric field exceeds a critical value. For such electrons with energies typically reaching tens of MeV, the Abraham–Lorentz–Dirac (ALD) radiation force, in reaction to the synchrotron emission, is significant and can be the dominant process limiting electron acceleration. The effect of the ALD force on runaway electron dynamics in a homogeneous plasma is investigated using the relativistic finite-difference Fokker–Planck codes LUKE (Decker and Peysson 2004 Report EUR-CEA-FC-1736, Euratom-CEA), and CODE (Landreman et al 2014 Comput. Phys. Commun. 185 847). The time evolution of the distribution function is analyzed as a function of the relevant parameters: parallel electric field, background magnetic field, and effective charge. Under the action of the ALD force, we find that runaway electrons are subject to an energy limit, and that the electron distribution evolves towards a steady-state. In addition, a bump is formed in the tail of the electron distribution function if the electric field is sufficiently strong. The mechanisms leading to the bump formation and energy limit involve both the parallel and perpendicular momentum dynamics; they are described in detail. An estimate for the bump location in momentum space is derived. We observe that the energy of runaway electrons in the bump increases with the electric field amplitude, while the population increases with the bulk electron temperature. The presence of the bump divides the electron distribution into a runaway beam and a bulk population. This mechanism may give rise to beam-plasma types of instabilities that could, in turn, pump energy from runaway electrons and alter their confinement.

SOFT: A synthetic synchrotron diagnostic for runaway electrons

Improved understanding of the dynamics of runaway electrons can be obtained by measurement and interpretation of their synchrotron radiation emission. Models for synchrotron radiation emitted by relativistic electrons are well established, but the question of how various geometric effects -- such as magnetic field inhomogeneity and camera placement -- influence the synchrotron measurements and their interpretation remains open. In this paper we address this issue by simulating synchrotron images and spectra using the new synthetic synchrotron diagnostic tool SOFT (Synchrotron-detecting Orbit Following Toolkit). We identify the key parameters influencing the synchrotron radiation spot and present scans in those parameters. Using a runaway electron distribution function obtained by Fokker-Planck simulations for parameters from an Alcator C-Mod discharge, we demonstrate that the corresponding synchrotron image is well-reproduced by SOFT simulations, and we explain how it can be understood in terms of the parameter scans. Geometric effects are shown to significantly influence the synchrotron spectrum, and we show that inherent inconsistencies in a simple emission model (i.e. not modeling detection) can lead to incorrect interpretation of the images.

Kinetic modelling of runaway electrons in dynamic scenarios

Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concern. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.

Quasi-linear analysis of the extraordinary electron wave destabilized by runaway electrons

Runaway electrons with strongly anisotropic distributions present in post-disruption tokamak plasmas can destabilize the extraordinary electron (EXEL) wave. The present work investigates the dynamics of the quasi-linear evolution of the EXEL instability for a range of different plasma parameters using a model runaway distribution function valid for highly relativistic runaway electron beams produced primarily by the avalanche process. Simulations show a rapid pitch-angle scattering of the runaway electrons in the high energy tail on the 100–1000 μs time scale. Due to the wave-particle interaction, a modification to the synchrotron radiation spectrum emitted by the runaway electron population is foreseen, exposing a possible experimental detection method for such an interaction.

Effect of bremsstrahlung radiation emission on fast electrons in plasmas

Bremsstrahlung radiation emission is an important energy loss mechanism for energetic electrons in plasmas. In this paper we investigate the effect of spontaneous bremsstrahlung emission on the momentum-space structure of the electron distribution, fully accounting for the emission of finite-energy photons. We find that electrons accelerated by electric fields can reach significantly higher energies than what is expected from energy-loss considerations. Furthermore, we show that the emission of soft photons can contribute significantly to the dynamics of electrons with an anisotropic distribution.