Search

Accurately quantifying the ionizing photon production efficiency ($\xi_\rm{ion}$) of $z>6$ star-forming galaxies (SFGs) is necessary to understand their contribution to reionization. We investigate the ionizing properties of N=289 SFGs selected at $z=6.9-7.6$ from the JWST Cycle-1 imaging programmes; PRIMER and JADES. We use BAGPIPES to consistently infer the equivalent widths of their [OIII]+$\rm{H\beta}$ emission lines ($W_\lambda$) and their physical properties. To supplement our high-redshift galaxies, we measure $W_\lambda$ photometrically for a sample of N=253 $z=3.2-3.6$ SFGs selected from the VANDELS spectroscopic survey. Comparing these samples, we find a strong apparent redshift evolution in their median $W_\lambda$, increasing from $W_\lambda =310\pm25\,\r{A}$ in VANDELS to $W_\lambda =630\pm30\,\r{A}$ in our JWST-based sample. In the JWST sample at $z>7$, we find that $W_\lambda$ correlates with both stellar mass and UV luminosity, with high-mass, $M_{ UV}-$faint galaxies producing systematically weaker emission lines. Moreover, we discover a departure from the standard log-normal shape of the $W_\lambda$ distribution, characterised by a more pronounced tail at lower $W_\lambda$, consistent with increasingly bursty star formation. Using $W_\lambda$ as a proxy for $\xi_\rm{ion}$, and UV spectral slope as a proxy for LyC escape fraction ($f_\rm{esc}$), we find a minority of galaxies with high $\xi_\rm{ion}$ and $f_\rm{esc}$ (e.g., $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.6$ and $f_\rm{esc}\sim0.15$). However, we find that the LyC photon budget at $z>7$ is dominated by galaxies with more moderate output, close to the sample median of $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.3$ and $f_\rm{esc}\sim0.05$. This is consistent with estimates for the number of LyC photons required to power reionization at $z>7$, with no evidence for over or under-production., Comment: 18 pages, 9 figures + 1 appendix. Submitted to MNRAS

We present the JWST Emission Line Survey (JELS), a Cycle 1 JWST imaging programme exploiting the wavelength coverage and sensitivity of NIRCam to extend narrow-band rest-optical emission line selection into the epoch of reionization (EoR) for the first time, and to enable unique studies of the resolved ionised gas morphology in individual galaxies across cosmic history. The primary JELS observations comprise $\sim4.7\mu$m narrow-band imaging over $\sim63$ arcmin$^{2}$ designed to enable selection of H$\alpha$ emitters at $z\sim6.1$, as well as the selection of a host of novel emission-line samples including [OIII] at $z\sim8.3$ and Pa $\alpha/\beta$ at $z\sim1.5/2.8$. For the prime F466N and F470N narrow-band observations, the emission-line sensitivities achieved are up to $\sim2\times$ more sensitive than current slitless spectroscopy surveys (5$\sigma$ limits of 1.1-1.6$\times10^{-18}\text{erg s}^{-1}\text{cm}^{-2}$), corresponding to unobscured H$\alpha$ star-formation rates (SFRs) of 1-1.6 $\text{M}_{\odot}\,\text{yr}^{-1}$ at $z\sim6.1$ and extending emission-line selections in the EoR to fainter populations. Simultaneously, JELS also obtained F200W broadband and F212N narrow-band imaging (H$\alpha$ at $z\sim2.23$) that probes SFRs $\gtrsim5\times$ fainter than previous ground-based narrow-band studies ($\sim0.2 \text{M}_{\odot}\text{yr}^{-1}$), offering an unprecedented resolved view of star formation at cosmic noon. In this paper we describe the detailed JELS survey design, key data processing steps specific to the survey observations, and demonstrate the exceptional data quality and imaging sensitivity achieved. We then summarise the key scientific goals of JELS and present some early science results, including examples of spectroscopically confirmed H$\alpha$ and [OIII] emitters discovered by JELS that illustrate the novel parameter space probed., Comment: Submitted to MRNAS; 13 pages, 8 figures. Accompanying data will be publicly released upon acceptance of this and companion paper (Pirie et al., arXiv:2410.11808). V2 updated to reflect revised author list and cross-citation to Pirie et al

