Search

We analyze a sample of 25 [Ne v] (λ3426) emission-line galaxies at 1.4 < z < 2.3 using Hubble Space Telescope/Wide Field Camera 3 G102 and G141 grism observations from the CANDELS Lyα Emission at Reionization (CLEAR) survey. [Ne v] emission probes extremely energetic photoionization (creation potential of 97.11 eV) and is often attributed to energetic radiation from active galactic nuclei (AGNs), shocks from supernovae, or an otherwise very hard ionizing spectrum from the stellar continuum. In this work, we use [Ne v] in conjunction with other rest-frame UV/optical emission lines ([O ii] λ λ3726, 3729, [Ne iii] λ3869, Hβ, [O iii] λ λ4959, 5007, Hα+[N ii] λ λ6548, 6583, [S ii] λ λ6716, 6731), deep (2–7 Ms) X-ray observations (from Chandra), and mid-infrared imaging (from Spitzer) to study the origin of this emission and to place constraints on the nature of the ionizing engine. The majority of the [Ne v]-detected galaxies have properties consistent with ionization from AGNs. However, for our [Ne v]-selected sample, the X-ray luminosities are consistent with local (z ≲ 0.1) X-ray-selected Seyferts, but the [Ne v] luminosities are more consistent with those from z ∼ 1 X-ray-selected QSOs. The excess [Ne v] emission requires either reduced hard X-rays or a ∼0.1 keV excess. We discuss possible origins of the apparent [Ne v] excess, which could be related to the “soft (X-ray) excess” observed in some QSOs and Seyferts and/or be a consequence of a complex/anisotropic geometry for the narrow-line region, combined with absorption from a warm, relativistic wind ejected from the accretion disk. We also consider implications for future studies of extreme high-ionization systems in the epoch of reionization (z ≳ 6) with the James Webb Space Telescope.

We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at $z = 2.696$. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at $z > 2$. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 $\mu$m image ($0.13''$ resolution) and subtracted, in the visibility plane, the central starburst disk from ALMA rest-frame 240 $\mu$m observations ($0.03''$ resolution). The residual images revealed substructures at rest-frame 1.2 $\mu$m co-located with those found at rest-frame 240 $\mu$m, $\simeq 2$ kpc away from the galactic center. The largest substructure contains $\simeq20$% of the total stellar mass and $\simeq5$% of the total SFR of the galaxy. While UDF2 exhibits a kinematically-ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region ($N \geqslant 4$ massive galaxies within 70 kpc projected distance at $z=2.690-2.697$) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy., Comment: 7 pages, 5 figures; ApJL, accepted

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.

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

We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$\lambda$5008/H$\beta$ ratios, similar [OIII]$\lambda$4364/H$\gamma$ ratios, and higher ($\sim$0.5 dex) [NeIII]$\lambda$3870/[OII]$\lambda$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$\lambda$3870/[OII]$\lambda$3728, [OIII]$\lambda$4364/H$\gamma$, and [OIII]$\lambda$5008/H$\beta$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$\lambda$4364/H$\gamma$ and [OIII]$\lambda$(4960+5008)/H$\beta$ line ratios indicate very high electron temperatures of $4.1, Comment: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

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

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

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3, Comment: Published in CEERS ApJL Focus Issue, ApJL 946, L16

We use the SED-fitting code Prospector to reconstruct the nonparametric star formation history (SFH) of massive ($\log M_*>10.3$) star-forming galaxies (SFGs) and quiescent galaxies (QGs) at redshift $z_{\rm{obs}}\sim2$ to investigate the joint evolution of star-formation activity and structural properties. We find significant correlations between the SFH of the galaxies and their morphology. Compared to extended SFGs, compact SFGs are more likely to have experienced multiple star-formation episodes, with the fractional mass formed during the older ($\ge1$ Gyr) episode being larger, suggesting that high-redshift SFGs assembled their central regions earlier and then kept growing in central mass as they become more compact. The SFH of compact QGs does not significantly differ from the average for this category, and shows an early burst followed by a gradual decline of the star formation rate. The SFH of extended QGs, however, is similar to that of post-starburst galaxies and their morphology is also frequently disturbed. Knowledge of the SFH also enables us to empirically reconstruct the structural evolution of individual galaxies. While the progenitor effect is clearly observed and accounted for in our analysis, it alone is insufficient to explain the observed structural evolution. We show that, as they evolve from star-forming phase to quiescence, galaxies grow massive dense stellar cores. Quenching begins at the center and then propagates outward to the rest of the structure. We discuss possible physical scenarios for the observed evolution and find that our empirical constraints are in good quantitative agreement with the model predictions from dissipative accretion of gas to the center followed by massive starbursts before final quiescence (wet compaction)., 31 pages, 20 figures; Accepted for publication in ApJ

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

