High-z stellar masses can be recovered robustly with JWST photometry

TL;DR

This study demonstrates that JWST photometry allows for robust recovery of high-redshift galaxy stellar masses within a factor of three, despite some systematic biases caused by emission lines affecting broadband measurements.

Contribution

It provides a validation of SED fitting methods for high-z galaxy mass estimation using JWST data, highlighting systematic trends and their effects on the stellar mass function.

Findings

Stellar masses are generally recovered within a factor of 3.

Systematic overestimation for low-mass galaxies and underestimation for high-mass galaxies.

Emission lines significantly impact broadband photometry and mass estimates.

Abstract

Robust inference of galaxy stellar masses from photometry is crucial for constraints on galaxy assembly across cosmic time. Here, we test a commonly-used Spectral Energy Distribution (SED) fitting code, using simulated galaxies from the SPHINX20 cosmological radiation hydrodynamics simulation, with JWST NIRCam photometry forward-modelled with radiative transfer. Fitting the synthetic photometry with various star formation history models, we show that recovered stellar masses are, encouragingly, generally robust to within a factor of ~3 for galaxies in the range M*~10^7-10^9M_sol at z=5-10. These results are in stark contrast to recent work claiming that stellar masses can be underestimated by as much as an order of magnitude in these mass and redshift ranges. However, while >90% of masses are recovered to within 0.5dex, there are notable systematic trends, with stellar masses typically overestimated for low-mass galaxies (M*<~10^8M_sol) and slightly underestimated for high-mass galaxies (M*>~10^9M_sol). We demonstrate that these trends arise due to the SED fitting code poorly modelling the impact of strong emission lines on broadband photometry. These systematic trends, which exist for all star formation history parametrisations tested, have a tilting effect on the inferred stellar mass function, with number densities of massive galaxies underestimated (particularly at the lowest redshifts studied) and number densities of lower-mass galaxies typically overestimated. Overall, this work suggests that we should be optimistic about our ability to infer the masses of high-z galaxies observed with JWST (notwithstanding contamination from AGN) but careful when modelling the impact of strong emission lines on broadband photometry.