Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal powerlaw distribution functions in the momentum (or energy) space,
Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for particle acceleration within the plasma’s coherent structures. Here, we investigate how intermittency modeling connects statistical fluctuations in turbulence to regions of highenergy dissipation. This connection is established by employing a generalized She–Leveque model to characterize the exponents
 Award ID(s):
 2308090
 NSFPAR ID:
 10495104
 Publisher / Repository:
 DOI PREFIX: 10.3847
 Date Published:
 Journal Name:
 The Astrophysical Journal
 Volume:
 964
 Issue:
 1
 ISSN:
 0004637X
 Format(s):
 Medium: X Size: Article No. 14
 Size(s):
 Article No. 14
 Sponsoring Org:
 National Science Foundation
More Like this

Abstract f (γ )d γ ∝γ ^{−α}d γ , whereγ is the Lorenz factor. We argue that in addition to exhibiting nonGaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particleincell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles,n , has a powerlaw distribution function, . We argue that the scaling exponents are related as $P(n)\mathit{dn}\phantom{\rule{0.25em}{0ex}}\propto \phantom{\rule{0.25em}{0ex}}{n}^{\beta}\mathit{dn}$β ≈α + 1, which is broadly consistent with numerical simulations. Nonspacefilling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence. 
Abstract We perform particleincell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electronpositron pairs. Various external magnetizations
σ ≲ 10^{−4}and pairloading factorsZ _{±}≲ 10 are studied, whereZ _{±}is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field whenσ exceeds a critical valueσ _{L}that decreases withZ _{±}. Atσ ≲σ _{L}the shock is mediated by particle scattering in the selfgenerated microturbulent fields, the strength and scale of which decrease withZ _{±}, leading to lowerσ _{L}. (2) The energy fraction carried by the postshock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per postshock electron scales as . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as ${\overline{E}}_{\mathrm{e}}\propto {\left({Z}_{\pm}+1\right)}^{1}$σ ≈ 5 × 10^{−6}. The ions then become essentially thermal with mean energy , while electrons form a nonthermal tail, extending from ${\overline{E}}_{\mathrm{i}}$ to $E\sim {\left({Z}_{\pm}+1\right)}^{1}{\overline{E}}_{\mathrm{i}}$ . When ${\overline{E}}_{\mathrm{i}}$σ = 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gammaray burst afterglows. 
Abstract Magnetic reconnection is often invoked as a source of highenergy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with largescale threedimensional kinetic particleincell simulations in electron–positron plasmas—that elucidates the physics governing the generation of powerlaw energy spectra in relativistic reconnection. Particles with Lorentz factor
γ ≳ 3σ (here,σ is the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is , where ${t}_{\mathrm{acc}}\sim \gamma \phantom{\rule{0.25em}{0ex}}{\eta}_{\mathrm{rec}}^{1}{\omega}_{\mathrm{c}}^{1}\simeq 20\phantom{\rule{0.25em}{0ex}}\gamma \phantom{\rule{0.25em}{0ex}}{\omega}_{\mathrm{c}}^{1}$η _{rec}≃ 0.06 is the inflow speed in units of the speed of light andω _{c}=eB _{0}/mc is the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert _{esc}, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret _{esc}in our simulations and find thatt _{esc}∼t _{acc}forσ ≳ few. This leads to a universal (i.e.,σ independent) powerlaw spectrum for the particles undergoing active acceleration, and ${\mathit{dN}}_{\mathrm{free}}/d\gamma \propto {\gamma}^{1}$ for the overall particle population. Our results help to shed light on the ubiquitous presence of powerlaw particle and photon spectra in astrophysical nonthermal sources. $\mathit{dN}/d\gamma \propto {\gamma}^{2}$ 
Abstract We present a new suite of numerical simulations of the starforming interstellar medium (ISM) in galactic disks using the TIGRESSNCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, selfgravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct raytracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over
, gas surface density Σ_{gas}∼ 5–150 $Z\prime \equiv Z/{Z}_{\odot}\sim 0.13$M _{⊙}pc^{−2}, and stellar surface density Σ_{star}∼ 1–50M _{⊙}pc^{−2}. The range of emergent star formation rate surface density is Σ_{SFR}∼ 10^{−4}–0.5M _{⊙}kpc^{−2}yr^{−1}, and ISM total midplane pressure isP _{tot}/k _{B}= 10^{3}–10^{6}cm^{−3}K, withP _{tot}equal to the ISM weight . For given Σ_{gas}and Σ_{star}, we find $\mathcal{W}$ . We provide an interpretation based on the pressureregulated feedbackmodulated (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 ${\mathrm{\Sigma}}_{\mathrm{SFR}}\propto Z{\prime}^{0.3}$ , while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon}}_{\mathrm{th}}\propto {\mathcal{W}}^{0.46}Z{\prime}^{0.53}$ . The reduction in Σ_{SFR}at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicitydependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. ${\mathrm{\Upsilon}}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{0.22}Z{\prime}^{0.18}$ 
Abstract While it is well known that cosmic rays (CRs) can gain energy from turbulence via secondorder Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time
becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gasCR fluid with diffusive CR transport show clear imprints of CRinduced damping, saturating at ${t}_{\mathrm{damp}}\sim \rho {v}^{2}/{\stackrel{\u0307}{E}}_{\mathrm{CR}}\propto {E}_{\mathrm{CR}}^{1}$ , where ${\stackrel{\u0307}{E}}_{\mathrm{CR}}\sim \tilde{\u03f5}$ is the turbulent energy input rate. In that case, almost all of the energy in largescale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering smallscale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering of $\tilde{\u03f5}$E ≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp largescale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticE _{CR}∼E _{g}scenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such largescale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e., . $\tilde{\u03f5}\gg \rho {v}^{3}/L$