The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape., ApJS, accepted. 28 pages, 19 figures, 8 tables. Full machine readable tables will be made available by the publisher at the time of publication

We study the relationship between the morphology and star formation history (SFH) of 361 quiescent galaxies (QGs) at redshift $\langle z_{obs}\rangle\approx 2$, with stellar mass $\log M_*\ge10.3$, selected with the UVJ technique. Taking advantage of panchromatic photometry covering the rest-frame UV-to-NIR spectral range ($\approx40$ bands), we reconstruct the non-parametric SFH of the galaxies with the fully Bayesian SED fitting code Prospector. We find that the half-light radius $R_e$, observed at $z_{obs}$, depends on the formation redshift of the galaxies, $z_{form}$, and that this relationship depends on stellar mass. At $\log M_*11$, the relationship between $R_e$ and $z_{form}$ flattens, suggesting that mergers become increasingly important for the size growth of more massive galaxies after they quenched. We also find that the relationship between $z_{form}$ and galaxy compactness similarly depends on stellar mass. While no clear trend is observed for QGs with $\log M_*>11$, lower-mass QGs that formed earlier, i.e. with larger $z_{form}$, have larger central stellar mass surface densities, both within the $R_e$ ($\Sigma_e$) and central 1 kpc ($\Sigma_{1kpc}$), and also larger $M_{1kpc}/M_*$, the fractional mass within the central 1 kpc. These trends between $z_{form}$ and compactness, however, essentially disappear, if the progenitor effect is removed by normalizing the stellar density with the cosmic density at $z_{form}$. Our findings highlight the importance of reconstructing the SFH of galaxies before attempting to infer their intrinsic structural evolution., Comment: 28 pages, 29 figures; published in ApJ

A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1, Comment: Accepted for publication in MNRAS - This is the final version

Using analytic modeling and simulations, we address the origin of an abundance of star-forming, clumpy, extended gas rings about massive central bodies in massive galaxies at $z \!\! 10^{9.5-10} M_\odot$ where merger-driven spin flips and supernova feedback are ineffective. The rings survive after events of compaction to central nuggets. Ring longevity was unexpected based on inward mass transport driven by torques from violent disc instability. However, evaluating the torques from a tightly wound spiral structure, we find that the timescale for transport per orbital time is long and $\propto\! \delta_{\rm d}^{-3}$, with $\delta_{\rm d}$ the cold-to-total mass ratio interior to the ring. A long-lived ring forms when the ring transport is slower than its replenishment by accretion and the interior depletion by SFR, both valid for $\delta_{\rm d} \!, Comment: 40 pages, 30 figures

We present the results of a new study investigating the relationship between observed Ly α equivalent width (Wλ(Ly α)) and the metallicity of the ionizing stellar population ( Z⋆) for a sample of 768 star-forming galaxies at 3 ≤ z ≤ 5 drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame Wλ(Ly α) across the range $-58 \,\rm {\mathring{\rm A}} \lesssim$ Wλ(Ly α) $\lesssim 110 \,\rm {\mathring{\rm A}}$ , we determine Z⋆ from full spectral fitting of composite far-ultraviolet spectra and find a clear anticorrelation between Wλ(Ly α) and Z⋆. Our results indicate that Z⋆ decreases by a factor ≳ 3 between the lowest Wλ(Ly α) quartile (λ(Ly α) $\rangle =-18\,\rm {\mathring{\rm A}}$ ) and the highest Wλ(Ly α) quartile (λ(Ly α) $\rangle =24\,\rm {\mathring{\rm A}}$ ). Similarly, galaxies typically defined as Lyman alpha emitters (LAEs; Wλ(Ly α) $\gt 20\,\rm {\mathring{\rm A}}$ ) are, on average, metal poor with respect to the non-LAE galaxy population (Wλ(Ly α) $\le 20\,\rm {\mathring{\rm A}}$ ) with Z⋆non-LAE ≳ 2 × Z⋆LAE. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower Z⋆ is responsible for ${\simeq}15{-}25{{\ \rm per\ cent}}$ of the observed variation in Wλ(Ly α) across our sample, with the remaining contribution ( ${\simeq}75{-}85{{\ \rm per\ cent}}$ ) being due to a decrease in the H I/dust covering fractions in low- Z⋆ galaxies.

