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.


This content will become publicly available on November 1, 2025

Title: Electron Acceleration at Quasi-parallel Nonrelativistic Shocks: A 1D Kinetic Survey
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 M A = 5 40 , sonic Mach numbers M s = 5 160 , 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 M A 10 , 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 p i , e / m i v sh 3 m i , e / m i ), 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
Award ID(s):
2308021 2009326 2010240 2206607
PAR ID:
10555552
Author(s) / Creator(s):
; ;
Publisher / Repository:
ApJ
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
976
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
10
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The polarization of the cosmic microwave background is rich in information but obscured by foreground emission from the Milky Way’s interstellar medium (ISM). To uncover relationships between the underlying turbulent ISM and the foreground power spectra, we simulated a suite of driven, magnetized, turbulent models of the ISM, varying the fluid properties via the sonic Mach number, M S , and magnetic (Alfvén) Mach number, M A . We measure the power spectra of density (ρ), velocity (v), magnetic field (H), total projected intensity (T), parity-even polarization (E), and parity-odd polarization (B). We find that the slopes of all six quantities increase with M S . Most increase with M A , while the magnetic field spectrum steepens with M A . By comparing spectral slopes ofEandBto those measured by Planck, we infer typical values of M S and M A for the ISM. As the fluid velocity increases, M S > 4 , the ratio of BB power to EE power increases to approach a constant value near the Planck-observed value of ∼0.5, regardless of the magnetic field strength. We also examine correlation coefficients between projected quantities, and find thatrTE≈ 0.3, in agreement with Planck, for appropriate combinations of M S and M A . Finally, we consider parity-violating correlationsrTBandrEB
    more » « less
  2. Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z Z / Z 0.1 - 3 , gas surface density Σgas∼ 5–150Mpc−2, and stellar surface density Σstar∼ 1–50Mpc−2. The range of emergent star formation rate surface density is ΣSFR∼ 10−4–0.5Mkpc−2yr−1, and ISM total midplane pressure isPtot/kB= 103–106cm−3K, withPtotequal to the ISM weight W . For given Σgasand Σstar, we find Σ SFR Z 0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒ th W 0.46 Z 0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒ turb + mag W 0.22 Z 0.18 . The reduction in ΣSFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 
    more » « less
  3. Abstract The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of M 50 = 10 8 M at the 1σlevel; however, a tail toward low M 50 prevents a 2σmeasurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on M 50 at the 1σlevel, but the 2σtail toward low M 50 persists. We project that observations of one (two) complete satellite populations can constrain warm DM models withmWDM≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010M, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics. 
    more » « less
  4. Abstract The Kruskal–Szekeres coordinate construction for the Schwarzschild spacetime could be interpreted simply as a squeezing of thet-line into a single point, at the event horizon r = 2 M . Starting from this perspective, we extend the Kruskal charting to spacetimes with two horizons, in particular the Reissner–Nordström manifold, M RN . We develop a new method to construct Kruskal-like coordinates through casting the metric in new null coordinates, and find two algebraically distinct ways to chart M RN , referred to as classes: type-I and type-II within this work. We pedagogically illustrate our method by crafting two compact, conformal, and global coordinate systems labeled GK I and GK II as an example for each class respectively, and plot the corresponding Penrose diagrams. In both coordinates, the metric differentiability can be promoted to C in a straightforward way. Finally, the conformal metric factor can be written explicitly in terms of thetandrfunctions for both types of charts. We also argued that the chart recently reported in Soltani (2023 arXiv:2307.11026) could be viewed as another example for the type-II classification, similar to GK II
    more » « less
  5. Abstract We survey 20 reconnection outflow events observed by Magnetospheric MultiScale in the low-βand high-Alfvén-speed regime of the Earth’s magnetotail to investigate the scaling of ion bulk heating produced by reconnection. The range of inflow Alfvén speeds (800–4000 km s−1) and inflow ionβ(0.002–1) covered by this study is in a plasma regime that could be applicable to the solar corona and flare environments. We find that the observed ion heating increases with increasing inflow (upstream) Alfvén speed,VA, based on the reconnecting magnetic field and the upstream plasma density. However, ion heating does not increase linearly as a function of available magnetic energy per particle, m i V A 2 . Instead, the heating increases progressively less as m i V A 2 rises. This is in contrast to a previous study using the same data set, which found that electron heating in this high-Alfvén-speed and low-βregime scales linearly with m i V A 2 , with a scaling factor nearly identical to that found for the low-VAand high-βmagnetopause. Consequently, the ion-to-electron heating ratio in reconnection exhausts decreases with increasing upstreamVA, suggesting that the energy partition between ions and electrons in reconnection exhausts could be a function of the available magnetic energy per particle. Finally, we find that the observed difference in ion and electron heating scaling may be consistent with the predicted effects of a trapping potential in the exhaust, which enhances electron heating, but reduces ion heating. 
    more » « less