Search for: All records

Abstract This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on howimbalancesbetween heating and cooling cause either expansion or contraction of the CGM. It does this by trackingallof the mass and energy associated with a halo’s baryons, including their gravitational potential energy, even if feedback has pushed some of those baryons beyond the halo’s virial radius. We show how a star-forming galaxy’s equilibrium state can be algebraically derived within the context of this framework, and we analyze how the equilibrium star formation rate depends on supernova feedback. We consider the consequences of varying the mass loading parameter ${\eta }_M\equiv {\dot{M}}_{\mathrm{wind}}/{\dot{M}}_{*}$relating a galaxy’s gas mass outflow rate ( ${\dot{M}}_{\mathrm{wind}}$) to its star formation rate ( ${\dot{M}}_{*}$) and obtain results that challenge common assumptions. In particular, we find that equilibrium star formation rates in low-mass galaxies are generally insensitive to mass loading, and when mass loading does matter, increasing it actually results inmorestar formation because more supernova energy is needed to resist atmospheric contraction. 

Abstract The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This follow-up paper interprets that model and discusses its implications. The model itself is an accounting system that tracks all of the mass and energy associated with a halo’s circumgalactic baryons—the central galaxy’s atmosphere. Algebraic solutions for the equilibrium states of that model reveal that star formation in low-mass halos self-regulates primarily by expanding the atmospheres of those halos, ultimately resulting in stellar masses that are insensitive to the mass-loading properties of galactic winds. What matters most is the proportion of supernova energy that couples with circumgalactic gas. However, supernova feedback alone fails to expand galactic atmospheres in higher-mass halos. According to the minimalist regulator model, an atmospheric contraction crisis ensues, which may be what triggers strong black hole feedback. The model also predicts that circumgalactic medium properties emerging from cosmological simulations should depend largely on the specific energy of the outflows they produce, and we interpret the qualitative properties of several numerical simulations in light of that prediction. 

Abstract The baryonic physics shaping galaxy formation and evolution are complex, spanning a vast range of scales and making them challenging to model. Cosmological simulations rely on subgrid models that produce significantly different predictions. Understanding how models of stellar and active galactic nucleus (AGN) feedback affect baryon behavior across different halo masses and redshifts is essential. Using the SIMBA and IllustrisTNG suites from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, we explore the effect of parameters governing the subgrid implementation of stellar and AGN feedback. We find that while IllustrisTNG shows higher cumulative feedback energy across all halos, SIMBA demonstrates a greater spread of baryons, quantified by the closure radius and circumgalactic medium (CGM) gas fraction. This suggests that feedback in SIMBA couples more effectively to baryons and drives them more efficiently within the host halo. There is evidence that the different feedback modes are highly interrelated in these subgrid models. The parameters controlling the stellar feedback efficiency significantly impact AGN feedback, as seen in the suppression of black hole mass growth and delayed activation of AGN feedback to higher-mass halos with increasing stellar feedback efficiency in both simulations. Additionally, the AGN feedback efficiency parameters affect the CGM gas fraction at low halo masses in SIMBA, hinting at complex, nonlinear interactions between the AGN and supernova feedback modes. Overall, we demonstrate that stellar and AGN feedback are intimately interwoven, especially at low redshift, due to subgrid implementation, resulting in halo property effects that might initially seem counterintuitive. 

Abstract The circumgalactic medium (CGM) around massive galaxies plays a crucial role in regulating star formation and feedback. Using the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) suite, we develop emulators for the X-ray surface brightness profile and the X-ray luminosity–stellar mass scaling relation, to investigate how stellar and active galactic nucleus (AGN) feedback shape the X-ray properties of the hot CGM. Our analysis shows that at CGM scales (10¹²≲M_halo/M_⊙≲ 10¹³, 10 ≲rkpc⁻¹≲ 400), stellar feedback more significantly impacts the X-ray properties than AGN feedback within the parameters studied. Comparing the emulators to recent eROSITA All Sky Survey (eRASS) observations, it is found that stronger feedback than is currently implemented in the IllustrisTNG, SIMBA, and Astrid simulations is required to match the observed CGM properties. However, adopting these enhanced feedback parameters causes deviations in the stellar mass–halo mass relations from observational constraints below the group-mass scale. This tension suggests possible unaccounted-for systematics in X-ray CGM observations or inadequacies in the feedback models of cosmological simulations. 

ABSTRACT The circum-galactic medium (CGM) can feasibly be mapped by multiwavelength surveys covering broad swaths of the sky. With multiple large data sets becoming available in the near future, we develop a likelihood-free Deep Learning technique using convolutional neural networks (CNNs) to infer broad-scale physical properties of a galaxy’s CGM and its halo mass for the first time. Using CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations) data, including IllustrisTNG, SIMBA, and Astrid models, we train CNNs on Soft X-ray and 21-cm (H i) radio two-dimensional maps to trace hot and cool gas, respectively, around galaxies, groups, and clusters. Our CNNs offer the unique ability to train and test on ‘multifield’ data sets comprised of both H i and X-ray maps, providing complementary information about physical CGM properties and improved inferences. Applying eRASS:4 survey limits shows that X-ray is not powerful enough to infer individual haloes with masses log (Mhalo/M⊙) < 12.5. The multifield improves the inference for all halo masses. Generally, the CNN trained and tested on Astrid (SIMBA) can most (least) accurately infer CGM properties. Cross-simulation analysis – training on one galaxy formation model and testing on another – highlights the challenges of developing CNNs trained on a single model to marginalize over astrophysical uncertainties and perform robust inferences on real data. The next crucial step in improving the resulting inferences on the physical properties of CGM depends on our ability to interpret these deep-learning models. 