We present a spatially resolved analysis of four star-forming galaxies at $z = 4.44-5.64$ using data from the JWST PRIMER and ALMA-CRISTAL surveys to probe the stellar and inter-stellar medium properties on the sub-kpc scale. In the $1-5\,\mu{\rm m}$ JWST NIRCam imaging we find that the galaxies are composed of multiple clumps (between $2$ and $\sim 8$) separated by $\simeq 5\,{\rm kpc}$, with comparable morphologies and sizes in the rest-frame UV and optical. Using BAGPIPES to perform pixel-by-pixel SED fitting to the JWST data we show that the SFR ($\simeq 25\,{\rm M}_{\odot}/{\rm yr}$) and stellar mass (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot}) \simeq 9.5$) derived from the resolved analysis are in close ($ \lesssim 0.3\,{\rm dex}$) agreement with those obtained by fitting the integrated photometry. In contrast to studies of lower-mass sources, we thus find a reduced impact of outshining of the older (more massive) stellar populations in these normal $z \simeq 5$ galaxies. Our JWST analysis recovers bluer rest-frame UV slopes ($\beta \simeq -2.1$) and younger ages ($\simeq 100\,{\rm Myr}$) than archival values. We find that the dust continuum from ALMA-CRISTAL seen in two of these galaxies correlates, as expected, with regions of redder rest-frame UV slopes and the SED-derived $A_{\rm V}$, as well as the peak in the stellar mass map. We compute the resolved IRX-$\beta$ relation, showing that the IRX is consistent with the local starburst attenuation curve and further demonstrating the presence of an inhomogeneous dust distribution within the galaxies. A comparison of the CRISTAL sources to those from the FirstLight zoom-in simulation of galaxies with the same $M_{\star}$ and SFR reveals similar age and colour gradients, suggesting that major mergers may be important in the formation of clumpy galaxies at this epoch., Comment: 16 pages, 8 figures, 3 tables, plus 4 page appendix. Submitted to MNRAS

We present a set of photometric catalogs primarily aimed at providing the community with a comprehensive database for the study of galaxy populations in the high redshift Universe. The set gathers data from eight JWST NIRCam observational programs, targeting the Abell 2744 (GLASS-JWST, UNCOVER, DDT2756 and GO3990), EGS (CEERS), COSMOS and UDS (PRIMER), and GOODS North and South (JADES and NGDEEP) deep fields, for a total area of $\sim$0.2 sq. degrees. Photometric estimates are obtained by means of well-established techniques, including tailored improvements designed to enhance the performance on the specific dataset. We also include new measurements from HST archival data, thus collecting 16 bands spanning from 0.44 to 4.44 $\mu$m. A grand total of $\sim$530 thousand sources is detected on stacks of NIRCam 3.56 and 4.44 $\mu$m mosaics. We assess the photometric accuracy by comparing fluxes and colors against archival catalogs. We also provide photometric redshift estimates, statistically validated against a large set of robust spectroscopic data. The catalogs are publicly available on the Astrodeep website., Comment: 28 pages, 16 figures; accepted for publication on Astronomy & Astrophysics

We present a detailed JWST/NIRSpec and NIRCam analysis of a gravitationally-lensed galaxy ($\rm \mu=17-21$) at redshift 6.14 magnified by the Hubble Frontier Field galaxy cluster MACS J0416. The target galaxy is overall a typical compact and UV-faint ($\rm M_{UV}=-17.8$) Lyman-$\alpha$ emitter; yet, the large magnification allows the detailed characterisation of structures on sub-galactic (down to few parsec) scales. Prominent optical $\rm H\alpha$, $\rm H\beta$ and [OIII]$\lambda\lambda4959,5007$ lines are spatially resolved with the high spectral resolution grating (G395H, R~2700), with large equivalent widths, EW($\rm H\beta$+[OIII])$\gtrsim1000$ \AA, and elevated ionising photon production efficiencies $\rm log(\xi_{ion}/erg^{-1}Hz)=25.2-25.7$. NIRCam deep imaging reveals the presence of compact rest-UV bright regions along with individual star clusters of sizes $\rm R_{eff}=3-8~pc$ and masses $\rm M\sim2\cdot10^5-5\cdot10^{6}~M_\odot$ These clusters are characterised by steep UV slopes, $\rm\beta_{UV}\lesssim-2.5$, in some cases associated with a dearth of line emission, indicating possible leaking of the ionising radiation, as also supported by a Lyman-$\rm \alpha$ emission peaking at $\rm \sim100~km~s^{-1}$ from the systemic redshift. While the entire system is characterised by low-metallicity, $\sim0.1~Z_\odot$, the NIRSpec-IFU map also reveals the presence of a low-luminosity, metal-poor region with $\rm Z\lesssim2\%~Z_\odot$, barely detected in NIRCam imaging; this region is displaced by $\rm >200~pc$ from one of the UV brightest structures of the system, and it would have been too faint to detect if not for the large magnification of the system., Comment: 19 pages (11 figures, 2 tables) + appendix (3 pages, 4 figures, 1 table). Submitted to A&A; comments are welcome

