skip to main content


Title: Whistler-regulated Magnetohydrodynamics: Transport Equations for Electron Thermal Conduction in the High-β Intracluster Medium of Galaxy Clusters
Abstract Transport equations for electron thermal energy in the high- β e intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds V w ∼ v te / β e ≪ v te dominate the turbulent spectrum, with v te the electron thermal speed and β e ≫ 1 the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when L c > L / β e , with L c the classical electron mean free path and L the electron temperature scale length, and (2) in the whistler-dominated regime the electron thermal flux is controlled by both advection at V w and a comparable diffusive term. The findings suggest whistlers limit electron heat flux over large regions of the ICM, including locations unstable to isobaric condensation. Consequences include: (1) the Field length decreases, extending the domain of thermal instability to smaller length scales, (2) the heat flux temperature dependence changes from T e 7 / 2 / L to V w nT e ∼ T e 1 / 2 , (3) the magneto-thermal- and heat-flux-driven buoyancy instabilities are impaired or completely inhibited, and (4) sound waves in the ICM propagate greater distances, as inferred from observations. This description of thermal transport can be used in macroscale ICM models.  more » « less
Award ID(s):
1715140 2109083
NSF-PAR ID:
10312802
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
923
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A heat flux in a high- $\unicode[STIX]{x1D6FD}$ plasma with low collisionality triggers the whistler instability. Quasilinear theory predicts saturation of the instability in a marginal state characterized by a heat flux that is fully controlled by electron scattering off magnetic perturbations. This marginal heat flux does not depend on the temperature gradient and scales as $1/\unicode[STIX]{x1D6FD}$ . We confirm this theoretical prediction by performing numerical particle-in-cell simulations of the instability. We further calculate the saturation level of magnetic perturbations and the electron scattering rate as functions of $\unicode[STIX]{x1D6FD}$ and the temperature gradient to identify the saturation mechanism as quasilinear. Suppression of the heat flux is caused by oblique whistlers with magnetic-energy density distributed over a wide range of propagation angles. This result can be applied to high- $\unicode[STIX]{x1D6FD}$ astrophysical plasmas, such as the intracluster medium, where thermal conduction at sharp temperature gradients along magnetic-field lines can be significantly suppressed. We provide a convenient expression for the amount of suppression of the heat flux relative to the classical Spitzer value as a function of the temperature gradient and $\unicode[STIX]{x1D6FD}$ . For a turbulent plasma, the additional independent suppression by the mirror instability is capable of producing large total suppression factors (several tens in galaxy clusters) in regions with strong temperature gradients. 
    more » « less
  2. 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 fieldB(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-Bregions, quasi-parallel toB(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.

     
    more » « less
  3. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–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: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis 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 theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

     
    more » « less
  4. Abstract

    We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with aT>Ttemperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel instabilities. Electron velocity distributions taken after the system has reached a steady state show that these instabilities inhibit the heat flux and drive the total distributions toward isotropy. Electron trajectories in velocity space show a circular-like diffusion along constant energy surfaces in the wave frame. The key parameter controlling the scattering rate is the average speed, or drift speedvd, of the heat flux compared with the electron Alfvén speedvAe, with higher drift speeds producing stronger fluctuations and a more significant reduction of the heat flux. Reducing the density of the electrons carrying the heat flux by 50% does not significantly affect the scattering rate. A scaling law for the electron scattering rate versusvd/vAeis deduced from the simulations. The implications of these results for understanding energetic electron transport during energy release in solar flares are discussed.

     
    more » « less
  5. Abstract

    We present particle-in-cell simulations of a combined whistler heat flux and temperature anisotropy instability that is potentially operating in the solar wind. The simulations are performed in a uniform plasma and initialized with core and halo electron populations typical of the solar wind beyond about 0.3 au. We demonstrate that the instability produces whistler-mode waves propagating both along (anti-sunward) and opposite (sunward) to the electron heat flux. The saturated amplitudes of both sunward and anti-sunward whistler waves are strongly correlated with their initial linear growth rates,Bw/B0(γ/ωce)ν, where for typical electron betas we have 0.6 ≲ν≲ 0.9. We show that because of the relatively large spectral width of the whistler waves, the instability saturates through the formation of quasi-linear plateaus around the resonant velocities. The revealed correlations of whistler wave amplitudes and spectral widths with electron beta and temperature anisotropy are consistent with solar wind observations. We show that anti-sunward whistler waves result in an electron heat flux decrease, while sunward whistler waves actually lead to an electron heat flux increase. The net effect is the electron heat flux suppression, whose efficiency is larger for larger electron betas and temperature anisotropies. The electron heat flux suppression can be up to 10%–60% provided that the saturated whistler wave amplitudes exceed about 1% of the background magnetic field. The experimental applications of the presented results are discussed.

     
    more » « less