Showing 1–50 of 362 results for author: Ferguson, H C

Multivariate Predictors of LyC Escape II: Predicting LyC Escape Fractions for High-Redshift Galaxies

Authors: Anne E. Jaskot, Anneliese C. Silveyra, Anna Plantinga, Sophia R. Flury, Matthew Hayes, John Chisholm, Timothy Heckman, Laura Pentericci, Daniel Schaerer, Maxime Trebitsch, Anne Verhamme, Cody Carr, Henry C. Ferguson, Zhiyuan Ji, Mauro Giavalisco, Alaina Henry, Rui Marques-Chaves, Göran Östlin, Alberto Saldana-Lopez, Claudia Scarlata, Gábor Worseck, Xinfeng Xu

Abstract: JWST is uncovering the properties of ever increasing numbers of galaxies at z>6, during the epoch of reionization. Connecting these observed populations to the process of reionization requires understanding how efficiently they produce Lyman continuum (LyC) photons and what fraction (fesc) of these photons escape into the intergalactic medium. By applying the Cox proportional hazards model, a surv… ▽ More JWST is uncovering the properties of ever increasing numbers of galaxies at z>6, during the epoch of reionization. Connecting these observed populations to the process of reionization requires understanding how efficiently they produce Lyman continuum (LyC) photons and what fraction (fesc) of these photons escape into the intergalactic medium. By applying the Cox proportional hazards model, a survival analysis technique, to the Low-redshift Lyman Continuum Survey (LzLCS), we develop new, empirical, multivariate predictions for fesc. The models developed from the LzLCS reproduce the observed fesc for z~3 samples, which suggests that LyC emitters may share similar properties at low and high redshift. Our best-performing models for the z~3 galaxies include information about dust attenuation, ionization, and/or morphology. We then apply these models to z$\gtrsim$6 galaxies. For large photometric samples, we find a median predicted fesc=0.047-0.14. For smaller spectroscopic samples, which may include stronger emission line galaxies, we find that $\geq$33% of the galaxies have fesc >0.2, and we identify several candidate extreme leakers with fesc $\geq$0.5. The current samples show no strong trend between predicted fesc and UV magnitude, but limited spectroscopic information makes this result uncertain. Multivariate predictions can give significantly different results from single variable predictions, and the predicted fesc for high-redshift galaxies can differ significantly depending on whether star formation rate surface density or radius is used as a measure of galaxy morphology. We provide all parameters necessary to predict fesc for additional samples of high-redshift galaxies using these models. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 33 pages, 9 figures, 10 tables, plus appendix

Multivariate Predictors of LyC Escape I: A Survival Analysis of the Low-redshift Lyman Continuum Survey

Authors: Anne E. Jaskot, Anneliese C. Silveyra, Anna Plantinga, Sophia R. Flury, Matthew Hayes, John Chisholm, Timothy Heckman, Laura Pentericci, Daniel Schaerer, Maxime Trebitsch, Anne Verhamme, Cody Carr, Henry C. Ferguson, Zhiyuan Ji, Mauro Giavalisco, Alaina Henry, Rui Marques-Chaves, Göran Östlin, Alberto Saldana-Lopez, Claudia Scarlata, Gábor Worseck, Xinfeng Xu

Abstract: To understand how galaxies reionized the universe, we must determine how the escape fraction of Lyman Continuum (LyC) photons (fesc) depends on galaxy properties. Using the z~0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of fesc. These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detec… ▽ More To understand how galaxies reionized the universe, we must determine how the escape fraction of Lyman Continuum (LyC) photons (fesc) depends on galaxy properties. Using the z~0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of fesc. These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detections and upper limits. Our best model predicts the LzLCS fesc detections with a root-mean-square (RMS) scatter of 0.31 dex, better than single-variable correlations. According to ranking techniques, the most important predictors of fesc are the equivalent width (EW) of Lyman-series absorption lines and the UV dust attenuation, which track line-of-sight absorption due to HI and dust. The HI absorption EW is uniquely crucial for predicting fesc for the strongest LyC emitters, which show properties similar to weaker LyC emitters and whose high fesc may therefore result from favorable orientation. In the absence of HI information, star formation rate surface density ($Σ_{\rm SFR}$) and [O III]/[O II] ratio are the most predictive variables and highlight the connection between feedback and fesc. We generate a model suitable for z>6, which uses only the UV slope, $Σ_{\rm SFR}$, and [O III]/[O II]. We find that $Σ_{\rm SFR}$ is more important in predicting fesc at higher stellar masses, whereas [O III]/[O II] plays a greater role at lower masses. We also analyze predictions for other parameters, such as the ionizing-to-non ionizing flux ratio and Ly=alpha escape fraction. These multivariate models represent a promising tool for predicting fesc at high redshift. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 34 pages + appendix, 12 figures

The Rise of Faint, Red AGN at $z>4$: A Sample of Little Red Dots in the JWST Extragalactic Legacy Fields

Authors: Dale D. Kocevski, Steven L. Finkelstein, Guillermo Barro, Anthony J. Taylor, Antonello Calabrò, Brivael Laloux, Johannes Buchner, Jonathan R. Trump, Gene C. K. Leung, Guang Yang, Mark Dickinson, Pablo G. Pérez-González, Fabio Pacucci, Kohei Inayoshi, Rachel S. Somerville, Elizabeth J. McGrath, Hollis B. Akins, Micaela B. Bagley, Laura Bisigello, Rebecca A. A. Bowler, Adam Carnall, Caitlin M. Casey, Yingjie Cheng, Nikko J. Cleri, Luca Costantin , et al. (32 additional authors not shown)

Abstract: We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifti… ▽ More We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN. △ Less

Submitted 19 April, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: 23 pages, 17 figures, submitted to ApJ

Constraints on the Lyman Continuum Escape from Low-mass Lensed Galaxies at 1.3 $\leq$ z $\leq$ 3.0

Authors: Intae Jung, Henry C. Ferguson, Matthew J. Hayes, Alaina Henry, Anne E. Jaskot, Daniel Schaerer, Keren Sharon, Ricardo O. Amorín, Hakim Atek, Matthew B. Bayliss, Håkon Dahle, Steven L. Finkelstein, Andrea Grazian, Lucia Guaita, Göran Östlin, Laura Pentericci, Swara Ravindranath, Claudia Scarlata, Harry I. Teplitz, Anne Verhamme

Abstract: Low-mass galaxies can significantly contribute to reionization due to their potentially high Lyman continuum (LyC) escape fraction and relatively high space density. We present a constraint on the LyC escape fraction from low-mass galaxies at z = 1.3 - 3.0. We obtained rest-frame UV continuum imaging with the ACS/SBC and the WFC3/UVIS from the Hubble Space Telescope for eight strongly-lensed galax… ▽ More Low-mass galaxies can significantly contribute to reionization due to their potentially high Lyman continuum (LyC) escape fraction and relatively high space density. We present a constraint on the LyC escape fraction from low-mass galaxies at z = 1.3 - 3.0. We obtained rest-frame UV continuum imaging with the ACS/SBC and the WFC3/UVIS from the Hubble Space Telescope for eight strongly-lensed galaxies that were identified in the Sloan Giant Arc Survey (SGAS) and the Cluster Lensing And Supernova survey with Hubble (CLASH). The targeted galaxies were selected to be spectroscopically confirmed, highly magnified, and blue in their UV spectral shapes ($β<-1.7$). Our targets include intrinsically low luminosity galaxies down to a magnification-corrected absolute UV magnitude of $M_{\rm UV}\sim-14$. We perform custom-defined aperture photometry to place the most reliable upper limits of LyC escape from our sample. From our observations, we report no significant ($>$$2σ$) detections of LyC fluxes, placing 1$σ$ upper limits on the absolute LyC escape fractions of 3 - 15%. Our observations do not support the expected increased escape fractions of LyC photons from intrinsically UV faint sources. Considering the highly anisotropic geometry of LyC escape, increasing the sample size of faint galaxies in future LyC observations is crucial. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ (16 pages, 6 figures, 5 tables)

Linking Mg II and [O II] spatial distribution to ionizing photon escape in confirmed LyC leakers and non-leakers