We report the serendipitous discovery of a faint (M_UV > -12.2), low-metallicity (Z ~ 0.02 Zsun), ionizing source (dubbed T2c) with a spectroscopic redshift of z=6.146. T2c is part of a larger structure amplified by the Hubble Frontier Field galaxy cluster MACSJ0416, and was observed with JWST/NIRSpec IFU. Stacking the short-wavelength NIRCam data reveals no stellar continuum detection down to a magnitude limit of m_UV ~ 31.0 (3 sigma). However, prominent Hb, [OIII]4959,5007, and Ha emissions are detected, with equivalent widths exceeding 200A, 800A, and 1300A (3 sigma), respectively. The corresponding intrinsic (magnification-corrected x23 +/- 3) ultraviolet and optical rest-frame magnitudes exceed 34.4 and 33.9 (corresponding to M_uv and M_opt fainter than -12.2 and -12.8, at lambda_rest ~ 2000A and ~5000A, respectively), suggesting a stellar mass lower than a few 10^4 Msun under an instantaneous burst scenario. The inferred ionizing photon production efficiency (xi_ion) is high, xi_ion >~ 26.08(25.86) 3(5)sigma, assuming no dust attenuation and no Lyman continuum leakage, indicating the presence of massive stars despite the low mass of the object. The very poor sampling of the initial mass function at such low mass star-forming complex suggests that the formation of very massive stars might be favored in very low metallicity environments. T2c is surrounded by Balmer and weak oxygen emission on a spatial scale of a few hundred parsecs after correcting for lensing effects. This system resembles an HII region potentially powered by currently undetected, extremely efficient, low-metallicity star complexes or clusters. We propose that massive O-type stars populate this low-mass and metallicity high-redshift satellites, likely caught in an early and short formation phase, contributing to the ionization of the surrounding medium., Comment: 9 pages, 5 figures, 1 table. Submitted to A&A. Comments are welcome

Dusty star-forming galaxies (DSFGs) contribute significantly to the stellar buildup at cosmic noon. Major mergers and gas accretion are often invoked to explain DSFGs' prodigious star-formation rates (SFRs) and large stellar masses. We conducted a spatially-resolved morphological analysis of the rest-frame UV/NIR emission in three DSFGs at z~2.5. Initially discovered as CO emitters by NOEMA observations of a bright Herschel source, we observed them with the JWST/NIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies' stellar populations and dust distributions on scales of 250 pc. Spatial variations in stellar mass, SFR, and dust extinction are determined in resolved maps obtained through pixel-based SED fitting. The CO emitters are massive, dusty starburst galaxies with SFRs=340-2500 Msun/yr, positioning them among the most active SFGs at 2

We report ultra-deep, medium-resolution spectroscopic observations for 4 quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot})>11$ at $3 < z < 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by $860$ pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and $3.19$, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha-$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these 4 galaxies are consistent with the $\Lambda$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $\Lambda$-CDM cosmology., Comment: 25 pages, 10 figures, accepted for publication in MNRAS

Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at $z=1.783$, and is the first SN Ia with sufficient light curve sampling and long enough time delays for an $H_0$ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry we measure time delays of $\Delta t_{ab}=-116.6^{+10.8}_{-9.3}$ and $\Delta t_{cb}=-48.6^{+3.6}_{-4.0}$ observer-frame days relative to the last image to arrive (image 2b; all uncertainties are $1\sigma$), which corresponds to a $\sim5.6\%$ uncertainty contribution for $H_0$ assuming $70 \rm{km s^{-1} Mpc^{-1}}$. We also constrain the absolute magnification of each image to $\mu_{a}=4.3^{+1.6}_{-1.8}$, $\mu_{b}=7.6^{+3.6}_{-2.6}$, $\mu_{c}=6.4^{+1.6}_{-1.5}$ by comparing the observed peak near-IR magnitude of SN H0pe to the non-lensed population of SNe Ia., Comment: Accepted in ApJ

