Search

The ultraviolet continuum traces young stars while the near-infrared unveils older stellar populations and dust-obscured regions. Balmer emission lines provide insights on gas properties and young stellar objects but are highly affected by dust attenuation. The near-infrared Paschen lines suffer less dust attenuation and can be used to measure star formation rates (SFRs) in star-forming regions obscured by dust clouds. We select 13 sources between redshifts 1 and 3 observed with HST, JWST/NIRCam and NIRSpec based on the availability of at least one Balmer and one Paschen line with S/N > 5. With a newly-developed version of CIGALE, we fit their hydrogen line equivalent widths (EWs) and photometric data. We assess the impacts of the removal of spectroscopic data by comparing the quality of the fits of the spectro-photometric data to those with photometric data only. We compare the single (BC03) vs binary (BPASS) stellar populations models in the fitting process of spectro-photometric data. We derive the differential attenuation and explore different attenuation recipes by fitting spectro-photometric data with BC03. For each stellar model and for each input dataset (with and without EWs), we quantify the deviation on the SFRs and stellar masses from the "standard" choice. On average, the SFRs are overestimated and the stellar masses are underestimated when EWs are not included as input data. We find a major contribution of the H${\alpha}$ emission line to the broadband photometric measurements of our sources, and a trend of increasing contribution with specific SFR. Using the BPASS models has a significant impact on the derived SFRs and stellar masses. We show that a flexible attenuation recipe provides more accurate estimates of the dust attenuation parameters, especially the differential attenuation which agrees with the original value of Charlot & Fall (2000)., Comment: Submitted to Astronomy & Astrophysics, 19 pages, 15 figures

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., Comment: 23 pages, 17 figures, submitted to ApJ

The James Webb Space Telescope (JWST) has discovered a surprising population of bright galaxies in the very early universe (< 500 Myrs after the Big Bang) that is hard to explain with conventional galaxy formation models and whose physical properties remain to be fully understood. Insight into the internal physics of galaxies is captured best via observations of excited-state atomic transitions of ionized gas, but beyond z~7-9, the brightest spectral signatures are redshifted into the mid-infrared regime, where observations are increasingly more difficult. Here, we present the first detection of a hydrogen recombination line (H{\alpha}) and doubly-ionized oxygen ([OIII]4959,5007{\AA}) at z>10 using the JWST Mid-Infrared Instrument, MIRI. These detections place the bright galaxy GHZ2/GLASS-z12 at z=12.33+/-0.02, making it the most distant astronomical object with direct spectroscopic detection of these lines and the brightest confirmed object at this epoch. These observations provide key insights into the conditions of this primeval galaxy, which shows hard ionizing conditions rarely seen in the local Universe and likely driven by compact, young (<30 Myr) star formation. Its oxygen-to-hydrogen abundance is close to a tenth of the solar value, indicating a rapid metal enrichment during the earliest phases of galaxy formation. This study confirms the unique conditions of the brightest and most distant galaxies recently discovered by JWST and the huge potential of mid-IR observations to characterize these systems, opening a range of new possibilities in the study of the very early Universe., Comment: 14 pages, 6 figures. Submitted for publication. Comments are most welcome

Age-dating and weighting stellar populations in galaxies at various cosmic epochs are essential steps to study galaxy formation through cosmic times. Evolutionary population synthesis models with different input physics are used towards this aim. In particular, the contribution from the thermally pulsing asymptotic-giant-branch (TP-AGB) stellar phase, which peaks for intermediate-age 0.6-2 Gyr systems, has been debated upon for decades. Here we report the detection of strong cool star signatures in the rest-frame near-infrared spectra of three young (~1 Gyr), massive (~10^10 Msun) quiescent galaxies at large look-back time, z=1-2, using JWST/NIRSpec. The co-existence of oxygen- and carbon-type absorption features, spectral edges and features from rare species such as Vanadium, and possibly Zirconium, reveal a strong contribution from TP-AGB stars. Population synthesis models with significant TP-AGB contribution reproduce the observations considerably better than those with weak TP-AGB, which are those commonly used. These findings call for revisions of published stellar population fitting results, pointing to lower masses and younger ages, with additional implications on cosmic dust production and chemical enrichment. These results will stimulate new generations of improved models informed by these and future observations., Comment: submitted

