The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM coevolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (“thermalization”) from a cool, turbulence-supported phase to a virial-temperature, thermally supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high redshift in ultrafaint progenitors and at low redshift in classical
We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind with large fluctuations in density and temperature. In most of our simulations, the heating eventually overwhelms cooling, leading to a rapid increase in temperature and a thermally driven wind; the exception to this is in galaxies with the shallowest potentials, which produce nearly isothermal $T \approx 10^4\,$ K winds driven by CR pressure. Many of the time-averaged wind solutions found here have a remarkable critical point structure, with two critical points. Scaled to real galaxies, we find mass outflow rates $\dot{M}$ somewhat larger than the observed star-formation rate in low-mass galaxies, and an approximately ‘energy-like’ scaling $\dot{M} \propto v_{\rm esc}^{-2}$. The winds accelerate slowly and reach asymptotic wind speeds of only ∼0.4vesc. The total wind power is $\sim 1~{{\ \rm per\ cent}}$ of the power from supernovae, suggesting inefficient preventive CR feedback for the physical conditions modelled here. We predict significant spatially extended emission and absorption lines from 104–105.5 K gas; this may correspond to extraplanar diffuse ionized gas seen in star-forming galaxies.
more » « less- Award ID(s):
- 2107872
- NSF-PAR ID:
- 10439945
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 524
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 6374-6391
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract M vir∼ 1011M ⊙dwarfs, while Milky Way–mass halos undergo the transition atz ≈ 0.5. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multiphase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. -
ABSTRACT In several models of galaxy formation feedback occurs in cycles or mainly at high redshift. At times and in regions where feedback heating is ineffective, hot gas in the galaxy halo is expected to form a cooling flow, where the gas advects inward on a cooling timescale. Cooling flow solutions can thus be used as a benchmark for observations and simulations to constrain the timing and extent of feedback heating. Using analytic calculations and idealized 3D hydrodynamic simulations, we show that for a given halo mass and cooling function, steady-state cooling flows form a single-parameter family of solutions, while initially hydrostatic gaseous haloes converge on one of these solutions within a cooling time. The solution is thus fully determined once either the mass inflow rate ${\dot{M}}$ or the total halo gas mass are known. In the Milky Way halo, a cooling flow with ${\dot{M}}$ equal to the star formation rate predicts a ratio of the cooling time to the free-fall time of ∼10, similar to some feedback-regulated models. This solution also correctly predicts observed $\rm{O\,{\small VII}}$ and $\rm{O\,{\small VIII}}$ absorption columns, and the gas density profile implied by $\rm{O\,{\small VII}}$ and $\rm{O\,{\small VIII}}$ emission. These results suggest ongoing heating by feedback may be negligible in the inner Milky-Way halo. Extending similar solutions out to the cooling radius however underpredicts observed $\rm{O\,{\small VI}}$ columns around the Milky-Way and around other low-redshift star-forming galaxies. This can be reconciled with the successes of the cooling flow model with either a mechanism which preferentially heats the $\rm{O\,{\small VI}}$-bearing outer halo, or alternatively if $\rm{O\,{\small VI}}$ traces cool photoionized gas beyond the accretion shock. We also demonstrate that the entropy profiles of some of the most relaxed clusters are reasonably well described by a cooling flow solution.more » « less
-
more » « less
ABSTRACT We use FIRE simulations to study disc formation in z ∼ 0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar discs is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation prior to joining the galaxy. Among galaxies with a high fraction ($\gt 70{{\ \rm per\ cent}}$) of their young stars in a thin disc (h/R ∼ 0.1), we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (≳20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to $\lesssim 10^4\, {\rm K}$; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disc, and while cooling transitions from a quasi-spherical spatial configuration to a more-flattened, disc-like configuration. We show that the existence of this ‘rotating cooling flow’ accretion mode is strongly correlated with the fraction of stars forming in a thin disc, using a sample of 17 z ∼ 0 galaxies spanning a halo mass range of 1010.5 M⊙ ≲ Mh ≲ 1012 M⊙ and stellar mass range of 108 M⊙ ≲ M⋆ ≲ 1011 M⊙. Notably, galaxies with a thick disc or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion on to the galaxy are likely a necessary condition for the formation of thin, star-forming disc galaxies in a ΛCDM universe.
-
more » « less
ABSTRACT Cosmic rays (CRs) are thought to be an important feedback mechanism in star-forming galaxies. They can provide an important source of pressure support and possibly drive outflows. We perform multidimensional CR magnetohydrodynamic simulations including transport by streaming and diffusion to investigate wind launching from an initially hydrostatic atmosphere by CRs. We estimate a characteristic Eddington limit on the CR flux for which the CR force exceeds gravity and compare it to simulated systems. Scaling our results to conditions in star-forming galaxies, we find that CRs are likely to contribute to driving outflows for a broad range of star formation environments. We quantify the momentum and energy transfer between CRs and gas, along with the associated mass outflow rates under different assumptions about the relative importance of streaming and diffusion for transport. In simulations with streaming, we observe the growth and saturation of the CR acoustic instability, but the CRs and gas remain well coupled, with CR momentum transferred efficiently to the gas even when this instability is present. Higher CR fluxes transfer more energy to the gas and drive stronger outflows. When streaming is present, most of the transferred energy takes the form of Alfvén wave heating of the gas, raising its pressure and internal energy, with a lower fractional contribution to the kinetic energy of the outflow. We also consider runs with radiative cooling, which modifies gas temperature and pressure profiles but does not seem to have a large impact on the mass outflow for super-Eddington CR fluxes.
-
more » « less
ABSTRACT Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔPc is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures, on ∼kpc scales.
