Showing 1–50 of 109 results for author: Anglés-Alcázar, D

Deciphering the imprint of AGN feedback in Seyfert galaxies: Nuclear-scale molecular gas deficits

Authors: S. García-Burillo, E. K. S. Hicks, A. Alonso-Herrero, M. Pereira-Santaella, A. Usero, M. Querejeta, O. González-Martin, D. Delaney, C. Ramos Almeida, F. Combes, D. Anglés-Alcázar, A. Audibert, E. Bellocchi, R. I. Davies, T. A. Davis, J. S. Elford, I. García-Bernete, S. Hönig, A. Labiano, M. T. Leist, N. A. Levenson, E. López-Rodríguez, J. Mercedes-Feliz, C. Packham, C. Ricci , et al. (4 additional authors not shown)

Abstract: We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV ran… ▽ More We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV range ($L_{\rm X}$). We also study the concentration of the hot molecular gas using NIR data obtained for the H2 1-0S(1) line. We find a turnover in the distribution of the cold molecular gas concentration as a function of $L_{\rm X}$ with a breakpoint which divides the sample into two branches: the AGN build-up branch ($L_{\rm X}\leq10^{41.5\pm0.3}$erg/s) and the AGN feedback branch ($L_{\rm X}\geq10^{41.5\pm0.3}$erg/s) . Lower luminosity AGN and non-AGN of the AGN build-up branch show high cold molecular gas concentrations and centrally peaked radial profiles on nuclear ($r\leq50$~pc) scales. Higher luminosity AGN of the AGN feedback branch, show a sharp decrease in the concentration of molecular gas and flat or inverted radial profiles. The cold molecular gas concentration index ($CCI$), defined as the ratio of surface densities at $r\leq50$~pc and $r\leq200$~pc , namely $CCI \equiv$~log$_{\rm 10}(Σ^{\rm gas}_{\rm 50}/Σ^{\rm gas}_{\rm 200}$), spans a factor ~4-5 between the galaxies lying at the high end of the AGN build-up branch and the galaxies of the AGN feedback branch. The concentration and radial distributions of the hot molecular gas in our sample follow less extreme trends as a function of the X-ray luminosity. These observations confirm, on a three times larger sample, previous evidence found by the GATOS survey that the imprint of AGN feedback on the CND-scale distribution of molecular gas is more extreme in higher luminosity Seyfert galaxies of the local universe. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted by Astronomy & Astrophysics (A&A) (14/06/2024), 26 pages, 16 figures

Host-galaxy stars can dominate the ionizing radiation field of the circumgalactic medium in galaxies at Cosmic Noon

Authors: Francisco Holguin, Christopher C. Hayward, Xiangcheng Ma, Daniel Anglés-Alcázar, Rachel K. Cochrane

Abstract: Elucidating the processes that shape the circumgalactic medium (CGM) is crucial for understanding galaxy evolution. Absorption and emission diagnostics can be interpreted using photoionization calculations to obtain information about the phase and ionization structure of the CGM. For simplicity, typically only the metagalactic background is considered in photoionization calculations, and local sou… ▽ More Elucidating the processes that shape the circumgalactic medium (CGM) is crucial for understanding galaxy evolution. Absorption and emission diagnostics can be interpreted using photoionization calculations to obtain information about the phase and ionization structure of the CGM. For simplicity, typically only the metagalactic background is considered in photoionization calculations, and local sources are ignored. To test this simplification, we perform Monte Carlo radiation transfer on 12 cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project with halo masses $10^{10.5}-10^{13} \mathrm{M}_{\odot}$ in the redshift range $z = 0-3.5$ to determine the spatial extent over which local sources appreciably contribute to the ionizing radiation field in the CGM. We find that on average, the contribution of stars within the galaxy is small beyond one-tenth of the virial radius, $R_{\mathrm{vir}}$, for $z < 1$. For $1<z<2$ and $M_{\mathrm{vir}} \sim 10^{11.5}$, the radius at which the contribution to the ionizing radiation field from stars within the galaxy and that from the UV background are equal is roughly 0.2 $R_{\mathrm{vir}}$. For $M_{\mathrm{vir}} > 10^{12} \mathrm{M}_{\odot}$ at $z \sim 1.5-2.5$ and for all $M_{\mathrm{vir}}$ considered at $z>3$ , this transition radius can sometimes exceed 0.5 $R_{\mathrm{vir}}$. We also compute the escape fraction at $R_{\mathrm{vir}}$, finding typical values of less than $0.1$, except in higher-mass halos ($M_{\mathrm{halo}} \gtrsim 10^{12} \mathrm{M}_{\odot}$), which have consistently high values of $\sim 0.5-0.6$. Our results indicate that at low redshift, it is reasonable to ignore the ionizing radiation from host-galaxy stars outside of 0.2 $R_{\mathrm{vir}}$, while at Cosmic Noon, local stellar ionizing radiation likely extends further into the CGM and thus should be included in photoionization calculations. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Submitted to OJA

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine learning and Simulations

Authors: Jonah C. Rose, Paul Torrey, Francisco Villaescusa-Navarro, Mariangela Lisanti, Tri Nguyen, Sandip Roy, Kassidy E. Kollmann, Mark Vogelsberger, Francis-Yan Cyr-Racine, Mikhail V. Medvedev, Shy Genel, Daniel Anglés-Alcázar, Nitya Kallivayalil, Bonny Y. Wang, Belén Costanza, Stephanie O'Neil, Cian Roche, Soumyodipta Karmakar, Alex M. Garcia, Ryan Low, Shurui Lin, Olivia Mostow, Akaxia Cruz, Andrea Caputo, Arya Farahi , et al. (5 additional authors not shown)

Abstract: We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- f… ▽ More We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of Warm Dark Matter, each comprised of 1024 simulations generated using the Arepo code. One suite consists of uniform-box simulations covering a $(25~h^{-1}~{\rm M}_\odot)^3$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the Warm Dark Matter particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and Convolutional Neural Networks, to demonstrate how such simulation suites can be used to disentangle the effects of dark matter and baryonic physics on galactic properties. The DREAMS project can be extended further to include different dark matter models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of dark matter on these scales. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 28 pages, 8 figures, DREAMS website: https://www.dreams-project.org

Zooming by in the CARPoolGP lane: new CAMELS-TNG simulations of zoomed-in massive halos

Authors: Max E. Lee, Shy Genel, Benjamin D. Wandelt, Benjamin Zhang, Ana Maria Delgado, Shivam Pandey, Erwin T. Lau, Christopher Carr, Harrison Cook, Daisuke Nagai, Daniel Angles-Alcazar, Francisco Villaescusa-Navarro, Greg L. Bryan