The Balmer decrement (H$\alpha$/H$\beta$) provides a constraint on attenuation, the cumulative effects of dust grains in the ISM. The ratio is a reliable spectroscopic tool for deriving the dust properties of galaxies that determine many different quantities such as star formation rate, metallicity, and SED models. Here we measure independently both the attenuation and H$\alpha$/H$\beta$ of an occulting galaxy pair: VV 191. Attenuation measurements in the visible spectrum (A$_{V,stars}$) from dust maps derived from the F606W filter of HST and the F090W filter of JWST are matched with spaxel-by-spaxel H$\alpha$/H$\beta$ observations from the George and Cynthia Mitchell Spectrograph (GCMS) of the McDonald Observatory. The 0.5 to 0.7 micron bandpass covers the Balmer lines for VV 191. The dust maps of JWST and HST provide the high sensitivity necessary for comparisons and tracking trends of the geometrically favorable galaxy. We present maps and plots of the Balmer lines for the VV 191 galaxy pair and for a specific region highlighting dust lanes for VV 191b in the overlap region. We compute A$_{V, HII}$ from H$\alpha$/H$\beta$ and plot both quantities against A$_{V, stars}$. Our results show that regions with higher dust content, residing closer to the spiral center, dominate ionized gas attenuation, leading to an overestimation of A$_{V, HII}$ by a factor or 2. Further out in the spiral arms, the lower dust content leads to more agreement between the attenuations, indicating lower SFR and larger contribution from older stars to the stellar continuum outside the Petrosian radius., Comment: 21 pages, 14 figures, 2 tables, accepted for publication by The Astronomical Journal

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: Accepted for publication in A&A. We updated the manuscript following referee's suggestions

