Dr Harley Katz

Abstract:

The Hα emission line in galaxies is a powerful tracer of their recent star formation activity. With the advent of JWST, we are now able to routinely observe Hα in galaxies at high redshift (z ≳ 3) and thus measure their star formation rates (SFRs). However, using classical SFR(Hα) calibrations to derive the SFRs leads to biased results because high-redshift galaxies are commonly characterized by low metallicities and bursty star formation histories, affecting the conversion factor between the Hα luminosity (LHα ) and the SFR. We developed a set of new SFR(Hα) calibrations that allowed us to predict the SFRs of Hα-emitters at z ≳ 3 with very little error. We used the SPHINX cosmological simulations to select a sample of star-forming galaxies representative of the Hα-emitter population observed with JWST. We then derived linear corrections to the classical SFR(Hα) calibrations that took variations in the physical properties (e.g., stellar metallicities) among individual galaxies into account. We obtained two new SFR(Hα) calibrations that compared to the classical calibrations reduce the root mean squared error (RMSE) in the predicted SFRs by ΔRMSE ≈ 0.04 dex and ΔRMSE ≈ 0.06 dex, respectively. Using the recent JWST NIRCam/grism observations of Hα-emitters at z ∼ 6, we show that the new calibrations affect the high-redshift galaxy population statistics: (i) the estimated cosmic SFR density decreases by ΔρSFR ≈ 12%, and (ii) the observed slope of the star formation main sequence increases by Δ∂ log SFR/∂ log M★ = 0.08 ± 0.02.

Abstract:

Numerous high-z galaxies have recently been observed with the James Webb Space Telescope (JWST), providing new insights into early galaxy evolution. Their physical properties are typically derived through spectral energy distribution (SED) fitting, but the reliability of this approach remains uncertain owing to limited constraints on star formation histories (SFHs) and on the contribution from emission for such early systems. Applying BAGPIPES on simulated SEDs with SFR10 > 0.3 M ⊙ yr⁻¹ at z = 6 from the SPHINX ²⁰ cosmological simulation, we examine the uncertainties related to the recovery of stellar masses, star formation rates (SFR10), and stellar metallicities from mock JWST/Near-Infrared Camera photometry, spanning F115W–F444W. Even without dust or emission lines, fitting the intrinsic stellar continuum overestimates the stellar mass by about 60%, on average (and by up to a factor of five for low-mass galaxies with recent starbursts). It also underestimates the SFR10 by a factor of 2, due to inaccurate SFHs and age–metallicity degeneracies. In full SED-fitting models that include dust attenuation and nebular emission, stellar mass estimates are primarily affected by age–metallicity degeneracy and emission lines. Short-term SFRs are most sensitive to dust attenuation and nebular emission, while long-term SFRs additionally depend on the assumed SFHs. Incorporating bands that are free of strong emission lines, such as F410M, helps mitigate stellar mass overestimation by disentangling line emission from older stellar populations. We also find that best fit or likelihood-weighted estimates are generally more accurate than median posterior values. Although stellar mass functions are reproduced reasonably well (particularly when the minimum-χ ² estimates are used), the slope of the main sequence of star formation acutely depends on the adopted fitting model. Overall, these results underscore the importance of careful modelling when interpreting high-z photometry, particularly for galaxies with recent star formation burst and/or strong emission lines, to minimise systematic biases in derived physical properties.

Abstract:

Cosmic ray (CR) feedback has been proposed as a powerful mechanism for driving warm gas outflows in galaxies. We use cosmological magnetohydrodynamic simulations to investigate the impact of CR feedback on neutral hydrogen (H i) in a 1011 M⊙ dark matter halo at 2 < z < 4. To this end, we postprocess the simulations with ionizing radiative transfer and perform Monte Carlo Lyman-α (Lyα) transfer calculations. CR feedback reduces H i column densities around young stars, thereby allowing more Lyα photons to escape and consequently offering a better match to the Lyα luminosities of observed Lyα emitters. Although galaxies with CR-driven outflows have more extended H i in the circumgalactic medium, two Lyα line properties sensitive to optical depth and gas kinematic—the location of the red peak relative to the Lyα line center in velocity space (vred) and relative strength of the blue-to-red peaks (B/R)—cannot distinguish between the CR-driven and non-CR simulations. This is because Lyα photons propagate preferentially along low H i density channels created by the ionizing radiation, thereby limiting the scattering with volume-filling H i. In contrast, the observed low flux ratios between the valley and peak and the surface brightness profiles are better reproduced in the model with CR-driven outflows because the Lyα photons interact more before escaping, rather than being destroyed by dust as is the case in the non-CR simulation. We discuss the potential cause of the paucity of sightlines in simulations that exhibit prominent red peaks and large vred, which may require the presence of more volume-filling H i.