We use deep spectroscopy from the Hubble Space Telescope (HST) Wide-Field-Camera 3 (WFC3) IR grisms combined with broad-band photometry to study the stellar populations, gas ionization and chemical abundances in star-forming galaxies at $z\sim 1.1-2.3$. The data stem from the CANDELS Lyman-$\alpha$ Emission At Reionization (CLEAR) survey. At these redshifts the grism spectroscopy measure the [OII] 3727, 3729, [OIII] 4959, 5008, H-$\beta$ strong emission features, which constrain the ionization parameter and oxygen abundance of the nebular gas. We compare the line flux measurements to predictions from updated photoionization models (MAPPINGS (Kewley et al. 2019), which include an updated treatment of nebular gas pressure, log P/k = $n_e T_e$. Compared to low-redshift samples ($z\sim 0.2$) at fixed stellar mass, llog M / M$_\odot$ = 9.4-9.8, the CLEAR galaxies at z=1.35 (z=1.90) have lower gas-phase metallicity, $\Delta$(log Z) = 0.25 (0.35) dex, and higher ionization parameters, $\Delta$(log q) = 0.25 (0.35) dex, where U = q/c. We provide updated analytic calibrations between the [OIII], [OII], and H-$\beta$ emission line ratios, metallicity, and ionization parameter. The CLEAR galaxies show that at fixed stellar mass, the gas ionization parameter is correlated with the galaxy specific star-formation rates (sSFRs), where $\Delta$ log q = 0.4 $\Delta$(log sSFR), derived from changes in the strength of galaxy H-$\beta$ equivalent width. We interpret this as a consequence of higher gas densities, lower gas covering fractions, combined with higher escape fraction of H-ionizing photons. We discuss both tests to confirm these assertions and implications this has for future observations of galaxies at higher redshifts., Comment: 37 pages, plethora of figures, accepted for publication in ApJ. Updated with accepted version

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 4.91 − 1.47 + 0.91 Gyr, consistent with the slow processes of strangulation and delayed-then-rapid quenching although more data are needed to confirm this result.

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

We present a study of spatially-resolved star formation histories (SFHs) for 60 $z\sim2.3$ main-sequence, star-forming galaxies selected from the MOSDEF spectroscopic survey in the GOODS-N field, with median stellar mass $\log(M_\star/M_\odot)=9.75$ and spanning the range $8.6, 28 pages, 21 figures. Accepted for publication in ApJ

We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M$_{\ast}>10^{10}\,$M$_{\odot}$) simulated galaxies at $110^{11}\,$M$_{\odot}$ galaxy. In particular, we explore the potential of HST and JWST datasets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this dataset for characterizing the first episodes in the SFH of local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies by studying their progenitors at $z>1$. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1% to 50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%-30%. According to our measurements on Illustris data, we are able to recover that local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies would have started their formation by $z=16$, forming the first 5% of their stellar mass present at $z \sim 1$ by $z=4.5$, 10% by $z=3.7$, and 25% by $z=2.7$., Comment: 28 pages, 13 figures, 4 tables. ApJ in press. Summary of changes from original submission: the major change is that we now include in Sec. 6 the comparison of the results obtained for our sample of massive 1 < z < 4 progenitors with those obtained by considering all massive galaxies at 1 < z < 4 in the simulated images. Several figures and sections have been updated

We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample of $i, Comment: 15 pages, 14 figures, accepted for publication in ApJ

