Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsvsh≈ 0.067–0.267c, Alfvén Mach numbers , sonic Mach numbers , as well as the proton-to-electron mass ratiosmi/me= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when , the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around ), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation.
more »
« less
This content will become publicly available on May 26, 2027
Efficient Particle Acceleration in 2.5-dimensional Hybrid-kinetic Simulations of Decaying Supersonic Plasma Turbulence
Abstract Collisionless, turbulent plasmas surround the Earth, from the magnetosphere to the intergalactic medium, and the fluctuations within them affect nearly every field in the space sciences, from space weather forecasts to theories of galaxy formation. Where turbulent motions become supersonic, their interactions can lead to the formation of shocks, which are known to efficiently energize ions to cosmic-ray energies. We present 2.5-dimensional hybrid-kinetic simulations of decaying, supersonic, nonrelativistic turbulence in a collisionless plasma using the codedHybridR. The ions are initially cold ( ), although they heat rapidly at the onset of turbulence. Turbulence within these simulations is highly compressible; after accounting for this compression by taking the omnidirectional power spectrum of thedensity-weightedvelocity field, we find turbulent spectra with power-law slopes of for low Mach numbers, in the inertial range, andα ≈ −2 for high Mach numbers. Ions embedded in the highly supersonic simulations are accelerated to high energies at efficiencies similar to those seen in shocks, despite being in a nonrelativistic regime and lacking the large-scale structure of a shock. We observe that particles are accelerated into a power-law spectrum, with a slope ofq ≈ 2.5 in (nonrelativistic) energy. We compare these results to those obtained from the theory and simulations of diffusive shock acceleration and discuss the astrophysical implications of this theoretical work.
more »
« less
- PAR ID:
- 10700681
- Publisher / Repository:
- IOP
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 1003
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 212
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The hot plasma in galaxy clusters, the intracluster medium, is expected to be shaped by subsonic turbulent motions, which are key for heating, cooling, and transport mechanisms. The turbulent motions contribute to the nonthermal pressure, which, if not accounted for, consequently imparts a hydrostatic mass bias. Accessing information about turbulent motions is thus of major astrophysical and cosmological interest. Characteristics of turbulent motions can be indirectly accessed through surface brightness fluctuations. This study expands on our pilot investigations of surface brightness fluctuations in the Sunyaev–Zel’dovich and in X-ray data by examining, for the first time, a large sample of 60 clusters using both SPT-SZ and XMM-Newton data and spans the redshift range 0.2 < z < 1.5, thus constraining the respective pressure and density fluctuations within 0.6R500. We deem density fluctuations to be of sufficient quality for 32 clusters, finding mild correlations between the peak of the amplitude spectra of density fluctuations and various dynamical parameters. We infer turbulent velocities from density fluctuations with an average Mach number , in agreement with numerical simulations. For clusters with inferred turbulent Mach numbers from fluctuations in both pressure, , and density, , we find broad agreement between and . Our results suggest either a bimodal or a skewed unimodal Mach number distribution, with the majority of clusters being turbulence-dominated (subsonic) while the remainder are shock-dominated (supersonic).more » « less
-
Abstract Understanding the nature of compressible fluctuations in a broad range of turbulent plasmas, from the intracluster medium to the solar wind, has been an active field of research in the past decades. Theoretical frameworks for weakly compressible magnetohydrodynamical turbulence in an inhomogeneous background magnetic field predict a linear scaling of the normalized mass density fluctuation (δρ/ρ0), as a function of the turbulent Mach number ( ), . However, so far, the scaling relation has been tested only using moderate to low plasma beta (β ≲ 1) solar wind observational data, where the compressibility is weakδρ/ρ0 ∼ 0.1. Here, we combine NASA’s Magnetospheric Multiscale Mission data in Earth’s magnetosheath, whereβ ∼ 10 is high, andβ ∼ 1/8 highly compressible magnetohydrodynamic turbulence simulations at unprecedented resolutions. Both show that holds across a broad range ofδρ/ρ0, , andβ, demonstrating that is a robust compressible turbulence relation, going beyond the asymptotics of the weakly compressible theory. We discuss the findings in the context of understanding the nature of strongly compressible turbulent fluctuations and the driving parameter in astrophysical and space plasmas.more » « less
-
Abstract We investigate the physics of quasi-parallel transrelativistic shocks propagating in weakly magnetized plasmas by means of long-duration two-dimensional particle-in-cell simulations. The structure of the shock precursor is shaped by a competition between the Bell instability and the Weibel (filamentation) instability. The Bell instability is dominant at relatively high magnetizations (σ≳ 10−3), whereas the Weibel instability prevails at lower magnetizations (σ≲ 10−4). Shocks with precursors shaped by Bell modes efficiently accelerate ions, converting a fractionεi∼ 0.2 of the upstream flow energy into downstream nonthermal ion energy. The maximum energy of nonthermal ions exhibits a Bohm scaling in time, as . A much smaller fractionεe≪ 0.1 of the upstream flow energy goes into downstream nonthermal electrons in the Bell regime. On the other hand, when the precursor is dominated by Weibel modes, the shock efficiently generates both nonthermal ions and electrons withεi∼εe∼ 0.1, albeit with a slower scaling for the maximum energy, . Our results are applicable to a wide range of transrelativistic shocks, including the termination shocks of extragalactic jets, the late stages of gamma-ray burst afterglows, and shocks in fast blue optical transients.more » « less
-
Abstract In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, ΩL, is the local Keplerian frequency ΩK. In terms of the standard dimensionless Shakura–Sunyaevαparameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by and , in whichHandcare the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where ΩL≥ ΩKand show that, under these conditions, the characteristic length and velocity scales are respectively and , where is twice the Rossby number. It follows that , where is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For , the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation.more » « less
An official website of the United States government