Extreme emission line galaxies (EELGs) are typically characterized by high equivalent widths (EWs) which are driven by elevated specific star formation rates (sSFR) in low-mass galaxies with subsolar metallicities and little dust. Such extreme systems are rare in the local universe, but the number density of EELGs increases with redshift. Such starburst galaxies are currently presumed to be the main drivers of hydrogen reionization over 5.50.4-1 which indicates that star formation may be the dominant source of ionization. Based on the photometric fluxes, we find an increase of EW([OIII]+H$\beta$) with sSFR and $\Sigma_{SFR}$, and a decrease with age and stellar mass. The sample of EELGs can reach $\Sigma_{SFR}>$10Msun yr$^{-1}$kpc$^{-2}$ which indicate they are strong candidates of LyC leakers. Another indirect indicator is the high values of O32>5 that can be reached for some galaxies in the sample., Comment: Submitted to A&A

The galaxy integrated star-formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) has been proposed as a valuable diagnostic of the mass accumulation in galaxies as being more tightly related to the physics of star-formation (SF) and stellar feedback than other SF indicators. In this paper, we assemble a statistical sample of 230 galaxies observed with JWST in the GLASS and CEERS spectroscopic surveys to estimate Balmer line based dust attenuations and SFRs, and UV rest-frame effective radii. We study the evolution of galaxy SFR and $\Sigma_{\rm SFR}$ in the first 1.5 Billion years of our Universe, finding that $\Sigma_{\rm SFR}$ is mildly increasing with redshift with a linear slope of $0.16 \pm 0.06$. We also explore the dependence of SFR and $\Sigma_{\rm SFR}$ on stellar mass, showing that a SF 'Main-Sequence' and a $\Sigma_{\rm SFR}$ `Main-Sequence' are in place out to z=10, with a similar slope compared to the same relations at lower redshifts. We find that the specific SFR (sSFR) and $\Sigma_{\rm SFR}$ are correlated with the [OIII]5007/[OII]3727 ratio and with indirect estimates of the escape fraction of Lyman continuum photons, hence they likely play an important role in the evolution of ionization conditions and in the escape of ionizing radiation. We also search for spectral outflow signatures in a subset of galaxies observed at high resolution, finding an outflow incidence of $2/11$ ($=20\%^{32\%}_{9\%}$) at $z<6$, but no evidence at $z>6$ ($<26\%$). Finally, we find a positive correlation between A$_V$ and $\Sigma_{\rm SFR}$, and a flat trend as a function of sSFR, indicating that there is no evidence of a drop of A$_V$ in extremely star-forming galaxies between z=4 and 10. This might be at odds with a dust-clearing outflow scenario, which might instead take place at redshifts $z\geq 10$, as suggested by some theoretical models., Comment: Submitted to A&A ; 17 pages, 9 figures, 1 table

The Cosmological Advanced Survey Telescope for Optical and UV Research (CASTOR) is a planned flagship space telescope, covering the blue-optical and UV part of the spectrum. Here we introduce the CASTOR image simulator, a Python GalSim package-based script capable of generating mock CASTOR images from an input catalogue. We generate example images from the CASTOR Wide, Deep, and Ultra-Deep surveys using simulated light-cones from the Santa Cruz Semi-Analytic Model. We make predictions for the performance of these surveys by comparing galaxies that are extracted from each image using Source Extractor to the input catalogue. We find that the Wide, Deep, and Ultra-Deep surveys will be complete to ~27, 29 and 30 mag, respectively, in the UV, u, and g filters, with the UV-split and u-split filters reaching a shallower depth. With a large area of ~2200 deg$^2$, the Wide survey will detect hundreds of millions of galaxies out to z~4, mostly with $M_\ast \gtrsim 10^9 M_\odot$. The Ultra-Deep survey will probe to z~5, detecting a large fraction of $M_\ast \simeq 10^8 M_\odot$ galaxies. These powerful samples will enable precision measurements of the distribution of star formation in the cosmic web, connecting the growth of stellar mass to the assembly of dark matter halos over two thirds of the history of the Universe, and other core goals of CASTOR's legacy surveys. These image simulations and the tools developed to generate them will be a vital planning tool to estimate CASTOR's performance and iterate the telescope and survey designs prior to launch., Comment: Submitted to AJ

The resonant scattering interaction between Ly$\alpha$ photons and neutral hydrogen implies that a partially neutral IGM can significantly impact the detectability of Ly$\alpha$ emission in galaxies. The redshift evolution of the Ly$\alpha$ equivalent width distribution of galaxies thus offers a key probe of the degree of ionization during the Epoch of Reionization (EoR). Previous in-depth investigations at $z$ $\geq$ 7 were limited by ground-based instrument capabilities. We present an extensive study of Ly$\alpha$ emission from galaxies at 4 < $z$ < 8.5, observed from the CEERS and JADES surveys in the JWST NIRSpec/PRISM configuration. The sample consists of 235 galaxies, among which we identify 65 as Ly$\alpha$ emitters. We first measure Ly$\alpha$ escape fractions from Balmer lines, and explore the correlations with the inferred galaxies' physical properties, which are similar to those found at lower redshift. We also investigate the possible connection between the escape of Ly$\alpha$ photons and the inferred escape fractions of LyC photons obtained from indirect indicators. We then analyze the redshift evolution of the Ly$\alpha$ emitter fraction, finding lower average values at $z$ = 5 and 6 compared to ground-based observations. At $z$ = 7 we find a very large difference in Ly$\alpha$ visibility between the EGS and GOODS-South fields, possibly due to the presence of early reionized regions in the EGS. Such large variance is also expected in the Cosmic Dawn II radiation-hydrodynamical simulation. Our findings suggest a scenario in which the ending phase of the EoR is characterized by $\sim$ 1 pMpc ionized bubbles around a high fraction of moderately bright galaxies. Finally, we characterize such two ionized regions found in the EGS at $z$ = 7.18 and $z$ = 7.49 by estimating the radius of the ionized bubble that each of the spectroscopically-confirmed members could have created., Comment: Submitted to A&A

Ultraviolet (UV; rest-frame $\sim1200-2000$ A) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at $5.6

We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3$\lesssim$z$\lesssim$4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Of these 51, we find 12 AGN had $\geq$4$\sigma$ detected LyC flux in the WFC3/UVIS images. Using space- and ground-based data from X-ray to radio, we fit the multi-wavelength photometric data of each AGN to a CIGALE SED and correlate various SED parameters to the LyC flux. KS-tests of the SED parameter distributions for the LyC-detected and non-detected AGN showed they are likely not distinct samples. However, we find that X-ray luminosity, star-formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlation of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission, 6 $\mu m$ luminosity, and anti-correlation with metallicity and $A_{FUV}$. We simulate the LyC escape fraction ($f_{esc}$) using the CIGALE and IGM transmission models for the LyC-detected AGN and find an average $f_{esc}$$\simeq$18%, weighted by uncertainties. We stack the LyC flux of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our $sub-mm\ detected$ F336W sample shows the brightest stacked LyC flux. These findings indicate that LyC-production and -escape in AGN is more complicated than the simple assumption of thermal emission and a 100% escape fraction. Further testing of AGN models with larger samples than presented here is needed., Comment: 21 pages, 6 figures, 3 tables. Accepted for publication in The Astrophysical Journal

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, the basic demographics of clumps and their host galaxies remain uncertain. Using the Hubble Space Telescope/Wide Field Camera 3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we detect and analyze giant star-forming clumps in galaxies at 0.5 ≤ z ≤ 1, connecting two epochs when clumps are common (at cosmic high noon, z ∼ 2) and rare (in the local Universe). We construct a clump sample whose rest-frame 1600 Å luminosity is 3 times higher than the most luminous local H ii regions (MUV ≤ −16 AB). In our sample, 35% ± 3% of low-mass galaxies (log[M∗/M⊙] < 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22% ± 3% and 22% ± 4% for intermediate (10 ≤ log[M∗/M⊙] ≤ 10.5) and high-mass (log[M∗/M⊙] > 10.5) galaxies, in agreement with previous studies. When compared to similar-mass nonclumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher star formation rates (SFRs) and bluer rest-frame U − V colors, while high-mass clumpy SFGs tend to be larger than nonclumpy SFGs. However, clumpy and nonclumpy SFGs have similar Sérsic index, indicating a similar underlying density profile. Furthermore, we investigate how the UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher specific SFR galaxies and are found closer to their hosts' galactic centers.

High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high-redshift galaxies. Studying these substructures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W ≤ 25 AB mag sample of 6767 galaxies at 0.5 ≤ z ≤ 3 in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-passband filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10% of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies (fclumpy) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that fclumpy increases with redshift, reaching ∼65% at z ∼ 1.5. We also find that fclumpy in low-mass galaxies () is 10% higher compared to that of their high-mass counterparts (

We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in 0.2 < z < 0.8 groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20 k group catalog. Limiting our analysis to a complete sample of log (M*/M⊙) > 10.56 group galaxies, we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow (∼Gyr) timescales. We find a quenching time of approximately Gyr, consistent with the slow processes of strangulation and delayed-then-rapid quenching although more data are needed to confirm this result.

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 (AV'., >'., 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™ ) > 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.

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. Each galaxy has been imaged with the JWST Near-Infrared Camera, which provides rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host by fitting its spectral energy distribution using a combination of galaxy and AGN templates. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk, and one is dominated by pointlike emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. When compared to a sample of mass-matched inactive galaxies, we find that the AGN hosts have morphologies that are less disturbed and more bulge-dominated. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or poststarburst stellar population. The presence of AGN in passively evolving galaxies at z > 3 is significant because a rapid feedback mechanism is required in most semianalytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings show that AGN can continue to inject energy into these systems after their star formation is curtailed, potentially heating their halos and preventing renewed star formation. Additional observations will be needed to determine what role this feedback may play in helping to quench these systems and/or maintain their quiescent state.

We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 ${\mu}$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 ${\mu}$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 ${\mu}$m CEERS data to derive rest-frame 7.7 ${\mu}$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${{\sigma}_{{\Delta}SFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 ${\mu}$m luminosity, which has a scatter of ${{\sigma}_{{\Delta}SFR}}$ = 0.21 dex. The relatively small decrease in ${\sigma}$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 ${\mu}$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity., Comment: 20 pages, 11 figures, 2 tables, submitted to ApJ

Dark matter halo demographics and assembly histories are a manifestation of cosmological structure formation and have profound implications for the formation and evolution of galaxies. In particular, merger trees provide fundamental input for several modelling techniques, such as semi-analytic models (SAMs), sub-halo abundance matching (SHAM), and decorated halo occupation distribution models (HODs). Motivated by the new ultra-high-redshift (z > 10) frontier enabled by JWST, we present a new suite of Gadget at Ultrahigh Redshift with Extra-Fine Timesteps (GUREFT) dark matter-only cosmological simulations that are carefully designed to capture halo merger histories and structural properties in the ultra-z universe. The simulation suite consists of four 1024^3-particle simulations with box sizes of 5, 15, 35, and 90 Mpc h-1, each with 170 snapshots stored between 40 > z > 6. With the unprecedented number of available snapshots and strategically chosen dynamic range covered by these boxes, gureft uncovers the emerging dark matter halo populations and their assembly histories in the earliest epochs of cosmic history. In this work, we present the halo mass functions between z ~ 20 to 6 down to log(Mvir/Msun) ~ 5, and show that at high redshift, these robust halo mass functions can differ substantially from commonly used analytic approximations or older fitting functions in the literature. We also present key physical properties of the ultra-z halo population, such as concentration and spin, as well as their mass growth and merger rates, and again provide updated fitting functions., Comment: 18 pages, 17 figures, submitted to MNRAS (v2: corrected for a bug that caused the Tinker+08 line in Fig.5 to evolve more than it should)

The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction ($f_{esc}$). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their $f_{esc}$ can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at $6 \leq z \leq 9$ from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average $6

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$\mu$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$\mu$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., Comment: 21 pages, 10 figures. Resubmitted to ApJS after revision

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 {\AA} luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar S\'ersic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center., Comment: 21 pages, 13 figures, accepted for publication in ApJ

The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and Visible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic survey fields: GOODS-N, GOODS-S, EGS, and COSMOS. We introduce this survey by presenting a thorough search for galaxies at $z\gtrsim2.4$ that leak significant Lyman continuum (LyC) radiation, as well as a stringent constraint on the LyC escape fraction ($f_{\rm esc}$) from stacking the UV images of a population of star-forming galaxies with secure redshifts. Our extensive search for LyC emission and stacking analysis benefit from the catalogs of high-quality spectroscopic redshifts compiled from archival ground-based data and HST slitless spectroscopy, carefully vetted by dedicated visual inspection efforts. We report a sample of five galaxies as individual LyC leaker candidates, showing $f_{\rm esc}^{\rm rel}\gtrsim60\%$ estimated using detailed Monte Carlo analysis of intergalactic medium attenuation. We develop a robust stacking method to apply to five samples of in total 85 non-detection galaxies in the redshift range of $z\in[2.4,3.7]$. Most stacks give tight 2-$\sigma$ upper limits below $f_{\rm esc}^{\rm rel}<6\%$. A stack for a subset of 32 emission-line galaxies shows tentative LyC leakage detected at 2.9-$\sigma$, indicating $f_{\rm esc}^{\rm rel}=5.7\%$ at $z\sim2.65$, supporting the key role of such galaxies in contributing to the cosmic reionization and maintaining the UV ionization background. These new F275W and F435W imaging mosaics from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes., Comment: 33 pages, 21 figures, and 5 tables. Resubmitted after addressing the referee report

The mass assembly history (MAH) of dark matter halos plays a crucial role in shaping the formation and evolution of galaxies. MAHs are used extensively in semi-analytic and empirical models of galaxy formation, yet current analytic methods to generate them are inaccurate and unable to capture their relationship with the halo internal structure and large-scale environment. This paper introduces FLORAH, a machine-learning framework for generating assembly histories of ensembles of dark matter halos. We train FLORAH on the assembly histories from the GUREFT and VSMDPL N-body simulations and demonstrate its ability to recover key properties such as the time evolution of mass and concentration. We obtain similar results for the galaxy stellar mass versus halo mass relation and its residuals when we run the Santa Cruz semi-analytic model on FLORAH-generated assembly histories and halo formation histories extracted from an N-body simulation. We further show that FLORAH also reproduces the dependence of clustering on properties other than mass (assembly bias), which is not captured by other analytic methods. By combining multiple networks trained on a suite of simulations with different redshift ranges and mass resolutions, we are able to construct accurate main progenitor branches (MPBs) with a wide dynamic mass range from $z = 0$ up to an ultra-high redshift $z \approx 20$, currently far beyond that of a single N-body simulation. FLORAH is the first step towards a machine learning-based framework for planting full merger trees; this will enable the exploration of different galaxy formation scenarios with great computational efficiency at unprecedented accuracy., Comment: Submitted to MNRAS; 18 pages, 13 + 3 figures, 1 table

The Cosmic Evolution Early Release Science Survey (CEERS), targeting the Extended Groth Strip extragalactic field, is one of the JWST Director's Discretionary Early Release Science programs. To date, all observations have been executed and include NIRCam/MIRI imaging and NIRSpec/NIRCam spectroscopic exposures. Here, we discuss the MIRI imaging, which includes eight pointings, four of which provide deep imaging with the bluer bands (F560W, F770W) and four with contiguous wavelength coverage in F1000W, F1280W, F1500W, and F1800W, where two of these also include coverage in F770W and F2100W. We present a summary of the data, the data quality, and data reduction. The data reduction is based on the JWST Calibration Pipeline combined with custom modifications and additional steps designed to enhance the output quality, including improvements in astrometry and the removal of detector artifacts. We estimate the image depth of the reduced mosaics, and show that these generally agree with expectations from the Exposure Time Calculator. We compare the MIRI F560W and F770W flux densities for bright sources to measurements from Spitzer/IRAC Ch3 (5.8 $\mu$m) and Ch4 (8.0 $\mu$m), and we find that they agree with systematic differences of $<0.1$ mag. For the redder MIRI bands, we assess their quality by studying the spectral energy distributions (SEDs) of Galactic stars. The SEDs are consistent with the expected Rayleigh-Jeans law with a deviation $\sim 0.03$ mag, indicating that the MIRI colors are reliable. We also discuss all publicly released data products (images and source catalogs), which are available on the CEERS website (https://ceers.github.io/)., Comment: 17 pages, 11 figures, and 4 tables. ApJL in press

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$\sigma$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$\sigma$) 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

We investigate the morphology and physical properties of a sample of 22 IR-selected dusty star-forming galaxies at Cosmic Noon (z ~ 2), using James Webb Space Telescope NIRCam images obtained in the EGS field for the Cosmic Evolution Early Release Science survey. The exceptional resolution of the NIRCam images allows us to spatially resolve these galaxies up to 4.4um and identify their bulge/core even when very extinguished by dust. Based on red-green-blue images using the F115W, F200W and F444W filters, we divide each galaxy in several uniformly colored regions, fit their respective Spectral Energy Distribution and measure dust attenuations, stellar masses, star formation rates and ages. After classifying each region as star-forming or quiescent, we assign galaxies to three classes, depending on whether active star-formation is located in the core, in the disk or in both. (i) ~70% of our DSFGs have a compact highly dust attenuated star-forming core that can contain up to 80% of the star-formation of the galaxy but only 20-30% of its stellar mass, and is always surrounded by a larger, less attenuated massive disk (no blue nuggets); (ii) 64% (27%) of disks are significantly (strongly) lopsided, likely due to asymmetric cold gas accretion, major mergers and/or large scale instabilities; (iii) 23% of galaxies have a star-forming core embedded in a quiescent disk, they are undergoing outside-in quenching, often facilitated by their strong lopsidedness inducing small and large scale instabilities; (iv) some galaxies host highly heterogeneous disks in term of RGB colors: these are driven by in-homogeneous dust attenuation; and (v) we find surprising evidence for clump-like substructures being quiescent and/or residing in quiescent regions. This work demonstrates the major impact JWST/NIRCam has on understanding the complexity of the evolution of distant massive galaxies., Comment: 22 pages, 21 figures, submitted to A&A

We present a spectroscopic survey of Ly$\alpha$ emitters in the Extended Groth Strip (EGS) field, targeting the regime near the Epoch of Reionization. Using Keck/DEIMOS, we observed 947 high-$z$ candidates with photometric redshifts from 3 $< z_\text{phot} <$ 7 and down to an $H$-band (HST/WFC3 F160W) magnitude limit of < 27.5. Observations were taken over the course of 8 nights, with integration times ranging from 4 to 7.8 hours. Our survey secured 137 unique redshifts, 126 of which are Ly$\alpha$ emitters at 2.8 $< z <$ 6.5 with a mean redshift of $\overline{z} = 4.3$. We provide a comprehensive redshift catalog for our targets, as well as the reduced one- and two- dimensional spectra for each object. These observations will provide an important auxiliary dataset for the JWST Directors Discretionary Early Release Science (DD-ERS) program the Cosmic Evolution Early Release Science Survey (CEERS), which recently completed near- and mid-IR imaging and spectroscopy of galaxies in the EGS field., Comment: 8 pages, 7 figures, 2 tables, Accepted in MNRAS

Optical rest-frame spectroscopic diagnostics are usually employed to distinguish between star formation and AGN-powered emission. However, this method is biased against dusty sources, hampering a complete census of the AGN population across cosmic epochs. To mitigate this effect, it is crucial to observe at longer wavelengths in the rest-frame near-infrared (near-IR), which is less affected by dust attenuation and can thus provide a better description of the intrinsic properties of galaxies. AGN diagnostics in this regime have not been fully exploited so far, due to the scarcity of near-IR observations of both AGNs and star-forming galaxies, especially at redshifts higher than 0.5. Using Cloudy photoionization models, we identify new AGN - star formation diagnostics based on the ratio of bright near-infrared emission lines, namely [SIII] 9530 Angstrom, [CI] 9850 Angstrom, [PII] 1.188 $\mu m$, [FeII] $1.257 \mu m$, and [FeII] $1.64 \mu m$ to Paschen lines (either Pa$\gamma$ or Pa$\beta$), providing simple, analytical classification criteria. We apply these diagnostics to a sample of 64 star-forming galaxies and AGNs at 0 < z < 1, and 65 sources at 1 < z < 3 recently observed with JWST-NIRSpec in CEERS. We find that the classification inferred from the near-infrared is broadly consistent with the optical one based on the BPT and the [SII]/H$\alpha$ ratio. However, in the near-infrared, we find $\sim 60 \%$ more AGNs than in the optical (13 instead of 8), with 5 sources classified as 'hidden' AGNs, showing a larger AGN contribution at longer wavelengths, possibly due to the presence of optically thick dust. The diagnostics we present provide a promising tool to find and characterize AGNs from z=0 to z=3 with low and medium-resolution near-IR spectrographs in future surveys., Comment: Paper accepted for publication in A&A on 05/09/2023. Three public Github repositories include: (1) a table with emission line measurements for the paper sample : https://github.com/Anthony96/Line_measurements_nearIR , Cloudy emission line predictions for star-forming galaxies and AGN models : https://github.com/Anthony96/star-forming_models , https://github.com/Anthony96/AGN_models

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$\sigma$, 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., Comment: 21 pages, 10 figures. Submitted to ApJL

We study a new population of extremely red objects (EROs) recently discovered by JWST based on their NIRCam colors F277W$-$F444W $>1.5$ mag. We find 37 EROs in the CEERS field with F444W $<28$ mag and photometric redshifts between $57$ by up to a factor $\sim$60. Similarly, if they are OSOs with luminosities in the L$_{\rm bol}>10^{46-47}$ erg s$^{-1}$ range, their number would exceed that of bright blue QSOs by more than two orders of magnitude. Additional photometry at mid-IR wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution., Comment: 26 pages, 10 figures, submitted to ApJ

[Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.175$; 35 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and S\'ersic index become unreliable. We extended our study to fainter sources ($S/N>10$; 69 galaxies) by fixing their S\'ersic index to unity. The S\'ersic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to radial color gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component, the so-called blue nuggets., Comment: Published in A&A. 18 pages, 14 figures

High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W$\leq 25$ AB mag sample of 6767 galaxies at $0.5\leq z\leq 3$ in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-pass band filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10$\%$ of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies ($\rm f_{clumpy}$) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that $\rm f_{clumpy}$ increases with redshift, reaching $\sim 65\%$ at $z\sim 1.5$. We also find that $\rm f_{clumpy}$ in low-mass galaxies ($\rm 9.5\leq log(M_*/M_\odot)\leq 10$) is 10$\%$ higher compared to that of their high-mass counterparts ($\rm log(M_*/M_\odot)>10.5$). Moreover, we find no evidence of significant environmental dependence of $\rm f_{clumpy}$ for galaxies at the redshift range of this study. Our results suggest that the fragmentation of gas clouds under violent disk instability remains the primary driving mechanism for clump formation, and incidents common in dense environments, such as mergers, are not the dominant processes., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in ApJ

We analyze the Near Infrared ($\sim0.8-1\mu$m) rest-frame morphologies of galaxies with $\log M_*/M_\odot>9$ in the redshift range $010.5$) at $z\sim5$, and bulge-dominated galaxies also exist at these early epochs, confirming a rich and evolved morphological diversity of galaxies $\sim1$ Gyr after the Big Bang. Finally, we find that the morphology-quenching relation is already in place for massive galaxies at $z>3$, with massive quiescent galaxies ($\log M_*/M_\odot>10.5$) being predominantly bulge-dominated., Comment: Submitted to A&A, comments welcome

Substantial populations of massive quiescent galaxies at $z\ge3$ challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at $z>3$. As a result, studies report that as much as 70\% of quiescent galaxies at $z>3$ may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at $3\lesssim z \lesssim 6$ using JWST NIRCam imaging data. We use empirically-constrained galaxy SED templates to define a region in the $F277W-F444W$ vs. $F150W-F277W$ color plane that captures quiescent galaxies at $z>3$. We apply this color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and identify 44 candidate $z\gtrsim3$ quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of post-starburst galaxies. We derive volume density estimates of $n\sim1-4\times10^{-5}$\,Mpc$^{-3}$ at $3< z <5$, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam's wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies., Comment: 12 pages, 4 figures, 2 tables, submitted to ApJ

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

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., Comment: Accepted by A&A. 13 pages, 9 figures

We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$\alpha$-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$\alpha$ emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal poor. We investigate Ly$\alpha$ 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$\alpha$ emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly$\alpha$ 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$\alpha$ 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$\alpha$ from these galaxies suggests diverse scenarios on escape of Ly$\alpha$ during the epoch of reionization. High spectral resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization., Comment: Submitted to ApJ (18 pages, 7 figures, 2 tables)

We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of $\sim 20 \ M_\odot \textrm{yr}^{-1}$ at $z\sim4$ at a $1\sigma$ level, allowing us to quantify the amount of contamination from dusty star-forming sources in our quiescent sample. Starting with a parent sample of candidate massive quiescent galaxies from the Stevans et al. 2021 v1 SHELA catalog, we use the Bayesian \textsc{Bagpipes} spectral energy distribution fitting code to derive robust stellar masses ($M_*$) and star-formation rates (SFRs) for these sources, and select a conservative sample of 36 candidate massive ($M_* > 10^{11}M_\odot$) quiescent galaxies, with specific SFRs at $>2\sigma$ below the star-forming main sequence at $z\sim4$. Based on ALMA imaging, six of these candidate quiescent galaxies have the presence of significant dust-obscured star-formation, thus were removed from our final sample. This implies a $\sim 17\%$ contamination rate from dusty star-forming galaxies with our selection criteria using the v1 SHELA catalog. This conservatively-selected quiescent galaxy sample at $z=3-5$ will provide excellent targets for future observations to better constrain how massive galaxies can both grow and shut-down their star-formation in a relatively short time period., Comment: 16 pages, 13 figures

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.38 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., Comment: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

A substantial number of ultra-high redshift (8 < z < 17) galaxy candidates have been detected with JWST, posing the question: are these observational results surprising in the context of current galaxy formation models? We address this question using the well-established Santa Cruz semi-analytic models, implemented within merger trees from the new suite of cosmological N-body simulations GUREFT, which were carefully designed for ultra-high redshift studies. Using our fiducial models calibrated at z=0, we present predictions for stellar mass functions, rest-frame UV luminosity functions, and various scaling relations. We find that our (dust-free) models predict galaxy number densities at z~11 (z~13) that are an order of magnitude (a factor of ~30) lower than the observational estimates. We estimate the uncertainty in the observed number densities due to cosmic variance, and find that it leads to a fractional error of ~20-30% at z=11 (~30-80% at z=14) for a 100 sq arcmin field. We explore which processes in our models are most likely to be rate-limiting for the formation of luminous galaxies at these early epochs, considering the halo formation rate, gas cooling, star formation, and stellar feedback, and conclude that it is mainly efficient stellar-driven winds. We find that a modest boost of a factor of ~4 to the UV luminosities, which could arise from a top-heavy stellar initial mass function, would bring our current models into agreement with the observations. Adding a stochastic component to the UV luminosity can also reconcile our results with the observations., Comment: 20 pages, 16 figures, accepted for publication in MNRAS

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., Comment: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

We take advantage of the NIRCam photometric observations available as part of the Cosmic Evolution Early Release Science survey (CEERS) to identify and analyse very red sources in an effort to discover very dusty star forming galaxies. We select red galaxies as objects with a S/N>3 at 4.4 $\mu$m and a S/N<2 in all JWST and HST filters at $\lambda\leq2\mu$m, which corresponds to [F200W]-[F444W]>1.2 considering CEERS depths. This selection is ideal to identify very dusty (Av>1 mag) galaxies with stellar masses between $10^6$ to $10^{10}\, \rm M_{\odot}$ at z<5, more massive dusty galaxies at z=5-18 and galaxies at z>18 due to the Lyman absorption, independently of their dust extinction. Our sample of F200W-dropouts contains no strong candidates at z>6.5, instead it consists almost completely (~81%) of z<2 low-mass galaxies, with a median stellar mass of $10^{7.3} \rm M_{\odot}$. These galaxies show an exceptional dust extinction with median value of Av=4.9 mag, completely unexpected given their low stellar mass. The remaining galaxies, which are at z<6.5, show similar large dust extinction (Av>1), but they are generally more massive $>10^{7.5}\rm M_{\odot}$., Comment: 30 pages, 1 table, 19 figures, accepted for publication in A&A