skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: simba-EoR : early galaxy formation in the simba simulation including a new sub-grid interstellar medium model
ABSTRACT We update the dust model present within the simba galaxy simulations with a self-consistent framework for the co-evolution of dust and molecular hydrogen populations in the interstellar medium, and use this to explore $$z \ge 6$$ galaxy evolution. In addition to tracking the evolution of dust and molecular hydrogen abundances, our model fully integrates these species into the simba simulation, explicitly modelling their impact on physical processes such as star formation and cooling through the inclusion of a novel two-phase sub-grid model for interstellar gas. Running two cosmological simulations down to $$z \sim 6$$ we find that our simba-EoR model displays a generally tighter concordance with observational data than fiducial simba. Additionally we observe that our simba-EoR models increase star formation activity at early epochs, producing larger dust-to-gas ratios consequently. Finally, we discover a significant population of hot dust at $$\sim 100$$ K, aligning with contemporaneous observations of high-redshift dusty galaxies, alongside the large $$\sim 20$$ K population typically identified.  more » « less
Award ID(s):
1909153
PAR ID:
10554344
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
535
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 1293-1314
Size(s):
p. 1293-1314
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We study specific star formation rate (sSFR) and gas profiles of star-forming (SF) and green valley (GV) galaxies in the simba cosmological hydrodynamic simulation. SF galaxy half-light radii (Rhalf) at z = 0 and their evolution (∝(1 + z)−0.78) agree with observations. Passive galaxy Rhalf agree with observations at high redshift, but by z = 0 are too large, owing to numerical heating. We compare simbaz = 0 sSFR radial profiles for SF and GV galaxies to observations. simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in GV galaxies and massive SF systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, simba yields GV galaxies with strongly dropping sSFR profiles beyond ≳Rhalf, regardless of active galactic nucleus feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. simba’s satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At z = 2, simba does not show central depressions in massive SF galaxies, suggesting simba’s X-ray feedback should be more active at high-z. High-resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in z = 0 GV galaxies represents a unique success of simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms. 
    more » « less
  2. Abstract We explore the cosmic evolution of the fraction of dust-obscured star formation predicted by thesimbacosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up toz= 3, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, atz> 3 we find strong evolution at fixed galaxy stellar mass toward greater amounts of obscured star formation, in tension with high-redshift observations. We explain the trend of increasing obscuration at higher redshifts by evolving star-dust geometry, as the dust-to-stellar mass ratios remain relatively constant across cosmic time. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust-to-stellar mass ratios at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust-obscured star formation history peaks aroundz∼ 2−3, and becomes subdominant atz≳ 5. The dominance of obscured star formation at redshiftsz≲ 4 is consistent with our results for the evolution of the obscured star formation fraction at fixed stellar mass to higher values at higher redshift because there exist fewer massive, heavily obscured galaxies at high redshift. 
    more » « less
  3. null (Ed.)
    ABSTRACT We examine the 1.4 GHz radio luminosities of galaxies arising from star formation and active galactic nuclei (AGNs) within the state-of-the-art cosmological hydrodynamic simulation Simba. Simba grows black holes via gravitational torque limited accretion from cold gas and Bondi accretion from hot gas, and employs AGN feedback including jets at low Eddington ratios. We define a population of radio loud AGNs (RLAGNs) based on the presence of ongoing jet feedback. Within RLAGN, we define high and low excitation radio galaxies (HERGs and LERGs) based on their dominant mode of black hole accretion: torque limited accretion representing feeding from a cold disc, or Bondi representing advection-dominated accretion from a hot medium. Simba predicts good agreement with the observed radio luminosity function (RLF) and its evolution, overall as well as separately for HERGs and LERGs. Quiescent galaxies with AGN-dominated radio flux dominate the RLF at $$\gtrsim 10^{22-23}$$ W Hz−1, while star formation dominates at lower radio powers. Overall, RLAGNs have higher black hole accretion rates and lower star formation rates than non-RLAGN at a given stellar mass or velocity dispersion, but have similar black hole masses. Simba predicts an LERG number density of 8.53 Mpc−3, ∼10× higher than for HERGs, broadly as observed. While LERGs dominate among most massive galaxies with the largest black holes and HERGs dominate at high specific star formation rates, they otherwise largely populate similar-sized dark matter haloes and have similar host galaxy properties. Simba thus predicts that deeper radio surveys will reveal an increasing overlap between the host galaxy demographics of HERGs and LERGs. 
    more » « less
  4. ABSTRACT We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties. 
    more » « less
  5. Abstract Over the past decade, rest-frame color–color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color–color selection tool—the rest-frameU−VversusV−Jdiagram—by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties inUVJspace and the color–color evolution of massive galaxies as they quench at redshiftsz∼ 1–2. We find that our models broadly reproduce the observedUVJdiagram atz= 1–2, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of theUVJdiagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley insimba. When investigating the time evolution of galaxies on theUVJdiagram, we find that the quenching pathway on theUVJdiagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies. 
    more » « less