We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $8.58.5)>0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq4$ magnitudes in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($\rho_{\rm{SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the observed early evolution of the galaxy UV LF (and $\rho_{\rm{SFR}}$) can be reproduced in a ${\rm \Lambda}$CDM Universe, with no change in dust properties or star-formation efficiency required out to $z \simeq 12$. Instead, a progressive trend towards younger stellar population ages can reproduce the observations, and the typical ages required at $z \simeq$ 8, 9, 10, and 11 all converge on $\simeq 380-330$ Myr after the Big Bang, indicative of a rapid emergence of early galaxies at $z \simeq 12 - 13$. This is consistent with the first indications of a steeper drop-off in $\rho_{\rm{SFR}}$ we find beyond $z \simeq 13$, possibly reflecting the rapid evolution of the halo mass function at earlier times., Comment: 15 pages, 9 figures, accepted for publication in MNRAS

We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distributions of luminous sources. We fit S\'ersic profiles to eight NIRCam bands and also measure a non-parametric half-light radius. We find that the size distributions with both measurements are well-fit by a log-normal distribution at all redshifts, consistent with disk formation models where size is governed by host dark-matter halo angular momentum. We find a size-redshift evolution of $R_{e} = 3.51(1+z)^{-0.60\pm0.22}$ kpc, in agreement with JWST studies. When considering the typical (modal) size over $z=3-5$, we find little evolution with bright LBGs remaining compact at $R_{e}\simeq0.7-0.9$ kpc. Simultaneously, we find evidence for a build-up of large ($R_{e} > 2$ kpc) galaxies by $z=3$. We find some evidence for a negatively sloped size-mass relation at $z=5$ when S\'ersic profiles are used to fit the data in F200W. The intrinsic scatter in our size-mass relations increases at higher redshifts. Additionally, measurements probing the rest-UV (F200W) show larger intrinsic scatter than those probing the rest-optical (F356W). Finally, we leverage rest-UV and rest-optical photometry to show that disky galaxies are well established by $z=5$, but are beginning to undergo dissipative processes, such as mergers, by $z=3$. The agreement of our size-mass and size-luminosity relations with simulations provides tentative evidence for centrally concentrated star formation at high-redshift., Comment: 15 pages, 9 figures, 3 tables. Accepted to MNRAS. 2 figures added in appendix

The Time Domain Field (TDF) near the North Ecliptic Pole in JWST's continuous-viewing zone will become a premier "blank field" for extragalactic science. JWST/NIRCam data in a 16 arcmin$^2$ portion of the TDF identify 4.4 $\mu$m counterparts for 62 of 63 3 GHz sources with S(3 GHz) > 5 {\mu}Jy. The one unidentified radio source may be a lobe of a nearby Seyfert galaxy, or it may be an infrared-faint radio source. The bulk properties of the radio-host galaxies are consistent with those found by previous work: redshifts range from 0.14 to 4.4 with a median redshift of 1.33. The radio emission arises primarily from star formation in $\sim 2/3$ of the sample and from an active galactic nucleus in $\sim 1/3$, but just over half the sample shows evidence for an AGN either in the spectral energy distribution or by radio excess. All but three counterparts are brighter than magnitude 23 AB at 4.4 $\mu$m, and the exquisite resolution of JWST identifies correct counterparts for sources for which observations with lower angular resolution would mis-identify a nearby bright source as the counterpart when the correct one is faint and red. Up to 11% of counterparts might have been unidentified or misidentified absent NIRCam observations., Comment: Accepted for publication in ApJ. V2 adds an author and some acknowledgments inadvertently omitted

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

Supernova (SN) H0pe was discovered as a new transient in James Webb Space Telescope (JWST) NIRCam images of the galaxy cluster PLCK G165.7+67.0 taken as part of the "Prime Extragalactic Areas for Reionization and Lensing Science" (PEARLS) JWST GTO program (# 1176) on 2023 March 30 (AstroNote 2023-96; Frye et al. 2023). The transient is a compact source associated with a background galaxy that is stretched and triply-imaged by the cluster's strong gravitational lensing. This paper reports spectra in the 950-1370 nm observer frame of two of the galaxy's images obtained with Large Binocular Telescope (LBT) Utility Camera in the Infrared (LUCI) in longslit mode two weeks after the \JWST\ observations. The individual average spectra show the [OII] doublet and the Balmer and 4000 Angstrom breaks at redshift z=1.783+/-0.002. The CIGALE best-fit model of the spectral energy distribution indicates that SN H0pe's host galaxy is massive (Mstar~6x10^10 Msun after correcting for a magnification factor ~7) with a predominant intermediate age (~2 Gyr) stellar population, moderate extinction, and a magnification-corrected star formation rate ~13 Msun/yr, consistent with being below the main sequence of star formation. These properties suggest that H0pe might be a type Ia SN. Additional observations of SN H0pe and its host recently carried out with JWST (JWST-DD-4446; PI: B. Frye) will be able to both determine the SN classification and confirm its association with the galaxy analyzed in this work., Comment: 6 pages, 4 figures, Letter accepted for publication in Astronomy & Astrophysics

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

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 present the $NuSTAR$ extragalactic survey of the $James$ $Webb$ $Space$ $Telescope$ ($JWST$) North Ecliptic Pole (NEP) Time-Domain Field. The survey covers a $\sim$0.16 deg$^2$ area with a total exposure of 681 ks acquired in a total of nine observations from three epochs. The survey sensitivities at 20% of the area are 2.39, 1.14, 2.76, 1.52, and 5.20 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 3-24, 3-8, 8-24, 8-16, and 16-24 keV bands, respectively. The NEP survey is one of the most sensitive extragalactic surveys with $NuSTAR$ so far. A total of 33 sources were detected above 95% reliability in at least one of the five bands. We present the number counts, log$N$-log$S$, measured in the hard X-ray 8-24 and 8-16 keV bands, uniquely accessible by $NuSTAR$ down to such faint fluxes. We performed source detection on the XMM-$Newton$ and $Chandra$ observations of the same field to search for soft X-ray counterparts of each $NuSTAR$ detection. The soft band positions were used to identify optical and infrared associations. We present the X-ray properties (hardness ratio and luminosity) and optical-to-X-ray properties of the detected sources. The measured fraction of candidate Compton-thick (N$\rm _H\ge10^{24} cm^{-2}$) active galactic nuclei, derived from the hardness ratio, is between 3% to 27%. As this survey was designed to have variability as its primary focus, we present preliminary results on multi-epoch flux variability in the 3-24 keV band., Comment: 22 pages, 17 figures

Stellar masses in future James Webb Space Telescope (JWST) deep blank-fields will be mainly derived fitting the spectral energy distribution with theoretical galaxy templates. We investigate the uncertainties and biases of the stellar masses derived by using the LePhare code for spectral energy distribution fitting and the Yggdrasil theoretical templates. We consider a sample of mock galaxies at z$=$7-10 with mock JWST observations with S/N$_{F150W}\geqslant$ 10. Our goal is to provide a list of statistical stellar mass corrections to include on the stellar mass derivation for different output galaxy properties and JWST filter combinations to correct for template degeneracies. Median statistical stellar mass corrections vary from -0.83 dex to 0.87 dex, while 25$\%$ (75$\%$) quartiles range from -0.83 dex (-0.67 dex) to 0.51 dex (0.88 dex), depending on filter combinations and galaxy models. The most challenging cases are galaxies with nebular emission lines, especially the ones that are wrongly identified as galaxies without, relative dust-free galaxies and galaxies with small metallicities (i.e. Z$=1/50Z_{\odot}$). The stellar mass estimation of galaxies correctly identified without emission lines is generally fine, except at z$=$10 when considering only the 8 NIRCam bands, which make the MIRI bands very valuable. We have tested our stellar mass corrections using the public JAGUAR galaxy catalogue, deriving that the average discrepancy in the recovered stellar mass distribution decreases by 20-50$\%$ at z$>$7 after the correction. We found that without the stellar-mass corrections the number of low-mass galaxies (M$^{*}<10^{7}M_{\odot}$) is overestimated, which can potentially lead to systematic errors in the calculation of the galaxy stellar mass function faint-end slope at high z., Comment: 22 pages, 10 figures, 7 tables, accepted for publication at the ApJS

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role., Comment: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584

While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN., Comment: 13 pages, 7 figures, 1 table, Accepted for publication in the Astrophysical Journal

The analysis of galaxies on the star formation rate - stellar mass (SFR-M*) plane is a powerful diagnostic for galaxy evolution at different cosmic times. We consider a sample of 24463 galaxies from the CANDELS/GOODS-S survey to conduct a detailed analysis of the SFR-M* relation at redshifts 0.5$\leqslant z<$3 over more than three dex in stellar mass. To obtain SFR estimates, we utilise mid- and far-IR photometry when available, and rest-UV fluxes for all the other galaxies. We perform our analysis in different redshift bins, with two different methods: 1) a linear regression fitting of all star-forming galaxies, defined as those with specific star formation rates $\rm log_{10}(sSFR/yr^{-1}) > -9.8$, similarly to what is typically done in the literature; 2) a multi-Gaussian decomposition to identify the galaxy main sequence (MS), the starburst sequence and the quenched galaxy cloud. We find that the MS slope becomes flatter when higher stellar mass cuts are adopted, and that the apparent slope change observed at high masses depends on the SFR estimation method. In addition, the multi-Gaussian decomposition reveals the presence of a starburst population which increases towards low stellar masses and high redshifts. We find that starbursts make up ~5% of all galaxies at z=0.5-1.0, while they account for ~16% of galaxies at 2$

Imaging with the James Webb Space Telescope (JWST) will allow for observing the bulk of distant galaxies at the epoch of reionisation. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR) and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broad-band filters, from 0.6 {\mu}m to 7.7 {\mu}m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z=7-10. We found that, with only 8 NIRCam broad-bands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Besides, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z=7-9. Instead, at z=10, no NIRCam band traces purely the {\lambda}> 4000 {\AA} regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z=9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z=10. When nebular emission lines are present, deriving correct galaxy properties is challenging, at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations improves only marginally the estimation., Comment: 21 pages, 11 figures, 4 tables. Accepted for publication at the ApJS

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields., Comment: Accepted for publication in A&A, 16 pages, 6 figures

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates., Comment: 28 Pages, 22 Figures and 10 Tables. Accepted for Publication in the ApJS

We present quantified visual morphologies of approximately 48,000 galaxies observed in three Hubble Space Telescope legacy fields by the Cosmic And Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and classified by participants in the Galaxy Zoo project. 90% of galaxies have z < 3 and are observed in rest-frame optical wavelengths by CANDELS. Each galaxy received an average of 40 independent classifications, which we combine into detailed morphological information on galaxy features such as clumpiness, bar instabilities, spiral structure, and merger and tidal signatures. We apply a consensus-based classifier weighting method that preserves classifier independence while effectively down-weighting significantly outlying classifications. After analysing the effect of varying image depth on reported classifications, we also provide depth-corrected classifications which both preserve the information in the deepest observations and also enable the use of classifications at comparable depths across the full survey. Comparing the Galaxy Zoo classifications to previous classifications of the same galaxies shows very good agreement; for some applications the high number of independent classifications provided by Galaxy Zoo provides an advantage in selecting galaxies with a particular morphological profile, while in others the combination of Galaxy Zoo with other classifications is a more promising approach than using any one method alone. We combine the Galaxy Zoo classifications of "smooth" galaxies with parametric morphologies to select a sample of featureless disks at 1 < z < 3, which may represent a dynamically warmer progenitor population to the settled disk galaxies seen at later epochs., Comment: 30 pages, 14 figures, MNRAS accepted. Galaxy classifications available at data.galaxyzoo.org

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events., Comment: submitted to ApJ; 18 pages; see http://www.stsci.edu/hst/campaigns/frontier-fields/ and http://ssc.spitzer.caltech.edu/warmmission/scheduling/approvedprograms/ddt/frontier/ for data and more information

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization., Comment: Published in ApJL; matches published version

The determination of galaxy redshifts in James Webb Space Telescope (JWST)'s blank-field surveys will mostly rely on photometric estimates, based on the data provided by JWST's Near-Infrared Camera (NIRCam) at 0.6-5.0 {\mu}m and Mid Infrared Instrument (MIRI) at {\lambda}>5.0 {\mu}m. In this work we analyse the impact of choosing different combinations of NIRCam and MIRI broad-band filters (F070W to F770W), as well as having ancillary data at {\lambda}<0.6 {\mu}m, on the derived photometric redshifts (zphot) of a total of 5921 real and simulated galaxies, with known input redshifts z=0-10. We found that observations at {\lambda}<0.6 {\mu}m are necessary to control the contamination of high-z samples by low-z interlopers. Adding MIRI (F560W and F770W) photometry to the NIRCam data mitigates the absence of ancillary observations at {\lambda}<0.6 {\mu}m and improves the redshift estimation. At z=7-10, accurate zphot can be obtained with the NIRCam broad bands alone when S/N>=10, but the zphot quality significantly degrades at S/N<=5. Adding MIRI photometry with one magnitude brighter depth than the NIRCam depth allows for a redshift recovery of 83-99%, depending on SED type, and its effect is particularly noteworthy for galaxies with nebular emission. The vast majority of NIRCam galaxies with [F150W]=29 AB mag at z=7-10 will be detected with MIRI at [F560W, F770W]<28 mag if these sources are at least mildly evolved or have spectra with emission lines boosting the mid-infrared fluxes., Comment: 29 pages, 27 figures, 12 tables. Accepted for publication at the ApJS

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content., Comment: Accepted for publication in MNRAS. 22 pages, 18 figures. Data products are available at http://www.eso.org/~gpopping/Gergo_Poppings_Homepage/Data.html

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved., Comment: 23 pages, 19 figures, accepted by ApJS

We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$)., Comment: 15 pages, 14 figures

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored., Comment: Accepted for publication by ApJ. Replaced to match the final version (two tables and one figure added in the appendix)

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.54. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts., Comment: 26 pages, 23 figures. Accepted for publication on A&A

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged)., Comment: Accepted for publication on Astronomy and Astrophysics

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5., Comment: 25 pages (double-column), 17 Figures, under review at A&A