Authors: Floriane Leclercq, John Chisholm, Wichahpi King, Greg Zeimann, Anne E. Jaskot, Alaina Henry, Matthew Hayes, Sophia R. Flury, Yuri Izotov, Xavier J. Prochaska, Anne Verhamme, Ricardo O. Amorín, Hakim Atek, Omkar Bait, Jérémy Blaizot, Cody Carr, Zhiyuan Ji, Alexandra Le Reste, Harry C. Ferguson, Simon Gazagnes, Timothy Heckman, Lena Komarova, Rui Marques-Chaves, Göran Östlin, Alberto Saldana-Lopez , et al. (7 additional authors not shown)

Abstract: 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 S… ▽ More 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. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 31 pages, 24 figures, 4 tables, submitted to A&A

Journal ref: A&A 687, A73 (2024)

The Next Generation Deep Extragalactic Exploratory Public Near-Infrared Slitless Survey Epoch 1 (NGDEEP-NISS1): Extra-Galactic Star-formation and Active Galactic Nuclei at 0.5 < z < 3.6

Authors: Nor Pirzkal, Barry Rothberg, Casey Papovich, Lu Shen, Gene C. K. Leung, Micaela B. Bagley, Steven L. Finkelstein, Brittany N. Vanderhoof, Jennifer M. Lotz, Anton M. Koekemoer, Nimish P. Hathi, Yingjie Cheng, Nikko J. Cleri, Norman A. Grogin, L. Y. Aaron Yung, Mark Dickinson, Henry C. Ferguson, Jonathan P. Gardner, Intae Jung, Jeyhan S. Kartaltepe, Russell Ryan, Raymond C. Simons, Swara Ravindranath, Danielle A. Berg, Bren E. Backhaus , et al. (26 additional authors not shown)

Abstract: The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the firs… ▽ More The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$σ$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$β$), and dust content, of galaxies at $1<z<3.5$. The latter showing intriguing initial results on the applicability and assumptions made regarding the use of Case B recombination. Further, we identify the presence of active galactic nuclei (AGN) and infer the mass of their supermassive black holes (SMBHs) using broadened restframe MgII and H$β$ emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution (SED) models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts. △ Less

Submitted 20 April, 2024; v1 submitted 15 December, 2023; originally announced December 2023.

Comments: 33 pages, 28 Figures, Accepted (ApJ)

Rest-Frame UV Colors for Faint Galaxies at $z \sim 9-16$ with the \textit{JWST} NGDEEP Survey

Authors: Alexa M. Morales, Steven L. Finkelstein, Gene C. K. Leung, Micaela B. Bagley, Nikko J. Cleri, Romeel Dave, Mark Dickinson, Henry C. Ferguson, Nimish P. Hathi, Ewan Jones, Anton M. Koekemoer, Casey Papovich, Pablo G. Perez-Gonzalez, Nor Pirzkal, Britton Smith, Stephen M. Wilkins, L. Y. Aaron Yung

Abstract: We present measurements of the rest-frame UV spectral slope, $β$, for a sample of 36 faint star-forming galaxies at z ~ 9-16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey. We use robust photometric measurements for UV-faint galaxies (down to $M_{UV}$ ~ -16), originally published in Leung+23, and measure value… ▽ More We present measurements of the rest-frame UV spectral slope, $β$, for a sample of 36 faint star-forming galaxies at z ~ 9-16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey. We use robust photometric measurements for UV-faint galaxies (down to $M_{UV}$ ~ -16), originally published in Leung+23, and measure values of the UV spectral slope via photometric power-law fitting to both the observed photometry and to stellar population models obtained through spectral energy distribution (SED) fitting with Bagpipes. We obtain a median and 68% confidence interval for $β$ from photometric power-law fitting of $β_{PL} = -2.7^{+0.5}_{-0.5}$ and from SED-fitting, $β_{SED} = -2.3^{+0.2}_{-0.1}$ for the full sample. We show that when only 2-3 photometric detections are available, SED-fitting has a lower scatter and reduced biases than photometric power-law fitting. We quantify this bias and find that after correction, the median $β_{SED,corr} = -2.5^{+0.2}_{-0.2}$. We measure physical properties for our galaxies with Bagpipes and find that our faint ($M_{UV} = -18.1^{+0.7}_{-0.9}$) sample is low mass (${log}[M_{\ast}/M_\odot] = 7.7^{+0.5}_{-0.5}$), fairly dust-poor ($A_{v} = 0.1^{+0.2}_{-0.1}$ mag), and modestly young (${log[age]} = 7.8^{+0.2}_{-0.8}$ yr) with a median star formation rate of $\mathrm{log(SFR)} = -0.3^{+0.4}_{-0.4} M_\odot{/yr}$. We find no strong evidence for ultra-blue UV spectral slopes ($β$ ~ -3) within our sample, as would be expected for exotically metal-poor ($Z/Z_{\odot}$ < 10$^{-3}$) stellar populations with very high LyC escape fractions. Our observations are consistent with model predictions that galaxies of these stellar masses at z~9-16 should have only modestly low metallicities ($Z/Z_{\odot}$ ~ 0.1--0.2). △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 21 pages, 8 figures, 3 tables, submitted to ApjL

The Complete CEERS Early Universe Galaxy Sample: A Surprisingly Slow Evolution of the Space Density of Bright Galaxies at z ~ 8.5-14.5

Authors: Steven L. Finkelstein, Gene C. K. Leung, Micaela B. Bagley, Mark Dickinson, Henry C. Ferguson, Casey Papovich, Hollis B. Akins, Pablo Arrabal Haro, Romeel Dave, Avishai Dekel, Jeyhan S. Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Norbert Pirzkal, Rachel S. Somerville, L. Y. Aaron Yung, Ricardo Amorin, Bren E. Backhaus, Peter Behroozi, Laura Bisigello, Volker Bromm, Caitlin M. Casey, Oscar A. Chavez Ortiz, Yingjie Cheng, Katherine Chworowsky , et al. (30 additional authors not shown)

Abstract: We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than p… ▽ More We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z~9, 11 and 14, and show that the space density of bright (M_UV=-20) galaxies changes only modestly from z~14 to z~9, compared to a steeper increase from z~8 to z~4. While our candidates are photometrically selected, spectroscopic followup has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z>10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models which naturally have enhanced star-formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic followup of a large sample of early galaxies can distinguish between these competing scenarios. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Submitted to ApJL. Main paper is 33 pages, 13 figures, 5 tables. Two appendices with additional figures and tables

Low-redshift Lyman Continuum Survey: Radio continuum properties of low-$z$ Lyman continuum emitters

Authors: Omkar Bait, Sanchayeeta Borthakur, Daniel Schaerer, Emmanuel Momjian, Biny Sebastian, Alberto Saldana-Lopez, Sophia R. Flury, John Chisholm, Rui Marques-Chaves, Anne E. Jaskot, Harry C. Ferguson, Gabor Worseck, Zhiyuan Ji, Lena Komarova, Maxime Trebitsch, Matthew J. Hayes, Laura Pentericci, Goran Ostlin, Trinh Thuan, Ricardo O. Amorín, Bingjie Wang, Xinfeng Xu, Mark T. Sargent

Abstract: Sources that leak Lyman-continuum (LyC) photons and lead to the reionisation of the universe are intensely studied using multiple observing facilities. Recently, the Low-redshift LyC Survey (LzLCS) has found the first large sample of LyC emitting galaxies at low redshift ($z\sim 0.3$) with the Hubble Space Telescope/Cosmic Origins Spectrograph. The LzLCS sample contains a robust estimate of the Ly… ▽ More Sources that leak Lyman-continuum (LyC) photons and lead to the reionisation of the universe are intensely studied using multiple observing facilities. Recently, the Low-redshift LyC Survey (LzLCS) has found the first large sample of LyC emitting galaxies at low redshift ($z\sim 0.3$) with the Hubble Space Telescope/Cosmic Origins Spectrograph. The LzLCS sample contains a robust estimate of the LyC escape fraction ($f_\mathrm{esc}^\mathrm{LyC}$) for 66 galaxies spanning a wide range of $f_\mathrm{esc}^\mathrm{LyC}$. Here we, for the first time, aim to study the radio continuum (RC) properties of LzLCS sources and their dependence on $f_\mathrm{esc}^\mathrm{LyC}$. We present Karl G. Jansky Very Large Array RC observations at C (4-8 GHz), S (2-4 GHz) and L (1-2 GHz) bands for a sub-sample of the LzLCS sources. The radio spectral index ($α^{\mathrm{3GHz}}_\mathrm{6GHz}$) spans a wide range from being flat ( $\geq -0.1$) to very steep ($\leq -1.0$). The strongest leakers in our sample show flat $α^{\mathrm{3GHz}}_\mathrm{6GHz}$, weak leakers have $α^{\mathrm{3GHz}}_\mathrm{6GHz}$ close to normal star-forming galaxies, and non-leakers are characterized by steep $α^{\mathrm{3GHz}}_\mathrm{6GHz}$. We argue that a combination of young ages, free-free absorption, and a flat cosmic-ray energy spectrum can altogether lead to a flat $α^{\mathrm{3GHz}}_\mathrm{6GHz}$ for strong leakers. Non-leakers are characterized by steep spectra which can arise due to break/cutoff at high frequencies. Such a cutoff in the spectrum can arise in a single injection model of CRs characteristic of galaxies which have recently stopped star formation. Such a relation between $α^{\mathrm{3GHz}}_\mathrm{6GHz}$ and $f_\mathrm{esc}^\mathrm{LyC}$ hints at the interesting role of supernovae, CRs, and magnetic fields in facilitating the escape ( and/or the lack) of LyC photons. (Abridged) △ Less

