skip to main content


Title: Cooling flow solutions for the circumgalactic medium
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
Award ID(s):
1652522
NSF-PAR ID:
10184328
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
488
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2549 to 2572
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    A systematic comparison of the models of the circumgalactic medium (CGM) and their observables is crucial to understanding the predictive power of the models and constraining physical processes that affect the thermodynamics of CGM. This paper compares four analytic CGM models: precipitation, isentropic, cooling flow, and baryon pasting models for the hot, volume-filling CGM phase, all assuming hydrostatic or quasi-hydrostatic equilibrium. We show that for fiducial parameters of the CGM of a Milky Way (MW)-like galaxy ($M_{\rm vir} \sim 10^{12}~{\rm M}_{\odot }$ at $z\sim 0$), the thermodynamic profiles – entropy, density, temperature, and pressure – show most significant differences between different models at small ($r\lesssim 30$ kpc) and large scales ($r\gtrsim 100$ kpc) while converging at intermediate scales. The slope of the entropy profile, which is one of the most important differentiators between models, is $\approx 0.8$ for the precipitation and cooling flow models, while it is $\approx 0.6$ and 0 for the baryon pasting and isentropic models, respectively. We make predictions for various observational quantities for an MW mass halo for the different models, including the projected Sunyaev–Zeldovich effect, soft X-ray emission (0.5–2 keV), dispersion measure, and column densities of oxygen ions (O vi, O vii, and O viii) observable in absorption. We provide Python packages to compute the thermodynamic and observable quantities for the different CGM models.

     
    more » « less
  2. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $\sim 10^{12}\, {\rm M_{\odot }}$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of ${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $\gtrsim {\dot{M}}_{\rm crit}$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare ${\dot{M}}_{\rm crit}$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that ${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $\sim 10^{12}\, {\rm M_{\odot }}$. 
    more » « less
  3. ABSTRACT

    Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterize the gas in the vicinity of Milky Way and Andromeda analogues in the hestia (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The hestia suite comprise of a set of three high-resolution arepo-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the z = 0 simulation data sets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (H i and Si iii) are more clumpy in comparison to warmer gas tracers (O vi, O vii, and O viii). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. hestia underproduces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.

     
    more » « less
  4. ABSTRACT

    Cold, non-self-gravitating clumps occur in various astrophysical systems, ranging from the interstellar and circumgalactic medium (CGM), to active galactic nucleus outflows and solar coronal loops. Cold gas has diverse origins such as turbulent mixing or precipitation from hotter phases. We obtain the analytical solution for a steady pressure-driven 1D cooling flow around cold, local overdensities, irrespective of their origin. Our solutions describe the slow and steady radiative cooling-driven gas inflow in the saturated regime of non-linear thermal instability in clouds, sheets, and filaments. Such a cooling flow develops when the gas around small clumps undergoes radiative cooling. These small-scale, cold ‘seeds’ are embedded in a large volume-filling hot CGM maintained by feedback. We use a simple two-fluid treatment to include magnetic fields as an additional polytropic fluid. To test the limits of applicability of these analytical solutions, we compare with the gas structure found in and around small-scale cold clouds in the CGM of massive haloes in the TNG50 cosmological magnetohydrodynamic simulation from the IllustrisTNG suite. Despite qualitative resemblance of the gas structure, we find deviations from steady-state profiles generated by our model. Complex geometries and turbulence all add complexity beyond our analytical solutions. We derive an exact relation between the mass cooling rate ($\dot{\rm M}_{\rm cool}$) and the radiative cooling rate ($\dot{\rm E}_{\rm cool}$) for a steady cooling flow. A comparison with the TNG50 clouds shows that this cooling flow relation only applies in a narrow temperature range around $\rm \sim 10^{4.5}$ K where the isobaric cooling time is the shortest. In general, turbulence and mixing, instead of radiative cooling, may dominate the transition of gas between different temperature phases.

     
    more » « less
  5. ABSTRACT

    We study the intrinsic large-scale distribution and evolution of seven ionized metals in the IllustrisTNG magnetohydrodynamical cosmological simulation. We focus on the fractions of C ii, C iv, Mg ii, N v, Ne viii, O vi, and Si iv in different cosmic web structures (filaments, haloes, and voids) and gas phases (warm–hot intergalactic medium, hot, diffuse, and condensed gas) from z = 6 to z = 0. Our analysis provides a new perspective to the study of the distribution and evolution of baryons across cosmic time while offering new hints in the context of the well-known missing baryons problem. The cosmic web components are here identified using the local comoving dark matter density, which provides a simple but effective way of mapping baryons on large scales. Our results show that C ii and Mg ii are mostly located in condensed gas inside haloes in high-density and low-temperature star-forming regions ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\gtrsim 10^3$, and T ≲ 105 K). C iv and Si iv present similar evolution of their mass fractions in haloes and filaments across cosmic time. In particular, their mass budgets in haloes in condensed phase ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\gtrsim 10^3$, and T ≲ 105 K) are driven by gas cooling and star formation with a peak at z ∼ 2. Finally, our results confirm that O vi, Ne viii, and N v are good tracers of warm/hot and low-density gas at low redshift ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\lesssim 10^3$, and T ≳ 105 K), regions that are likely to contain most of the missing baryons in the local Universe.

     
    more » « less