We present the results of a study investigating the sizes and morphologies of redshift 4 < z < 8 galaxies in the CANDELS GOODS-S, HUDF and HUDF parallel fields. Based on non-parametric measurements and incorporating a careful treatment of measurement biases, we quantify the typical size of galaxies at each redshift as the peak of the log-normal size distribution, rather than the arithmetic mean size. Parameterizing the evolution of galaxy half-light radius as $r_{50} \propto (1+z)^n$, we find $n = -0.20 \pm 0.26$ at bright UV-luminosities ($0.3L_{*(z=3)} < L < L_*$) and $n = -0.47 \pm 0.62$ at faint luminosities ($0.12L_* < L < 0.3L_*$). Furthermore, simulations based on artificially redshifting our z~4 galaxy sample show that we cannot reject the null hypothesis of no size evolution. We show that this result is caused by a combination of the size-dependent completeness of high-redshift galaxy samples and the underestimation of the sizes of the largest galaxies at a given epoch. To explore the evolution of galaxy morphology we first compare asymmetry measurements to those from a large sample of simulated single S\'ersic profiles, in order to robustly categorise galaxies as either `smooth' or `disturbed'. Comparing the disturbed fraction amongst bright ($M_{UV} \leq -20$) galaxies at each redshift to that obtained by artificially redshifting our z~4 galaxy sample, while carefully matching the size and UV-luminosity distributions, we find no clear evidence for evolution in galaxy morphology over the redshift interval 4 < z < 8. Therefore, based on our results, a bright ($M_{UV} \leq -20$) galaxy at z~6 is no more likely to be measured as `disturbed' than a comparable galaxy at z~4, given the current observational constraints., Comment: 29 pages, 25 figures, 4 tables, published in MNRAS