Abstract Quenching of star formation in the central galaxies of cosmological halos is thought to result from energy released as gas accretes onto a supermassive black hole. The same energy source also appears to lower the central density and raise the cooling time of baryonic atmospheres in massive halos, thereby limiting both star formation and black hole growth, by lifting the baryons in those halos to greater altitudes. One predicted signature of that feedback mechanism is a nearly linear relationship between the central black hole’s mass (M_BH) and the original binding energy of the halo’s baryons. We present the increasingly strong observational evidence supporting a such a relationship, showing that it extends up to halos of massM_halo∼ 10¹⁴M_⊙. We then compare current observational constraints on theM_BH–M_halorelation with numerical simulations, finding that black hole masses in IllustrisTNG appear to exceed those constraints atM_halo< 10¹³M_⊙and that black hole masses in EAGLE fall short of observations atM_halo∼ 10¹⁴M_⊙. A closer look at IllustrisTNG shows that quenching of star formation and suppression of black hole growth do indeed coincide with black hole energy input that lifts the halo’s baryons. However, IllustrisTNG does not reproduce the observedM_BH–M_halorelation because its black holes gain mass primarily through accretion that does not contribute to baryon lifting. We suggest adjustments to some of the parameters in the IllustrisTNG feedback algorithm that may allow the resulting black hole masses to reflect the inherent links between black hole growth, baryon lifting, and star formation among the massive galaxies in those simulations. 

Abstract We present an analysis of Hubble Space Telescope COS/G160M observations of CIVin the inner circumgalactic medium (CGM) of a novel sample of eightz∼ 0,L≈L^⋆galaxies, paired with UV-bright QSOs at impact parameters (R_proj) between 25 and 130 kpc. The galaxies in this stellar-mass-controlled sample (log₁₀M_⋆/M_⊙∼ 10.2–10.9M_⊙) host supermassive black holes (SMBHs) with dynamically measured masses spanning log₁₀M_BH/M_⊙∼ 6.8–8.4; this allows us to compare our results with models of galaxy formation where the integrated feedback history from the SMBH alters the CGM over long timescales. We find that the CIVcolumn density measurements (N_{C IV}; average log₁₀N_{C IV,CH}= 13.94 ± 0.09 cm⁻²) are largely consistent with existing measurements from other surveys ofN_{C IV}in the CGM (average log₁₀N_{C IV,Lit}= 13.90 ± 0.08 cm⁻²), but do not show obvious variation as a function of the SMBH mass. By contrast, specific star formation rate (sSFR) is highly correlated with the ionized content of the CGM. We find a large spread in sSFR for galaxies with log₁₀M_BH/M_⊙> 7.0, where the CGM CIVcontent shows a clear dependence on galaxy sSFR but notM_BH. Our results do not indicate an obvious causal link between CGM CIVand the mass of the galaxy’s SMBH; however, through comparisons to the EAGLE, Romulus25, and IllustrisTNG simulations, we find that our sample is likely too small to constrain such causality. 

Abstract Dwarf galaxies are found to have lost most of their metals via feedback processes; however, there still lacks consistent assessment on the retention rate of metals in their circumgalactic medium (CGM). Here we investigate the metal content in the CGM of 45 isolated dwarf galaxies withM_*= 10^6.5–9.5M_⊙(M_200m= 10^10.0–11.5M_⊙) using the Hubble Space Telescope/Cosmic Origins Spectrograph. While Hi(Lyα) is ubiquitously detected (89%) within the CGM, we find low detection rates (≈5%–22%) in Cii, Civ, Siii, Siiii, and Siiv, largely consistent with literature values. Assuming these ions form in the cool (T≈ 10⁴K) CGM with photoionization equilibrium, the observed Hiand metal column density profiles can be best explained by an empirical model with low gas density and high volume filling factor. For a typical galaxy withM_200m= 10^10.9M_⊙(median of the sample), our model predicts a cool gas mass ofM_CGM,cool∼ 10^8.4M_⊙, corresponding to ∼2% of the galaxy’s baryonic budget. Assuming a metallicity of 0.3 Z_⊙, we estimate that the dwarf galaxy’s cool CGM likely harbors ∼10% of the metals ever produced, with the rest either in more ionized states in the CGM or transported to the intergalactic medium. We further examine the EAGLE simulation and show that Hiand low ions may arise from a dense cool medium, while Civarises from a diffuse warmer medium. Our work provides the community with a uniform data set on dwarf galaxies’ CGM that combines our recent observations, additional archival data and literature compilation, which can be used to test various theoretical models of dwarf galaxies. 

ABSTRACT Known as the ‘Missing Baryon Problem’, about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm–hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an O vii absorption line towards the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the O vii absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGNs) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated O vii column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the O vii distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the O vii column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated O vii column densities cannot be explained by the wide range of feedback models implemented in CAMELS. 

Abstract The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 H i -selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields at z ≲ 1. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) at z < 1. Here we present Keck/KCWI and Very Large Telescope/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with Hubble Space Telescope/Advanced Camera for Surveys imaging to identify and characterize the absorber-associated galaxies at 0.16 ≲ z ≲ 0.84. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with <10% solar metallicities, we find no associated galaxies with log M ⋆ ≳ 9.0 within ρ / R vir and ∣Δ v ∣/ v esc ≤ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or H i column densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] ≤ −1.4) systems have a probability of 0.39 − 0.15 + 0.16 for having a host galaxy with log M ⋆ ≥ 9.0 within ρ / R vir ≤ 1.5, while the higher metallicity absorbers have a probability of 0.78 − 0.13 + 0.10 . This implies metal-enriched pLLSs/LLSs at z < 1 are typically associated with the CGM of galaxies with log M ⋆ > 9.0 , whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] ≲ −2.1 at z < 1, which is lower than previously estimated.