Explore the vast range of titles available.

We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10 m Keck telescope to search for Lyα emission from candidate galaxies at z ∼ 9–10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Lyα from one object in ∼8.1 hr of integration, at z = 8.665 ± 0.001 with an integrated signal-to-noise ratio > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Lyα detected at z = 8.683. The detection of Lyα suggests the existence of large (∼1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV luminosity function down to M UV = −13). The Lyα we detect would be able to escape the predominantly neutral intergalactic medium at this epoch if our detected galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 that is powering one or more ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field.

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z ≳ 6) owing to intergalactic H i absorption. To identify LyC emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O iii]/[O ii] flux ratio and star formation rate surface density, predict whether a galaxy is an LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H i column density, ionization parameter, and stellar feedback. We find that LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the Epoch of Reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with the Hubble Space Telescope Cosmic Origins Spectrograph for a z = 0.2–0.4 sample of 66 galaxies. After careful processing of the far-UV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape.

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg ii 2796 Å, 2804 Å absorption and emission lines in 29 galaxies taken from the Low-z LyC Survey to investigate the impact of (radiation and mechanical) feedback on LyC escape. Using constraints on Mg+ and photoionization models, we map the outflows’ neutral hydrogen content and predict fescLyC with a multiphase wind model. We measure mass-, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in star formation rate (SFR)-based deposition rates. We use spectral energy distribution template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high-SFR surface density and compactness. Our findings indicate that the strongest leakers are radiation feedback dominant, lack deep Mg ii absorption features, but have extended broad components in higher-ionization lines like [O iii] 5007 Å, as observed by Amorín et al. In contrast, galaxies experiencing supernovae feedback typically exhibit weaker fescLyC and show evidence of outflows in both Mg ii and higher-ionization lines. We attribute these findings to enhanced LyC escape facilitated by turbulence and cloud fragmentation in intense radiation fields, prolonged in low-metallicity environments experiencing delayed supernova feedback.

We investigate broad emission-line wings, reaching ≤800 km s−1, observed in 26 galaxies with Lyman continuum (LyC) observations, primarily from the Low-redshift Lyman Continuum Survey. Using Magellan/MIKE, Very Large Telescope/X-shooter, and WHT/ISIS high-resolution spectroscopy, we show that this fast gas appears to probe the dominant feedback mechanisms linked to LyC escape. We find that in 14 galaxies, the wings are best-fit with power laws of slope α ∼ −3.5 to −1.6, with four others best fit by Gaussians of width σBW ∼ 300 km s−1; the remaining eight show ambiguous wing morphologies. Gaussian wings are found only at low O32 = [O III]λ5007/[O II]λ3726, 3729 and high metallicity, while power-law wings span the full range of these parameters. The general evidence suggests a dual-mode paradigm for LyC escape: radiation-driven superwinds traced by power-law wings and supernova-driven feedback traced by Gaussian wings. For the former, the <3 Myr-old, pre-supernova stellar population correlates with more luminous, faster winds. The data also show that radiation-driven wind parameters like wind luminosity and power-law slope α depend on the UV luminosity more than the optically thick covering fraction, consistent with “picket-fence” radiative transfer. Observed α values flatten with both escaping LyC luminosity and higher extinction, while still preserving the anticorrelation between these two quantities. Additionally, the differential between red and blue slopes implies that extinction and dense gas are centrally concentrated relative to the wind emission. Overall, our results show that power-law emission-line wings probe LyC-driven winds and LyC escape in metal-poor starbursts.

We use the Illustris-1 simulation to explore the capabilities of the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broadband ACS, WFC3, and NIRCam data and two-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M * > 1010 M ⊙) simulated galaxies at 1 < z < 4 that evolve into a local M * > 1011 M ⊙ galaxy. In particular, we explore the potential of HST and JWST data sets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this data set for characterizing the first episodes in the SFH of local M * > 1011 M ⊙ galaxies by studying their progenitors at z > 1. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1%–50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%–30%. According to our measurements on Illustris data, we are able to recover that local M * > 1011 M ⊙ galaxies would have started their formation by z = 16, forming the first 5% of their stellar mass present at z ∼ 1 by z = 4.5, 10% by z = 3.7, and 25% by z = 2.7.

The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at z ∼ 1 and z ∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and A V values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.

We investigate broad emission-line wings, reaching \\( 800~km~s^-1\\), observed in 26 galaxies with Lyman continuum (LyC) observations, primarily from the Low-redshift Lyman Continuum Survey (LzLCS). Using Magellan/MIKE, VLT/X-shooter, and WHT/ISIS high-resolution spectroscopy, we show that this fast gas appears to probe the dominant feedback mechanisms linked to LyC escape. We find that in 14 galaxies, the wings are best fit with power laws of slope \\( -3.5 to -1.6\\), with four others best fit by Gaussians of width \\(_ BW 300~ km~s^-1\\); the remaining eight show ambiguous wing morphologies. Gaussian wings are found only at low \\(O_32\\) = \\([ O~III]5007/[O~II]3726,3729\\) and high metallicity, while power-law wings span the full range of these parameters. The general evidence suggests a dual-mode paradigm for LyC escape: radiation-driven superwinds traced by power-law wings and supernova-driven feedback traced by Gaussian wings. For the former, the \\(<3\\) Myr-old, pre-supernova stellar population correlates with more luminous, faster winds. The data also show that radiation-driven wind parameters like wind luminosity and power-law slope \\(\\) depend on the UV luminosity more than the optically thick covering fraction, consistent with ``picket-fence\" radiative transfer. Observed \\(\\) values flatten with both escaping LyC luminosity and higher extinction, while still preserving the anticorrelation between these two quantities. Additionally, the differential between red and blue slopes implies that extinction and dense gas are centrally concentrated relative to the wind emission. Overall, our results show that power-law emission-line wings probe LyC-driven winds and LyC escape in metal-poor starbursts.

The geometry of the neutral gas in and around galaxies is a key regulator of the escape of ionizing photons. We present the first statistical study aiming at linking the neutral and ionized gas distributions to the Lyman continuum (LyC) escape fraction (fesc(LyC)) in a sample of 22 confirmed LyC leakers and non-leakers at z~0.35 using the Keck Cosmic Web Imager (Keck/KCWI) and the Low Resolution Spectrograph 2 (HET/LRS2). Our integral field unit data enable the detection of neutral and low-ionization gas, as traced by Mg II, and ionized gas, as traced by [O II], extending beyond the stellar continuum for 7 and 10 objects, respectively. All but one object with extended Mg II emission also shows extended [O II] emission; in this case, Mg II emission is always more extended than [O II] by a factor 1.3 on average. Most of the galaxies with extended emission are non or weak LyC leakers (fesc(LyC) < 5%), but we find a large diversity of neutral gas configurations around these weakly LyC-emitting galaxies. Conversely, the strongest leakers (fesc(LyC) > 10%) appear uniformly compact in both Mg II and [O II] with exponential scale lengths <1 kpc. We also find a trend between fesc(LyC) and the spatial offsets of the nebular gas and the stellar continuum emission. Moreover, we find significant anti-correlations between the spatial extent of the neutral gas and the [O III]/[O II] ratio, and H\\(\\) equivalent width, as well as positive correlations with metallicity and UV size, suggesting that galaxies with more compact neutral gas sizes are more highly ionized. The observations suggest that strong LyC emitters do not have extended neutral gas halos and ionizing photons may be emitted in many directions. Combined with high ionization diagnostics, we propose the Mg II, and potentially [O II], spatial compactness are indirect indicators of LyC emitting galaxies at high-redshift.