The formation of bars in disk galaxies is a tracer of the dynamical maturity of the population. Previous studies have found that the incidence of bars in disks decreases from the local Universe to z ~ 1, and by z > 1 simulations predict that bar features in dynamically mature disks should be extremely rare. Here we report the discovery of strong barred structures in massive disk galaxies at z ~ 1.5 in deep rest-frame optical images from CANDELS. From within a sample of 876 disk galaxies identified by visual classification in Galaxy Zoo, we identify 123 barred galaxies. Selecting a sub-sample within the same region of the evolving galaxy luminosity function (brighter than L*), we find that the bar fraction across the redshift range 0.5< z < 2 (f_bar = 10.7 +6.3 -3.5% after correcting for incompleteness) does not significantly evolve. We discuss the implications of this discovery in the context of existing simulations and our current understanding of the way disk galaxies have evolved over the last 11 billion years., Comment: 11 pages, 6 figures; accepted for publication in MNRAS

We have constructed a mass-selected sample of Mstar>10^11Msolar galaxies at 1

We present the results of a new and improved study of the morphological and spectral evolution of massive galaxies over the redshift range 110^11 Msolar from the CANDELS WFC3/IR imaging within the COSMOS and UKIDSS UDS survey fields. We find that, by modelling the H(160) image of each galaxy with a combination of a de Vaucouleurs bulge (Sersic index n=4) and an exponential disk (n=1), we can then lock all derived morphological parameters for the bulge and disk components, and successfully reproduce the shorter-wavelength J(125), i(814), v(606) HST images simply by floating the magnitudes of the two components. This then yields sub-divided 4-band HST photometry for the bulge and disk components which, with no additional priors, is well described by spectrophotometric models of galaxy evolution. Armed with this information we are able to properly determine the masses and star-formation rates for the bulge and disk components, and find that: i) from z=3 to z=1 the galaxies move from disk-dominated to increasingly bulge-dominated, but very few galaxies are pure bulges/ellipticals by z=1; ii) while most passive galaxies are bulge-dominated, and most star-forming galaxies disk-dominated, 18+/-5% of passive galaxies are disk-dominated, and 11+/-3% of star-forming galaxies are bulge-dominated, a result which needs to be explained by any model purporting to connect star-formation quenching with morphological transformations; iii) there exists a small but significant population of pure passive disks, which are generally flatter than their star-forming counterparts (whose axial ratio distribution peaks at b/a~0.7); iv) flatter/larger disks re-emerge at the highest star-formation rates, consistent with recent studies of sub-mm galaxies, and with the concept of a maximum surface-density for star-formation activity., Comment: 38 pages, 35 figures. Submitted to MNRAS

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 03x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times., Comment: Published in ApJ. Structural parameter measurements publicly available at http://www.mpia-hd.mpg.de/homes/vdwel/3dhstcandels.html