Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published February 2021 | Accepted Version + Published
Journal Article Open

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∗∼10¹¹M⊙ at z = 2, hosted by haloes of mass M_(halo)∼10^(12.5)M⊙⁠) 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 (R_e), as well as the stellar mass surface density within this radius and within 1kpc (Σ_e and Σ₁, 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 Σ₁−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.

Additional Information

© 2020 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society. This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model). Accepted 2020 December 1. Received 2020 November 28; in original form 2020 September 21. Published: 05 December 2020. We thank the anonymous reviewer for detailed suggestions that helped improve the content and clarity of the paper. This work was initiated as a project for the Kavli Summer Program in Astrophysics held at the Center for Computational Astrophysics of the Flatiron Institute in 2018. The programme was co-funded by the Kavli Foundation and the Simons Foundation. We thank them for their generous support. TP would like to thank all of the participants of the summer programme including Marius Ramsoy, Ulrich Steinwandel, and Daisy Leung for stimulating discussion. TP acknowledges funding from the Future Investigators in NASA Earth and Space Science and Technology (FINESST) Fellowship, NASA grant 80NSSC19K1610. RKC acknowledges funding from the John Harvard Distinguished Science Fellowship and thanks Sandro Tacchella for helpful discussions. The Flatiron Institute is supported by the Simons Foundation. DAA was supported in part by NSF grant AST-2009687. RF acknowledges financial support from the Swiss National Science Foundation (grant number 157591). CAFG was supported by the NSF through grants AST-1517491, AST-1715216 and CAREER award AST-1652522, by NASA through grant 17-ATP17-0067, by STScI through grant HST-AR-14562.001, and by a Cottrell Scholar Award from the Research Corporation for Science Advancement. SW was supported by the CIERA Postdoctoral Fellowship Program (Center for Interdisciplinary Exploration and Research in Astrophysics, Northwestern University) and by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-2001905. The simulations were run using XSEDE (TG-AST160048, TG-AST140023), supported by NSF grant ACI-1053575, and Northwestern University's compute cluster 'Quest', as well as using NASA HEC allocations SMD-16-7561 and SMD-17-1204. Support for PFH was provided by NSF Collaborative Research Grants 1715847 and 1911233, NSF CAREER grant 1455342, NASA grants 80NSSC18K0562, JPL 1589742. Numerical calculations were run on the Caltech compute cluster 'Wheeler,' allocations FTA-Hopkins supported by the NSF and TACC, and NASA HEC SMD-16-7592. The data used in this work were, in part, hosted on facilities supported by the Scientific Computing Core at the Flatiron Institute. This research has made use of the SVO Filter Profile Service, supported from the Spanish MINECO through grant AyA2014-55216. Data Availability: The data underlying this article will be shared on reasonable request to the corresponding author. Additional data including simulation snapshots, initial conditions, and derived data products are available at https://fire.northwestern.edu/data/.

Attached Files

Published - staa3765.pdf

Accepted Version - 2009.10161.pdf

Files

2009.10161.pdf
Files (3.1 MB)
Name Size Download all
md5:be0407c3c8fb11805fa1cd6f3f9fc3ad
1.9 MB Preview Download
md5:61432d028c3a3ceb2db85647f39098e4
1.2 MB Preview Download

Additional details

Created:
August 20, 2023
Modified:
October 23, 2023