Abstract: Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we dev… ▽ More Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method, CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) to include a set of 768 zoom-in simulations of halos in the mass range of $10^{13} - 10^{14.5} M_\odot\,h^{-1}$ that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the Compton $Y-M$, black hole mass-halo mass, and metallicity-mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of the $Y-M$ relation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev-Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: The manuscript was submitted to arxiv after receiving and responding to comments from the first referee report

The baryon cycle in modern cosmological hydrodynamical simulations

Authors: Ruby J. Wright, Rachel S. Somerville, Claudia del P. Lagos, Matthieu Schaller, Romeel Davé, Daniel Anglés-Alcázar, Shy Genel

Abstract: In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predict… ▽ More In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, M_{\odot}$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2-3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the ISM, recycle within the CGM (within $R_{\rm 200c}$). We find that AGN-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx10^{13.5}\, M_{\odot}$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2-3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform sub-grid feedback models in the next generation of simulations. △ Less

Submitted 9 July, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures. Accepted for publication in MNRAS 03/07/24

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Boon Kiat Oh, Christopher C. Hayward, Rachel K. Cochrane, Alexander J. Richings, Claude-André Faucher-Giguère, Sarah Wellons, Bryan A. Terrazas, Jorge Moreno, Kung Yi Su, Philip F. Hopkins

Abstract: We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerfu… ▽ More We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $Σ_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $Σ_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters. △ Less

Submitted 17 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS; doi:10.1093/mnras/stae1021

Beware the recent past: a bias in spectral energy distribution modelling due to bursty star formation

Authors: P. Haskell, S. Das, D. J. B. Smith, R. K. Cochrane, C. C. Hayward, D. Anglés-Alcázar

Abstract: We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We id… ▽ More We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We identify a previously-unknown systematic error in the Magphys results due to bursty star formation: the derived SFRs can differ from the truth by as much as 1 dex, at large statistical significance ($>5σ$), depending on the details of their recent SFH. SFRs inferred using Prospector with non-parametric SFHs do not exhibit this trend. We show that using parametric SFHs (pSFHs) causes SFR uncertainties to be underestimated by a factor of up to $5\times$. Although this undoubtedly contributes to the significance of the systematic, it cannot explain the largest biases in the SFRs of the starbursting galaxies, which could be caused by details of the stochastic prior sampling or the burst implementation in the Magphys libraries. We advise against using pSFHs and urge careful consideration of starbursts when SED modelling galaxies where the SFR may have changed significantly over the last ~100 Myr, such as recently quenched galaxies, or those experiencing a burst. This concern is especially relevant, e.g. when fitting JWST observations of very high-redshift galaxies. △ Less

Submitted 8 March, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS Letters

Effects of multi-channel AGN feedback in FIRE cosmological simulations of massive galaxies

Authors: Lindsey Byrne, Claude-André Faucher-Giguère, Sarah Wellons, Philip F. Hopkins, Daniel Anglés-Alcázar, Imran Sultan, Nastasha Wijers, Jorge Moreno, Sam Ponnada

Abstract: Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. AGN feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of FIRE-3 cosmological zoom-in simulations w… ▽ More Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. AGN feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of FIRE-3 cosmological zoom-in simulations with AGN feedback evolved to $z\sim0$, examining the impact of AGN feedback on a set of galaxies with halos in the mass range $10^{12}-10^{13} M_{\odot}$. These simulations combine detailed stellar and ISM physics with multi-channel AGN feedback including radiative feedback, mechanical outflows, and in some simulations, cosmic rays (CRs). We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass-halo mass relation and the $M_{\rm BH}$-$σ$ relation; in the stronger model with CRs, they also match the size-mass relation and the Faber-Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber-Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: while the stronger model produces the most realistic massive galaxies, it tends to over-quench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed. △ Less

Submitted 7 May, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 18 pages, 8 figures, submitted to ApJ

Disappearing galaxies: the orientation dependence of JWST-bright, HST-dark, star-forming galaxy selection

Authors: R. K. Cochrane, D. Anglés-Alcázar, F. Cullen, C. C. Hayward

