skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Interaction between Multiple Current Sheets and a Shock Wave: 2D Hybrid Kinetic Simulations
Abstract Particle acceleration behind a shock wave due to interactions between magnetic islands in the heliosphere has attracted attention in recent years. The downstream acceleration may yield a continuous increase of particle flux downstream of the shock wave. Although it is not obvious how the downstream magnetic islands are produced, it has been suggested that current sheets are involved in the generation of magnetic islands due to their interaction with a shock wave. We perform 2D hybrid kinetic simulations to investigate the interaction between multiple current sheets and a shock wave. In the simulation, current sheets are compressed by the shock wave and a tearing instability develops at the compressed current sheets downstream of the shock. As the result of this instability, the electromagnetic fields become turbulent and magnetic islands form well downstream of the shock wave. We find a “post-cursor” region in which the downstream flow speed normal to the shock wave in the downstream rest frame is decelerated to ∼ 1 V A immediately behind the shock wave, where V A is the upstream Alfvén speed. The flow speed then gradually decelerates to 0 accompanied by the development of the tearing instability. We also observe an efficient production of energetic particles above 100 E 0 during the development of the instability some distance downstream of the shock wave, where E 0 = m p V A 2 and m p is the proton mass. This feature corresponds to Voyager observations showing that the anomalous cosmic-ray intensity increase begins some distance downstream of the heliospheric termination shock.  more » « less
Award ID(s):
1655280
PAR ID:
10308614
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
922
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The possibility that charged particles are accelerated statistically in a “sea” of small-scale interacting magnetic flux ropes in the supersonic solar wind is gaining credence. In this Letter, we extend the Zank et al. statistical transport theory for a nearly isotopic particle distribution by including an escape term corresponding to particle loss from a finite acceleration region. Steady-state 1D solutions for both the accelerated particle velocity distribution function and differential intensity are derived. We show Ulysses observations of an energetic particle flux enhancement event downstream of a shock near 5 au that is inconsistent with the predictions of classical diffusive shock acceleration (DSA) but may be explained by local acceleration associated with magnetic islands. An automated Grad-Shafranov reconstruction approach is employed to identify small-scale magnetic flux ropes behind the shock. For the first time, the observed energetic particle “time-intensity” profile and spectra are quantitatively compared with theoretical predictions. The results show that stochastic acceleration by interacting magnetic islands accounts successfully for the observed (i) peaking of particle intensities behind the shock instead of at the shock front as standard DSA predicts; (ii) increase in the particle flux amplification factor with increasing particle energy; (ii) increase in distance between the particle intensity peak and the shock front with increasing energy; and (iv) hardening of particle power-law spectra with increasing distance downstream of the shock. 
    more » « less
  2. The possibility that charged particles are accelerated statistically in a “sea” of small-scale interacting magnetic flux ropes in the supersonic solar wind is gaining credence. In this Letter, we extend the Zank et al. statistical transport theory for a nearly isotopic particle distribution by including an escape term corresponding to particle loss from a finite acceleration region. Steady-state 1D solutions for both the accelerated particle velocity distribution function and differential intensity are derived. We show Ulysses observations of an energetic particle flux enhancement event downstream of a shock near 5 au that is inconsistent with the predictions of classical diffusive shock acceleration (DSA) but may be explained by local acceleration associated with magnetic islands. An automated Grad-Shafranov reconstruction approach is employed to identify small-scale magnetic flux ropes behind the shock. For the first time, the observed energetic particle “time-intensity” profile and spectra are quantitatively compared with theoretical predictions. The results show that stochastic acceleration by interacting magnetic islands accounts successfully for the observed (i) peaking of particle intensities behind the shock instead of at the shock front as standard DSA predicts; (ii) increase in the particle flux amplification factor with increasing particle energy; (ii) increase in distance between the particle intensity peak and the shock front with increasing energy; and (iv) hardening of particle power-law spectra with increasing distance downstream of the shock. 
    more » « less
  3. We investigate particle acceleration in an MHD-scale system of multiple current sheets by performing 2D and 3D MHD simulations combined with a test particle simulation. The system is unstable for the tearing-mode instability, and magnetic islands are produced by magnetic reconnection. Due to the interaction of magnetic islands, the system relaxes to a turbulent state. The 2D (3D) case both yield −5/3 (− 11/3 and −7/3) power-law spectra for magnetic and velocity fluctuations. Particles are efficiently energized by the generated turbulence, and form a power-law tail with an index of −2.2 and −4.2 in the energy distribution function for the 2D and 3D case, respectively. We find more energetic particles outside magnetic islands than inside. We observe super-diffusion in the 2D (∼ t 2.27 ) and 3D (∼ t 1.2 ) case in the energy space of energetic particles. 
    more » « less
  4. Abstract We investigate the interaction of turbulence with shock waves by performing 2D hybrid kinetic simulations. We inject force-free magnetic fields upstream that are unstable to the tearing-mode instability. The magnetic fields evolve into turbulence and interact with a shock wave whose sonic Mach number is 2.4. Turbulence properties, the total and normalized residual energy and the normalized cross helicity, change across the shock wave. While the energy of velocity and magnetic fluctuations is mostly distributed equally upstream, the velocity fluctuations are amplified dominantly downstream of the shock wave. The amplitude of turbulence spectra for magnetic, velocity, and density fluctuations are also increased at the shock wave while their spectral index remains unchanged. We compare our results with the Zank et al. model of turbulence transmission across a shock, and find that it provides a reasonable explanation for the spectral change across the shock wave. We find that particles are efficiently accelerated at the shock front, and a power-law spectrum forms downstream. This can be explained by diffusive shock acceleration, in which particles gain energy by being scattered upstream and downstream of a shock wave. The trajectory of an accelerated particle suggests that upstream turbulence plays a role scattering of particles. 
    more » « less
  5. Abstract The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extendingupstream of the shockfor up to seven fluxe-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy asϵ−1, consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivityκon the particle fluxF(nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ3/2that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration. 
    more » « less