We present a detailed stellar population analysis of 11 bright ($H10$., Comment: 27 pages + Appendix; accepted by ApJ

Using spatially resolved Hα emission line maps of star-forming galaxies, we study the spatial distribution of star formation over a wide range in redshift (0.5 ≲ z ≲ 1.7). Our z ∼ 0.5 measurements come from deep Hubble Space Telescope (HST) Wide Field Camera 3 G102 grism spectroscopy obtained as part of the CANDELS Lyα Emission at Reionization Experiment. For star-forming galaxies with log(M */M ⊙) ≥ 8.96, the mean Hα effective radius is 1.2 ± 0.1 times larger than that of the stellar continuum, implying inside-out growth via star formation. This measurement agrees within 1σ with those measured at z ∼ 1 and z ∼ 1.7 from the 3D-HST and KMOS3D surveys, respectively, implying no redshift evolution. However, we observe redshift evolution in the stellar mass surface density within 1 kpc (Σ1kpc). Star-forming galaxies at z ∼ 0.5 with a stellar mass of log(M */M ⊙) = 9.5 have a ratio of Σ1kpc in Hα relative to their stellar continuum that is lower by (19 ± 2)% compared to z ∼ 1 galaxies. Σ1kpc,Hα /Σ1kpc,Cont decreases toward higher stellar masses. The majority of the redshift evolution in Σ1kpc,Hα /Σ1kpc,Cont versus stellar mass stems from the fact that log(Σ1kpc,Hα ) declines twice as much as log(Σ1kpc,Cont) from z ∼ 1 to 0.5 (at a fixed stellar mass of log(M */M ⊙) = 9.5). By comparing our results to the TNG50 cosmological magneto-hydrodynamical simulation, we rule out dust as the driver of this evolution. Our results are consistent with inside-out quenching following in the wake of inside-out growth, the former of which drives the significant drop in Σ1kpc,Hα from z ∼ 1 to z ∼ 0.5.

The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a data-driven, flexible theoretical framework to probe the structural evolution of MGs. We assign galaxies to dark matter haloes via stellar mass-halo mass (SMHM) relations with varying high-mass slopes and scatters $\sigma_{\rm SMHM}$ in stellar mass at fixed halo mass, and assign sizes to galaxies using an empirically-motivated, constant and linear relationship between $R_e$ and the host dark matter halo radius $R_h$. We find that: 1) the fast mean size growth of MGs is well reproduced independently of the shape of the input SMHM relation; 2) the numbers of compact MGs grow steadily until $z\gtrsim2$ and fall off at lower redshifts, suggesting a lesser role of progenitor bias at later epochs; 3) a time-independent scatter $\sigma_{\rm SMHM}$ is consistent with a scenario in which compact starforming MGs transition into quiescent MGs in a few $10^8$yr with a negligible structural evolution during the compact phase, while a scatter increasing at high redshift implies significant size growth during the starforming phase. A robust measurement of the size function of MGs at high redshift can set strong constraints on the scatter of the SMHM relation and, by extension, on models of galaxy evolution., Comment: Accepted for publication in MNRAS. Comments still welcome

