An official website of the United States government
Abstract The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.
more »
« less
A Comparison of Particle-in-cell and Hybrid Simulations of the Heliospheric Termination Shock
Abstract We compare hybrid (kinetic proton, fluid electron) and particle-in-cell (kinetic proton, kinetic electron) simulations of the solar wind termination shock with parameters similar to those observed by Voyager 2 during its crossing. The steady-state results show excellent agreement between the downstream variations in the density, plasma velocity, and magnetic field. The quasi-perpendicular shock accelerates interstellar pickup ions to a maximum energy limited by the size of the computational domain, with somewhat higher fluxes and maximal energies observed in the particle-in-cell simulation, likely due to differences in the cross-shock electric field arising from electron kinetic-scale effects. The higher fluxes may help address recent discrepancies noted between observations and large-scale hybrid simulations.
more »
« less
- Award ID(s):
- 2148653
- PAR ID:
- 10476372
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 959
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 4
- Size(s):
- Article No. 4
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration.more » « less
-
Large-amplitude electrostatic fluctuations are routinely observed by spacecraft upon traversal of collisionless shocks in the heliosphere. Kinetic simulations of shocks have struggled to reproduce the amplitude of such fluctuations, complicating efforts to un- derstand their influence on energy dissipation and shock structure. In this paper, 1D particle-in-cell simulations with realistic proton-to-electron mass ratio are used to show that in cases with upstream electron temperature Te exceeding the ion temperature Ti, the magnitude of the fluctuations increases with the electron plasma-to-cyclotron frequency ratio ωpe/Ωce, reaching realistic values at ωpe/Ωce ≳ 30. The large-amplitude fluctuations in the simulations are shown to be associated with electrostatic solitary structures, such as ion phase-space holes. In the cases where upstream temperature ratio is reversed, the magnitude of the fluctuations remains small.more » « less
-
Abstract The acceleration of charged particles by interplanetary shocks (IPs) can drain a nonnegligible fraction of the plasma pressure. In this study, we have selected 17 IPs observed in situ at 1 au by the Advanced Composition Explorer and the Wind spacecraft, and 1 shock at 0.8 au observed by Parker Solar Probe. We have calculated the time-dependent partial pressure of suprathermal and energetic particles (smaller and greater than 50 keV for protons and 30 keV for electrons, respectively) in both the upstream and downstream regions. The particle fluxes were averaged for 1 hr before and 1 hr after the shock time to remove short timescale effects. Using the MHD Rankine–Hugoniot jump conditions, we find that the fraction of the total upstream energy flux transferred to suprathermal and energetic downstream particles is typically ≲16%, in agreement with previous observations and simulations. Notably, by accounting for errors on all measured shock parameters, we have found that for any given fast magnetosonic Mach number,Mf< 7, the angle between the shock normal and average upstream magnetic field,θBn, is not correlated with the energetic particle pressure; in particular, the partial pressure of energized particles does not decrease forθBn≳ 45°. The downstream electron-to-proton energy ratio in the range ≳ 140 eV for electrons and ≳ 70 keV for protons exceeds the expected ∼1% and nears equipartition (>0.1) for the Wind events.more » « less
-
Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (vsh≈ (5500–22000) km s−1) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (vsh≈ 300 km s−1), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations.more » « less