Submitted 19 June, 2024; v1 submitted 28 October, 2023; originally announced October 2023.

Comments: 20 pages, 17 figures, 4 tables, 8 appendix pages, Astronomy & Astrophysics in press

Galaxies Going Bananas: Inferring the 3D Geometry of High-Redshift Galaxies with JWST-CEERS

Authors: Viraj Pandya, Haowen Zhang, Marc Huertas-Company, Kartheik G. Iyer, Elizabeth McGrath, Guillermo Barro, Steven L. Finkelstein, Martin Kuemmel, William G. Hartley, Henry C. Ferguson, Jeyhan S. Kartaltepe, Joel Primack, Avishai Dekel, Sandra M. Faber, David C. Koo, Greg L. Bryan, Rachel S. Somerville, Ricardo O. Amorin, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Emmanuel Bertin, Luca Costantin, Romeel Dave, Mark Dickinson , et al. (31 additional authors not shown)

Abstract: The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean e… ▽ More The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications. △ Less

Submitted 15 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

CEERS Key Paper VII: JWST/MIRI Reveals a Faint Population of Galaxies at Cosmic Noon Unseen by Spitzer

Authors: Allison Kirkpatrick, Guang Yang, Aurelien Le Bail, Greg Troiani, Eric F. Bell, Nikko J. Cleri, David Elbaz, Steven L. Finkelstein, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Dale D. Kocevski, Ray A. Lucas, Jed McKinney, Casey Papovich, Pablo G. Perez-Gonzalez, Alexander de la Vega, Micaela B. Bagley, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Adriano Fontana, Andrea Grazian, Norman A. Grogin, Pablo Arrabal Haro , et al. (11 additional authors not shown)

Abstract: The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24$μ$m po… ▽ More The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24$μ$m population in the EGS field. We find that MIRI can observe an order of magnitude deeper than MIPS in significantly shorter integration times, attributable to JWST's much larger aperture and MIRI's improved sensitivity. MIRI is exceptionally good at finding faint ($L_{\rm IR}<10^{10} L_\odot$) galaxies at $z\sim1-2$. We find that a significant portion of MIRI galaxies are "mid-IR weak"--they have strong near-IR emission and relatively weaker mid-IR emission, and most of the star formation is unobscured. We present new IR templates that capture how the mid-IR to near-IR emission changes with increasing infrared luminosity. We present two color-color diagrams to separate mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find that these color diagrams are most effective when used in conjunction with each other. We present the first number counts of 10$μ$m sources and find that there are $\lesssim10$ IR AGN per MIRI pointing, possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to low metallicities or low dust content). We conclude that MIRI is most effective at observing moderate luminosity ($L_{\rm IR}=10^9-10^{10}L_\odot$) galaxies at $z=1-2$, and that photometry alone is not effective at identifying AGN within this faint population. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 21 pages, 10 figures. Resubmitted to ApJS after revision

CEERS Key Paper VIII: Emission Line Ratios from NIRSpec and NIRCam Wide-Field Slitless Spectroscopy at z>2

Authors: Bren E. Backhaus, Jonathan R. Trump, Nor Pirzkal, Guillermo Barro, Steven L. Finkelstein, Pablo Arrabal Haro, Raymond C. Simons, Jessica Wessner, Nikko J. Cleri, Michaela Hirschmann, Micaela B. Bagley, David C. Nicholls, Mark Dickinson, Jeyhan S. Kartaltepe, Casey Papovich, Dale D. Kocevski, Anton M. Koekemoer, Laura Bisigello, Anne E. Jaskot, Ray A. Lucas, Intae Jung, Stephen M. Wilkins, L. Y. Aaron Yung, Henry C. Ferguson, Adriano Fontana , et al. (15 additional authors not shown)

Abstract: We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey (CEERS) to measure rest-frame optical emission-line of 155 galaxies at z>2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks.… ▽ More We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey (CEERS) to measure rest-frame optical emission-line of 155 galaxies at z>2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the Ha, [OIII]/Hb, and [NeIII]/[OII] emission lines in terms of redshift by comparing to lower redshift SDSS and CLEAR samples. We find a significant (>3$σ$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$σ$) correlation with redshift. We compare [OIII]/Hb and [NeIII]/[OII] to stellar mass and Hb SFR. We find that both emission-line ratios have a correlation with Hb SFR and an anti-correlation with stellar mass across the redshifts 0<z<9. Comparison with MAPPINGS~V models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We additionally compare to IllustriousTNG predictions and find that they effectively describe the highest [OIII]/Hb ratios observed in our sample, without the need to invoke MAPPINGS models with significant shock ionizionation components. △ Less

Submitted 7 September, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 16 pages, 11 figures

NGDEEP Epoch 1: The Faint-End of the Luminosity Function at $z \sim$ 9-12 from Ultra-Deep JWST Imaging

Authors: Gene C. K. Leung, Micaela B. Bagley, Steven L. Finkelstein, Henry C. Ferguson, Anton M. Koekemoer, Pablo G. Perez-Gonzalez, Alexa Morales, Dale D. Kocevski, Guang Yang, Rachel S. Somerville, Stephen M. Wilkins, L. Y. Aaron Yung, Seiji Fujimoto, Rebecca L. Larson, Casey Papovich, Nor Pirzkal, Danielle A. Berg, Jennifer M. Lotz, Marco Castellano, Oscar A. Chavez Ortiz, Yingjie Cheng, Mark Dickinson, Mauro Giavalisco, Nimish P. Hathi, Taylor A. Hutchison , et al. (4 additional authors not shown)

Abstract: We present a robust sample of very high-redshift galaxy candidates from the first epoch of {\it JWST}/NIRCam imaging from the Next Generation Extragalactic Exploratory Deep (NGDEEP) Survey. The NGDEEP NIRCam imaging in the Hubble Ultra Deep Field Parallel Field 2 (HUDF-Par2) reaches $m=30.4$ (5$σ$, point-source) in F277W, making it the deepest public {\it JWST} GO imaging dataset to date. We descr… ▽ More We present a robust sample of very high-redshift galaxy candidates from the first epoch of {\it JWST}/NIRCam imaging from the Next Generation Extragalactic Exploratory Deep (NGDEEP) Survey. The NGDEEP NIRCam imaging in the Hubble Ultra Deep Field Parallel Field 2 (HUDF-Par2) reaches $m=30.4$ (5$σ$, point-source) in F277W, making it the deepest public {\it JWST} GO imaging dataset to date. We describe our detailed data reduction process of the six-filter broad-band {\it JWST}/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 $z \gtrsim 9$ galaxy candidates. These objects span a redshift range of $z=8.5-15.8$, and apparent magnitudes of $m_\mathrm{F277W} = 27-30.5$ AB mag, reaching $\sim 1.5$ mag deeper than previous public {\it JWST} imaging surveys. We calculate the rest-frame ultraviolet (UV) luminosity function at $z \sim$ 9 and 11, and present a new measurement of the luminosity function faint-end slope at $z \sim 11$. There is no significant evolution in the faint-end slope and number density from $z=9$ to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs. △ Less

Submitted 9 June, 2023; originally announced June 2023.

Comments: 21 pages, 10 figures. Submitted to ApJL

JWST CEERS probes the role of stellar mass and morphology in obscuring galaxies

