skip to main content


Title: Long-term Evolution of Relativistic Unmagnetized Collisionless Shocks
Abstract

We study a relativistic collisionless electron–positron shock propagating into an unmagnetized ambient medium using 2D particle-in-cell simulations of unprecedented duration and size. The shock generates intermittent magnetic structures of increasingly larger size as the simulation progresses. Toward the end of our simulation, at around 26,000 plasma times, the magnetic coherence scale approachesλ∼ 100 plasma skin depths, both ahead and behind the shock front. We anticipate a continued growth ofλbeyond the time span of our simulation, as long as the shock accelerates particles to increasingly higher energies. The post-shock field is concentrated in localized patches, which maintain a local magnetic energy fractionεB∼ 0.1. Particles randomly sampling the downstream fields spend most of their time in low field regions (εB≪ 0.1) but emit a large fraction of the synchrotron power in the localized patches with strong fields (εB∼ 0.1). Our results have important implications for models of gamma-ray burst afterglows.

 
more » « less
NSF-PAR ID:
10494965
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
963
Issue:
2
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L44
Size(s):
["Article No. L44"]
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–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: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis 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 theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

     
    more » « less
  2. Abstract

    We present a statistical analysis of >2,100 bipolar electrostatic solitary waves (ESWs) collected from 10 quasi‐perpendicular Earth's bow shock crossings by Magnetospheric Multiscale spacecraft. We developed and implemented a correction procedure for reconstruction of actual electric fields, velocities, and other properties of ESW, whose spatial scales are typically comparable with or smaller than spatial distance between voltage‐sensitive probes. We found that more than 95% of the ESW are of negative polarity with amplitudes typically below a few Volts and 0.1Te(5–30 V or 0.1–0.3Tefor a few percent of ESW), spatial scales of 10–100 m orλD–10λD, and velocities from a few tens to a few hundred km/s that is on the order of local ion‐acoustic speed. The spatial scales of ESW are correlated with local Debye lengthλD. The ESW have electric fields generally oblique to magnetic field and they propagate highly oblique to shock normalN; more than 80% of ESW propagate within 30° of the shock planeLM. In the shock plane, ESW typically propagates within a few tens of degrees of local magnetic field projectionBLMand preferentially opposite toN × BLM. We argue that the ESW of negative polarity are ion holes produced by ion‐ion streaming instabilities. We estimate ion hole lifetimes to be 10–100 ms, or 1–10 km in terms of traveling distance. The revealed statistical properties will be useful for quantitative studies of electron thermalization in the Earth's bow shock.

     
    more » « less
  3. 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
  4. Abstract

    Long periods of strong southward magnetic fields are known to be the primary cause of intense geomagnetic storms. The majority of such events are caused by the passage over Earth of a magnetic ejecta. Irrespective of the interplanetary cause, fast‐forward shocks often precede such strong southwardBzperiods. Here we first look at all long periods of strong southward magnetic fields as well as fast‐forward shocks measured by theWindspacecraft in a 22.4‐year span. We find that 76% of strong southwardBzperiods are preceded within 48 hr by at least a fast‐forward shock, but only about 23% of all shocks are followed within 48 hr by strong southwardBzperiods. Then, we devise a threshold‐based probabilistic forecasting method based on the shock properties and the pre‐shock near‐Earth solar wind plasma and interplanetary magnetic field characteristics adopting asuperposed epoch analysis‐like approach. Our analysis shows that the solar wind conditions in the 30‐min interval around the arrival of fast‐forward shocks have a significant contribution to the prediction of long‐duration southwardBzperiods. This probabilistic model may provide on average a 14‐hr warning time for an intense and long‐duration southwardBzperiod. Evaluating the forecast capability of the model through a statistical and skill score‐based approach reveals that it outperforms a coin‐flipping forecast. By using the information provided by the arrival of a fast‐forward shock at L1, this model represents a marked improvement over similar forecasting methods. We outline a number of future potential improvements.

     
    more » « less
  5. Abstract

    This paper presents a stochastic three-dimensional focused transport simulation of solar energetic particles (SEPs) produced by a data-driven coronal mass ejection (CME) shock propagating through a data-driven model of coronal and heliospheric magnetic fields. The injection of SEPs at the CME shock is treated using diffusive shock acceleration of post-shock suprathermal solar wind ions. A time-backward stochastic simulation is employed to solve the transport equation to obtain the SEP time–intensity profile at any location, energy, and pitch angle. The model is applied to a SEP event on 2020 May 29, observed by STEREO-A close to ∼1 au and by Parker Solar Probe (PSP) when it was about 0.33 au away from the Sun. The SEP event was associated with a very slow CME with a plane-of-sky speed of 337 km s−1at a height below 6RSas reported in the SOHO/LASCO CME catalog. We compute the time profiles of particle flux at PSP and STEREO-A locations, and estimate both the spectral index of the proton energy spectrum for energies between ∼2 and 16 MeV and the equivalent path length of the magnetic field lines experienced by the first arriving SEPs. We find that the simulation results are well correlated with observations. The SEP event could be explained by the acceleration of particles by a weak CME shock in the low solar corona that is not magnetically connected to the observers.

     
    more » « less