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: 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
Award ID(s):
2206609 2206607
PAR ID:
10539873
Author(s) / Creator(s):
; ;
Publisher / Repository:
The Astrophysical Journal Letters
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
963
Issue:
2
ISSN:
2041-8205
Page Range / eLocation ID:
L44
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Relativistic magnetic turbulence has been proposed as a process for producing nonthermal particles in high-energy astrophysics. The particle energization may be contributed by both magnetic reconnection and turbulent fluctuations, but their interplay is poorly understood. It has been suggested that during magnetic reconnection the parallel electric field dominates the particle acceleration up to the lower bound of the power-law particle spectrum, but recent studies show that electric fields perpendicular to the magnetic field can play an important, if not dominant role. In this study, we carry out two-dimensional fully kinetic particle-in-cell simulations of magnetically dominated decaying turbulence in a relativistic pair plasma. For a fixed magnetization parameterσ0 = 20, we find that the injection energyεinjconverges with increasing domain size toεinj ≃ 10mec2. In contrast, the power-law index, the cut-off energy, and the power-law extent increase steadily with domain size. We trace a large number of particles and evaluate the contributions of the work done by the parallel (W) and perpendicular (W) electric fields during both the injection phase and the postinjection phase. We find that during the injection phase, theWcontribution increases with domain size, suggesting that it may eventually dominate injection for a sufficiently large domain. In contrast, on average, both components contribute equally during the postinjection phase, insensitive to the domain size. For high energy (ε ≫ εinj) particles,Wdominates the subsequent energization. These findings may improve our understanding of nonthermal particles and their emissions in astrophysical plasmas. 
    more » « less
  2. Abstract In the present study we investigate the response of the dayside ground magnetic field to the sequence of interplanetary magnetic field (IMF)BYchanges during the May 2024 geomagnetic storm. We pay particular attention to its extraordinarily large (>120 nT) and abrupt flip, and use GOES‐18 (G18) magnetic field measurements in the dayside magnetosheath as a time reference. In the dayside auroral zone, the northward magnetic component changed by as much as 4,300 nT from negative to positive indicating that the direction of the auroral electrojet changed from westward to eastward. The overall sequence was consistent with the conventional understanding of the IMFBYdriving of zonal ionospheric flows and Hall currents, which is also confirmed by a global simulation conducted for this storm. Surprisingly, however, the time delay from G18 to the ground increased significantly in time. The delay was 2–3 min for a sharpBYreduction ∼30 min prior to theBYflip, but it became as long as 10 min for the zero‐crossing of theBYflip. It is suggested that the prolonged time delay reflected the travel time from G18 to the reconnection site, which sensitively depends on the final velocity at the magnetopause, that is, the inflow velocity of the magnetic reconnection. Around theBYflip, the solar wind number density transiently exceeded 100 cm−3, and should have increased further through the bow shock crossing. It is suggested that this unusually dense plasma reduced the reconnection rate, and therefore, the solar wind‐magnetosphere energy coupling due to the extraordinary IMF. 
    more » « less
  3. Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure isβpi1∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed isMs∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperatureTe2exceeds the adiabatic expectationTe2,adby an amount T e 2 / T e 2 , ad 1 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratioTe1/Ti1, which we vary from 0.1 to unity. We investigate the heating physics forMs∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping atTe1/Ti1= 1, whereas heating byB-parallel electric fields (i.e., parallel to the local magnetic field) dominates atTe1/Ti1= 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows. 
    more » « less
  4. Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫σ, and form a hard power-law energy distribution,fi∝γ−1. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to eitherσin the strong cooling regime or the radiation burnoff limit,γsyn, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit,εc≈ 16 MeV, is only efficiently produced in the strong cooling regime,γsyn<σ. In this regime, we find that the spectral cutoff scales asεcut≈εc(σ/γsyn) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula. 
    more » « less
  5. Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫σ, and form a hard power-law energy distribution,fi∝γ−1. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to eitherσin the strong cooling regime or the radiation burnoff limit,γsyn, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit,εc≈ 16 MeV, is only efficiently produced in the strong cooling regime,γsyn<σ. In this regime, we find that the spectral cutoff scales asεcut≈εc(σ/γsyn) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula. 
    more » « less