skip to main content


Title: Numerical Modeling of Suprathermal Electron Transport in the Solar Wind: Effects of Whistler Turbulence with a Full Diffusion Tensor
Abstract The electron VDF in the solar wind consists of a Maxwellian core, a suprathermal halo, a field-aligned component strahl, and an energetic superhalo that deviates from the equilibrium. Whistler wave turbulence is thought to resonantly scatter the observed electron velocity distribution. Wave–particle interactions that contribute to Whistler wave turbulence are introduced into a Fokker–Planck kinetic transport equation that describes the interaction between the suprathermal electrons and the Whistler waves. A recent numerical approach for solving the Fokker–Planck kinetic transport equation has been extended to include a full diffusion tensor. Application of the extended numerical approach to the transport of solar wind suprathermal electrons influenced by Whistler wave turbulence is presented. Comparison and analysis of the numerical results with observations and diagonal-only model results are made. The off-diagonal terms in the diffusion tensor act to depress effects caused by the diagonal terms. The role of the diffusion coefficient on the electron heat flux is discussed.  more » « less
Award ID(s):
1655280
NSF-PAR ID:
10313999
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
924
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Many spacecraft fly within or through a natural and variable particle accelerator powered by the coupling between the magnetosphere and the solar wind: the Earth’s radiation belts. Determining the dominant pathways to plasma energization is a central challenge for radiation belt science and space weather alike. Inward radial transport from an external source was originally thought to be the most important acceleration process occurring in the radiation belts. Yet, when modeling relied on a radial diffusion equation including electron lifetimes, notable discrepancies in model-observation comparisons highlighted a need for improvement. Works by Professor Richard M. Thorne and others showed that energetic (hundreds of keV) electrons interacting with whistler-mode chorus waves could be efficiently accelerated to very high energies. The same principles were soon transposed to understand radiation belt dynamics at Jupiter and Saturn. These results led to a paradigm shift in our understanding of radiation belt acceleration, supported by observations of a growing peak in the radial profile of the phase space density for the most energetic electrons of the Earth’s outer belt. Yet, quantifying the importance of local acceleration at the gyroscale, versus large-scale acceleration associated with radial transport, remains controversial due to various sources of uncertainty. The objective of this review is to provide context to understand the variety of challenges associated with differentiating between the two main radiation belt acceleration processes: radial transport and local acceleration. Challenges range from electron flux measurement analysis to radiation belt modeling based on a three-dimensional Fokker-Planck equation. We also provide recommendations to inform future research on radiation belt radial transport and local acceleration. 
    more » « less
  2. Abstract A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space (r,u) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker-Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε=u-eφ(r) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in Argon gas at low discharge currents using a two-level excitation-ionization model and a “full-chemistry” model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the 2-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures. 
    more » « less
  3. Abstract

    Electron diffusion by whistler‐mode chorus waves is one of the key processes controlling the dynamics of relativistic electron fluxes in the Earth's radiation belts. It is responsible for the acceleration of sub‐relativistic electrons injected from the plasma sheet to relativistic energies as well as for their precipitation and loss into the atmosphere. Based on analytical estimates of chorus wave‐driven quasi‐linear electron energy and pitch‐angle diffusion rates, we provide analytical steady‐state solutions to the corresponding Fokker‐Planck equation for the relativistic electron distribution and flux. The impact on these steady‐state solutions of additional electromagnetic ion cyclotron waves, and of ultralow frequency waves are examined. Such steady‐state solutions correspond to hard energy spectra at 1–4 MeV, dangerous for satellite electronics, and represent attractors for the system dynamics in the presence of sufficiently strong driving by continuous injections of 10–300 keV electrons. Therefore, these analytical steady‐state solutions provide a simple means for estimating the most extreme electron energy spectra potentially encountered in the outer radiation belt, despite the great variability of injections and plasma conditions. These analytical steady‐state solutions are compared with numerical simulations based on the full Fokker‐Planck equation and with relativistic electron flux spectra measured by satellites during one extreme event and three strong events of high time‐integrated geomagnetic activity, demonstrating a good agreement.

     
    more » « less
  4. Abstract

     We present a stochastic field line mapping model where the interplanetary magnetic field lines are described by a density distribution function satisfying a Fokker–Planck equation that is solved numerically. Due to the spiral geometry of the nominal Parker field and to the evolving nature of solar wind turbulence, the heliospheric diffusion of the magnetic field lines is both heterogeneous and anisotropic, including a radial component. The longitudinal distributions of the magnetic field lines are shown to be close to circular Gaussian distributions, although they develop a noticeable skewness. The magnetic field lines emanating from the Sun are found to differ, on average, from the spirals predicted by Parker. Although the spirals remain close to Archimedean, they are here underwound, on average. Our model predicts a spiral angle that is smaller by ∼5° than the Parker spiral angle at Earth’s orbit for the same solar wind speed ofVsw= 400 km s−1. It also predicts an angular position on the solar disk of the best magnetically connected footpoint to an observer at 1 au that is shifted westward by ∼10° with respect to the Parker’s field model. This significantly changes the angle of the most probable magnetic connection between possible sources on the Sun and observers in the inner heliosphere. The results have direct implications for the heliospheric transport of “scatter-free” electrons accelerated in the aftermath of solar eruptions.

     
    more » « less
  5. null (Ed.)
    ABSTRACT We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl. 
    more » « less