We examine the properties of damped Lyman-α absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: simba and technicolor dawn. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multifrequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly α forest mean transmission at z > 5. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at z ≥ 3, and both reproduce the observed neutral hydrogen cosmological mass density, $\Omega _{\rm H\, \small{I}}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows’ feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows underproduces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an H2-regulated star formation recipe includes a tail extending to [M/H] ≪ −3, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an important role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions.
Reproducing submillimetre galaxy number counts with cosmological hydrodynamic simulations
ABSTRACT Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity as a function of SFR and dust mass. Our results demonstrate that solutions to the discrepancy between submm counts in simulations and observations, such as a top-heavy initial mass function, are unnecessary, and that submillimetre-bright phases are a natural consequence of massive galaxy evolution.
more »
« less
- NSF-PAR ID:
- 10390151
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 502
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 772 to 793
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
Abstract Recent discoveries of a significant population of bright galaxies at cosmic dawn
z ≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright (M UV≲ −20) galaxies atz ≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models. -
ABSTRACT In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simulations, at any redshift.more » « less
-
null (Ed.)ABSTRACT The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$ in galaxies of $M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$. While the time evolution of the mean MBH − M⋆ relation is mild ($\rm \Delta M_{\rm BH}\leqslant 1\, dex$ for 0 $\leqslant z \leqslant$ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs ($\log _{10}\, f_{\rm Edd}\geqslant -2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.more » « less
-
The Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) targets the [CII] 158 μ m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z = 4.4 and z = 5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details regarding the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates (SFRs), we measured the conversion factor from the ALMA 158 μ m rest-frame dust continuum luminosity to the total infrared luminosity ( L IR ) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L IR of 4.4 × 10 11 L ⊙ . We also detected 57 additional continuum sources in our ALMA pointings. They are at a lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5 ± 0.2. We measured the 850 μ m number counts between 0.35 and 3.5 mJy, thus improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40% of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8 × 10 8 L ⊙ and their median full width at half maximum is 252 km s −1 . After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR– L [CII] relations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa4043
