Particle-in-cell (PIC) simulations have shown that relativistic collisionless shocks mediated by the Weibel instability accelerate ∼1% of incoming particles, while the majority are transmitted through the shock and become thermalized. The microphysical processes that determine whether an incoming particle will be transmitted or reflected are poorly understood. We study the microphysics of particle reflection in Weibel-mediated shocks by tracking a shell of test particles in a PIC simulation of a shock in pair plasma. We find that electrons in positron-dominated filaments and positrons in electron-dominated filaments efficiently reflect off of strong magnetic structures at the shock. To participate in diffusive shock acceleration, however, these reflected particles headed toward the upstream must avoid getting advected downstream. This is enabled by incoming filaments, which trap reflected particles carrying the same sign of current as the filaments. The final injection efficiency on the order of ∼1% thus results from the effectiveness of the initial reflection at the shock and the reflected particles’ probability of survival in the upstream postreflection. We develop a model that predicts the fraction of high-energy particles as a function of the properties of Weibel filamentation.
This content will become publicly available on June 1, 2025
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- PAR ID:
- 10542827
- Publisher / Repository:
- ApJ
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 968
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 17
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
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,
M f < 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. -
We have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We find that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can eventually trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq en_{\text{CR}}v_{\text{sh}}$ , where $en_{\text{CR}}$ is the seed charge density and $v_{\text{sh}}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than ${\sim}10~\text{yr}$ , thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.more » « less
-
Abstract The evolution of plasma entropy and the process of plasma energy redistribution at the collisionless plasma shock front are evaluated based on the high temporal resolution data from the four Magnetospheric Multiscale spacecraft during the crossing of the terrestrial bow shock. The ion distribution function has been separated into the populations with different characteristic behaviors in the vicinity of the shock: the upstream core population, the reflected ions, the gyrating ions, the ions trapped in the vicinity of the shock, and the downstream core population. The values of ion and electron moments (density, bulk velocity, and temperature) have been determined separately for these populations. It is shown that the solar wind core population bulk velocity slows down mainly in the ramp with the electrostatic potential increase but not in the foot region as it was supposed. The reflected ion population determines the foot region properties, so the proton temperature peak in the foot region is an effect of the relative motion of the different ion populations, rather than an actual increase in the thermal speed of any of the ion population. The ion entropy evaluated showed a significant increase across the shock: the enhancement of the ion entropy occurs in the foot of the shock front and at the ramp, where the reflected ions are emerging in addition to the upstream solar wind ions, the anisotropy growing to generate the bursts of ion-scale electrostatic waves. The entropy of electrons across the shock does not show a significant change: electron heating goes almost adiabatically.
-
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.