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,
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
- NSF-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 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. -
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 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 We investigate a secondary proton beam instability coexisting with the ambient solar wind turbulence at 50
R ☉. Three-dimensional hybrid numerical simulations (particle ions and a quasi-neutralizing electron fluid) are carried out with the plasma parameters in the observed range. In the turbulent background, the particle distribution function, in particular the slope of the “bump-on-tail” responsible for the instability, is time-dependent and inhomogeneous. The presence of the turbulence substantially reduces the growth rate and saturation level of the instability. We derive magnetic power spectra from the observational data and perform a statistical analysis to evaluate the average turbulence intensity at 50R ☉. This information is used to link the observed frequency spectrum to the wavenumber spectrum in the simulations. We verify that Taylor’s frozen-in hypothesis is valid for this purpose to a sufficient extent. To reproduce the typical magnetic power spectrum of the instability observed concurrently with the background turbulence, an artificial spacecraft probe is run through the simulation box. The thermal-ion instabilities are often seen as power elevations in the kinetic range of scales above an extrapolation of the turbulence spectrum from larger scales. We show that the elevated power in the simulations is much higher than the background level. Therefore, the turbulence at the average intensity does not obscure the secondary proton beam instability, as opposed to the solar wind at 1 au, in which the ambient turbulence typically obscures thermal-ion instabilities. -
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.
-
Abstract Identifying the accelerators of Galactic cosmic ray (CR) protons with energies up to a few PeV (1015eV) remains a theoretical and observational challenge. Supernova remnants (SNRs) represent strong candidates because they provide sufficient energetics to reproduce the CR flux observed at Earth. However, it remains unclear whether they can accelerate particles to PeV energies, particularly after the very early stages of their evolution. This uncertainty has prompted searches for other source classes and necessitates comprehensive theoretical modeling of the maximum proton energy,
, accelerated by an arbitrary shock. While analytic estimates of have been put forward in the literature, they do not fully account for the complex interplay between particle acceleration, magnetic field amplification, and shock evolution. This paper uses a multizone, semianalytic model of particle acceleration based on kinetic simulations to place constraints on for a wide range of astrophysical shocks. In particular, we develop relationships between , shock velocity, size, and ambient medium. We find that SNRs can only accelerate PeV particles under a select set of circumstances, namely, if the shock velocity exceeds ∼104km s−1and escaping particles drive magnetic field amplification. However, older and slower SNRs may still produce observational signatures of PeV particles due to populations accelerated when the shock was younger. Our results serve as a reference for modelers seeking to quickly produce a self-consistent estimate of the maximum energy accelerated by an arbitrary astrophysical shock.1Presented as a thesis to the Department of Astronomy and Astrophysics, The University of Chicago, in partial fulfillment of the requirements for a Ph.D. degree.