How robustly can we constrain the low-mass end of the z ∼ 6−7 stellar mass function? The limits of lensing models and stellar population assumptions in the Hubble Frontier Fields

Abstract

We present new measurements of the very low mass end of the galaxy stellar mass function (GSMF) at z ∼ 6−7 computed from a rest-frame ultraviolet selected sample of dropout galaxies. These galaxies lie behind the six Hubble Frontier Field clusters and are all gravitationally magnified. Using deep Spitzer/IRAC and Hubble Space Telescope imaging, we derive stellar masses by fitting galaxy spectral energy distributions and explore the impact of different model assumptions and parameter degeneracies on the resulting GSMF. Our sample probes stellar masses down to $M_{\star }\gt 10^{6}\, \text{M}_{\odot}$ and we find the z ∼ 6−7 GSMF to be best parametrized by a modified Schechter function that allows for a turnover at very low masses. Using a Monte Carlo Markov chain analysis of the GSMF, including accurate treatment of lensing uncertainties, we obtain a relatively steep low-mass end slope $\alpha \simeq -1.96_{-0.08}^{+0.09}$ and a turnover at $\log (M_T/\text{M}_{\odot})\simeq 7.10_{-0.56}^{+0.17}$ with a curvature of $\beta \simeq 1.00_{-0.73}^{+0.87}$ for our minimum assumption model with constant star formation history (SFH) and low dust attenuation, AV ≤ 0.2. We find that the z ∼ 6−7 GSMF, in particular its very low mass end, is significantly affected by the assumed functional form of the star formation history and the degeneracy between stellar mass and dust attenuation. For example, the low-mass end slope ranges from $\alpha \simeq -1.82_{-0.07}^{+0.08}$ for an exponentially rising SFH to $\alpha \simeq -2.34_{-0.10}^{+0.11}$ when allowing AV of up to 3.25. Future observations at longer wavelengths and higher angular resolution with the James Webb Space Telescope are required to break these degeneracies and to robustly constrain the stellar mass of galaxies on the extreme low-mass end of the GSMF.

Implications of a Temperature-dependent Initial Mass Function. II. An Updated View of the Star-forming Main Sequence

“…M M are rare in COSMOS2015 (Furtak et al 2021), this is likely a reflection of a large population of poorly constrained galaxies being fit with these low masses, perhaps combined with selection effects.…”

Revealing Galaxy Candidates out to $z \sim 16$ with JWST Observations of the Lensing Cluster SMACS0723

“…Moreover, Spitzer Space Telescope observations allowed to probe the rest-frame optical emission of these high-redshift galaxies, which enabled us to constrain their stellar mass (e.g. Song et al 2016;Bhatawdekar et al 2019;Kikuchihara et al 2020;Stefanon et al 2021), albeit with important uncertainties due to source confusion and nebular emission contamination (Grazian et al 2015;Furtak et al 2021).…”

The effects of binary stars on galaxies and metal-enriched gas during reionization

“…The relatively small simulation volume could provide an explanation for the discrepancies seen between observations and the simulated metals, at least for some ions. The small volume results in only a subset of the true galaxy population at z ≥ 5, with TD's simulated stellar masses only reaching log M * = 8.8 at z = 6, while observations detect galaxies up to log M * = 11 (Furtak et al 2021). Thus, if strong absorption systems are associated with higher stellar mass galaxies then the dearth of strong absorbers in our simulations would be expected.…”