skip to main content


Title: Cosmic-Ray Drag and Damping of Compressive Turbulence
Abstract

While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order 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 timetdampρv2/ĖCRECR1becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating atĖCRϵ˜, whereϵ˜is the turbulent energy input rate. In that case, almost all of the energy in large-scale 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 small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticECREgscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,ϵ˜ρv3/L.

 
more » « less
Award ID(s):
1911198
NSF-PAR ID:
10463840
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
955
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 64
Size(s):
["Article No. 64"]
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 report on a search for electron antineutrinos (ν¯e) from astrophysical sources in the neutrino energy range 8.3–30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we find no significant excess over background model predictions. Assuming several supernova relic neutrino spectra, we give upper flux limits of 60–110 cm−2s−1(90% confidence level, CL) in the analysis range and present a model-independent flux. We also set limits on the annihilation rates for light dark matter pairs to neutrino pairs. These data improve on the upper probability limit of8B solar neutrinos converting intoν¯e,Pνeν¯e<3.5×105(90% CL) assuming an undistortedν¯eshape. This corresponds to a solarν¯eflux of 60 cm−2s−1(90% CL) in the analysis energy range.

     
    more » « less
  3. Abstract

    We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10−4and pair-loading 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σLthat decreases withZ±. AtσσLthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ±, leading to lowerσL. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales asE¯eZ±+11. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10−6. The ions then become essentially thermal with mean energyE¯i, while electrons form a nonthermal tail, extending fromEZ±+11E¯itoE¯i. Whenσ= 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 gamma-ray burst afterglows.

     
    more » « less
  4. Abstract

    We use ALMA observations of CO(2–1) in 13 massive (M*≳ 1011M) poststarburst galaxies atz∼ 0.6 to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The poststarburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (Data Release 14) based on their spectral shapes, as part of the Studying QUenching at Intermediate-z Galaxies: Gas, anguLarmomentum, and Evolution (SQuIGGLE) program. Early results showed that two poststarburst galaxies host large H2reservoirs despite their low inferred star formation rates (SFRs). Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, withMH2109M. Given their high stellar masses, this mass limit corresponds to an average gas fraction offH2MH2/M*7%or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theDn4000 spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched ≲150 Myr prior to observation host detectable CO(2–1) emission, while older poststarburst galaxies are undetected. The large H2reservoirs and low SFRs in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100–200 Myr, theSQuIGGLEgalaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H2scaling relations.

     
    more » « less
  5. Abstract

    Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofRB≈ 2 × 105GM/c2, whereMis the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r−1rather thanr−3/2and the accretion rateṀis consequently suppressed by over 2 orders of magnitude below the Bondi rateṀB. We find continuous energy feedback from the accretion flow to the external medium at a level of102Ṁc25×105ṀBc2. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

     
    more » « less