skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Realistic mock observations of the sizes and stellar mass surface densities of massive galaxies in FIRE-2 zoom-in simulations
Title: Realistic mock observations of the sizes and stellar mass surface densities of massive galaxies in FIRE-2 zoom-in simulations
ABSTRACT The galaxy size–stellar mass and central surface density–stellar mass relationships are fundamental observational constraints on galaxy formation models. However, inferring the physical size of a galaxy from observed stellar emission is non-trivial due to various observational effects, such as the mass-to-light ratio variations that can be caused by non-uniform stellar ages, metallicities, and dust attenuation. Consequently, forward-modelling light-based sizes from simulations is desirable. In this work, we use the skirt  dust radiative transfer code to generate synthetic observations of massive galaxies ($M_{*}\sim 10^{11}\, \rm {M_{\odot }}$ at z = 2, hosted by haloes of mass $M_{\rm {halo}}\sim 10^{12.5}\, \rm {M_{\odot }}$) from high-resolution cosmological zoom-in simulations that form part of the Feedback In Realistic Environments project. The simulations used in this paper include explicit stellar feedback but no active galactic nucleus (AGN) feedback. From each mock observation, we infer the effective radius (Re), as well as the stellar mass surface density within this radius and within $1\, \rm {kpc}$ (Σe and Σ1, respectively). We first investigate how well the intrinsic half-mass radius and stellar mass surface density can be inferred from observables. The majority of predicted sizes and surface densities are within a factor of 2 of the intrinsic values. We then compare our predictions to the observed size–mass relationship and the Σ1−M⋆ and Σe−M⋆ relationships. At z ≳ 2, the simulated massive galaxies are in general agreement with observational scaling relations. At z ≲ 2, they evolve to become too compact but still star forming, in the stellar mass and redshift regime where many of them should be quenched. Our results suggest that some additional source of feedback, such as AGN-driven outflows, is necessary in order to decrease the central densities of the simulated massive galaxies to bring them into agreement with observations at z ≲ 2.  more » « less
Award ID(s):
2009687 2001905 1715216
NSF-PAR ID:
10237289
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
501
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
1591 to 1602
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    Without active galactic nucleus (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 }\sim 10^{11}\, \rm {{\rm M}_{\odot }}$) 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 multiwavelength 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 to $\sim 0.1\, \rm {kpc}$ in $\lesssim 40\, \rm {Myr}$ when AGN winds are not included, but increases to $\sim 2\, \rm {kpc}$ when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over $35\, \rm {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.

     
    more » « less
  2. ; ; ; ; ; ; ; ; ; ; et al ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (∼1000 galaxies with $M_{\rm star}\gt 10^8\, M_\odot$ at z  = 0) at a resolution ($\Delta {}x\sim {}20\, {\rm pc}$ , $m_{\rm b}\sim {}6\times {}10^4\, M_\odot$ ) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L = 22.1 Mpc) thanks to its exceptionally high dynamic range (≳106) and the inclusion of multichannel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that star formation rates, gas masses, and metallicities of simulated galaxies with $M_{\rm star}\lt 10^{10.5-11}\, M_\odot$ broadly agree with observations. These galaxy scaling relations extend to low masses ($M_{\rm star}\sim {}10^7\, M_\odot$ ) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z ∼ 0–5. At low z , FIREbox predicts a peak in the stellar-mass–halo-mass relation but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a ΛCDM Universe while also highlighting modelling challenges to be addressed in next-generation galaxy simulations. 
    more » « less
  3. ; ; ; ; ; ; ; ; ; ; et al ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We perform a consistent comparison of the mass and mass profiles of massive (M⋆ > 1011.4 M⊙) central galaxies at z ∼ 0.4 from deep Hyper Suprime-Cam (HSC) observations and from the Illustris, TNG100, and Ponos simulations. Weak lensing measurements from HSC enable measurements at fixed halo mass and provide constraints on the strength and impact of feedback at different halo mass scales. We compare the stellar mass function (SMF) and the Stellar-to-Halo Mass Relation (SHMR) at various radii and show that the radius at which the comparison is performed is important. In general, Illustris and TNG100 display steeper values of α where $M_{\star } \propto M_{\rm vir}^{\alpha }$. These differences are more pronounced for Illustris than for TNG100 and in the inner rather than outer regions of galaxies. Differences in the inner regions may suggest that TNG100 is too efficient at quenching in situ star formation at Mvir ≃ 1013 M⊙ but not efficient enough at Mvir ≃ 1014 M⊙. The outer stellar masses are in excellent agreement with our observations at Mvir ≃ 1013 M⊙, but both Illustris and TNG100 display excess outer mass as Mvir ≃ 1014 M⊙ (by ∼0.25 and ∼0.12 dex, respectively). We argue that reducing stellar growth at early times in $M_\star \sim 10^{9-10} \, \mathrm{M}_{\odot }$ galaxies would help to prevent excess ex-situ growth at this mass scale. The Ponos simulations do not implement AGN feedback and display an excess mass of ∼0.5 dex at r < 30 kpc compared to HSC which is indicative of overcooling and excess star formation in the central regions. The comparison of the inner profiles of Ponos and HSC suggests that the physical scale over which the central AGN limits star formation is r ≲ 20 kpc. Joint comparisons between weak lensing and galaxy stellar profiles are a direct test of whether simulations build and deposit galaxy mass in the correct dark matter haloes and thereby provide powerful constraints on the physics of feedback and galaxy growth. Our galaxy and weak lensing profiles are publicly available to facilitate comparisons with other simulations. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z = 6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified ‘delayed cooling’ feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n ∼ 103–105 cm−3). The sites where star formation ignites are characterized by maximum densities of the order of 105 cm−3 and gravitational pressures Pgrav/k >107 K cm−3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z = 6.14 is 3.4$\times 10^7~\rm M_{\odot }$, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of $\sim 10^6~\rm M_{\odot }$ and half-mass sizes of ∼100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one orders of magnitude lower than the observed one. In the size–mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

     
    more » « less
  5. ; ; ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy–halo size relation (GHSR) for low-mass ($M_\star \sim 10^{7-9}\, {\rm M}_\odot$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size R1/2 and the parent dark matter halo virial radius Rvir. This relation evolves only weakly since redshift z = 5: $R_{1/2}\, [{\rm kpc}] = (0.053\pm 0.002)(R_{\rm vir}/35\, {\rm kpc})^{0.934\pm 0.054}$, with a nearly constant scatter $\langle \sigma \rangle = 0.084\, [{\rm dex}]$. While this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios vrot/vcirc ∼ 0.15. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate – including spin and concentration – suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the GHSR. Given the relatively small scatter and the weak dependence of the GHSR on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email