We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedure, including screening for persistence, stellar contamination, inclusion of ground-based imaging, and followup space-based imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M_UV=-22.3) is consistent with predictions from simulations which assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z=4-8, or an accelerate decline at z > 8. We calculate that the presence of seven galaxies in a single field (EGS) is an outlier at the 2-sigma significance level, implying the discovery of a significant overdensity. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope., 47 figures, 25 pages, 10 tables, Accepted to the Astrophysical Journal. Photometric catalogs available via email request to lead author

Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$\alpha$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$\alpha$ from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at $6.06$. Furthermore, we compare the redshift evolution of $W_0$ between galaxies at different UV luminosities. The UV-bright ($M_{\text{UV}}L^{*}$) galaxies show weaker evolution with a decrease of 0.4 ($\pm$0.2) dex in $W_0$ at $z>6$ while UV-faint ($M_{\text{UV}}>-21$, or $L_{\text{UV}}6$. Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Ly$\alpha$ transmission in the IGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity., Comment: 21 pages, 9 figures, 2 tables, 1st revision submitted to ApJ

Star-forming galaxies are considered the likeliest source of the H i ionizing Lyman continuum (LyC) photons that reionized the intergalactic medium at high redshifts. However, above z ≳ 6, the neutral intergalactic medium prevents direct observations of LyC. Therefore, recent years have seen the development of indirect indicators for LyC that can be calibrated at lower redshifts and applied in the epoch of reionization. Emission from the Mg ii λλ2796, 2803 doublet has been proposed as a promising LyC proxy. In this paper, we present new Hubble Space Telescope/Cosmic Origins Spectrograph observations for eight LyC emitter candidates, selected to have strong Mg ii emission lines. We securely detect LyC emission in 50% (4/8) of the galaxies with 2σ significance. This high detection rate suggests that strong Mg ii emitters might be more likely to leak LyC than similar galaxies without strong Mg ii. Using photoionization models, we constrain the escape fraction of Mg ii as ∼15%–60%. We confirm that the escape fraction of Mg ii correlates tightly with that of Lyα, which we interpret as an indication that the escape fraction of both species is controlled by resonant scattering in the same low column density gas. Furthermore, we show that the combination of the Mg ii emission and dust attenuation can be used to estimate the escape fraction of LyC statistically. These findings confirm that Mg ii emission can be adopted to estimate the escape fraction of Lyα and LyC in local star-forming galaxies and may serve as a useful indirect indicator at the epoch of reionization.

Many X-ray bright active galactic nuclei (AGNs) are predicted to follow an extended stage of obscured black hole growth. In support of this picture we examine the X-ray undetected AGNs in the COSMOS field and compare their host galaxies with X-ray bright AGNs. We examine galaxies with M * > 109.5 M ⊙ for the presence of AGNs at redshifts z = 0.5–3. We select AGNs in the infrared using Spitzer and Herschel detections and use color selection techniques to select AGNs within strongly star-forming hosts. We stack Chandra X-ray data of galaxies with an infrared (IR) detection but lacking an X-ray detection to obtain soft and hard fluxes, allowing us to measure the energetics of these AGNs. We find a clear correlation between X-ray luminosity and IR AGN luminosity in the stacked galaxies. We also find that X-ray undetected AGNs all lie on the main sequence—the tight correlation between the star formation rate and M * that holds for the majority of galaxies, regardless of mass or redshift. This work demonstrates that there is a higher population of obscured AGNs than previously thought.

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5�� allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar), 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

We use Hubble Space Telescope WFC3 G102 and G141 grism spectroscopy to measure rest-optical emission-line ratios of 533 galaxies at $z\sim1.5$ in the CANDELS Ly$\alpha$ Emission at Reionization (CLEAR) survey. We compare $\frac{[OIII]}{H\beta}$ vs. $\frac{[SII]}{(H\alpha +[NII])}$ as an "unVO87" diagram for 461 galaxies and $\frac{[OIII]}{Hb}$ vs. $\frac{[NeIII]}{[OII]}$ as an "OHNO" diagram for 91 galaxies. The unVO87 diagram does not effectively separate active galactic nuclei (AGN) and $[NeV]$ sources from star-forming galaxies, indicating that the unVO87 properties of star-forming galaxies evolve with redshift and overlap with AGN emission-line signatures at $z>1$. The OHNO diagram effectively separates X-ray AGN and $[NeV]$-emitting galaxies from the rest of the population. We find that the $\frac{[OIII]}{H\beta}$ line ratios are significantly anti-correlated with stellar mass and significantly correlated with $\log(L_{H\beta})$, while $\frac{[SII]}{(H\alpha +[NII])}$ is significantly anti-correlated with $\log(L_{H\beta})$. Comparison with MAPPINGS~V photoionization models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We do not find evidence for redshift evolution of the emission-line ratios outside of the correlations with mass and SFR.Our results suggest that the OHNO diagram of $\frac{[OIII]}{Hb}$ vs. $\frac{[NeIII]}{[OII]}$ will be a useful indicator of AGN content and gas conditions in very high-redshift galaxies to be observed by the James Webb Space Telescope., Comment: 20 pages, 22 figures

We use \Pab\ (1282~nm) observations from the Hubble Space Telescope ($\HST$) G141 grism to study the star-formation and dust attenuation properties of a sample of 29 low-redshift ($z < 0.287$) galaxies in the CANDELS Ly$\alpha$ Emission at Reionization (CLEAR) survey. We first compare the nebular attenuation from $\Pab/\Ha$ with the stellar attenuation inferred from the spectral energy distribution, finding that the galaxies in our sample are consistent with an average ratio of the continuum attenuation to the nebular gas of 0.44, but with a large amount of excess scatter beyond the observational uncertainties. Much of this scatter is linked to a large variation between the nebular dust attenuation as measured by (space-based) $\Pab$ to (ground-based) $\Ha$ to that from (ground-based) $\Ha/\Hb$. This implies there are important differences between attenuation measured from grism-based / wide-aperture $\Pab$ fluxes and the ground-based / slit-measured Balmer decrement. We next compare star-formation rates (SFRs) from $\Pab$ to those from dust-corrected UV. We perform a survival analysis to infer a census of \Pab\ emission implied by both detections and non-detections. We find evidence that galaxies with lower stellar mass have more scatter in their ratio of \Pab\ to attenuation-corrected UV SFRs. When considering our \Pab\ detection limits, this observation supports the idea that lower mass galaxies experience "burstier" star-formation histories. Together, these results show that \Pab\ is a valuable tracer of a galaxy's SFR, probing different timescales of star-formation and potentially revealing star-formation that is otherwise missed by UV and optical tracers., Comment: 20 pages, 15 figures, 2 tables, submitted to ApJ

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues., Comment: 18 pages, 13 figures. Accepted for publication in A&A

