An official website of the United States government
Abstract We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-L⋆ starburst, and L⋆ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of γ-ray emission from nearby and starburst galaxies. We reproduce the γ-ray observations of dwarf and L⋆ galaxies with constant isotropic diffusion coefficient κ ∼ 3 × 1029 cm2 s−1. Advection-only and streaming-only models produce order-of-magnitude too large γ-ray luminosities in dwarf and L⋆ galaxies. We show that in models that match the γ-ray observations, most CRs escape low-gas-density galaxies (e.g. dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for γ-ray emissivities. Models where CRs are “trapped” in the star-forming disk have lower star formation efficiency, but these models are ruled out by γ-ray observations. For models with constant κ that match the γ-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc.
more »
« less
Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport
Abstract Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000M⊙MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low (ζ≲ 2 × 10−19s−1) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated (ζ≲ 6 × 10−18), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.
more »
« less
- Award ID(s):
- 2108318
- PAR ID:
- 10541519
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 973
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 16
- Size(s):
- Article No. 16
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 ratevA· ∇PCR/PCR∼vA/L0, 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
-
ABSTRACT We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass haloes hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, and collisional and streaming losses, with constant parallel diffusivity $$\kappa \sim 3\times 10^{29}\, \mathrm{cm^2\ s^{-1}}$$ chosen to match γ-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass haloes at z ≲ 1–2. The gas in these ‘CR-dominated’ haloes differs significantly from runs without CRs: the gas is primarily cool (a few $${\sim}10^{4}\,$$ K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ‘low’ and ‘mid’ ions in this diffuse cool gas is dominated by photoionization, with O vi columns $${\gtrsim}10^{14.5}\, \mathrm{cm^{-2}}$$ at distances $${\gtrsim}150\, \mathrm{kpc}$$. CR and thermal gas pressure are locally anticorrelated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same haloes are primarily warm/hot ($$T\gtrsim 10^{5}\,$$K) with thermal pressure balancing gravity, collisional ionization dominates, O vi columns are lower and Ne viii higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.more » « less
-
Abstract We study the propagation of mildly relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta et al. to postprocess three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the postprocessed simulations to make predictions regarding the potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extraplanar region (∣ z ∣ > 0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR.more » « less
-
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCHframework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCHby implementing the N-body code PETAR, thereby enabling TORCHto handle massive clusters forming from 106M⊙clouds with ≥105individual stars. We present results from TORCHsimulations of star clusters forming from 104, 105, and 106M⊙turbulent spherical gas clouds (named M4, M5, M6) of radiusR= 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 101,2,3M⊙pc−2, after a free-fall time oftff= 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 areϵ⋆ = M⋆/Mcloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies ofϵinst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 102M⊙pc−2) and massive (Mcloud≥ 105M⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars.more » « less
