skip to main content


This content will become publicly available on October 1, 2024

Title: A Unified Model for the Coevolution of Galaxies and Their Circumgalactic Medium: The Relative Roles of Turbulence and Atomic Cooling Physics
Abstract

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 classicalMvir∼ 1011Mdwarfs, 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.

 
more » « less
Award ID(s):
1652522 1715216 2307327 1835509
NSF-PAR ID:
10475777
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
IOP
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
956
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (ηM,ηE, andηZ, here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 1010to 1012M, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (ηM∼ 0.1–10) but high-energy-loaded (ηE∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringηE∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.

     
    more » « less
  2. ABSTRACT

    We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma β ≡ Pg/PB to fall, the CGM becomes globally unstable. At the cool disc-to-hot−halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower β haloes, CR heating drives a hot, single-phase diffuse wind with velocities v ∝ (theat/tff)−1, which exceeds the escape velocity when theat/tff ≲ 0.4. In higher β haloes, where the Alfven Mach number is higher, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to ‘trapping’ of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.

     
    more » « less
  3. Abstract

    Massive elliptical galaxies harbor large amounts of hot gas (T≳ 106K) in their interstellar medium (ISM) but are typically quiescent in star formation. The jets of active galactic nuclei (AGNs) and Type Ia supernovae (SNe Ia) inject energy into the ISM, which offsets its radiative losses and keeps it hot. SNe Ia deposit their energy locally within the galaxy compared to the larger few ×10 kiloparsec-scale AGN jets. In this study, we perform high-resolution (5123) hydrodynamic simulations of a local (1 kpc3) density-stratified patch of the ISM of massive galaxies. We include radiative cooling and shell-averaged volume heating, as well as randomly exploding SN Ia. We study the effect of different fractions of supernova (SN) heating (with respect to the net cooling rate), different initial ISM density/entropy (which controls the growth timettiof the thermal instability), and different degrees of stratification (which affect the freefall timetff). We find that SNe Ia drive predominantly compressive turbulence in the ISM with a velocity dispersion ofσvup to 40 km s−1and logarithmic density dispersion ofσs∼ 0.2–0.4. These fluctuations trigger multiphase condensation in regions of the ISM, wheremin(tti)/tff0.6exp(6σs), in agreement with theoretical expectations that large density fluctuations efficiently trigger multiphase gas formation. Since the SN Ia rate is not self-adjusting, when the net cooling drops below the net heating rate, SNe Ia drive a hot wind which sweeps out most of the mass in our local model. Global simulations are required to assess the ultimate fate of this gas.

     
    more » « less
  4. Abstract

    Absorption-line measurements of the circumgalactic medium (CGM) display a highly nonuniform distribution of lower ionization state species accompanied by more widespread higher ionization state material. This suggests that the CGM is a dynamic, multiphase medium, such as arises in the presence of turbulence. To better understand this evolution, we perform hydrodynamic and magnetohydrodynamic (MHD) simulations of the CGM surrounding Milky Way–like galaxies. In both cases, the CGM is initially in hydrostatic balance in a 1012Mdark matter gravitational potential, and the simulations include rotation in the inner halo and turbulence that decreases radially. They also track ionizations, recombinations, and species-by-species radiative cooling in the presence of the redshift-zero UV background, employing the MAIHEM nonequilibrium chemistry package. We find that after 9 Gyr of evolution, the presence of a magnetic field leads to an overall hotter CGM, with cool gas in the center where magnetic pressure dominates. While the non-MHD run produces more cold clouds overall, we find similar Siiv/Oviand Nv/Oviratios between the MHD and non-MHD runs, which are both very different from their equilibrium values. The non-MHD halo develops cool, low angular momentum filaments above the central disk, in comparison to the MHD run that has more efficient angular momentum transport, especially for the cold gas, which forms a more ordered and extended disk late into its evolution.

     
    more » « less
  5. Abstract

    This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshiftz≈ 0.4–1. Through spectrally resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small aslcl∼ 1 pc and density lower thannH= 10−3cm−3are identified in galaxy halos. In addition, comparing the line widths between different elements for kinematically matched components provides robust empirical constraints on the thermal temperatureTand the nonthermal motionsbNT, independent of the ionization models. On average,bNTis found to increase withlclfollowingbNTlcl0.3over three decades in spatial scale fromlcl≈ 1 pc tolcl≈ 1 kpc. Attributing the observedbNTto turbulent motions internal to the clumps, the best-fitbNTlclrelation shows that the turbulence is consistent with Kolmogorov at <1 kpc with a roughly constant energy transfer rate per unit mass ofϵ≈ 0.003 cm2s−3and a dissipation timescale of ≲100 Myr. No significant difference is found between massive quiescent and star-forming halos in the sample on scales less than 1 kpc. While the inferredϵis comparable to what is found in Civabsorbers at high redshift, it is considerably smaller than observed in star-forming gas or in extended line-emitting nebulae around distant quasars. A brief discussion of possible sources to drive the observed turbulence in the cool CGM is presented.

     
    more » « less