In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (��_{obs} =151.44GHz with a S/N=8.58) and AGS17 (��_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG)., 17 pages, 13 figures, 3 tables, accepted by A&A

The origin of the correlations between mass, morphology, quenched fraction, and formation history in galaxies is difficult to define, primarily due to the uncertainties in galaxy star-formation histories. Star-formation histories are better constrained for higher redshift galaxies, observed closer to their formation and quenching epochs. Here we use "non-parametric" star-formation histories and a nested sampling method to derive constraints on the formation and quenching timescales of quiescent galaxies at $0.7$10.25, where $\Sigma_1$ is the stellar mass within 1~pkpc (proper kpc). Quiescent galaxies with the highest stellar-mass surface density, $\log\Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) > 10.25$, } show a \textit{minimum} formation redshift: all such objects in our sample have $z_{50} > 2.9$. Quiescent galaxies with lower surface density, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) = 9.5 - 10.25$, show a range of formation epochs ($z_{50} \simeq 1.5 - 8$), implying these galaxies experienced a range of formation and assembly histories. We argue that the surface density threshold $\log\Sigma_1/(M_\odot\ \mathrm{kpc}^{-2})>10.25$ uniquely identifies galaxies that formed in the first few Gyr after the Big Bang, and we discuss the implications this has for galaxy formation models., Comment: 13 pages, 7 figures, accepted for publication in ApJ. Includes an interactive online appendix (https://vince-ec.github.io/appendix/appendix)

UNSTRUCTURED We discovered that the time evolution of the inverse fractional daily growth of new infections, N/dN, in the current outbreak of COVID-19 is accurately described by a universal function, namely the two-parameter Gumbel cumulative function, in all countries that we have investigated. While the two Gumbel parameters, as determined bit fits to the data, vary from country to country (and even within different regions of the same country), reflecting the diversity and efficacy of the adopted containment measures, the functional form of the evolution of N/dN appears to be universal. The result of the fit in a given region or country appears to be stable against variations of the selected time interval. This makes it possible to robustly estimate the two parameters from the data data even over relatively small time periods. In turn, this allows one to predict with large advance and well-controlled confidence levels, the time of the peak in the daily new infections, its magnitude and duration (hence the total infections), as well as the time when the daily new infections decrease to a pre-set value (e.g. less than about 2 new infections per day per million people), which can be very useful for planning the reopening of economic and social activities. We use this formalism to predict and compare these key features of the evolution of the COVID-19 disease in a number of countries and provide a quantitative assessment of the degree of success in in their efforts to countain the outbreak.

The Lya line in the UV and the [CII] line in the FIR are widely used tools to identify galaxies and to obtain insights into ISM properties in the early Universe. By combining data obtained with ALMA in band 7 at ~ 320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at Keck, VIMOS and FORS2 at the VLT, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 < z < 6 in which we detect both the Lya line and the [CII]. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we find that 90% of the selected objects have Lya-[CII] velocity offsets in the range 0 < Dv_Lya-[CII] < 400 km/s, in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower z. We observe ISM-[CII] offsets in the range -500 < Dv_ISM-[CII] < 0 km/s, in line with values at all redshifts. We find significant anticorrelations between Dv_Lya-[CII] and the Lya rest-frame equivalent width EW0(Lya) (or equivalently, the Lya escape fraction f_esc(Lya)). According to available models for the radiative transfer of Lya photons, the escape of Lya photons would be favored in galaxies with high outflow velocities, in agreement with our observations. The uniform shell model would also predict that the Lya escape in galaxies with slow outflows (0 < v_out < 300 km/s) is mainly determined by the neutral hydrogen column density (NHI), while the alternative model by Steidel+10 would favor a combination of NHI and covering fraction as driver of the Lya escape. We suggest that the observed increase in Lya escape that is observed between z~2 and z~6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction. [abridged], Comment: 22 pages, 6 figures, accepted for publication in A&A. This version replaces the previous one

The ALMA-ALPINE [CII] survey is aimed at characterizing the properties of a sample of normal star-forming galaxies (SFGs). The ALMA Large Program to INvestigate (ALPINE) features 118 galaxies observed in the [CII]-158$\mu$m line and far infrared (FIR) continuum emission during the period of rapid mass assembly, right after the end of the HI reionization, at redshifts of 45). Based on a visual inspection of the [CII] data cubes together with the large wealth of ancillary data, we find a surprisingly wide range of galaxy types, including 40% that are mergers, 20% extended and dispersion-dominated, 13% compact, and 11% rotating discs, with the remaining 16% too faint to be classified. This diversity indicates that a wide array of physical processes must be at work at this epoch, first and foremost, those of galaxy mergers. This paper sets a reference sample for the gas distribution in normal SFGs at 4, Comment: 20 pages, 13 figures, 1 table, accepted by A&A

