We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (
- Award ID(s):
- 1847905
- NSF-PAR ID:
- 10399403
- Date Published:
- Journal Name:
- Plasma Sources Science and Technology
- Volume:
- 31
- Issue:
- 10
- ISSN:
- 0963-0252
- Page Range / eLocation ID:
- 105001
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract β p= 0.25) collisionless ion–electron shocks with mass ratiom i/m e= 200, fast Mach number –4, and upstream magnetic field angleθ Bn= 55°–85° from the shock normal . It is known that shock electron heating is described by an ambipolar, -parallel electric potential jump, ΔB ϕ ∥, that scales roughly linearly with the electron temperature jump. Our simulations have –0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ ∥, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕ ∥in our low-β pshocks. We further focus on twoθ Bn= 65° shocks: a ( ) case with a long, 30d iprecursor of whistler waves along , and a ( ) case with a shorter, 5d iprecursor of whistlers oblique to both and ;B d iis the ion skin depth. Within the precursors,ϕ ∥has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the ,θ Bn= 65° case,ϕ ∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω pe/Ωce, andϕ ∥decreases by a factor of 2 asm i/m eis raised to the true proton–electron value of 1836. -
Abstract Electron cyclotron waves (whistlers) are commonly observed in plasmas near Earth and the solar wind. In the presence of nonlinear mirror modes, bursts of whistlers, usually called lion roars, have been observed within low magnetic field regions associated with these modes. In the intracluster medium (ICM) of galaxy clusters, the excitation of the mirror instability is expected, but it is not yet clear whether electron and ion cyclotron (IC) waves can also be present under conditions where gas pressure dominates over magnetic pressure (high
β ). In this work, we perform fully kinetic particle-in-cell simulations of a plasma subject to a continuous amplification of the mean magnetic field (B t ) to study the nonlinear stages of the mirror instability and the ensuing excitation of whistler and IC waves under ICM conditions. Once mirror modes reach nonlinear amplitudes, both whistler and IC waves start to emerge simultaneously, with subdominant amplitudes, propagating in low- regions, quasi-parallel toB (B t ). We show that the underlying source of excitation is the pressure anisotropy of electrons and ions trapped in mirror modes with loss-cone-type distributions. We also observe that IC waves play an essential role in regulating the ion pressure anisotropy at nonlinear stages. We argue that whistler and IC waves are a concomitant feature at late stages of the mirror instability even at highβ , and therefore, expected to be present in astrophysical environments like the ICM. We discuss the implications of our results for collisionless heating and dissipation of turbulence in the ICM. -
On electron kinetic scales, ions and electrons decouple, and electron velocity shear on electron inertial length ∼de can trigger electromagnetic (EM) electron Kelvin–Helmholtz instability (EKHI). In this paper, we present an analytic study of EM EKHI in an inviscid collisionless plasma with a step-function electron shear flow. We show that in incompressible collisionless plasma, the ideal electron frozen-in condition E+ve×B/c=0 must be broken for the EM EKHI to occur. In a step-function electron shear flow, the ideal electron frozen-in condition is replaced by magnetic flux conservation, i.e., ∇×(E+ve×B/c)=0, resulting in a dispersion relation similar to that of the standard ideal and incompressible magnetohydrodynamics KHI. The magnetic field parallel to the electron streaming suppresses the EM EKHI due to magnetic tension. The threshold for the EM mode of the EKHI is (k·ΔUe)2>ne1+ne2ne1ne2[ne1(vAe1·k)2+ne2(vAe2·k)2], where vAe=B/(4πmene)1/2, ΔUe, and ne are the electron streaming velocity shear and densities, respectively. The growth rate of the EM mode is γem∼Ωce, which is the electron gyro-frequency.more » « less
-
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.more » « less
-
ABSTRACT Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modelled as a collisional, fluid magnetohydrodynamic (MHD) discontinuity, despite the fact that many shocks occur in the collisionless regime. In light of this, using fully kinetic, 3D simulations of non-relativistic, parallel propagating collisionless shocks comprised of electron-positron plasma, we detail the deviation of collisionless shocks form MHD predictions for varying magnetization/Alfvénic Mach numbers, with particular focus on systems with Alfénic Mach numbers much smaller than sonic Mach numbers. We show that the shock compression ratio decreases for sufficiently large upstream magnetic fields, in agreement with theoretical predictions from previous works. Additionally, we examine the role of magnetic field strength on the shock front width. This work reinforces a growing body of work that suggest that modelling many astrophysical systems with only a fluid plasma description omits potentially important physics.more » « less