Search

We present a new method to determine the star formation rate (SFR) density of the Universe at $z \gtrsim 5$ that includes the contribution of dust-obscured star formation. For this purpose, we use a [CII] (158 $\mu$m) selected sample of galaxies serendipitously identified in the fields of known $z\gtrsim 4.5$ objects to characterize the fraction of obscured SFR. The advantage of a [CII] selection is that our sample is SFR-selected, in contrast to a UV-selection that would be biased towards unobscured star formation. We obtain a sample of 23 [CII] emitters near star-forming (SF) galaxies and QSOs -- three of which we identify for the first time -- using previous literature and archival ALMA data. 18 of these serendipitously identified galaxies have sufficiently deep rest-UV data and are used to characterize the obscured fraction of the star formation in galaxies with SFRs $\gtrsim 30\ \text{M}_{\odot} \ \text{yr}^{-1}$. We find that [CII] emitters identified around SF galaxies have $\approx$63\% of their SFR obscured, while [CII] emitters around QSOs have $\approx$93\% of their SFR obscured. By forward modeling existing wide-area UV luminosity function (LF) determinations, we derive the intrinsic UV LF using our characterization of the obscured SFR. Integrating the intrinsic LF to $M_{UV}$ = $-$20 we find that the obscured SFRD contributes to $>3\%$ and $>10\%$ of the total SFRD at $z \sim 5$ and $z \sim 6$ based on our sample of companions galaxies near SFGs and QSOs, respectively. Our results suggest that dust obscuration is not negligible at $z\gtrsim 5$, further underlining the importance of far-IR observations of the $z\gtrsim 5$ Universe., Comment: Accepted for publication in MNRAS. 24 pages, 14 figures, 8 tables (including appendices)

We present the census of H$\beta$+[OIII] 4960,5008 \AA emitters at 6.87$ galaxies. We report a rapid decline of the [OIII] luminosity density at $z\gtrsim 6-7$ which cannot be explained by the evolution of the cosmic star-formation rate density. Finally we find that FRESCO detects in only 2h galaxies likely accounting for $\sim 10-20\%$ of the ionising budget at $z=7-8$ (assuming an escape fraction of 10%), raising the prospect of directly detecting a significant fraction of the sources of reionisation with JWST., Comment: 20 pages + appendices. Accepted in MNRAS. Public catalogue release at https://github.com/rameyer/fresco. V3: matching accepted version

The advent of the JWST has revolutionised our understanding of high-redshift galaxies. In particular, the NIRCam instrument on-board JWST has revealed a population of Hubble Space Telescope (HST)-dark galaxies that had previously evaded optical detection, potentially due to significant dust obscuration, quiescence, or simply extreme redshift. Here, we present the first NIRSpec spectra of 23 HST-dark galaxies ($\mathrm{H-F444W>1.75}$), unveiling their nature and physical properties. This sample includes both dusty and quiescent galaxies with spectroscopic data from NIRSpec/PRISM, providing accurate spectroscopic redshifts with $\mathrm{\overline{z}_{spec} = 4.1 \pm 0.7}$. The spectral features demonstrate that, while the majority of HST-dark galaxies are dusty, a substantial fraction, $\mathrm{13^{+9}_{-6} \%}$, are quiescent. For the dusty galaxies, we have quantified the dust attenuation using the Balmer decrement ($\mathrm{H\alpha / H\beta}$), finding attenuations $\mathrm{A_{V} > 2\ mag}$. We find that HST-dark dusty galaxies are $\mathrm{H\alpha}$ emitters with equivalent widths spanning the range $\mathrm{ 68 A < EW_{H\alpha} < 550 A }$, indicative of a wide range of recent star-formation activity. Whether dusty or quiescent, we find that HST-dark galaxies are predominantly massive, with 85\% of the galaxies in the sample having masses $\mathrm{log(M_{*}/M_{\odot}) > 9.8}$. This pilot NIRSpec program reveals the diverse nature of HST-dark galaxies and highlights the effectiveness of NIRSpec/PRISM spectroscopic follow-up in distinguishing between dusty and quiescent galaxies and properly quantifying their physical properties. Upcoming research utilising higher-resolution NIRSpec data and combining JWST with ALMA observations will enhance our understanding of these enigmatic and challenging sources., Comment: 13 pages, 11 figures

