An official website of the United States government
Abstract Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma β . Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low- β environments like the interstellar medium but remain unchanged in high- β environments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss rate v A · ∇ P CR / P CR ∼ v A / L 0 , due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength over β ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent spectrum.
more »
« less
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 time t damp ∼ ρ v 2 / E ̇ CR ∝ E CR − 1 becomes 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 E ̇ 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 of E ≳ 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 realistic E CR ∼ E g scenarios, 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., ϵ ˜ ≫ ρ v 3 / L .
more »
« less
- Award ID(s):
- 1911198
- PAR ID:
- 10465360
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 955
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 64
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) – high-energy charged particles accelerated in supernovae and young stars – can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfvén Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfvénic streaming of CRs enhances galactic winds, in this paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disc and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases.more » « less
-
Abstract The structure of shocks and turbulence are strongly modified during the acceleration of cosmic rays (CRs) at a shock wave. The pressure and the collisionless viscous stress decelerate the incoming thermal gas and thus modify the shock structure. A CR streaming instability ahead of the shock generates the turbulence on which CRs scatter. The turbulent magnetic field in turn determines the CR diffusion coefficient and further affects the CR energy spectrum and pressure distribution. The dissipation of turbulence contributes to heating the thermal gas. Within a multicomponent fluid framework, CRs and thermal gas are treated as fluids and are closely coupled to the turbulence. The system equations comprise the gas dynamic equations, the CR pressure evolution equation, and the turbulence transport equations, and we adopt typical parameters for the hot ionized interstellar medium. It is shown that the shock has no discontinuity but possesses a narrow but smooth transition. The self-generated turbulent magnetic field is much stronger than both the large-scale magnetic field and the preexisting turbulent magnetic field. The resulting CR diffusion coefficient is substantially suppressed and is more than three orders smaller near the shock than it is far upstream. The results are qualitatively consistent with certain observations.more » « less
-
ABSTRACT Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ∼MeV-TeV, these fluctuations are usually attributed either to ‘extrinsic turbulence’ (ET) – a cascade from larger scales – or ‘self-confinement’ (SC) – self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed ‘live’ numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally required scattering rates disagree with ET predictions by orders of magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvénic streaming are unstable. In live simulations, they revert to either arbitrarily rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints.more » « less
-
Abstract Feedback processes in galaxies dictate their structure and evolution. Baryons can be cycled through stars, which inject energy into the interstellar medium in supernova explosions, fueling multiphase galactic winds. Cosmic rays (CRs) accelerated at supernova remnants are an important component of feedback. CRs can effectively contribute to wind driving; however, their impact heavily depends on the assumed CR transport model. We run high-resolution “tallbox” simulations of a patch of a galactic disk using the moving mesh magnetohydrodynamics code Arepo, including varied CR implementations and the Crispnonequilibrium thermochemistry model. We characterize the impact of CR feedback on star formation and multiphase outflows. While CR-driven winds are able to supply energy to a global-scale wind, a purely thermal wind loses most of its energy by the time it reaches 3 kpc above the disk midplane. We further find that the adopted CR transport model significantly affects the steady state of the wind. In the model with CR advection, streaming, diffusion, and nonlinear Landau damping, CRs provide very strong feedback. Additionally, accounting for ion-neutral damping (IND) decouples CRs from the cold ISM, which reduces the impact of CRs on the star formation rate. Nevertheless, CRs in this most realistic model are able to accelerate warm gas and levitate cool gas in the wind but have little effect on cold gas and hot gas. This model displays moderate mass loading and significant CR energy loading, demonstrating that IND does not prevent CRs from providing effective feedback.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3847/1538-4357/aceef9