We report the detection of [CII]158um emission from a system of three closely-separated sources in the COSMOS field at z~4.56, as part of the ALMA Large Program to INvestigate CII at Early times (ALPINE). The two dominant sources are closely associated, both spatially (1.6"~11kpc) and in velocity (~100km/s), while the third source is slightly more distant (2.8"~18kpc, ~300km/s). The second strongest source features a slight velocity gradient, while no significant velocity gradient is seen in the other two sources. Using the observed [CII] luminosities, we derive a total log(SFR_[CII]/[Msol/year])=2.8+/-0.2, which may be split into contributions of 59%, 31%, and 10% from the central, east, and west sources, respectively. Comparison of these [CII] detections to recent zoom-in cosmological simulations suggests an ongoing major merger. We are thus witnessing a system in a major phase of mass build-up by merging, including an on-going major merger and an upcoming minor merger, which is expected to end up in a single massive galaxy by z~2.5., 5 pages, 4 figures, submitted to MNRASL

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy 'Ion2' at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the 'Sunburst' arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]4959-5007 / [OII]3727-3729) of 9.18(-1.32/+1.82) is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII1640 appears broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and from VLT/MUSE spectroscopy the presence of spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 20), suggests that this exceptional object is a gravitationally-bound star cluster observed at a cosmological distance, with a stellar mass M, Comment: 13 pages, 9 figures and 1 table, MNRAS accepted. Some typos fixed

We report on the gas-phase metallicity gradients of a sample of 264 star-forming galaxies at 0.6 < z < 2.6, measured through deep near-infrared Hubble Space Telescope slitless spectroscopy. The observations include 12-orbit depth Hubble/WFC3 G102 grism spectra taken as a part of the CANDELS Lya Emission at Reionization (CLEAR) survey, and archival WFC3 G102+G141 grism spectra overlapping the CLEAR footprint. The majority of galaxies (84%) in this sample are consistent with a zero or slightly positive metallicity gradient across the full mass range probed (8.5 < log M_*/M_sun < 10.5). We measure the intrinsic population scatter of the metallicity gradients, and show that it increases with decreasing stellar mass---consistent with previous reports in the literature, but confirmed here with a much larger sample. To understand the physical mechanisms governing this scatter, we search for correlations between the observed gradient and various stellar population properties at fixed mass. However, we find no evidence for a correlation with the galaxy properties we consider---including star-formation rates, sizes, star-formation rate surface densities, and star-formation rates per gravitational potential energy. We use the observed weakness of these correlations to provide material constraints for predicted intrinsic correlations from theoretical models., 19 pages, 10 figures (v2: typo fixed in Figure 10 label); submitted to ApJ