Authors: Carlos Gómez-Guijarro, Benjamin Magnelli, David Elbaz, Stijn Wuyts, Emanuele Daddi, Aurélien Le Bail, Mauro Giavalisco, Mark Dickinson, Pablo G. Pérez-González, Pablo Arrabal Haro, Micaela B. Bagley, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Yingjie Cheng, Laure Ciesla, Avishai Dekel, Henry C. Ferguson, Steven L. Finkelstein, Maximilien Franco, Norman A. Grogin, Benne W. Holwerda, Shuowen Jin , et al. (16 additional authors not shown)

Abstract: In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JW… ▽ More In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs. △ Less

Submitted 4 September, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: Accepted by A&A. 13 pages, 9 figures

Journal ref: A&A 677, A34 (2023)

CEERS: Diversity of Lyman-Alpha Emitters during the Epoch of Reionization

Authors: Intae Jung, Steven L. Finkelstein, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Rebecca L. Larson, Raymond C. Simons, Casey Papovich, Hyunbae Park, Laura Pentericci, Jonathan R. Trump, Ricardo O. Amorin, Bren E. Backhaus, Caitlin M. Casey, Yingjie Cheng, Nikko J. Cleri, M. C. Cooper, Olivia R. Cooper, Jonathan P. Gardner, Eric Gawiser, Andrea Grazian, Nimish P. Hathi, Michaela Hirschmann , et al. (7 additional authors not shown)