Abstract: Galaxies that are invisible in deep optical-NIR imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at $z\gtrsim4$. The depths now achievable with JWST at the longest wavelengths probed by HST, coupled with the tra… ▽ More Galaxies that are invisible in deep optical-NIR imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at $z\gtrsim4$. The depths now achievable with JWST at the longest wavelengths probed by HST, coupled with the transformative resolution at longer wavelengths, are already enabling detailed, spatially-resolved characterisation of sources that were invisible to HST, often known as `HST-dark' galaxies. However, until now, there has been little theoretical work to compare against. We present the first simulation-based study of this population, using highly-resolved galaxies from the Feedback in Realistic Environments (FIRE) project, with multi-wavelength images along several lines of sight forward-modelled using radiative transfer. We naturally recover a population of modelled sources that meet commonly-used selection criteria ($H_{\rm{AB}}>27\,\rm{mag}$ and $H_{\rm{AB}}-\rm{F444W}>2.3$). These simulated HST-dark galaxies lie at high redshifts ($z=4-7$), have high levels of dust attenuation ($A_{V}=2-4$), and display compact recent star formation ($R_{1/2,\,\rm{4.4\,μ\rm{m}}}\lesssim1\,\rm{kpc}$). Orientation is very important: for all but one of the 17 simulated galaxy snapshots with HST-dark sightlines, there exist other sightlines that do not meet the criteria. This result has important implications for comparisons between observations and models that do not resolve the detailed star-dust geometry, such as semi-analytic models or coarsely-resolved hydrodynamical simulations. Critically, we demonstrate that HST-dark sources are not an unexpected or exotic population, but a subset of high-redshift, highly-dust-attenuated sources viewed along certain lines of sight. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures. Accepted for publication in ApJ

FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Jonathan Squire, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Yanlong Shi, Michael Y. Grudic, Sarah Wellons, Claude-Andre Faucher-Giguere, Daniel Angles-Alcazar, Norman Murray, Eliot Quataert

Abstract: In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc s… ▽ More In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $β\sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk. △ Less

Submitted 18 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

FORGE'd in FIRE: Resolving the End of Star Formation and Structure of AGN Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Michael Y. Grudic, Kung-Yi Su, Sarah Wellons, Daniel Angles-Alcazar, Ulrich P. Steinwandel, David Guszejnov, Norman Murray, Claude-Andre Faucher-Giguere, Eliot Quataert, Dusan Keres

Abstract: It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydro… ▽ More It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to <100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{\rm M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{\rm schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed. △ Less

Submitted 12 March, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 38 pages, 18 figures. Replaced with accepted version

Cosmology with multiple galaxies

Authors: Chaitanya Chawak, Francisco Villaescusa-Navarro, Nicolas Echeverri Rojas, Yueying Ni, ChangHoon Hahn, Daniel Angles-Alcazar

Abstract: Recent works have discovered a relatively tight correlation between $Ω_{\rm m}$ and properties of individual simulated galaxies. Because of this, it has been shown that constraints on $Ω_{\rm m}$ can be placed using the properties of individual galaxies while accounting for uncertainties on astrophysical processes such as feedback from supernova and active galactic nuclei. In this work, we quantif… ▽ More Recent works have discovered a relatively tight correlation between $Ω_{\rm m}$ and properties of individual simulated galaxies. Because of this, it has been shown that constraints on $Ω_{\rm m}$ can be placed using the properties of individual galaxies while accounting for uncertainties on astrophysical processes such as feedback from supernova and active galactic nuclei. In this work, we quantify whether using the properties of multiple galaxies simultaneously can tighten those constraints. For this, we train neural networks to perform likelihood-free inference on the value of two cosmological parameters ($Ω_{\rm m}$ and $σ_8$) and four astrophysical parameters using the properties of several galaxies from thousands of hydrodynamic simulations of the CAMELS project. We find that using properties of more than one galaxy increases the precision of the $Ω_{\rm m}$ inference. Furthermore, using multiple galaxies enables the inference of other parameters that were poorly constrained with one single galaxy. We show that the same subset of galaxy properties are responsible for the constraints on $Ω_{\rm m}$ from one and multiple galaxies. Finally, we quantify the robustness of the model and find that without identifying the model range of validity, the model does not perform well when tested on galaxies from other galaxy formation models. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures

Energy balance SED modelling can be effective at high redshifts regardless of UV-FIR offsets

Authors: P. Haskell, D. J. B. Smith, R. K. Cochrane, C. C. Hayward, D. Anglés-Alcázar

Abstract: Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such… ▽ More Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies' true properties. In this work, we test how well the SED fitting code Magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass and luminosity by fitting 6,706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4-500$μ$m, and spanning $1<z<8$) with fits based on either the starlight ($λ_\mathrm{eff} \le 2.2\,μ$m) or dust ($\ge 100\,μ$m) alone, we highlight the power of considering the full range of multi-wavelength data alongside an energy balance criterion. Overall, we obtain acceptable fits for 83 per cent of the synthetic SEDs, though the success rate falls rapidly beyond $z \approx 4$, in part due to the sparser sampling of the priors at earlier times since SFHs must be physically plausible (i.e. shorter than the age of the Universe). We use the ground truth from the simulations to show that when the quality of fit is acceptable, the fidelity of Magphys estimates is independent of the degree of UV\FIR offset, with performance very similar to that previously reported for local galaxies. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 18 pages, 9 figures, published in MNRAS

Journal ref: MNRAS, vol 525, p 1535 (2023)

An Observationally Driven Multifield Approach for Probing the Circum-Galactic Medium with Convolutional Neural Networks

Authors: Naomi Gluck, Benjamin D. Oppenheimer, Daisuke Nagai, Francisco Villaescusa-Navarro, Daniel Anglés-Alcázar

Abstract: The circum-galactic medium (CGM) can feasibly be mapped by multiwavelength surveys covering broad swaths of the sky. With multiple large datasets 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 (Cosmolog… ▽ More The circum-galactic medium (CGM) can feasibly be mapped by multiwavelength surveys covering broad swaths of the sky. With multiple large datasets 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 (HI) radio 2D 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'' datasets comprised of both HI 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 halos with masses $\log(M_{\rm{halo}}/M_{\odot}) < 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 physical CGM properties hinges on our ability to interpret these deep-learning models. △ Less

Submitted 16 January, 2024; v1 submitted 14 September, 2023; originally announced September 2023.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 4, February 2024, Pages 10038-10058

Cosmological baryon spread and impact on matter clustering in CAMELS

Authors: Matthew Gebhardt, Daniel Anglés-Alcázar, Josh Borrow, Shy Genel, Francisco Villaescusa-Navarro, Yueying Ni, Christopher Lovell, Daisuke Nagai, Romeel Davé, Federico Marinacci, Mark Vogelsberger, Lars Hernquist

Abstract: We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales com… ▽ More We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighboring dark matter owing to feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial \textsc{SIMBA} model spreading $\sim$40\% of baryons $>$1\,Mpc away compared to $\sim$10\% for the IllustrisTNG and \textsc{ASTRID} models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired $N$-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number ($k$) and baryonic spread up to $k \sim 10\,h$\,Mpc$^{-1}$ while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 17 pages, 15 figures

A Hierarchy of Normalizing Flows for Modelling the Galaxy-Halo Relationship

Authors: Christopher C. Lovell, Sultan Hassan, Daniel Anglés-Alcázar, Greg Bryan, Giulio Fabbian, Shy Genel, ChangHoon Hahn, Kartheik Iyer, James Kwon, Natalí de Santi, Francisco Villaescusa-Navarro

Abstract: Using a large sample of galaxies taken from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, a suite of hydrodynamic simulations varying both cosmological and astrophysical parameters, we train a normalizing flow (NF) to map the probability of various galaxy and halo properties conditioned on astrophysical and cosmological parameters. By leveraging the learnt cond… ▽ More Using a large sample of galaxies taken from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, a suite of hydrodynamic simulations varying both cosmological and astrophysical parameters, we train a normalizing flow (NF) to map the probability of various galaxy and halo properties conditioned on astrophysical and cosmological parameters. By leveraging the learnt conditional relationships we can explore a wide range of interesting questions, whilst enabling simple marginalisation over nuisance parameters. We demonstrate how the model can be used as a generative model for arbitrary values of our conditional parameters; we generate halo masses and matched galaxy properties, and produce realisations of the halo mass function as well as a number of galaxy scaling relations and distribution functions. The model represents a unique and flexible approach to modelling the galaxy-halo relationship. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 8 pages, 2 figures, accepted for ICML 2023 Workshop on Machine Learning for Astrophysics

An Exploration of AGN and Stellar Feedback Effects in the Intergalactic Medium via the Low Redshift Lyman-$α$ Forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Sultan Hassan, Rachel S. Somerville, Romeel Davé, Federico Marinacci, Lars Hernquist, Mark Vogelsberger

Abstract: We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both… ▽ More We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both AGN and stellar feedback in Simba play a role in setting the Ly$α$ forest column density distribution function (CDD) and the Doppler width ($b$-value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM causing changes in the shape and normalization of the CDD and a broadening of the $b$-value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Ly$α$ forest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low-$z$ Ly$α$ forest can be sensitive to changes in the ultraviolet background (UVB), stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM. △ Less

Submitted 1 November, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 27 pages, 12 figures, 2 tables, accepted to ApJ

Inferring Warm Dark Matter Masses with Deep Learning

Authors: Jonah C. Rose, Paul Torrey, Francisco Villaescusa-Navarro, Mark Vogelsberger, Stephanie O'Neil, Mikhail V. Medvedev, Ryan Low, Rakshak Adhikari, Daniel Angles-Alcazar

Abstract: We present a new suite of over 1,500 cosmological N-body simulations with varied Warm Dark Matter (WDM) models ranging from 2.5 to 30 keV. We use these simulations to train Convolutional Neural Networks (CNNs) to infer WDM particle masses from images of DM field data. Our fiducial setup can make accurate predictions of the WDM particle mass up to 7.5 keV at a 95% confidence level from small maps t… ▽ More We present a new suite of over 1,500 cosmological N-body simulations with varied Warm Dark Matter (WDM) models ranging from 2.5 to 30 keV. We use these simulations to train Convolutional Neural Networks (CNNs) to infer WDM particle masses from images of DM field data. Our fiducial setup can make accurate predictions of the WDM particle mass up to 7.5 keV at a 95% confidence level from small maps that cover an area of (25 h$^{-1}$ Mpc)$^2$. We vary the image resolution, simulation resolution, redshift, and cosmology of our fiducial setup to better understand how our model is making predictions. Using these variations, we find that our models are most dependent on simulation resolution, minimally dependent on image resolution, not systematically dependent on redshift, and robust to varied cosmologies. We also find that an important feature to distinguish between WDM models is present with a linear size between 100 and 200 h$^{-1}$ kpc. We compare our fiducial model to one trained on the power spectrum alone and find that our field-level model can make 2x more precise predictions and can make accurate predictions to 2x as massive WDM particle masses when used on the same data. Overall, we find that the field-level data can be used to accurately differentiate between WDM models and contain more information than is captured by the power spectrum. This technique can be extended to more complex DM models and opens up new opportunities to explore alternative DM models in a cosmological environment. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 16 pages, 12 figures

Massive black holes in galactic nuclei: Theory and Simulations

Authors: Tiziana Di Matteo, Daniel Angles-Alcazar, Francesco Shankar

Abstract: Massive black holes are fundamental constituents of our cosmos, from the Big Bang to today. Understanding their formation from cosmic dawn, their growth, and the emergence of the first, rare quasars in the early Universe remains one of our greatest theoretical and observational challenges. Hydrodynamic cosmological simulations self-consistently combine the processes of structure formation at cosmo… ▽ More Massive black holes are fundamental constituents of our cosmos, from the Big Bang to today. Understanding their formation from cosmic dawn, their growth, and the emergence of the first, rare quasars in the early Universe remains one of our greatest theoretical and observational challenges. Hydrodynamic cosmological simulations self-consistently combine the processes of structure formation at cosmological scales with the physics of smaller, galaxy scales. They capture our most realistic understanding of massive black holes and their connection to galaxy formation and have become the primary avenue for theoretical research in this field. The space-based gravitational wave interferometer, LISA, will open up new investigations into the dynamical processes involving massive black holes. Multi-messenger astrophysics brings new exciting prospects for tracing the origin, growth and merger history of massive black holes across cosmic ages. △ Less

Submitted 23 April, 2023; originally announced April 2023.

Comments: 77 pages, 20 figures, This chapter is the preprint of the version currently in production. Please cite this chapter as the following: T.DiMatteo, D. Angles-Alcazar, and F. Shankar. Massive black holes in galactic nuclei: Theory and simulations, in The Encyclopedia of Cosmology (Set 2): Black Holes, edited by Z. Haiman (World Scientific, New Jersey, 2023)

MSC Class: none

Cosmology with one galaxy? -- The ASTRID model and robustness

Authors: Nicolas Echeverri, Francisco Villaescusa-Navarro, Chaitanya Chawak, Yueying Ni, ChangHoon Hahn, Elena Hernandez-Martinez, Romain Teyssier, Daniel Angles-Alcazar, Klaus Dolag, Tiago Castro

Abstract: Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid… ▽ More Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid physics: Astrid. We find that also in this case, neural networks are able to infer the value of $Ω_{\rm m}$ with a $\sim10\%$ precision from the properties of individual galaxies while accounting for astrophysics uncertainties as modeled in CAMELS. This tight relationship is present at all considered redshifts, $z\leq3$, and the stellar mass, the stellar metallicity, and the maximum circular velocity are among the most important galaxy properties behind the relation. In order to use this method with real galaxies, one needs to quantify its robustness: the accuracy of the model when tested on galaxies generated by codes different from the one used for training. We quantify the robustness of the models by testing them on galaxies from four different codes: IllustrisTNG, SIMBA, Astrid, and Magneticum. We show that the models perform well on a large fraction of the galaxies, but fail dramatically on a small fraction of them. Removing these outliers significantly improves the accuracy of the models across simulation codes. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 16 pages, 12 figures

The CAMELS project: Expanding the galaxy formation model space with new ASTRID and 28-parameter TNG and SIMBA suites

Authors: Yueying Ni, Shy Genel, Daniel Anglés-Alcázar, Francisco Villaescusa-Navarro, Yongseok Jo, Simeon Bird, Tiziana Di Matteo, Rupert Croft, Nianyi Chen, Natalí S. M. de Santi, Matthew Gebhardt, Helen Shao, Shivam Pandey, Lars Hernquist, Romeel Dave

Abstract: We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies.… ▽ More We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2,124 hydrodynamic simulation runs that vary 3 cosmological parameters ($Ω_m$, $σ_8$, $Ω_b$) and 4 parameters controlling stellar and AGN feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex non-linear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set. △ Less

Submitted 4 April, 2023; originally announced April 2023.

The impact of AGN-driven winds on physical and observable galaxy sizes

Authors: R. K. Cochrane, D. Anglés-Alcázar, J. Mercedes-Feliz, C. C. Hayward, C. -A. Faucher-Giguère, S. Wellons, B. A. Terrazas, A. Wetzel, P. F. Hopkins, J. Moreno, K. -Y. Su, R. S. Somerville

Abstract: Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half… ▽ More Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multi-wavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ~1 kpc to ~0.1 kpc in <40 Myr when AGN winds are not included, but increases to ~2 kpc when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over 35 Myr, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations. △ Less

Submitted 17 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

A universal equation to predict $Ω_{\rm m}$ from halo and galaxy catalogues

Authors: Helen Shao, Natalí S. M de Santi, Francisco Villaescusa-Navarro, Romain Teyssier, Yueying Ni, Daniel Angles-Alcazar, Shy Genel, Lars Hernquist, Ulrich P. Steinwandel, Tiago Castro, Elena Hernandez-Martınez, Klaus Dolag, Christopher C. Lovell, Eli Visbal, Lehman H. Garrison, Mihir Kulkarni

Abstract: We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our… ▽ More We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $Ω_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $Ω_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $Ω_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $Ω_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 32 pages, 13 figures, summary video: https://youtu.be/STZHvDHkVgo

Robust Field-level Likelihood-free Inference with Galaxies

Authors: Natalí S. M. de Santi, Helen Shao, Francisco Villaescusa-Navarro, L. Raul Abramo, Romain Teyssier, Pablo Villanueva-Domingo, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Elena Hernandez-Martinez, Ulrich P. Steinwandel, Christopher C. Lovell, Klaus Dolag, Tiago Castro, Mark Vogelsberger

Abstract: We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny… ▽ More We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny $(25~h^{-1}{\rm Mpc})^3$ volumes our models can infer the value of $Ω_{\rm m}$ with approximately $12$ % precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and AGN feedback, run with five different codes and subgrid models - IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE -, we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on $1,024$ simulations that cover a vast region in parameter space - variations in $5$ cosmological and $23$ astrophysical parameters - finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network have likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than $\sim10~h^{-1}{\rm kpc}$. △ Less

Submitted 18 July, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 34 pages, 12 figures. For a video summarizing the results, see https://youtu.be/b59ep7cyPOs

Journal ref: Volume 952, Number 1, Year 2023, Pages 69

$\rm [C_{II}]$ 158 $\rm μm$ emission as an indicator of galaxy star formation rate

Authors: Lichen Liang, Robert Feldmann, Norman Murray, Desika Narayanan, Christopher C. Hayward, Daniel Anglés-Alcázar, Luigi Bassini, Alexander J. Richings, Claude-André Faucher-Giguère, Dongwoo T. Chung, Jennifer Y. H. Chan, Doǧa Tolgay, Onur Çatmabacak, Dušan Kereš, Philip F. Hopkins

Abstract: Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-lumin… ▽ More Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments. △ Less

Submitted 6 December, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

Predicting the impact of feedback on matter clustering with machine learning in CAMELS

Authors: Ana Maria Delgado, Daniel Angles-Alcazar, Leander Thiele, Shivam Pandey, Kai Lehman, Rachel S. Somerville, Michelle Ntampaka, Shy Genel, Francisco Villaescusa-Navarro, Lars Hernquist

Abstract: Extracting information from the total matter power spectrum with the precision needed for upcoming cosmological surveys requires unraveling the complex effects of galaxy formation processes on the distribution of matter. We investigate the impact of baryonic physics on matter clustering at $z=0$ using a library of power spectra from the Cosmology and Astrophysics with MachinE Learning Simulations… ▽ More Extracting information from the total matter power spectrum with the precision needed for upcoming cosmological surveys requires unraveling the complex effects of galaxy formation processes on the distribution of matter. We investigate the impact of baryonic physics on matter clustering at $z=0$ using a library of power spectra from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, containing thousands of $(25\,h^{-1}{\rm Mpc})^3$ volume realizations with varying cosmology, initial random field, stellar and AGN feedback strength and sub-grid model implementation methods. We show that baryonic physics affects matter clustering on scales $k \gtrsim 0.4\,h\,\mathrm{Mpc}^{-1}$ and the magnitude of this effect is dependent on the details of the galaxy formation implementation and variations of cosmological and astrophysical parameters. Increasing AGN feedback strength decreases halo baryon fractions and yields stronger suppression of power relative to N-body simulations, while stronger stellar feedback often results in weaker effects by suppressing black hole growth and therefore the impact of AGN feedback. We find a broad correlation between mean baryon fraction of massive halos ($M_{\rm 200c} > 10^{13.5}$\,\Msun) and suppression of matter clustering but with significant scatter compared to previous work owing to wider exploration of feedback parameters and cosmic variance effects. We show that a random forest regressor trained on the baryon content and abundance of halos across the full mass range $10^{10} \leq M_\mathrm{halo}/$\Msun$< 10^{15}$ can predict the effect of galaxy formation on the matter power spectrum on scales $k = 1.0$--20.0\,$h\,\mathrm{Mpc}^{-1}$. △ Less

Submitted 5 October, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Inferring the impact of feedback on the matter distribution using the Sunyaev Zel'dovich effect: Insights from CAMELS simulations and ACT+DES data

Authors: Shivam Pandey, Kai Lehman, Eric J. Baxter, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Francisco Villaescusa-Navarro, Ana Maria Delgado, Tiziana di Matteo

Abstract: Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Ex… ▽ More Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, $k$, we find that the baryon fraction in halos contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel'dovich (tSZ) signal from group-scale halos contains similar information. Using recent Dark Energy Survey (DES) weak lensing and Atacama Cosmology Telescope (ACT) tSZ cross-correlation measurements and models trained on CAMELS, we obtain $10\%$ constraints on feedback effects on the power spectrum at $k \sim 5\, h/{\rm Mpc}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal{O}(<1\%)$ at $k = 1\, h/{\rm Mpc}$ and $\mathcal{O}(3\%)$ at $k = 5\, h/{\rm Mpc}$ using the methods that we introduce here. Finally, we investigate the impact of feedback on the matter bispectrum, finding that tSZ observables are highly informative in this case. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 18 pages, 15 figures. Comments are welcome

Local positive feedback in the overall negative: the impact of quasar winds on star formation in the FIRE cosmological simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Christopher C. Hayward, Rachel K. Cochrane, Bryan A. Terrazas, Sarah Wellons, Alexander J. Richings, Claude-André Faucher-Giguère, Jorge Moreno, Kung Yi Su, Philip F. Hopkins, Eliot Quataert, Dušan Kereš

Abstract: Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using co… ▽ More Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environments (FIRE) project, including a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at $z \sim 2$ ($M_{\rm halo} \sim 10^{12.5} \, {\rm M}_{\odot}$), we show that strong quasar winds with kinetic power $\sim$10$^{46}$ erg/s acting for $>$20$\,$Myr drive the formation of a central gas cavity and can dramatically reduce the star formation rate surface density across the galaxy disc. The suppression of star formation is primarily driven by reducing the amount of gas that can become star-forming, compared to directly evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the global negative impact of quasar winds, we identify several plausible signatures of local positive feedback, including: (1) spatial anti-correlation of wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas near the edge of the cavity, and (3) increased local contribution of outflowing material to star formation. Stars forming under the presence of quasar winds tend to do so at larger radial distances. Our results suggest that positive and negative AGN feedback can coexist in galaxies, but local positive triggering of star formation plays a minor role in global galaxy growth. △ Less

Submitted 1 August, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 17 pages, 12 figures. Accepted for publication in MNRAS

Calibrating cosmological simulations with implicit likelihood inference using galaxy growth observables

Authors: Yongseok Jo, Shy Genel, Benjamin Wandelt, Rachel Somerville, Francisco Villaescusa-Navarro, Greg L. Bryan, Daniel Angles-Alcazar, Daniel Foreman-Mackey, Dylan Nelson, Ji-hoon Kim

Abstract: In a novel approach employing implicit likelihood inference (ILI), also known as likelihood-free inference, we calibrate the parameters of cosmological hydrodynamic simulations against observations, which has previously been unfeasible due to the high computational cost of these simulations. For computational efficiency, we train neural networks as emulators on ~1000 cosmological simulations from… ▽ More In a novel approach employing implicit likelihood inference (ILI), also known as likelihood-free inference, we calibrate the parameters of cosmological hydrodynamic simulations against observations, which has previously been unfeasible due to the high computational cost of these simulations. For computational efficiency, we train neural networks as emulators on ~1000 cosmological simulations from the CAMELS project to estimate simulated observables, taking as input the cosmological and astrophysical parameters, and use these emulators as surrogates to the cosmological simulations. Using the cosmic star formation rate density (SFRD) and, separately, stellar mass functions (SMFs) at different redshifts, we perform ILI on selected cosmological and astrophysical parameters (Omega_m, sigma_8, stellar wind feedback, and kinetic black hole feedback) and obtain full 6-dimensional posterior distributions. In the performance test, the ILI from the emulated SFRD (SMFs) can recover the target observables with a relative error of 0.17% (0.4%). We find that degeneracies exist between the parameters inferred from the emulated SFRD, confirmed with new full cosmological simulations. We also find that the SMFs can break the degeneracy in the SFRD, which indicates that the SMFs provide complementary constraints for the parameters. Further, we find that the parameter combination inferred from an observationally-inferred SFRD reproduces the target observed SFRD very well, whereas, in the case of the SMFs, the inferred and observed SMFs show significant discrepancies that indicate potential limitations of the current galaxy formation modeling and calibration framework, and/or systematic differences and inconsistencies between observations of the stellar mass function. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: This is the revised version from the reviewer's report (submitted to ApJ)

Predicting sub-millimeter flux densities from global galaxy properties

Authors: R. K. Cochrane, C. C. Hayward, D. Angles-Alcazar, R. S. Somerville

Abstract: Recent years have seen growing interest in post-processing cosmological simulations with radiative transfer codes to predict observable fluxes for simulated galaxies. However, this can be slow, and requires a number of assumptions in cases where simulations do not resolve the ISM. Zoom-in simulations better resolve the detailed structure of the ISM and the geometry of stars and gas, however statis… ▽ More Recent years have seen growing interest in post-processing cosmological simulations with radiative transfer codes to predict observable fluxes for simulated galaxies. However, this can be slow, and requires a number of assumptions in cases where simulations do not resolve the ISM. Zoom-in simulations better resolve the detailed structure of the ISM and the geometry of stars and gas, however statistics are limited due to the computational cost of simulating even a single halo. In this paper, we make use of a set of high resolution, cosmological zoom-in simulations of massive M_star>10^10.5M_sol at z=2), star-forming galaxies from the FIRE suite. We run the SKIRT radiative transfer code on hundreds of snapshots in the redshift range 1.5<z<5 and calibrate a power law scaling relation between dust mass, star formation rate and 870um flux density. The derived scaling relation shows encouraging consistency with observational results from the sub-millimeter-selected AS2UDS sample. We extend this to other wavelengths, deriving scaling relations between dust mass, stellar mass, star formation rate and redshift and sub-millimeter flux density at observed-frame wavelengths between 340um and 870um. We then apply the scaling relations to galaxies drawn from EAGLE, a large box cosmological simulation. We show that the scaling relations predict EAGLE sub-millimeter number counts that agree well with previous results that were derived using far more computationally expensive radiative transfer techniques. Our scaling relations can be applied to other simulations and semi-analytical or semi-empirical models to generate robust and fast predictions for sub-millimeter number counts. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 14 pages, 11 figures; resubmitted to MNRAS following minor comments from reviewer

A unified model for the co-evolution of galaxies and their circumgalactic medium: the relative roles of turbulence and atomic cooling physics

Authors: Viraj Pandya, Drummond B. Fielding, Greg L. Bryan, Christopher Carr, Rachel S. Somerville, Jonathan Stern, Claude-Andre Faucher-Giguere, Zachary Hafen, Daniel Angles-Alcazar, John C. Forbes

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 co-evolve 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… ▽ More 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 co-evolve 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 $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. △ Less

Submitted 27 August, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model

Dust temperature uncertainties hamper the inference of dust and molecular gas masses from the dust continuum emission of quiescent high-redshift galaxies

Authors: R. K. Cochrane, C. C. Hayward, D. Anglés-Alcázar

Abstract: Single flux density measurements at observed-frame sub-millimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive haloes from the Feedback in Realistic Environments (FIRE) project. Our starting p… ▽ More Single flux density measurements at observed-frame sub-millimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive haloes from the Feedback in Realistic Environments (FIRE) project. Our starting point is four star-forming, central galaxies at seven redshifts between z=1.5 and z=4.5. We generate modified quiescent galaxies that have been quenched for 100Myr, 500Myr, or 1Gyr prior to each of the studied redshifts by re-assigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature, T_dust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumed T_dust~25K, we show that the ratio between dust mass and 1.3mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non star-forming galaxies. This can lead to an underestimation of dust mass (and, when sub-mm flux density is used as a proxy for molecular gas content, gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3mm flux density probes rest-frame wavelengths far from the Rayleigh-Jeans regime, and hence depends super-linearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 10 pages, 4 figures; accepted for publication in ApJL

Stellar feedback-regulated black hole growth: driving factors from nuclear to halo scales

Authors: Lindsey Byrne, Claude-André Faucher-Giguère, Jonathan Stern, Daniel Anglés-Alcázar, Sarah Wellons, Alexander B. Gurvich, Philip F. Hopkins

Abstract: Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are under-massive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project,… ▽ More Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are under-massive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multi-channel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fueling: constant thresholds in these properties are typically crossed within ~0.1 Hubble time of accelerated BH fueling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold Sigma1 ~ 10^9.5 Msun/kpc^2, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived, thin gas disks, as well as with virialization of the inner circumgalactic medium. The halo mass Mh ~ 10^12 Msun and stellar mass Mstar ~ 10^10.5 Msun at which BH growth accelerates correspond to ~L* galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above ~L*. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 21 pages, 9 figures, submitted to MNRAS

Efficient long-range AGN feedback affects the low redshift Lyman-$α$ forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Romeel Davé, Shy Genel

Abstract: Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedbac… ▽ More Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedback. We show that TNG over-predicts the number density of absorbers at column densities $N_{\rm HI} < 10^{14}$ cm$^{-2}$ compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model which includes AGN jet feedback is the preferred fit to the observed CDD of the $z = 0.1$ Lyman-$α$ forest across five orders of magnitude in column density. We show that the Simba results with jets produce a quantitatively better fit to the observational data than the Simba results without jets, even when the UVB is left as a free parameter. AGN jets in Simba are high speed, collimated, weakly-interacting with the interstellar medium (via brief hydrodynamic decoupling) and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the IGM when compared to TNG's kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low redshift Lyman-$α$ forest provides plausible evidence for long-range AGN jet feedback. △ Less

Submitted 7 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 15 pages, 5 figures, 1 table, accepted to ApJL

Journal ref: ApJL 945 L17 (2023)

Robust field-level inference with dark matter halos

Authors: Helen Shao, Francisco Villaescusa-Navarro, Pablo Villanueva-Domingo, Romain Teyssier, Lehman H. Garrison, Marco Gatti, Derek Inman, Yueying Ni, Ulrich P. Steinwandel, Mihir Kulkarni, Eli Visbal, Greg L. Bryan, Daniel Angles-Alcazar, Tiago Castro, Elena Hernandez-Martinez, Klaus Dolag

Abstract: We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, co… ▽ More We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 25 pages, 11 figures, summary video: https://youtu.be/qkw92Z6owJU

The SZ flux-mass ($Y$-$M$) relation at low halo masses: improvements with symbolic regression and strong constraints on baryonic feedback

Authors: Digvijay Wadekar, Leander Thiele, J. Colin Hill, Shivam Pandey, Francisco Villaescusa-Navarro, David N. Spergel, Miles Cranmer, Daisuke Nagai, Daniel Anglés-Alcázar, Shirley Ho, Lars Hernquist

Abstract: Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with exten… ▽ More Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the $Y-M$ relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, h^{-1} \, M_\odot$); we find that simply replacing $Y\rightarrow Y(1+M_*/M_\mathrm{gas})$ in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations. We also forecast that measurements of the $Y-M$ relation could provide percent-level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g., SO, CMB-S4) and galaxy surveys (e.g., DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the an alternative relation, $Y-M_*$, provides complementary information on feedback than $Y-M$ △ Less

Submitted 28 April, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: Version appearing in MNRAS. Minor change to Fig.6 and added Fig. A5 compared to the previous version. 7+5 pages. The code and data associated with this paper are available at https://github.com/JayWadekar/ScalingRelations_ML

Journal ref: MNRAS, 522, 2628 (2023)

Constraining cosmology with machine learning and galaxy clustering: the CAMELS-SAM suite

Authors: Lucia A. Perez, Shy Genel, Francisco Villaescusa-Navarro, Rachel S. Somerville, Austen Gabrielpillai, Daniel Anglés-Alcázar, Benjamin D. Wandelt, L. Y. Aaron Yung

Abstract: As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophy… ▽ More As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project: CAMELS-SAM, encompassing one thousand dark-matter only simulations of (100 $h^{-1}$ cMpc)$^3$ with different cosmological parameters ($Ω_m$ and $σ_8$) and run through the Santa Cruz semi-analytic model for galaxy formation over a broad range of astrophysical parameters. As a proof-of-concept for the power of this vast suite of simulated galaxies in a large volume and broad parameter space, we probe the power of simple clustering summary statistics to marginalize over astrophysics and constrain cosmology using neural networks. We use the two-point correlation function, count-in-cells, and the Void Probability Function, and probe non-linear and linear scales across $0.68<$ R $<27\ h^{-1}$ cMpc. Our cosmological constraints cluster around 3-8$\%$ error on $Ω_{\text{M}}$ and $σ_8$, and we explore the effect of various galaxy selections, galaxy sampling, and choice of clustering statistics on these constraints. We additionally explore how these clustering statistics constrain and inform key stellar and galactic feedback parameters in the Santa Cruz SAM. CAMELS-SAM has been publicly released alongside the rest of CAMELS, and offers great potential to many applications of machine learning in astrophysics: https://camels-sam.readthedocs.io. △ Less

Submitted 22 May, 2023; v1 submitted 5 April, 2022; originally announced April 2022.

Comments: 28 pages main text, additional 12 in appendices, 23 figures, 9 tables

Exploring supermassive black hole physics and galaxy quenching across halo mass in FIRE cosmological zoom simulations

Authors: Sarah Wellons, Claude-André Faucher-Giguère, Philip F. Hopkins, Eliot Quataert, Daniel Anglés-Alcázar, Robert Feldmann, Christopher C. Hayward, Dušan Kereš, Kung-Yi Su, Andrew Wetzel

Abstract: Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosm… ▽ More Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosmological zoom-in simulations across a wide range of halo mass (10^10-10^13 Msun). Within the suite, we vary the numerical schemes for BH accretion and feedback, the accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models that satisfy the observational constraints, and for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many successful models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, over-massive BHs tend to over-quench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass halos, but by feedback efficiency in low-mass halos. The amount of star formation suppression by SMBH feedback in low-mass halos is determined primarily by the time-integrated feedback energy. For massive galaxies, the "responsiveness" of a model (i.e. how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 18 pages, 9 figures. Submitted to MNRAS, comments welcome

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Authors: Jorge Moreno, Shany Danieli, James S. Bullock, Robert Feldmann, Philip F. Hopkins, Onur Catmabacak, Alexander Gurvich, Alexandres Lazar, Courtney Klein, Cameron B. Hummels, Zachary Hafen, Francisco J. Mercado, Sijie Yu, Fangzhou Jiang, Coral Wheeler, Andrew Wetzel, Daniel Angles-Alcazar, Michael Boylan-Kolchin, Eliot Quataert, Claude-Andre Faucher-Giguere, Dusan Keres

Abstract: The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark ma… ▽ More The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 55 pages, 4 figures, 13 supplementary figures. 29 pages, 4 figures. Accepted for publication in Nature Astronomy. To appear on 14-February-2022. Published version: https://www.nature.com/articles/s41550-021-01598-4

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

Authors: Melanie Habouzit, Masafusa Onoue, Eduardo Banados, Marcel Neeleman, Daniel Angles-Alcazar, Fabian Walter, Annalisa Pillepich, Romeel Dave, Knud Jahnke, Yohan Dubois

Abstract: The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars w… ▽ More The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with L_bol= 10^45-10^46 erg/s opens a path to understand the build-up of more normal BHs at z>6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z>4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the M_BH-M_star scaling relation at z=0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z>4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with L_bol= 10^45-10^46 erg/s observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z=6, with the same diversity of sub-grid physics, in order to gain statistics on the most extreme objects at high redshift. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: published in MNRAS, 19 pages, 8 figures, key figures: Fig. 3, Fig.5, and Fig. 7

The black hole population in low-mass galaxies in large-scale cosmological simulations

Authors: Houda Haidar, Melanie Habouzit, Marta Volonteri, Mar Mezcua, Jenny Greene, Nadine Neumayer, Daniel Angles-Alcazar, Ignacio Martin-Navarro, Nils Hoyer, Yohan Dubois, Romeel Dave

Abstract: Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that ar… ▽ More Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction at z=0 is not inherited from the simulation seeding modeling. The ability of BHs and their host galaxies to power an AGN depends on BH and galaxy subgrid modeling. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus can be used to differentiate galaxy formation models. AGN fractions at z=0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in the masses of simulated and observed BHs, BH occupation fraction affected by numerical choices, and an unknown fraction of obscured AGN. Our work advocates for more thorough comparisons with observations to improve the modeling of cosmological simulations, and our understanding of BH and galaxy physics in the low-mass regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time. △ Less

Submitted 14 June, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted in MNRAS, 21 pages, 11 figures, 1 table

Hot-mode accretion and the physics of thin-disk galaxy formation

Authors: Zachary Hafen, Jonathan Stern, James Bullock, Alex B. Gurvich, Sijie Yu, Claude-Andre Faucher-Giguere, Drummond B. Fielding, Daniel Angles-Alcazar, Eliot Quataert, Andrew Wetzel, Tjitske Starkenburg, Michael Boylan-Kolchin, Jorge Moreno, Robert Feldmann, Kareem El-Badry, T. K. Chan, Cameron Trapp, Dusan Keres, Philip F. Hopkins

Abstract: We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks 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 (>70%) of their young stars in a thin disk (h/R~0.1) we find that:… ▽ More We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks 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 (>70%) of their young stars in a thin disk (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 T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-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 disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk 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 onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe. △ Less

Submitted 6 June, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 17 pages, 11 figures. Comments welcome. Accompanying interactive visualization at zhafen.github.io/rotating-cooling-flows

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Authors: Andrina Nicola, Francisco Villaescusa-Navarro, David N. Spergel, Jo Dunkley, Daniel Anglés-Alcázar, Romeel Davé, Shy Genel, Lars Hernquist, Daisuke Nagai, Rachel S. Somerville, Benjamin D. Wandelt

Abstract: Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, an… ▽ More Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 31 pages, 10 figures, to be submitted to JCAP

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

Authors: Emily Moser, Nicholas Battaglia, Daisuke Nagai, Erwin Lau, Luis Fernando Machado Poletti Valle, Francisco Villaescusa-Navarro, Stefania Amodeo, Daniel Angles-Alcazar, Greg L. Bryan, Romeel Dave, Lars Hernquist, Mark Vogelsberger

Abstract: The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Lea… ▽ More The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: 23 pages, 7 figures, comments welcome

Cosmology with one galaxy?

Authors: Francisco Villaescusa-Navarro, Jupiter Ding, Shy Genel, Stephanie Tonnesen, Valentina La Torre, David N. Spergel, Romain Teyssier, Yin Li, Caroline Heneka, Pablo Lemos, Daniel Anglés-Alcázar, Daisuke Nagai, Mark Vogelsberger

Abstract: Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulatio… ▽ More Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Comments: 20+6 pages, 15 figures, data and codes to reproduce the results publicly available at https://github.com/franciscovillaescusa/Cosmo1gal

Percent-level constraints on baryonic feedback with spectral distortion measurements

Authors: Leander Thiele, Digvijay Wadekar, J. Colin Hill, Nicholas Battaglia, Jens Chluba, Francisco Villaescusa-Navarro, Lars Hernquist, Mark Vogelsberger, Daniel Anglés-Alcázar, Federico Marinacci

Abstract: High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodyn… ▽ More High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models. △ Less

Submitted 27 April, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

Comments: 10+5 pages, 7 figures; matches PRD published version

Augmenting astrophysical scaling relations with machine learning: application to reducing the Sunyaev-Zeldovich flux-mass scatter

Authors: Digvijay Wadekar, Leander Thiele, Francisco Villaescusa-Navarro, J. Colin Hill, Miles Cranmer, David N. Spergel, Nicholas Battaglia, Daniel Anglés-Alcázar, Lars Hernquist, Shirley Ho

Abstract: Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple exten… ▽ More Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple extensions to existing relations) in abstract high-dimensional parameter spaces. We use a machine learning tool called symbolic regression (SR), which models patterns in a dataset in the form of analytic equations. We focus on the Sunyaev-Zeldovich flux$-$cluster mass relation ($Y_\mathrm{SZ}-M$), the scatter in which affects inference of cosmological parameters from cluster abundance data. Using SR on the data from the IllustrisTNG hydrodynamical simulation, we find a new proxy for cluster mass which combines $Y_\mathrm{SZ}$ and concentration of ionized gas ($c_\mathrm{gas}$): $M \propto Y_\mathrm{conc}^{3/5} \equiv Y_\mathrm{SZ}^{3/5} (1-A\, c_\mathrm{gas})$. $Y_\mathrm{conc}$ reduces the scatter in the predicted $M$ by $\sim 20-30$\% for large clusters ($M\gtrsim 10^{14}\, h^{-1} \, M_\odot$), as compared to using just $Y_\mathrm{SZ}$. We show that the dependence on $c_\mathrm{gas}$ is linked to cores of clusters exhibiting larger scatter than their outskirts. Finally, we test $Y_\mathrm{conc}$ on clusters from CAMELS simulations and show that $Y_\mathrm{conc}$ is robust against variations in cosmology, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from upcoming CMB and X-ray surveys like ACT, SO, eROSITA and CMB-S4. △ Less

Submitted 17 March, 2023; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: Version appearing in PNAS. Added additional tests but results unchanged compared to previous version. The code and data associated with this paper are available at https://github.com/JayWadekar/ScalingRelations_ML

Journal ref: PNAS March 17, 2023, 120 (12) e2202074120

The CAMELS project: public data release

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch , et al. (22 additional authors not shown)

Abstract: The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present… ▽ More The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$α$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org