An official website of the United States government
ABSTRACT Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($$f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$$), helium ionization ($$f_\text{He} \approx 8.7\,\rm{per\,\,cent}$$), and anisotropic escape fractions ($$f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$$), as these reduce the available budget for hydrogen line emission ($$f_\text{H} \approx 55.9\,\rm{per\,\,cent}$$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.
more »
« less
H α emission in local galaxies: star formation, time variability, and the diffuse ionized gas
ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$$10{{\ \rm per\ cent}}$$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $$\sim 40{{\ \rm per\ cent}}$$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $$f_{\rm abs}\approx 28{{\ \rm per\ cent}}$$ and $$f_{\rm He}\approx 9{{\ \rm per\ cent}}$$, respectively. Together with an escape fraction of $$f_{\rm esc}\approx 6{{\ \rm per\ cent}}$$, this reduces the available budget for hydrogen line emission by nearly half ($$f_{\rm H}\approx 57{{\ \rm per\ cent}}$$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.
more »
« less
- PAR ID:
- 10366765
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 513
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2904-2929
- Size(s):
- p. 2904-2929
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We measure the Lyman continuum (LyC) escape fraction in 54 faint Lyman-alpha emitters (LAEs) at $$z$$ ≃ 3.1 in the GOODS-South field. With the average magnitude of R = 26.7 AB (MUV = −18.8 and L ≃ 0.1L*), these galaxies represent a population of compact young dwarf galaxies. Their properties are likely to resemble those in the galaxies responsible for reionizing the Universe at $$z$$ > 6. We do not detect LyC emission in any individual LAEs in the deep HST F336W images, which covers the rest-frame 820 Å. We do not detect the LyC emission of these LAEs in the stacked F336W images either. The 3σ upper limit of LyC escape fractions is $$f_{\rm esc}\lt 14\!-\!32{{\ \rm per\ cent}}$$. However, the high Ly α rest-frame equivalent width (EW), low stellar mass, and UV luminosity of these LAEs suggest that they should have $$f_{\rm esc}\gt 50{{\ \rm per\ cent}}$$. The low LyC escape fraction from this work and other stacking analyses suggests that the LyC-leaking galaxies with $$f_{\rm esc}\gt 50{{\ \rm per\ cent}}$$ at $$z$$ = 2–3 do not follow the relation between fesc and UV luminosity and Ly α EW derived from typical galaxies at similar redshifts. Therefore, the UV luminosity and Ly α EW are not the best indicators for the LyC escape fraction.more » « less
-
We use a well-motivated galaxy formation framework to predict stellar masses, star formation rates (SFR), and ultraviolet (UV) luminosities of galaxy populations at redshifts $$z\in 5-16$$, taking into account stochasticity of SFR in a controlled manner. We demonstrate that the model can match observational estimates of UV luminosity functions (LFs) at $5<10$ with a modest level of SFR stochasticity, resulting in the scatter of absolute UV luminosity at a given halo mass of $$\sigma_{M_{\rm UV}}\approx 0.75$$. To match the observed UV LFs at $$z\approx 11-13$$ and $$z\approx 16$$ the SFR stochasticity should increase so that $$\sigma_{M_{\rm UV}}\approx 1-1.3$$ and $$\approx 2$$, respectively. Model galaxies at $$z\approx 11-13$$ have stellar masses and SFRs in good agreement with existing measurements. The median fraction of the baryon budget that was converted into stars, $$f_\star$$, is only $$f_\star\approx 0.005-0.05$$, but a small fraction of galaxies at $z=16$ have $$f_\star>1$$ indicating that SFR stochasticity cannot be higher. We discuss several testable consequences of the increased SFR stochasticity at $z>10$. The increase of SFR stochasticity with increasing $$z$$, for example, prevents steepening of UV LF and even results in some flattening of UV LF at $$z\gtrsim 13$$. The median stellar ages of model galaxies at $$z\approx 11-16$$ are predicted to decrease from $$\approx 20-30$$ Myr for $$M_{\rm UV}\gtrsim -21$$ galaxies to $$\approx 5-10$$ Myr for brighter ones. Likewise, the scatter in median stellar age is predicted to decrease with increasing luminosity. The scatter in the ratio of star formation rates averaged over 10 and 100 Myr should increase with redshift. Fluctuations of ionizing flux should increase at $z>10$ resulting in the increasing scatter in the line fluxes and their ratios for the lines sensitive to ionization parameter.more » « less
-
ABSTRACT The connection between the escape fraction of ionizing radiation (fesc) and the properties of galaxies, such as stellar mass ($$\rm M_{\rm *}$$), age, star-formation rate (SFR), and dust content, are key inputs for reionization models, but many of these relationships remain untested at high redshift. We present an analysis of a sample of 96 $$z$$ ∼ 3 galaxies from the Keck Lyman Continuum Spectroscopic Survey (KLCS). These galaxies have both sensitive Keck/LRIS spectroscopic measurements of the Lyman continuum (LyC) region, and multiband photometry that places constraints on stellar population parameters. We construct composite spectra from subsamples binned as a function of galaxy property and quantify the ionizing-photon escape for each composite. We find a significant anti-correlation between fesc and $$\rm M_{\rm *}$$, consistent with predictions from cosmological zoom-in simulations. We also find significant anti-correlation between fesc and E(B−V), encoding the underlying physics of LyC escape in our sample. We also find no significant correlation between fesc and either stellar age or specific SFR (= SFR/$$\rm M_{\rm *}$$), challenging interpretations that synchronize recent star formation and favorable conditions for ionizing escape. The galaxy properties now shown to correlate with fesc in the KLCS are Lyα equivalent width, UV Luminosity, $$\rm M_{\rm *}$$, SFR, and E(B−V), but not age or sSFR. This comprehensive analysis of galaxy properties and LyC escape at high redshift will be used to guide future models and observations of the reionization epoch.more » « less
-
ABSTRACT We present 10 main-sequence ALPINE galaxies (log (M/M⊙) = 9.2−11.1 and $${\rm SFR}=23-190\, {\rm M_{\odot }\, yr^{-1}}$$) at z ∼ 4.5 with optical [O ii] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multiwavelength sample at these redshifts, combining various measurements in the ultraviolet, optical, and far-infrared including [C ii]158 $$\mu$$m line emission and dust continuum from ALMA and H α emission from Spitzer photometry. For the first time, this unique sample allows us to analyse the relation between [O ii] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [O ii]/[C ii] and [O ii]/H α luminosity ratios at z ∼ 4.5. The [O ii]−SFR relation at z ∼ 4.5 cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $$50{{\ \rm per\ cent}}$$ solar. To explain the measured dust-corrected luminosity ratios of $$\log (L_{\rm [OII]}/L_{\rm [CII]}) \sim 0.98^{+0.21}_{-0.22}$$ and $$\log (L_{\rm [OII]}/L_{\rm H\alpha }) \sim -0.22^{+0.13}_{-0.15}$$ for our sample, ionization parameters log (U) < −2 and electron densities $$\log (\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$$ are required. The former is consistent with galaxies at z ∼ 2−3, however lower than at z > 6. The latter may be slightly higher than expected given the galaxies’ specific SFR. The analysis of this pilot sample suggests that typical log (M/M⊙) > 9 galaxies at z ∼ 4.5 to have broadly similar ISM properties as their descendants at z ∼ 2 and suggest a strong evolution of ISM properties since the epoch of reionization at z > 6.more » « less