We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization ($5.0\le{}z\le6.6$) from a large sample of 52 galaxies with spectroscopically confirmed redshifts from the VIMOS UltraDeep Survey. We identify galaxies with Ly$\alpha$ either in absorption or in emission, at variance with most spectroscopic samples in the literature where Ly$\alpha$ emitters dominate. We find that star-forming galaxies at these redshifts are distributed along a main sequence in the stellar mass vs. SFR plane. We report a flat evolution of the sSFR(z) in 32. The bright end of the Ly$\alpha$ luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest Ly$\alpha$ emitters (LAE) at these redshifts. LAE with EW>25\AA ~contribute to about 75\% of the total UV-derived SFRD. While our analysis favors a low dust content in 5.0, Comment: 27 pages, submitted to A&A

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v| 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs., Comment: 16 pages, 10 figures. Accepted for publication in A&A

The relationship between galaxy characteristics and the reionization of the universe remains elusive, mainly due to the observational difficulty in accessing the Lyman continuum (LyC) at these redshifts. It is thus important to identify low-redshift LyC-leaking galaxies that can be used as laboratories to investigate the physical processes that allow LyC photons to escape. The weakness of the [S II] nebular emission lines relative to typical star-forming galaxies has been proposed as a LyC predictor. In this paper, we show that the [S II]-deficiency is an effective method to select LyC-leaking candidates using data from the Low-redshift LyC Survey, which has detected flux below the Lyman edge in 35 out of 66 star-forming galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We show that LyC leakers tend to be more [S II]-deficient and that the fraction of their detections increases as [S II]-deficiency becomes more prominent. Correlational studies suggest that [S II]-deficiency complements other LyC diagnostics (such as strong Lyman-$\alpha$ emission and high [O III]/[O II]). Our results verify an additional technique by which reionization-era galaxies could be studied., Comment: 15 pages, 6 figures, 3 tables; published in ApJ

We use the Cosmic Assembly Deep Near-infrared Extragalactic Legacy Survey (CANDELS) data to study the relationship between quenching and the stellar mass surface density within the central radius of 1 kpc ($\Sigma_1$) of low-mass galaxies (stellar mass $M_* \lesssim 10^{9.5} M_\odot$) at $0.5 \leq z < 1.5$. Our sample is mass complete down to $\sim 10^9 M_\odot$ at $0.5 \leq z < 1.0$. We compare the mean $\Sigma_1$ of star-forming galaxies (SFGs) and quenched galaxies (QGs) at the same redshift and $M_*$. We find that low-mass QGs have higher $\Sigma_1$ than low-mass SFGs, similar to galaxies above $10^{10} M_\odot$. The difference of $\Sigma_1$ between QGs and SFGs increases slightly with $M_*$ at $M_* \lesssim 10^{10} M_\odot$ and decreases with $M_*$ at $M_* \gtrsim 10^{10} M_\odot$. The turnover mass is consistent with the mass where quenching mechanisms transition from internal to environmental quenching. At $0.5 \leq z < 1.0$, we find that the $\Sigma_1$ of galaxies increases by about 0.25 dex in the green valley (i.e., the transitioning region from star forming to fully quenched), regardless of their $M_*$. Using the observed specific star formation rate (sSFR) gradient in the literature as a constraint, we estimate that the quenching timescale (i.e., time spent in the transition) of low-mass galaxies is a few ($\sim4$) Gyrs at $0.5 \leq z < 1.0$. The mechanisms responsible for quenching need to gradually quench star formation in an outside-in way, i.e., preferentially ceasing star formation in outskirts of galaxies while maintaining their central star formation to increase $\Sigma_1$. An interesting and intriguing result is the similarity of the growth of $\Sigma_1$ in the green valley between low-mass and massive galaxies, which suggests that the role of internal processes in quenching low-mass galaxies is a question worthy of further investigation., Comment: 12 pages, 11 figures. Accepted by ApJ

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws., 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work., Comment: Re-submitted to MNRAS after minor revisions

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors., Comment: The catalogs, images and spectra are available in Rainbow-slicer (https://bit.ly/2OxKKx1 ), navigator (https://bit.ly/2GDS180 ) and MAST (https://bit.ly/2YtoBQ4 ). In addition to the photometry and other added-value catalogs we release UV+IR star formation rates based on Spitzer (MIPS) and Herschel (PACS and SPIRE) in the 5 CANDELS fields

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 \AA with HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 \AA with the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$\alpha$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$\alpha$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40, Comment: 24 pages, 13 figures, accepted for publication in ApJ

We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.068$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$., Comment: 32 pages, 23 figures, Resubmited to ApJ after revisions in response to referee report