Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological fire-2 CR-magnetohydrodynamic simulations of L* galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics-driven hysteresis. SC models show a higher tendency to undergo ‘ejective’ feedback events due to a runaway buildup of CR pressure in dense gas due to the behaviour of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, interstellar medium gas, and B properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L* galaxies, and indicating the prospect for non-thermal radio continuum observations to constrain CR transport physics.
Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relatively strong dynamical coupling to the gas. In galaxies, CRs can be close to equipartition with the thermal, magnetic, and turbulent energy density in the interstellar medium, and can be dynamically very important in driving large-scale galactic winds. Similarly, CRs may provide a significant contribution to the pressure in the circumgalactic medium. In galaxy clusters, CRs may play a key role in addressing the classic cooling flow problem by facilitating efficient heating of the intracluster medium and preventing excessive star formation. Overall, the underlying physics of CR interactions with plasmas exhibit broad parallels across the entire range of scales characteristic of the interstellar, circumgalactic, and intracluster media. Here we present a review of the state-of-the-art of this field and provide a pedagogical introduction to cosmic ray plasma physics, including the physics of wave–particle interactions, acceleration processes, CR spatial and spectral transport, and important cooling processes. The field is ripe for discovery and will remain the subject of intense theoretical, computational, and observational research over the next decade with profound implications for the interpretation of the observations of stellar and supermassive black hole feedback spanning the entire width of the electromagnetic spectrum and multi-messenger data.
more » « less- NSF-PAR ID:
- 10478002
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- The Astronomy and Astrophysics Review
- Volume:
- 31
- Issue:
- 1
- ISSN:
- 0935-4956
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
ABSTRACT Black hole feedback plays a central role in shaping the circumgalactic medium (CGM) of elliptical galaxies. We systematically study the impact of plasma physics on the evolution of ellipticals by performing three-dimensional non-ideal magnetohydrodynamic simulations of the interactions of active galactic nucleus (AGN) jets with the CGM including magnetic fields, and cosmic rays (CRs) and their transport processes. We find that the physics of feedback operating on large galactic scales depends very sensitively on plasma physics operating on small scales. Specifically, we demonstrate that (i) in the purely hydrodynamical case, the AGN jets initially maintain the atmospheres in global thermal balance. However, local thermal instability generically leads to the formation of massive cold discs in the vicinity of the central black hole in disagreement with observations; (ii) including weak magnetic fields prevents the formation of the discs because local B-field amplification in the precipitating cold gas leads to strong magnetic breaking, which quickly extracts angular momentum from the accreting clouds. The magnetic fields transform the cold clouds into narrow filaments that do not fall ballistically; (iii) when plasma composition in the AGN jets is dominated by CRs, and CR transport is neglected, the atmospheres exhibit cooling catastrophes due to inefficient heat transfer from the AGN to CGM despite Coulomb/hadronic CR losses being present; (iv) including CR streaming and heating restores agreement with the observations, i.e. cooling catastrophes are prevented and massive cold central discs do not form. The AGN power is reduced as its energy is utilized efficiently.more » « less
-
more » « less
ABSTRACT We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV–TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live kinetic-magnetohydrodynamics galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighbourhood [Local interstellar medium (LISM)] observations (spectra, B/C, e+/e−, $\mathrm{\bar{p}}/\mathrm{p}$, 10Be/9Be, ionization, and γ-rays). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circumgalactic medium in realistic galaxies necessarily imposes an $\sim 10\,$ kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by $\sim 0.1-1\,$ kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov–Taylor stages of SNe remnants.
-
more » « less
ABSTRACT Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGNs) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure-dominated circumgalactic medium (CGM). But explicit CR-magnetohydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether those results can be applied to simulations that do not explicitly treat magnetic fields or resolved interstellar medium phase structure. We therefore present an intentionally extremely simplified ‘sub-grid’ model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytical models including some approximate CR effects on galactic (≳ kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently developed sub-grid models for radiative feedback, and allows for a simple constant parametrization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star–particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty.
-
more » « less
ABSTRACT Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work, we showed that streaming cosmic rays (CRs) destabilize sound waves in dilute ICM plasmas. Here, we show that CR streaming in the presence of gravity also destabilizes a pressure-balanced wave. We term this new instability the CR buoyancy instability (CRBI). In stark contrast to standard results without CRs, the pressure-balanced mode is highly compressible at short wavelengths due to CR streaming. Maximal growth rates are of order (pc/pg)β1/2ωff, where pc/pg is the ratio of CR pressure to thermal gas pressure, β is the ratio of thermal to magnetic pressure, and ωff is the free-fall frequency. The CRBI operates alongside buoyancy instabilities driven by background heat fluxes, i.e. the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI). When the thermal mean free path lmfp is ≪ the gas scale height H, the HBI/MTI set the growth rate on large scales, while the CRBI sets the growth rate on small scales. Conversely, when lmfp ∼ H and (pc/pg)β1/2 ≳ 1, CRBI growth rates exceed HBI/MTI growth rates even on large scales. Our results suggest that CR-driven instabilities may be partially responsible for the sound waves/weak shocks and turbulence observed in galaxy clusters. CR-driven instabilities generated near radio bubbles may also play an important role redistributing AGN energy throughout clusters.
