skip to main content

Title: Cool and gusty, with a chance of rain: dynamics of multiphase CGM around massive galaxies in the Romulus simulations

Using high-resolution Romulus simulations, we explore the origin and evolution of the circumgalactic medium (CGM) in the region 0.1 ≤ R/R500 ≤ 1 around massive central galaxies in group-scale halos. We find that the CGM is multiphase and highly dynamic. Investigating the dynamics, we identify seven patterns of evolution. We show that these are robust and detected consistently across various conditions. The gas cools via two pathways: (1) filamentary cooling inflows and (2) condensations forming from rapidly cooling density perturbations. In our cosmological simulations, the perturbations are mainly seeded by orbiting substructures. The condensations can form even when the median tcool/tff of the X-ray emitting gas is above 10 or 20. Strong amplitude perturbations can provoke runaway cooling regardless of the state of the background gas. We also find perturbations whose local tcool/tff ratios drop below the threshold but which do not condense. Rather, the ratios fall to some minimum value and then bounce. These are weak perturbations that are temporarily swept up in satellite wakes and carried to larger radii. Their tcool/tff ratios decrease because tff is increasing, not because tcool is decreasing. For structures forming hierarchically, our study highlights the challenge of using a simple threshold argument to infer the CGM’s evolution. It also highlights that the median hot gas properties are suboptimal determinants of the CGM’s state and dynamics. Realistic CGM models must incorporate the impact of mergers and orbiting satellites, along with the CGM’s heating and cooling cycles.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 5677-5701
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We examine the thermodynamic state and cooling of the low-z circumgalactic medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely non-linear density perturbations sourced by the accretion of gas from the intergalactic medium (IGM) and outflows from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle-tracking analysis and find that most of the low-entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall time-scales tcool/tff in the hot component of the CGM spans a wide range of ∼1−100 at a given radius but exhibits approximately constant median values of ∼5−20 at all radii 0.1Rvir < r < Rvir. These are similar to the ≈10−20 value typically adopted as the thermal instability threshold in ‘precipitation’ models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant tcool/tff and hydrostatic equilibrium approximately reproduces the number density and entropy profiles of each simulation but only if it assumes the metallicity profile and temperature boundary condition taken directly from the simulation. We explicitly show that the tcool/tff value of a gas parcel in the hot component of the CGM does not predict its probability of subsequently accreting on to the central galaxy. This suggests that the value of tcool/tff is a poor predictor of thermal stability in gaseous haloes in which large-amplitude density perturbations are prevalent. 
    more » « less

    Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase ($\sim 10^6\, {\rm K}$) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support ($\sim 15\, {\rm kpc}$ for the Milky Way) the hot flow flattens into a disc geometry and then cools from $\sim 10^6$ to $\sim 10^4\, {\rm K}$ at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with ‘precipitation’ models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈tcool/tff radians of rotation, where tcool/tff is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio tcool/tff may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We show that if rotating hot inflows are common in Milky-Way-size disc galaxies, as predicted, then signatures of the expected hot gas rotation profile should be observable with X-ray telescopes and fast radio burst surveys.

    more » « less

    Understanding the survival, growth, and dynamics of cold gas is fundamental to galaxy formation. While there has been a plethora of work on ‘wind tunnel’ simulations that study such cold gas in winds, the infall of this gas under gravity is at least equally important, and fundamentally different since cold gas can never entrain. Instead, velocity shear increases and remains unrelenting. If these clouds are growing, they can experience a drag force due to the accretion of low-momentum gas, which dominates over ram pressure drag. This leads to subvirial terminal velocities, in line with observations. We develop simple analytic theory and predictions based on turbulent radiative mixing layers. We test these scalings in 3D hydrodynamic simulations, both for an artificial constant background and a more realistic stratified background. We find that the survival criterion for infalling gas is more stringent than in a wind, requiring that clouds grow faster than they are destroyed ($t_{\rm grow} \lt 4\, t_{\rm cc}$). This can be translated to a critical pressure, which for Milky Way-like conditions is $P \sim 3000 \, {k}_\mathrm{ B} \, {\rm K}\, {\rm cm}^{-3}$. Cold gas that forms via linear thermal instability (tcool/tff < 1) in planar geometry meets the survival threshold. In stratified environments, larger clouds need only survive infall until cooling becomes effective. We discuss applications to high-velocity clouds and filaments in galaxy clusters.

    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

    We present high-precision radial velocity observations of Gaia BH1, the nearest known black hole (BH). The system contains a solar-type G star orbiting a massive dark companion, which could be either a single BH or an inner BH + BH binary. A BH + BH binary is expected in some models where Gaia BH1 formed as a hierarchical triple, which is attractive because they avoid many of the difficulties associated with forming the system through isolated binary evolution. Our observations test the inner binary scenario. We have measured 115 precise RVs of the G star, including 40 from ESPRESSO with a precision of 3–5 m s−1, and 75 from other instruments with a typical precision of 30–100 m s−1. Our observations span 2.33 orbits of the G star and are concentrated near a periastron passage, when perturbations due to an inner binary would be largest. The RVs are well-fit by a Keplerian two-body orbit and show no convincing evidence of an inner binary. UsingREBOUNDsimulations of hierarchical triples with a range of inner periods, mass ratios, eccentricities, and orientations, we show that plausible inner binaries with periodsPinner≳ 1.5 days would have produced larger deviations from a Keplerian orbit than observed. Binaries withPinner≲ 1.5 days are consistent with the data, but these would merge within a Hubble time and would thus imply fine-tuning. We present updated parameters of Gaia BH1's orbit. The RVs yield a spectroscopic mass functionfMBH=3.9358±0.0002M—about 7000σabove the ∼2.5Mmaximum neutron star mass. Including the inclination constraint from Gaia astrometry, this implies a BH mass ofMBH= 9.27 ± 0.10M.

    more » « less