Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


Title: Kinetic simulations of strongly magnetized parallel shocks: deviations from MHD jump conditions
ABSTRACT Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modelled as a collisional, fluid magnetohydrodynamic (MHD) discontinuity, despite the fact that many shocks occur in the collisionless regime. In light of this, using fully kinetic, 3D simulations of non-relativistic, parallel propagating collisionless shocks comprised of electron-positron plasma, we detail the deviation of collisionless shocks form MHD predictions for varying magnetization/Alfvénic Mach numbers, with particular focus on systems with Alfénic Mach numbers much smaller than sonic Mach numbers. We show that the shock compression ratio decreases for sufficiently large upstream magnetic fields, in agreement with theoretical predictions from previous works. Additionally, we examine the role of magnetic field strength on the shock front width. This work reinforces a growing body of work that suggest that modelling many astrophysical systems with only a fluid plasma description omits potentially important physics.  more » « less
Award ID(s):
1936393
PAR ID:
10348637
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
509
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2084 to 2090
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. We have already devised a model for the kinetic history of the plasma through the shock front ( J. Plasma Phys. , vol. 84, issue 6, 2018, 905840604), allowing to self-consistently compute the downstream anisotropy, and hence the density jump, in terms of the upstream parameters. This model deals with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only one switch-on solution within a limited range of Alfvén Mach numbers, our model offers two kinds of solutions within a slightly different range of Alfvén Mach numbers. These two solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold. 
    more » « less
  2. Abstract Collisionless plasma shocks are a common feature of many space and astrophysical systems. They are sources of high-energy particles and nonthermal emission, channeling as much as 20% of the shock’s energy into nonthermal particles. The generation and acceleration of these nonthermal particles have been previously studied and shown to affect shock hydrodynamics to the zeroth order. In this work, we use self-consistent hybrid particle-in-cell simulations to examine the effect of self-generated nonthermal ion populations on the nature of collisionless, quasi-parallel shocks. Accelerated nonthermal particles downstream of the shock diffuse into the upstream region, taking energy away from the shock, which increases the compression ratio, slows the shock down, and flattens the nonthermal population’s spectral index for lower-Mach-number shocks. We show that this enhances shock compressibility when the heat flux is included in the Rankine–Hugoniot jump conditions, results that are roughly consistent with previous theories of “cosmic-ray-modified shocks.” Additionally, the simulation data show that heat flux and enthalpy flux cancels out in the upstream region, yielding a relatively simple, alternative closure for the jump conditions which accurately predict for the shock speed and compression ratio. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies in predictions of the arrival times of coronal mass ejections and the conflicting radio and X-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to predict shock evolution accurately. 
    more » « less
  3. The study of collisionless shocks and their role in cosmic ray acceleration has gained importance through observations and simulations, driving interest in reproducing these conditions in laboratory experiments using high-power lasers. In this work, we examine the role of three-dimensional (3D) effects in ion acceleration in quasi-perpendicular shocks under laboratory-relevant conditions. Using hybrid particle-in-cell (PIC) simulations (kinetic ions and fluid electrons), we explore how the Alfvénic and sonic Mach numbers, along with plasma beta, influence ion energization, unlocked only in 3D, and establish scaling criteria for when conducting 3D simulations is necessary. Our results show that efficient ion acceleration requires Alfvénic Mach numbers ≥25 and sonic Mach numbers ≥13, with plasma-β≤5. We theoretically found that, while two-dimensional (2D) simulations suffice for current laboratory-accessible shock conditions, 3D effects become crucial for shock velocities exceeding 1000 km/s and experiments sustaining the shock for at least 10 ns. We surveyed previous laboratory experiments on collisionless shocks and found that 3D effects are unimportant under those conditions, implying that one-dimensional and 2D simulations should be enough to model the accelerated ion spectra. However, we do find that the same experiments are realistically close to accessing the regime relevant to 3D effects, an exciting prospect for future laboratory efforts. We propose modifications to past experimental configurations to optimize and control 3D effects on ion acceleration. These proposed experiments could be used to benchmark plasma astrophysics kinetic codes and/or employed as controllable sources of energetic particles. 
    more » « less
  4. 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
  5. Shock waves in plasma are usually dealt with using magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing one to account for kinetic effects within a MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalogue, and a density jump that is notably less in the quasi-parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined. 
    more » « less