Search for: All records

Abstract We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω =v_circ/Rvia the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation. 

Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over $Z\prime \equiv Z/Z_{\odot }\sim 0.1-3$, gas surface density Σ_gas∼ 5–150M_⊙pc⁻², and stellar surface density Σ_star∼ 1–50M_⊙pc⁻². The range of emergent star formation rate surface density is Σ_SFR∼ 10⁻⁴–0.5M_⊙kpc⁻²yr⁻¹, and ISM total midplane pressure isP_tot/k_B= 10³–10⁶cm⁻³K, withP_totequal to the ISM weight $\mathcal{W}$. For given Σ_gasand Σ_star, we find ${\mathrm{\Sigma }}_{\mathrm{SFR}}\propto Z{\prime }^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to Σ_SFR. The thermal feedback yield varies sensitively with both weight and metallicity as ${\mathrm{\Upsilon }}_{\mathrm{th}}\propto {\mathcal{W}}^{-0.46}Z{\prime }^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon }}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{-0.22}Z{\prime }^{-0.18}$. The reduction in Σ_SFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 

Abstract Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of V. Springel & L. Hernquist (hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially resolved simulations of the turbulent, multiphase, star-forming interstellar medium (ISM) have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local-Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the pressure-regulated feedback-modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the “native” simulation outputs with postprocessed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step toward implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations. 

Abstract It has been claimed that traditional models struggle to explain the tentative detection of the 21 cm absorption trough centered atz∼ 17 measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power law of deep Hubble Space Telescope observations of 4 <z< 9 galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction atz∼ 6–8, Thomson-scattering optical depth, and ionizing emissivity atz∼ 5. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile existing high-redshift (cosmic dawn) and low-redshift (reionization) constraints. In particular, we find that high contributions from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with a gradual and longer duration, resulting in a higher optical depth. While the models dominated by faint galaxies successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically motivated semianalytical prescription produces results that are consistent with the (sophisticated) state-of-the-artTHESANradiation-magnetohydrodynamic simulation of the reionization. 

Abstract We explore the role of galactic feedback on the low-redshift Lyα(Lyα) forest (z≲ 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 subgrid 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 theb-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-zLyαforest can be sensitive to changes in the ultraviolet background, stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM. 

Abstract Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modeling reionization of using more detailed models for the ionizing sources as identified within the 100 h −1 Mpc cosmological hydrodynamic simulation S imba , coupled with postprocessed radiative transfer. Comparing with simple (one-to-one) models, the main difference with using S imba sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5%–10%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects. 

Abstract A wealth of cosmological and astrophysical information is expected from many ongoing and upcoming large-scale surveys. It is crucial to prepare for these surveys now and develop tools that can efficiently extract most information. We present HIF low : a fast generative model of the neutral hydrogen (H i ) maps that is conditioned only on cosmology (Ω m and σ 8 ) and designed using a class of normalizing flow models, the masked autoregressive flow. HIF low is trained on the state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. HIF low has the ability to generate realistic diverse maps without explicitly incorporating the expected two-dimensional maps structure into the flow as an inductive bias. We find that HIF low is able to reproduce the CAMELS average and standard deviation H i power spectrum within a factor of ≲2, scoring a very high R 2 > 90%. By inverting the flow, HIF low provides a tractable high-dimensional likelihood for efficient parameter inference. We show that the conditional HIF low on cosmology is successfully able to marginalize over astrophysics at the field level, regardless of the stellar and AGN feedback strengths. This new tool represents a first step toward a more powerful parameter inference, maximizing the scientific return of future H i surveys, and opening a new avenue to minimize the loss of complex information due to data compression down to summary statistics. 

Abstract We present the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) Multifield Data set (CMD), a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from more than 2000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span ∼100 million light-years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N -body simulations from the CAMELS project. Designed to train machine-learning models, CMD is the largest data set of its kind containing more than 70 TB of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io . 

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 4233 cosmological simulations, 2049 N -body simulations, and 2184 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 catalogs, power spectra, bispectra, Ly α spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over 1000 catalogs that contain billions of galaxies from CAMELS-SAM: a large collection of N -body simulations that have been combined with the Santa Cruz semianalytic 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 https://camels.readthedocs.io .