Abstract: We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$α$-emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) survey. From NIRSpec data, we confirm the systemic redshifts of these Ly$α$ emitters, and emission-line ratio diagnostics indicate… ▽ More We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$α$-emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) survey. From NIRSpec data, we confirm the systemic redshifts of these Ly$α$ emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal poor. We investigate Ly$α$ line properties, including the line flux, velocity offset, and spatial extension. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements ($\sim5\times$ greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly$α$ emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly$α$ in one galaxy is significantly more extended than the non-resonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly$α$ sources, discussing viable scenarios for the creation of sizable ionized bubbles ($>$1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, requiring additional sources of ionizing photons. Therefore, the fact that we detect Ly$α$ from these galaxies suggests diverse scenarios on escape of Ly$α$ during the epoch of reionization. High spectral resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ (18 pages, 7 figures, 2 tables)

Spectroscopic Confirmation of CEERS NIRCam-selected Galaxies at $\boldsymbol{z \simeq 8-10}$

Authors: Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Seiji Fujimoto, Vital Fernández, Jeyhan S. Kartaltepe, Intae Jung, Justin W. Cole, Denis Burgarella, Katherine Chworowsky, Taylor A. Hutchison, Alexa M. Morales, Casey Papovich, Raymond C. Simons, Ricardo O. Amorín, Bren E. Backhaus, Micaela B. Bagley, Laura Bisigello, Antonello Calabrò, Marco Castellano, Nikko J. Cleri, Romeel Davé, Avishai Dekel, Henry C. Ferguson, Adriano Fontana , et al. (23 additional authors not shown)

Abstract: We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a… ▽ More We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV<-20), with blue continua (-2.3<beta<-1.9) and low dust attenuation (A_V=0.15(+0.3,-0.1)); and at least one of them has high stellar mass for a galaxy at that redshift (log(M_*/M_sol)=9.3(+0.2,-0.3)). Considered together with spectroscopic observations of other CEERS NIRCam-selected high-z galaxy candidates in the literature, we find a high rate of redshift confirmation and low rate of confirmed interlopers (8.3%). Ten out of 34 z>8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8. △ Less

Submitted 6 July, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

Journal ref: Astrophys. J. Lett. 951, L22 (2023)

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Confirmation and refutation of very luminous galaxies in the early universe

Authors: Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Callum T. Donnan, Denis Burgarella, Adam Carnall, Fergus Cullen, James S. Dunlop, Vital Fernández, Seiji Fujimoto, Intae Jung, Melanie Krips, Rebecca L. Larson, Casey Papovich, Pablo G. Pérez-González, Ricardo O. Amorín, Micaela B. Bagley, Véronique Buat, Caitlin M. Casey, Katherine Chworowsky, Seth H. Cohen, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company , et al. (12 additional authors not shown)

Abstract: During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far… ▽ More During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times. △ Less

Submitted 15 August, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: Accelerated Article Preview (AAP) version published in Nature

CEERS Key Paper VI: JWST/MIRI Uncovers a Large Population of Obscured AGN at High Redshifts

Authors: G. Yang, K. I. Caputi, C. Papovich, P. Arrabal Haro, M. B. Bagley, P. Behroozi, E. F. Bell, L. Bisigello, V. Buat, D. Burgarella, Y. Cheng, N. J. Cleri, R. Dave, M. Dickinson, D. Elbaz, H. C. Ferguson, S. L. Finkelstein, N. A. Grogin, N. P. Hathi, M. Hirschmann, B. W. Holwerda, M. Huertas-Company, T. Hutchison, E. Iani, J. S. Kartaltepe , et al. (13 additional authors not shown)

Abstract: Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN populat… ▽ More Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 27 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 0.5$ dex at $z \approx 3$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts. △ Less

Submitted 15 May, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

CLEAR: Survey Overview, Data Analysis and Products

Authors: Raymond C. Simons, Casey Papovich, Ivelina G. Momcheva, Gabriel Brammer, Vicente Estrada-Carpenter, Steven L. Finkelstein, Catherine M. Gosmeyer, Jasleen Matharu, Jonathan R. Trump, Bren E. Backhaus, Yingjie Cheng, Nikko J. Cleri, Henry C. Ferguson, Kristian Finlator, Mauro Giavalisco, Zhiyuan Ji, Intae Jung, Jennifer M. Lotz, Rosalia O'Brien, Rosalind E. Skelton, Vithal Tilvi, Benjamin Weiner

Abstract: We present an overview of the CANDELS Lyman-a Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other… ▽ More We present an overview of the CANDELS Lyman-a Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other programs, the full CLEAR dataset includes spectroscopic imaging of these fields over 0.8-1.7 um. In this Paper, we describe the CLEAR survey, the survey strategy, the data acquisition, reduction, processing, and science products and catalogs released alongside this paper. The catalogs include emission line fluxes and redshifts derived from the combination of the photometry and grism spectroscopy for 6048 galaxies, primarily ranging from 0.2 < z < 3. We also provide an overview of CLEAR science goals and results. In conjunction with this Paper we provide links to electronic versions of the data products, including 1D + 2D extracted spectra and emission line maps. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: Accepted in ApJS. 29 pages, 14 figures. The data products (1D + 2D spectra, emission line maps) and catalogs (photometric and spectroscopic catalogs of GOODS-N and GOODS-S) are available at https://archive.stsci.edu/hlsp/clear

On the nature of disks at high redshift seen by JWST/CEERS with contrastive learning and cosmological simulations

Authors: J. Vega-Ferrero, M. Huertas-Company, L. Costantin, P. G. Pérez-González, R. Sarmiento, J. S. Kartaltepe, A. Pillepich, M. B. Bagley, S. L. Finkelstein, E. J. McGrath, J. H. Knapen, P. Arrabal Haro, E. F. Bell, F. Buitrago, A. Calabrò, A. Dekel, M. Dickinson, H. Domínguez Sánchez, D. Elbaz, H. C. Ferguson, M. Giavalisco, B. W. Holwerda, D. D. Kocesvski, A. M. Koekemoer, V. Pandya , et al. (4 additional authors not shown)

Abstract: Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to… ▽ More Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $ 3 \leq z \leq 6$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks -- either CNN-based or visually -- populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and non-oblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at $z \geq 3$ remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift. △ Less

Submitted 25 October, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 28 pages, 22 figures. Accepted for publication in ApJ

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey

Authors: Micaela B. Bagley, Nor Pirzkal, Steven L. Finkelstein, Casey Papovich, Danielle A. Berg, Jennifer M. Lotz, Gene C. K. Leung, Henry C. Ferguson, Anton M. Koekemoer, Mark Dickinson, Jeyhan S. Kartaltepe, Dale D. Kocevski, Rachel S. Somerville, L. Y. Aaron Yung, Bren E. Backhaus, Caitlin M. Casey, Marco Castellano, Óscar A. Chávez Ortiz, Katherine Chworowsky, Isabella G. Cox, Romeel Davé, Kelcey Davis, Vicente Estrada-Carpenter, Adriano Fontana, Seiji Fujimoto , et al. (23 additional authors not shown)

Abstract: We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation r… ▽ More We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5<z<4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam (m=30.6-30.9, 5sigma) to discover galaxies to z>12, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 21 pages, 9 figures, submitted to AAS Journals

CEERS Spectroscopic Confirmation of NIRCam-Selected z > 8 Galaxy Candidates with JWST/NIRSpec: Initial Characterization of their Properties

Authors: Seiji Fujimoto, Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Rebecca L. Larson, Denis Burgarella, Micaela B. Bagley, Peter Behroozi, Katherine Chworowsky, Michaela Hirschmann, Jonathan R. Trump, Stephen M. Wilkins, L. Y. Aaron Yung, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Henry C. Ferguson, Adriano Fontana, Norman A. Grogin, Andrea Grazian, Lisa J. Kewley, Dale D. Kocevski, Jennifer M. Lotz, Laura Pentericci , et al. (16 additional authors not shown)

Abstract: We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings… ▽ More We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$β$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$α$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$β$ equivalent widths ($\simeq$1100$Å$), and elevated production efficiency of ionizing photons ($\log(ξ_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$β$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models. △ Less

Submitted 26 May, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

Spatially Resolved Stellar Populations of $0.3<z<6.0$ Galaxies in WHL0137-08 and MACS0647+70 Clusters as Revealed by JWST: How do Galaxies Grow and Quench Over Cosmic Time?

Authors: Abdurro'uf, Dan Coe, Intae Jung, Henry C. Ferguson, Gabriel Brammer, Kartheik G. Iyer, Larry D. Bradley, Pratika Dayal, Rogier A. Windhorst, Adi Zitrin, Ashish Kumar Meena, Masamune Oguri, Jose M. Diego, Vasily Kokorev, Paola Dimauro, Angela Adamo, Christopher J. Conselice, Brian Welch, Eros Vanzella, Tiger Yu-Yang Hsiao, Xinfeng Xu, Namrata Roy, Celia R. Mulcahey

Abstract: We study the spatially resolved stellar populations of 444 galaxies at $0.3<z<6.0$ in two clusters (WHL0137-08 and MACS0647+70) and a blank field, combining imaging data from HST and JWST to perform spatially resolved spectral energy distribution (SED) modeling using pixedfit. The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fiel… ▽ More We study the spatially resolved stellar populations of 444 galaxies at $0.3<z<6.0$ in two clusters (WHL0137-08 and MACS0647+70) and a blank field, combining imaging data from HST and JWST to perform spatially resolved spectral energy distribution (SED) modeling using pixedfit. The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fields allows us to resolve some galaxies to sub-kpc scales (for 109 of our galaxies). At redshifts around cosmic noon and higher ($2.5\lesssim z\lesssim 6.0$), we find mass doubling times to be independent of radius, inferred from flat specific star formation rate (sSFR) radial profiles and similarities between the half-mass and half-SFR radii. At lower redshifts ($1.5\lesssim z\lesssim 2.5$), a significant fraction of our star-forming galaxies show evidence for nuclear starbursts, inferred from centrally elevated sSFR, and a much smaller half-SFR radius compared to the half-mass radius. At later epochs, we find more galaxies suppress star formation in their center but are still actively forming stars in the disk. Overall, these trends point toward a picture of inside-out galaxy growth consistent with theoretical models and simulations. We also observe a tight relationship between the central mass surface density and global stellar mass with $\sim 0.38$ dex scatter. Our analysis demonstrates the potential of spatially resolved SED analysis with JWST data. Future analysis with larger samples will be able to further explore the assembly of galaxy mass and the growth of their structures △ Less

Submitted 7 February, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 31 pages, 18 figures, accepted for publication in ApJ. Some examples and tutorials of spatially resolved SED analysis will be available at https://github.com/aabdurrouf/JWST-HST_resolvedSEDfits

CEERS Key Paper IV: Galaxies at $4 < z < 9$ are Bluer than They Appear -- Characterizing Galaxy Stellar Populations from Rest-Frame $\sim 1$ micron Imaging

Authors: Casey Papovich, Justin Cole, Guang Yang, Steven L. Finkelstein, Guillermo Barro, Véronique Buat, Denis Burgarella, Pablo G. Pérez-González, Paola Santini, Lise-Marie Seillé, Lu Shen, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, Caitlin M. Casey, Marco Castellano, Katherine Chworowsky, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Mark Dickinson, Henry C. Ferguson, Adriano Fontana , et al. (24 additional authors not shown)

Abstract: We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame… ▽ More We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame spectral-energy distribution (SED) to nearly 1 micron for galaxies in this redshift range. By modeling the galaxies' SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å) $-$ $I$-band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower (stellar)-mass-to-light ratios. The MIRI data reduce the stellar masses by $\langle Δ\log M_\ast\rangle=0.25$ dex at $4<z<6$ (a factor of 1.8) and 0.37 dex at $6<z<9$ (a factor of 2.3). This also reduces the star-formation rates (SFRs) by $\langle Δ\log\mathrm{SFR} \rangle=0.14$ dex at $4<z<6$ and 0.27 dex at $6<z<9$. The MIRI data also improve constraints on the allowable stellar mass formed in early star-formation. We model this using a star-formation history that includes both a "burst' at $z_f=100$ and a slowly varying ("delayed-$τ$") model. The MIRI data reduce the allowable stellar mass by 0.6 dex at $4<z< 6$ and by $\approx$1 dex at $6<z<9$. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early universe ($z\approx9$). Therefore, observations of rest-frame $\gtrsim$1 $μ$m are paramount for constraining the stellar-mass build-up in galaxies at very high-redshifts. △ Less

Submitted 25 March, 2023; v1 submitted 30 December, 2022; originally announced January 2023.

Comments: Updated with accepted ApJ version. Part of the CEERS Focus Issue. 27 pages, many figures (4 Figure Sets, available upon reasonable request)

New $z > 7$ Lyman-alpha Emitters in EGS: Evidence of an Extended Ionized Structure at $z \sim 7.7$

Authors: Intae Jung, Steven L. Finkelstein, Rebecca L. Larson, Taylor A. Hutchison, Amber N. Straughn, Micaela B. Bagley, Marco Castellano, Nikko J. Cleri, M. C. Cooper, Mark Dickinson, Henry C. Ferguson, Benne W. Holwerda, Jeyhan S. Kartaltepe, Seonwoo Kim, Anton M. Koekemoer, Casey Papovich, Hyunbae Park, Laura Pentericci, Pablo G. Perez-Gonzalez, Mimi Song, Sandro Tacchella, Benjamin J. Weiner, Christopher N. A. Willmer, Jorge A. Zavala

Abstract: We perform a ground-based near-infrared spectroscopic survey using the Keck/MOSFIRE spectrograph to target Ly$α$ emission at $7.0<z<8.2$ from 61 galaxies to trace the ionization state of the intergalactic medium (IGM). We cover a total effective sky area of $\sim10^\prime\times10^\prime$ in the Extended Groth Strip field of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. From o… ▽ More We perform a ground-based near-infrared spectroscopic survey using the Keck/MOSFIRE spectrograph to target Ly$α$ emission at $7.0<z<8.2$ from 61 galaxies to trace the ionization state of the intergalactic medium (IGM). We cover a total effective sky area of $\sim10^\prime\times10^\prime$ in the Extended Groth Strip field of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. From our observations, we detect Ly$α$ emission at a $>$4$σ$ level in eight $z>7$ galaxies, which include additional members of the known $z\sim7.7$ Ly$α$-emitter (LAE) cluster (Tilvi et al. 2020). With the addition of these newly-discovered $z\sim7.7$ LAEs, this is currently the largest measured LAE cluster at $z>7$. The unusually-high Ly$α$ detection rate at $z\sim7.7$ in this field suggests significantly stronger Ly$α$ emission from the clustered LAEs than from the rest of our targets. We estimate the ionized bubble sizes around these LAEs and conclude that the LAEs are clustered within an extended ionized structure created by overlapping ionized bubbles which allow the easier escape of Ly$α$ from galaxies. It is remarkable that the brightest object in the cluster has the lowest measured redshift of the Ly$α$ line, being placed in front of the other LAEs in the line-of-sight direction. This suggests that we are witnessing the enhanced IGM transmission of Ly$α$ from galaxies on the rear side of an ionized area. This could be a consequence of Ly$α$ radiative transfer: Ly$α$ close to the central velocity is substantially scattered by the IGM while Ly$α$ from the rear-side galaxies is significantly redshifted to where it has a clear path. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Submitted to ApJ (17 pages, 8 figures, 4 tables)

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength 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 the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength 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 the wide and deep multi-waveband 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 fourteen 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 AGN 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.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

CEERS Key Paper I: An Early Look into the First 500 Myr of Galaxy Formation with JWST

Authors: Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Stephen M. Wilkins, Jeyhan S. Kartaltepe, Casey Papovich, L. Y. Aaron Yung, Pablo Arrabal Haro, Peter Behroozi, Mark Dickinson, Dale D. Kocevski, Anton M. Koekemoer, Rebecca L. Larson, Aurelien Le Bail, Alexa M. Morales, Pablo G. Perez-Gonzalez, Denis Burgarella, Romeel Dave, Michaela Hirschmann, Rachel S. Somerville, Stijn Wuyts, Volker Bromm, Caitlin M. Casey, Adriano Fontana, Seiji Fujimoto , et al. (42 additional authors not shown)

Abstract: We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a det… ▽ More We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions. △ Less

Submitted 4 November, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: Replaced with published version

ALMA FIR View of Ultra High-redshift Galaxy Candidates at $z\sim$ 11-17: Blue Monsters or Low-$z$ Red Interlopers?

Authors: Seiji Fujimoto, Steven L. Finkelstein, Denis Burgarella, Chris L. Carilli, Véronique Buat, Caitlin M. Casey, Laure Ciesla, Sandro Tacchella, Jorge A. Zavala, Gabriel Brammer, Yoshinobu Fudamoto, Masami Ouchi, Francesco Valentino, M. C. Cooper, Mark Dickinson, Maximilien Franco, Mauro Giavalisco, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Takashi Kojima, Rebecca L. Larson, Eric J. Murphy, Casey Papovich, Pablo G. Pérez-González , et al. (28 additional authors not shown)

Abstract: We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr… ▽ More We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$σ$ line feature at $338.726\pm0.007$~GHz exactly coinciding with the JWST source position, with a 2\% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$μ$m at $z=16.01$ or [CII]158$μ$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11--13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11--17 is favored in every candidate as a very blue (UV continuum slope of $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($Δ$ $χ^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$--5, where such populations predicted from their luminosity functions and EW([OIII]+H$β$) distributions are abundant in survey volumes used for the identification of the $z\sim$11--17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search. △ Less

Submitted 26 July, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 25 pages, 8 figures, 7 tables. ApJ in press

CEERS Epoch 1 NIRCam Imaging: Reduction Methods and Simulations Enabling Early JWST Science Results

Authors: Micaela B. Bagley, Steven L. Finkelstein, Anton M. Koekemoer, Henry C. Ferguson, Pablo Arrabal Haro, Mark Dickinson, Jeyhan S. Kartaltepe, Casey Papovich, Pablo G. Pérez-González, Nor Pirzkal, Rachel S. Somerville, Christopher N. A. Willmer, Guang Yang, L. Y. Aaron Yung, Adriano Fontana, Andrea Grazian, Norman A. Grogin, Michaela Hirschmann, Lisa J. Kewley, Allison Kirkpatrick, Dale D. Kocevski, Jennifer M. Lotz, Aubrey Medrano, Alexa M. Morales, Laura Pentericci , et al. (12 additional authors not shown)

Abstract: We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations (Yang et al. in prep). W… ▽ More We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations (Yang et al. in prep). We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known pre-launch issues such as 1/f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with pre-launch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated dataset in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public datasets released from the Directors Discretionary Early Release Science (DD-ERS) program. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: 27 pages, 14 figures, submitted to ApJ. Accompanying CEERS public Data Release 0.5 available at ceers.github.io/releases.html

CEERS Key Paper V: A triality on the nature of HST-dark galaxies

Authors: Pablo G. Pérez-González, Guillermo Barro, Marianna Annunziatella, Luca Costantin, Ángela García-Argumánez, Elizabeth J. McGrath, Rosa M. Mérida, Jorge A. Zavala, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Peter Behroozi, Eric F. Bell, Laura Bisigello, Véronique Buat, Antonello Calabrò, Caitlin M. Casey, Nikko J. Cleri, Rosemary T. Coogan, M. C. Cooper, Asantha R. Cooray, Avishai Dekel, Mark Dickinson, David Elbaz, Henry C. Ferguson , et al. (30 additional authors not shown)

Abstract: The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometr… ▽ More The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2<z<6 with masses 9<log M/M_sun<11 and a variety of specific SFRs (<1 to >100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3<z<5, masses log M/M_sun~10 and post-starburst stellar mass-weighted ages (0.5-1 Gyr); and (3) 11% are strong young starbursts with indications of high-EW emission lines (typically, [OIII]+Hbeta) at 6<z<7 and log M/M_sun~9.5. The sample is dominated by disk-like galaxies with a remarkable compactness for XELG-z6 (effective radii smaller than 0.4 kpc). Large attenuations in SFGs, 2<A(V)<5 mag, are found within 1.5 times the effective radius, approximately 2 kpc, while QGs present A(V)~0.2 mag. Our SED-fitting technique reproduces the expected dust emission luminosities of IR-bright and sub-millimeter galaxies. This study implies high levels of star formation activity between z~20 and z~10, where virtually 100% of our galaxies had already formed 10^8 M_sun of their stellar content, 60% of them had assembled 10^9 M_sun, and 10% up to 10^10 M_sun (in situ or ex situ). (abridged) △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Published in CEERS ApJL Focus Issue, ApJL 946, L16

Journal ref: ApJL 946, L16 https://iopscience.iop.org/article/10.3847/2041-8213/acb3a5

First Look at z > 1 Bars in the Rest-Frame Near-Infrared with JWST Early CEERS Imaging

Authors: Yuchen Guo, Shardha Jogee, Steven L. Finkelstein, Zilei Chen, Eden Wise, Micaela B. Bagley, Guillermo Barro, Stijn Wuyts, Dale D. Kocevski, Jeyhan S. Kartaltepe, Elizabeth J. McGrath, Henry C. Ferguson, Bahram Mobasher, Mauro Giavalisco, Ray A. Lucas, Jorge A. Zavala, Jennifer M. Lotz, Norman A. Grogin, Marc Huertas-Company, Jesús Vega-Ferrero, Nimish P. Hathi, Pablo Arrabal Haro, Mark Dickinson, Anton M. Koekemoer, Casey Papovich , et al. (23 additional authors not shown)

Abstract: Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stell… ▽ More Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0. △ Less

Submitted 11 December, 2022; v1 submitted 16 October, 2022; originally announced October 2022.

Comments: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

Semi-analytic forecasts for Roman -- the beginning of a new era of deep-wide galaxy surveys

Authors: L. Y. Aaron Yung, Rachel S. Somerville, Steven L. Finkelstein, Peter Behroozi, Romeel Davé, Henry C. Ferguson, Jonathan P. Gardner, Gergö Popping, Sangeeta Malhotra, Casey Papovich, James E. Rhoads, Micaela B. Bagley, Michaela Hirschmann, Anton M. Koekemoer

Abstract: The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizatio… ▽ More The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizations of 2-deg^2 lightcones, containing a total of >25 million simulated galaxies with -16 < MUV < -25 spanning z ~ 0 to 10. This dataset enables a new set of experiments with the impacts of survey size on the derived galaxy formation and cosmological constraints. The intrinsic and observable galaxy properties are predicted using a well-established, physics-based semi-analytic modelling approach. We provide forecasts for number density, cosmic SFR, field-to-field variance, and angular two-point correlation functions, and demonstrate how the future wide-field surveys will be able to improve these measurements relative to current generation surveys. We also present a comparison between these lightcones and others that have been constructed with empirical models. The mock lightcones are designed to facilitate the exploration of multi-instrument synergies and connecting with current generation instruments and legacy surveys. In addition to Roman, we also provide photometry for a number of other instruments on upcoming facilities, including Euclid and Rubin, as well as the instruments that are part of many legacy surveys. Full object catalogues and data tables for the results presented in this work are made available through a web-based, interactive portal https://www.simonsfoundation.org/semi-analytic-forecasts. △ Less

Submitted 5 December, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures. Accepted for publication in MNRAS

The hidden side of cosmic star formation at z > 3: Bridging optically-dark and Lyman break galaxies with GOODS-ALMA

Authors: Mengyuan Xiao, David Elbaz, Carlos Gómez-Guijarro, Lucas Leroy, Longji Bing, Emanuele Daddi, Benjamin Magnelli, Maximilien Franco, Luwenjia Zhou, Mark Dickinson, Tao Wang, Wiphu Rujopakarn, Georgios E. Magdis, Ezequiel Treister, Hanae Inami, Ricardo Demarco, Mark T. Sargent, Xinwen Shu, Jeyhan S. Kartaltepe, David M. Alexander, Matthieu Béthermin, Frederic Bournaud, Laure Ciesla, Henry C. Ferguson, Steven L. Finkelstein , et al. (15 additional authors not shown)

Abstract: Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order… ▽ More Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift. △ Less

Submitted 10 February, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 672, A18 (2023)

A Machine Learning Approach to Predict Missing Flux Densities in Multi-band Galaxy Surveys

Authors: Nima Chartab, Bahram Mobasher, Asantha Cooray, Shoubaneh Hemmati, Zahra Sattari, Henry C. Ferguson, David B. Sanders, John R. Weaver, Daniel Stern, Henry J. McCracken, Daniel C. Masters, Sune Toft, Peter L. Capak, Iary Davidzon, Mark Dickinson, Jason Rhodes, Andrea Moneti, Olivier Ilbert, Lukas Zalesky, Conor McPartland, Istvan Szapudi, Anton M. Koekemoer, Harry I. Teplitz, Mauro Giavalisco

Abstract: We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample o… ▽ More We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample of $i<25$ AB mag galaxies. We find that with available $i$-band fluxes, $r$, $u$, IRAC/$ch2$ and $z$ bands provide most of the information regarding the redshift with importance decreasing from $r$-band to $z$-band. We also find that for the same sample, IRAC/$ch2$, $Y$, $r$ and $u$ bands are the most relevant bands in stellar mass measurements with decreasing order of importance. Investigating the inter-correlation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in $YJH$ bands can be simulated/predicted with an accuracy of $1σ$ mag scatter $\lesssim 0.2$ for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template-fitting. Such a machine learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template-fitting inevitable in the presence of a few bands. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 15 pages, 14 figures, accepted for publication in ApJ

CEERS Key Paper III: The Resolved Host Properties of AGN at 3 < z < 5 with JWST

Authors: Dale D. Kocevski, Guillermo Barro, Elizabeth J. McGrath, Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Shardha Jogee, Guang Yang, Mark Dickinson, Nimish P. Hathi, Bren E. Backhaus, Eric F. Bell, Laura Bisigello, Véronique Buat, Denis Burgarella, Caitlin M. Casey, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Darren Croton, Emanuele Daddi, Adriano Fontana, Seiji Fujimoto, Jonathan P. Gardner, Eric Gawiser , et al. (34 additional authors not shown)

Abstract: We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshif… ▽ More We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ; 12 Pages, 6 Figures

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may als… ▽ More Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

Probing the earliest phases in the formation of massive galaxies with simulated HST+JWST imaging data from Illustris

Authors: Ángela García-Argumánez, Pablo G. Pérez-González, Armando Gil de Paz, Gregory F. Snyder, Pablo Arrabal Haro, Micaela B. Bagley, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Anton Koekemoer, Casey Papovich, Nor Pirzkal, Harry C. Ferguson, L. Y. Aaron Yung, Marianna Annunziatella, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Benne W. Holwerda, Rosa María Mérida González, Caitlin Rose, Mauro Giavalisco, Norman A. Grogin, Dale D. Kocevski

Abstract: We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ 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 broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis… ▽ More We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ 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 broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M$_{\ast}>10^{10}\,$M$_{\odot}$) simulated galaxies at $1<z<4$ that evolve into a local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxy. In particular, we explore the potential of HST and JWST datasets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this dataset for characterizing the first episodes in the SFH of local M$_{\ast}>10^{11}\,$M$_{\odot}$ 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% to 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$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies would have started their formation by $z=16$, forming the first 5% of their stellar mass present at $z \sim 1$ by $z=4.5$, 10% by $z=3.7$, and 25% by $z=2.7$. △ Less

Submitted 19 December, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, 4 tables. ApJ in press. Summary of changes from original submission: the major change is that we now include in Sec. 6 the comparison of the results obtained for our sample of massive 1 < z < 4 progenitors with those obtained by considering all massive galaxies at 1 < z < 4 in the simulated images. Several figures and sections have been updated

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

The Physical Conditions of Emission-Line Galaxies at Cosmic Dawn from JWST/NIRSpec Spectroscopy in the SMACS 0723 Early Release Observations

Authors: Jonathan R. Trump, Pablo Arrabal Haro, Raymond C. Simons, Bren E. Backhaus, Ricardo O. Amorín, Mark Dickinson, Vital Fernández, Casey Papovich, David C. Nicholls, Lisa J. Kewley, Samantha W. Brunker, John J. Salzer, Stephen M. Wilkins, Omar Almaini, Micaela B. Bagley, Danielle A. Berg, Rachana Bhatawdekar, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Laure Ciesla, Nikko J. Cleri, Justin W. Cole , et al. (39 additional authors not shown)

Abstract: We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wav… ▽ More We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$λ$5008/H$β$ ratios, similar [OIII]$λ$4364/H$γ$ ratios, and higher ($\sim$0.5 dex) [NeIII]$λ$3870/[OII]$λ$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$λ$3870/[OII]$λ$3728, [OIII]$λ$4364/H$γ$, and [OIII]$λ$5008/H$β$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$λ$4364/H$γ$ and [OIII]$λ$(4960+5008)/H$β$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn. △ Less

Submitted 19 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

Semi-analytic forecasts for JWST -- VI. Simulated lightcones and galaxy clustering predictions

Authors: L. Y. Aaron Yung, Rachel S. Somerville, Henry C. Ferguson, Steven L. Finkelstein, Jonathan P. Gardner, Romeel Davé, Micaela B. Bagley, Gergö Popping, Peter Behroozi

Abstract: In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new s… ▽ More In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new set of 8 ultra-deep lightcones with 132 sq. arcmin footprints, containing galaxies up to z ~ 12 down to the magnitudes expected to be reached in the deepest JWST surveys. These mock lightcones are extracted from dissipationless N-body simulations and populated with galaxies using the well-established, computationally efficient Santa Cruz semi-analytic model for galaxy formation. We provide a wide range of predicted physical properties, and simulated photometry from NIRCam and many other instruments. We explore the predicted counts and luminosity functions and angular two-point correlation functions for galaxies in these simulated lightcones. We also explore the predicted field-to-field variance using multiple lightcone realizations. We find that these lightcones reproduce the available measurements of observed clustering from 0.2 < z < 7.5 very well. We provide predictions for galaxy clustering at high redshift that may be obtained from future JWST observations. All of the lightcones presented here are made available through a web-based, interactive data release portal https://flathub.flatironinstitute.org/group/sam-forecasts. △ Less

Submitted 27 July, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: Final paper of the Semi-analytic forecasts for JWST series. 24 pages, 14 figures, many years of exciting science to come, accepted for publication in MNRAS

The Low-Redshift Lyman Continuum Survey II: New Insights into LyC Diagnostics

Authors: Sophia R. Flury, Anne E. Jaskot, Harry C. Ferguson, Gabor Worseck, Kirill Makan, John Chisholm, Alberto Saldana-Lopez, Daniel Schaerer, Stephan McCandliss, Bingjie Wang, N. M. Ford, M. S. Oey, Timothy Heckman, Zhiyuan Ji, Mauro Giavalisco, Ricardo Amorin, Hakim Atek, Jeremy Blaizot, Sanchayeeta Borthakur, Cody Carr, Marco Castellano, Stefano Cristiani, Stephane de Barros, Mark Dickinson, Steven L. Finkelstein , et al. (22 additional authors not shown)

Abstract: The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum 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 t… ▽ More The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum 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 a 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 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. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: ApJ, accepted. 31 pages, 26 figures, 1 table

Searching for Islands of Reionization: A Potential Ionized Bubble Powered by a Spectroscopic Overdensity at z = 8.7

Authors: Rebecca L. Larson, Steven L. Finkelstein, Taylor A. Hutchison, Casey Papovich, Micaela B. Bagley, Mark E. Dickinson, Sofía Rojas-Ruiz, Henry C. Ferguson, Intae Jung, Mauro Giavalisco, A. Grazian, Laura Pentericci, Sandro Tacchella

Abstract: We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10m 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 hours of integration, at z = 8.665 +/- 0.001 with an integrated signal-to-noise… ▽ More We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10m 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 hours of integration, at z = 8.665 +/- 0.001 with an integrated signal-to-noise > 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 (Zitrin et al. 2015). 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 galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. 2021. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 which is powering one (or more) ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field. △ Less

Submitted 29 August, 2023; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 19 pages, 12 figures, 8 tables, Published in ApJ

Journal ref: The Astrophysical Journal, 2022, Volume 930, Issue 2, id.104, 19 pp

An Isolated Stellar-Mass Black Hole Detected Through Astrometric Microlensing

Authors: Kailash C. Sahu, Jay Anderson, Stefano Casertano, Howard E. Bond, Andrzej Udalski, Martin Dominik, Annalisa Calamida, Andrea Bellini, Thomas M. Brown, Marina Rejkuba, Varun Bajaj, Noe Kains, Henry C. Ferguson, Chris L. Fryer, Philip Yock, Przemek Mroz, Szymon Kozlowski, Pawel Pietrukowicz, Radek Poleski, Jan Skowron, Igor Soszynski, Michael K. Szymanski, Krzysztof Ulaczyk, Lukasz Wyrzykowski, Richard Barry , et al. (68 additional authors not shown)

Abstract: We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge.… ▽ More We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of six years, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of the Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 +/- 1.3 Msun and a distance of 1.58 +/- 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic-disk stars at similar distances by an amount corresponding to a transverse space velocity of ~45 km/s, suggesting that the BH received a 'natal kick' from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial-velocity measurements of Galactic X-ray binaries, and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique. △ Less

Submitted 22 July, 2022; v1 submitted 31 January, 2022; originally announced January 2022.

Comments: 37 pages, Published in ApJ

Journal ref: ApJ, 933, 83 (2022)

The Low-Redshift Lyman Continuum Survey. Unveiling the ISM properties of low-$z$ Lyman continuum emitters

Authors: Alberto Saldana-Lopez, Daniel Schaerer, John Chisholm, Sophia R. Flury, Anne E. Jaskot, Gábor Worseck, Kirill Makan, Simon Gazagnes, Valentin Mauerhofer, Anne Verhamme, Ricardo O. Amorín, Harry C. Ferguson, Mauro Giavalisco, Andrea Grazian, Matthew J. Hayes, Timothy M. Heckman, Alaina Henry, Zhiyuan Ji, Rui Marques-Chaves, Stephan R. McCandliss, M. S. Oey, Göran Östlin, Laura Pentericci, Trinh X. Thuan, Maxime Trebitsch , et al. (2 additional authors not shown)

Abstract: Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV at… ▽ More Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV attenuation, by fitting the stellar continuum with a combination of simple stellar population models. The models, together with accurate LyC flux measurements, allow to determine the absolute LyC photon escape fraction for each galaxy ($f_{\rm esc}^{\rm abs}$). We measure the equivalent widths and residual fluxes of multiple HI and low-ionization state (LIS) lines, and the geometrical covering fraction adopting the picket-fence model. The $f_{\rm esc}^{\rm abs}$ spans a wide range, with a median (0.16, 0.84 quantiles) of 0.04 (0.02, 0.20), and 50 out of the 89 galaxies detected in the LyC. The HI and LIS line equivalent widths scale with the UV luminosity and attenuation, and inversely with the residual flux of the lines. The HI and LIS residual fluxes are correlated, indicating that the neutral gas is spatially traced by the LIS transitions. We find the observed trends of the absorption lines and the UV attenuation are primarily driven by the covering fraction. The non-uniform gas coverage demonstrates that LyC photons escape through low-column density channels in the ISM. The equivalent widths and residual fluxes of the UV lines strongly correlate with $f_{\rm esc}^{\rm abs}$: strong LyC leakers show weak absorption lines, low UV attenuation, and large Ly$α$ equivalent widths. We finally show that simultaneous UV absorption line and dust attenuation measurements can predict, on average, the escape fraction of galaxies and the method can be applied to galaxies across a wide redshift range. △ Less

Submitted 31 March, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: 30 pages, 16 figures, 3 tables; accepted for publication in Astronomy and Astrophysics on December 16, 2021. Tables A1 to A4 are part of the LzLCS science products and will be publicly available in a dedicated website

Journal ref: A&A 663, A59 (2022)

The Low-Redshift Lyman Continuum Survey I: New, Diverse Local Lyman-Continuum Emitters

Authors: Sophia R. Flury, Anne E. Jaskot, Harry C. Ferguson, Gabor Worseck, Kirill Makan, John Chisholm, Alberto Saldana-Lopez, Daniel Schaerer, Stephan McCandless, Bingjie Wang, N. M. Ford, Timothy Heckman, Zhiyuan Ji, Mauro Giavalisco, Ricardo Amorin, Hakim Atek, Jeremy Blaizot, Sanchayeeta Borthakur, Cody Carr, Marco Castellano, Stefano Cristiani, Stephane de Barros, Mark Dickinson, Steven L. Finkelstein, Brian Fleming , et al. (22 additional authors not shown)

Abstract: 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 univ… ▽ More 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 HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV 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. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: ApJS, accepted. 28 pages, 19 figures, 8 tables. Full machine readable tables will be made available by the publisher at the time of publication

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

Authors: C. Gómez-Guijarro, D. Elbaz, M. Xiao, V. I. Kokorev, G. E. Magdis, B. Magnelli, E. Daddi, F. Valentino, M. T. Sargent, M. Dickinson, M. Béthermin, M. Franco, A. Pope, B. S. Kalita, L. Ciesla, R. Demarco, H. Inami, W. Rujopakarn, X. Shu, T. Wang, L. Zhou, D. M. Alexander, F. Bournaud, R. Chary, H. C. Ferguson , et al. (16 additional authors not shown)

Abstract: Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived… ▽ More Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: Accepted by A&A. 26 pages, 13 figures

Journal ref: A&A 659, A196 (2022)

On the Stellar Populations of Galaxies at z=9-11: The Growth of Metals and Stellar Mass at Early Times

Authors: Sandro Tacchella, Steven L. Finkelstein, Micaela Bagley, Mark Dickinson, Henry C. Ferguson, Mauro Giavalisco, Luca Graziani, Norman A. Grogin, Nimish Hathi, Taylor A. Hutchison, Intae Jung, Anton M. Koekemoer, Rebecca L. Larson, Casey Papovich, Norbert Pirzkal, Sofia Rojas-Ruiz, Mimi Song, Raffaella Schneider, Rachel S. Somerville, Stephen M. Wilkins, L. Y. Aaron Yung

Abstract: We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consis… ▽ More We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consistent modeling of emission lines. This approach allows us to assess how different priors affect our results, and how well we can break degeneracies between dust attenuation, stellar ages, metallicity and emission lines using data which probe only the rest-frame ultraviolet to optical wavelengths. We measure a median observed ultraviolet spectral slope $β=-1.87_{-0.43}^{+0.35}$ for relatively massive star-forming galaxies ($9<\log(M_{\star}/M_{\odot})<10$), consistent with no change from $z=4$ to $z=9-10$ at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sub-linear slope ($0.7\pm0.2$) and specific star-formation rates of $3-10~\mathrm{Gyr}^{-1}$. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of $\sim50-150$ Myr, which corresponds to 50\% formation redshifts of $z_{50}\sim10-12$ at $z\sim9$ and is of the order of the dynamical timescales of these systems. Importantly, the models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these $z=9-11$ galaxies, with our results consistent with SFHs implying both a shallow and steep increase of the cosmic SFR density with time at $z>10$. △ Less

Submitted 8 April, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 27 pages + Appendix; accepted by ApJ

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

Authors: Nicole E. Drakos, Bruno Villasenor, Brant E. Robertson, Ryan Hausen, Mark E. Dickinson, Henry C. Ferguson, Steven R. Furlanetto, Jenny E. Greene, Piero Madau, Alice E. Shapley, Daniel P. Stark, Risa H. Wechsler

Abstract: In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark… ▽ More In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 39 pages, 27 figures. Submitted to ApJ. Comments welcome

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Authors: C. Gómez-Guijarro, D. Elbaz, M. Xiao, M. Béthermin, M. Franco, B. Magnelli, E. Daddi, M. Dickinson, R. Demarco, H. Inami, W. Rujopakarn, G. E. Magdis, X. Shu, R. Chary, L. Zhou, D. M. Alexander, F. Bournaud, L. Ciesla, H. C. Ferguson, S. L. Finkelstein, M. Giavalisco, D. Iono, S. Juneau, J. S. Kartaltepe, G. Lagache , et al. (14 additional authors not shown)

Abstract: Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas… ▽ More Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks. △ Less

Submitted 21 December, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: Accepted by A&A. 30 pages, 18 figures

Journal ref: A&A 658, A43 (2022)