With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, making it possible to accurately characterise red, optically-faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour-selected using solely JWST bands. With deep JWST/NIRCam data from 1.15$\mu$m to 4.44$\mu$m and ancillary HST/ACS and WFC3 data, we determine the physical properties of our sample using spectral energy distribution fitting with BAGPIPES. We demonstrate that our selection method efficiently identifies massive ($\mathrm{\langle \log M_\star/M_\odot \rangle \sim 10}$) and dusty ($\mathrm{\langle A_V\rangle \sim 2.7\ mag}$) sources, with a majority at $z>3$ and predominantly lying on the galaxy main-sequence. The main results of this work are the stellar mass functions (SMF) of red, optically-faint galaxies from redshifts between $3

We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($\beta \simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $\beta_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$\beta$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$\beta$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$\beta$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs., Comment: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Submitted to MNRAS, comments welcome

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$., Comment: Accepted for publication in MNRAS. 21 pages, 8 figures

We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $\mu$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$\alpha$ and Pa$\beta$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$\alpha$ and [NII] at z~5-6.5, up to [OIII] and H$\beta$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $\sigma$ in 0.32" diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$\alpha$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations., Comment: 12 pages, 11 figures; MNRAS in press; for more information on the survey and data releases, see http://jwst-fresco.astro.unige.ch/ and https://archive.stsci.edu/hlsp/fresco

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations., Comment: 9 pages, 5 figures

Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies observed with very deep, high-resolution NIRCam imaging from the Early Release Science Program CEERS. 30 HST-dark sources are selected based on their red colours across 1.6 $\mu$m to 4.4 $\mu$m. Their physical properties are derived from 12-band multi-wavelength photometry, including ancillary HST imaging. We find that these galaxies are generally heavily dust-obscured ($A_{V}\sim2$ mag), massive ($\log (M/M_{\odot}) \sim10$), star-forming sources at $z\sim2-8$ with an observed surface density of $\sim0.8$ arcmin$^{-2}$. This suggests that an important fraction of massive galaxies may have been missing from our cosmic census at $z>3$ all the way into the Reionization epoch. The HST-dark sources lie on the main sequence of galaxies and add an obscured star formation rate density (SFRD) of $\mathrm{3.2^{+1.8}_{-1.3} \times 10^{-3} M_{\odot}/yr/Mpc^{3}}$ at $z\sim7$ showing likely presence of dust in the Epoch of Reionization. Our analysis shows the unique power of JWST to reveal this previously missing galaxy population and to provide a more complete census of galaxies at $z=2-8$ based on rest-frame optical imaging., Comment: 9 pages, 7 figures, accepted to MNRAS

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine HI 21-cm transition. Here we estimate the HI gas mass of a sample of main-sequence star-forming galaxies at $z\sim 6.5 - 7.8$ surveyed for [CII]$-158\mu$m emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [CII]-to-HI conversion factor. We find that the HI gas mass excess in galaxies increases as a function of redshift, with an average of $M_{\rm HI} / M_\star \approx 10$, corresponding to HI gas mass fractions of $f_{\rm HI} = M_{\rm HI} / (M_\star + M_{\rm HI}) = 90\%$, at $z\approx 7$. Based on the [CII]-$158\mu$m luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic HI gas mass density in galaxies ($\rho_{\rm HI}$) at this redshift, which we measure to be $\rho_{\rm HI} = 7.1^{+6.4}_{-3.0} \times 10^{6}\,M_{\odot}\,{\rm Mpc^{-3}}$. This estimate is substantially lower by a factor of $\approx 10$ than that inferred from an extrapolation of damped Lyman-$\alpha$ absorber (DLA) measurements, and largely depend on the exact [CII] LF adopted. However, we find this decrease in $\rho_{\rm HI}$ to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sightlines predominantly probing the substantial fraction of HI gas in high-$z$ galactic halos, whereas [CII] traces the HI in the ISM associated with star formation. We make predictions for this build-up of neutral gas in galaxies as a function of redshift, showing that at $z\gtrsim 5$ only $\approx 10\%$ of the cosmic HI gas content is confined in galaxies and associated with the star-forming ISM., Comment: The accepted version of the manuscript in press for ApJL

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4

We present new determinations of the rest-UV luminosity functions (LFs) at z=2-9 to extremely low luminosities (>-14 mag) from a sample of >2500 lensed galaxies found behind the HFF clusters. For the first time, we present faint-end slope results from lensed samples that are fully consistent with blank-field results over the redshift range z=2-9, while reaching to much lower luminosities than possible from the blank-field studies. Combining the deep lensed sample with the large blank-field samples allows us to set the tight constraints on the faint-end slope alpha of the z=2-9 UV LFs and its evolution. We find a smooth flattening in alpha from -2.28+/-0.10 (z=9) to -1.53+/-0.03 (z=2) with cosmic time (d(alpha)/dz=-0.11+/-0.01), fully consistent with dark matter halo buildup. We utilize these new results to present new measurements of the evolution in the UV luminosity density rho(UV) brightward of -13 mag from z~9 to z~2. Accounting for the SFR densities to faint luminosities implied by our LF results, we find that unobscured star formation dominates the SFR density at z>~4, with obscured star formation dominant thereafter. Having shown we can quantify the faint-end slope alpha of the LF accurately with our lensed HFF samples, we also quantify the apparent curvature in the shape of the UV LF through a curvature parameter delta. The constraints on the curvature delta strongly rule out the presence of a turn-over brightward of -13.1 mag at z~3, -14.3 mag at z~6, and -15.5 mag at all other redshifts between z~9 to z~2., Comment: 32 pages, 16 figures, 7 tables, ApJ, in press

We assemble a large comprehensive sample of 2534 z~2, 3, 4, 5, 6, 7, 8, and 9 galaxies lensed by the six clusters from the Hubble Frontier Fields (HFF) program. Making use of the availability of multiple independent magnification models for each of the HFF clusters and alternatively treating one of the models as the "truth," we show that the median magnification factors from the v4 parametric models are typically reliable to values of 30 to 50, and in one case to 100. Using the median magnification factor from the latest v4 models, we estimate the UV luminosities of the 2534 lensed z~2-9 galaxies, finding sources as faint as -12.4 mag at z~3 and -12.9 mag at z~7. We explicitly demonstrate the power of the surface density-magnification relations Sigma(z) vs. mu in the HFF clusters to constrain both distant galaxy properties and cluster lensing properties. Based on the Sigma(z) vs. mu relations, we show that the median magnification estimates from existing public models must be reliable predictors of the true magnification mu to mu<15 (95% confidence). We also use the observed Sigma(z) vs. mu relations to derive constraints on the evolution of the luminosity function faint-end slope from z~7 to z~2, showing that faint-end slope results can be consistent with blank-field studies if, and only if, the selection efficiency shows no strong dependence on the magnification factor mu. This can only be the case if very low luminosity galaxies are very small, being unresolved in deep lensing probes., Comment: 22 pages, 13 figures, 6 tables, accepted for publication in the Astrophysical Journal, updated to be consistent with the proofs

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures., Comment: Accepted for publication in MNRAS

Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$\alpha$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$\alpha$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$\alpha$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$\alpha$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$\alpha$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources., Comment: 13 pages, 14 figures, Submitted to MNRAS

We constrain the rest-UV size-luminosity relation for star-forming galaxies at z~4 and z~6, 7, and 8 identified behind clusters from the Hubble Frontier Fields (HFF) program. The size-luminosity relation is key to deriving accurate luminosity functions (LF) for faint galaxies. Making use of the latest lensing models and full data set for these clusters, lensing-corrected sizes and luminosities are derived for 68 z~4, 184 z~6, 93 z~7, and 53 z~8 galaxies. We show that size measurements can be reliably measured up to linear magnifications of 30x, where the lensing models are well calibrated. The sizes we measure span a >1-dex range, from <50 pc to >~500 pc. Uncertainties are based on both the formal fit errors and systematic differences between the public lensing models. These uncertainties range from ~20 pc for the smallest sources to 50 pc for the largest. Using a forward-modeling procedure to model the impact of incompleteness and magnification uncertainties, we characterize the size-luminosity relation at both z~4 and z~6-8. We find that the source sizes of star-forming galaxies at z~4 and z~6-8 scale with luminosity L as L^{0.54\pm0.08} and L^{0.40+/-0.04}, respectively, such that lower luminosity (>~-18 mag) galaxies are smaller than expected from extrapolating the size-luminosity relation at high luminosities (<~-18 mag). The new evidence for a steeper size-luminosity relation (3 sigma) adds to earlier evidence for small sizes based on the prevalence of highly magnified galaxies in high shear regions, theoretical arguments against upturns in the LFs, and other independent determinations of the size-luminosity relation from the HFF clusters., Comment: 22 pages, 18 figures, 2 tables, accepted for publication into ApJ. arXiv admin note: substantial text overlap with arXiv:1711.02090