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

A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$\pi$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies, and quasars. To construct training/validation datasets for the classification, we carefully built two different master stellar samples in addition to the spectroscopic galaxy and quasar samples. An ensemble classification model is obtained by averaging two XGBoost classifiers trained with different master stellar samples. Using a probability threshold of $p_{\mathrm{QSO\_mean}}>0.95$ in our ensemble classification model and an additional cut on the logarithmic probability density of zero proper motion, we retrieved 1,545,514 reliable quasar candidates from the parent Gaia DR3 quasar candidate catalog. We provide photometric redshifts for all candidates with an ensemble regression model. For a subset of 89,100 candidates, accurate spectroscopic redshifts are estimated with the Convolutional Neural Network from the Gaia BP/RP spectra. The CatNorth catalog has a high purity of ~ 90% while maintaining high completeness, which is an ideal sample to understand the quasar population and its statistical properties. The CatNorth catalog is used as the main source of input catalog for the LAMOST phase III quasar survey, which is expected to build a highly complete sample of bright quasars with $i < 19.5$., Comment: 27 pages, 15 figures, accepted by ApJS. Table 4 (The CatNorth quasar candidate catalog) is available at https://nadc.china-vo.org/res/r101313/

A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$\mu$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Bright Line Emission Survey, REBELS, with the Atacama Large Millimeter/submillimeter Array. We find substantial amounts of molecular gas ($\sim10^{10.5}\ M_\odot$) comparable to those found in lower redshift galaxies for similar stellar masses ($\sim10^{10}\ M_\odot$). The REBELS galaxies appear to follow the standard scaling relations of molecular gas to stellar mass ratio ($\mu_{\rm mol}$) and gas depletion timescale ($t_{\rm dep}$) with distance to the star-forming main-sequence expected from extrapolations of $z\sim1-4$ observations. We find median values at $z\sim7$ of $\mu_{\rm mol}=2.6_{-1.4}^{4.1}$ and $t_{\rm dep}=0.5_{-0.14}^{+0.26}$ Gyr, indicating that the baryonic content of these galaxies is gas-phase dominated and little evolution from $z\sim7$ to 4. Our measurements of the cosmic density of molecular gas, log$(\rho_{\rm mol}/(M_\odot {\rm Mpc}^{-3}))=6.34^{+0.34}_{-0.31}$, indicate a steady increase by an order of magnitude from $z\sim7$ to 4., Comment: Accepted for publication in Astronomy & Astrophysics. 11 pages, 7 figures. Typo fixed in an author's name

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

Observations of the rest-frame UV emission of high-redshift galaxies suggest that the early stages of galaxy formation involve disturbed structures. Imaging the cold interstellar medium can provide a unique view of the kinematics associated with the assembly of galaxies. In this paper, we analyzed the spatial distribution and kinematics of the cold ionized gas of the normal star-forming galaxy COS-2987030247 at z = 6.8076, based on new high-resolution observations of the [C II] 158um line emission obtained with the Atacama Large Millimeter/submillimeter Array. These observations allowed us to compare the spatial distribution and extension of the [C II] and rest-frame UV emission, model the [C II] line data-cube using 3DBarolo, and measure the [C II] luminosity and star formation rate (SFR) surface densities in the galaxy subregions. The system is found to be composed of a main central source, a fainter north extension, and candidate [C II] companions located 10-kpc away. We find similar rest-frame UV and [C II] spatial distributions, suggesting that the [C II] emission emerges from the star-forming regions. The agreement between the UV and [C II] surface brightness radial profiles rules out diffuse, extended [C II] emission in the main galaxy component. The [C II] velocity map reveals a velocity gradient in the north-south direction suggesting ordered motion, as commonly found in rotating-disk galaxies. But higher-resolution observations would be needed to rule out a compact merger scenario. Our model indicates a low average velocity dispersion, $\sigma$ < 30 km s$^{-1}$. This result implies a dispersion lower than the expected value from observations and semi-analytic models of high redshift galaxies. We argue that COS-2987030247 is a candidate rotating disk experiencing a short period of stability which will be possibly perturbed at later times by accreting sources., Comment: 13 pages (+ 8 Appendix), 15 Figures, Accepted for publication in A&A

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 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 present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman break galaxy in the Epoch of Reionisation, A1703-zD1, through the detection of [C II] at a redshift of z = 6.8269 +/- 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic L$_{UV}$ ~ L$^*$$_{z=7}$ luminosity and a very blue Spitzer/IRAC [3.6]-[4.5] colour, implying high equivalent width line emission of [O III]+H$\beta$. [C II] is reliably detected at 6.1$\sigma$ co-spatial with the rest-frame UV counterpart, showing similar spatial extent. Correcting for the lensing magnification, the [C II] luminosity in A1703-zD1 is broadly consistent with the local L$_{[CII]}$ - SFR relation. We find a clear velocity gradient of 103 +/- 22 km/s across the source which possibly indicates rotation or an ongoing merger. We furthermore present spectral scans with no detected [C II] above 4.6$\sigma$ in two unlensed Lyman break galaxies in the EGS-CANDELS field at z ~ 6.6 - 6.9. This is the first time that NOEMA has been successfully used to observe [C II] in a 'normal' star-forming galaxy at z > 6, and our results demonstrate its capability to complement ALMA in confirming galaxies in the Epoch of Reionisation., Comment: 9 pages, 15 figures. Accepted for publication in 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

Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density., Comment: 16 pages, 7 figures, 4 tables, published in Nature

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization., Comment: 30 pages, 4 tables, 18 figures, accepted for publication in ApJ

We compare the sizes and luminosities of faint $z=6$-8 galaxies magnified by the Hubble Frontier Fields (HFF) clusters with star-forming regions, as well as more evolved objects, in the nearby universe. Our high-redshift comparison sample includes 333 z=6-8 galaxies, for which size measurements were made as part of a companion study where lensing magnifications were estimated from various public models. Accurate size measurements for these sources are complicated by the lens model uncertainties, but other results and arguments suggest that faint galaxies are small, as discussed in a companion study. The measured sizes for sources in our comparison sample range from <50 pc to ~500 pc. For many of the lowest luminosity sources, extremely small sizes are inferred, reaching individual sizes as small as 10-30 pc, with several sources in the 10-15 pc range with our conservative magnification limits. The sizes and luminosities are similar to those of single star cluster complexes like 30 Doradus in the lower-redshift universe and -- in a few cases -- super star clusters. The identification of these compact, faint star-forming sources in the z~6-8 universe also allows us to set upper limits on the proto-globular cluster LF at z~6. By comparisons of the counts and sizes with recent models, we rule out (with some caveats) proto-globular cluster formation scenarios favoring substantial (xi=10) post-formation mass loss and set useful upper limits on others. Our size results suggest we may be very close to discovering a bona-fide population of forming globular clusters at high redshift., Comment: 14 pages, 5 figures, AJ, in press. arXiv admin note: substantial text overlap with arXiv:1711.02090

Here we provide the most comprehensive determinations of the rest-frame $UV$ LF available to date with HST at z~2, 3, 4, 5, 6, 7, 8, and 9. Essentially all of the non-cluster extragalactic legacy fields are utilized, including the Hubble Ultra Deep Field (HUDF), the Hubble Frontier Field parallel fields, and all five CANDELS fields, for a total survey area of 1136 arcmin^2. Our determinations include galaxies at z~2-3 leveraging the deep HDUV, UVUDF, and ERS WFC3/UVIS observations available over a ~150 arcmin^2 area in the GOODS North and GOODS South regions. All together, our collective samples include >24,000 sources, >2.3x larger than previous selections with HST. 5766, 6332, 7240, 3449, 1066, 601, 246, and 33 sources are identified at z~2, 3, 4, 5, 6, 7, 8, and 9, respectively. Combining our results with an earlier z~10 LF determination by Oesch+2018a, we quantify the evolution of the $UV$ LF. Our results indicate that there is (1) a smooth flattening of the faint-end slope alpha from alpha~-2.4 at z~10 to -1.5 at z~2, (2) minimal evolution in the characteristic luminosity M* at z>~2.5, and (3) a monotonic increase in the normalization log_10 phi* from z~10 to z~2, which can be well described by a simple second-order polynomial, consistent with an "accelerated" evolution scenario. We find that each of these trends (from z~10 to z~2.5 at least) can be readily explained on the basis of the evolution of the halo mass function and a simple constant star formation efficiency model., Comment: 18 pages, 10 figures, 5 tables, in submission to ApJ, figures 9 and 10 show the main result

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works., Comment: Accepted for publication in A&A, 28 pages and 20 figures. Full versions of Tables C.1, C.2, C.3, C.4, and C.5 are available for download at http://astro.puc.cl/~rcarvaja/docs/tables_ALMA_FF_V_RCP.tar.gz

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged), Comment: 32 pages, 22 figures

We report the results of an expanded search for z~9-10 candidates over the ~883 arcmin^2 CANDELS+ERS fields. This study adds 147 arcmin^2 to the search area we consider over the CANDELS COSMOS, UDS, and EGS fields, while expanding our selection to include sources with bluer J_{125}-H_{160} colors than our previous J_{125}-H_{160}>0.5 mag selection. In searching for new z~9-10 candidates, we make full use of all available HST, Spitzer/IRAC, and ground-based imaging data. As a result of our expanded search and use of broader color criteria, 3 new candidate z~9-10 galaxies are identified. We also find again the z=8.683 source previously confirmed by Zitrin+2015. This brings our sample of probable z~9-11 galaxy candidates over the CANDELS+ERS fields to 19 sources in total, equivalent to 1 candidate per 47 arcmin^2 (1 per 10 WFC3/IR fields). To be comprehensive, we also discuss 28 mostly lower likelihood z~9-10 candidates, including some sources that seem to be reliably at z>8 using the HST+IRAC data alone, but which the ground-based data show are much more likely at z<4. One case example is a bright z~9.4 candidate COS910-8 which seems instead to be at z~2. Based on this expanded sample, we obtain a more robust LF at z~9 and improved constraints on the volume density of bright z~9 and z~10 galaxies. Our improved z~9-10 results again reinforce previous findings for strong evolution in the UV LF at z>8, with a factor of ~10 evolution seen in the luminosity density from z~10 to z~8., Comment: 22 pages, 12 figures, 6 tables, accepted for publication in the Astrophysical Journal

The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$\beta$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$\beta$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(\xi_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$\beta$ and UV luminosity, which we use to derive for the first time the [OIII]+H$\beta$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$\beta$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$\beta$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$\beta$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST)., Comment: 11 pages, 10 figures, accepted for publication in ApJS

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event., Comment: 25 pages, 13 figures, 2 tables. Accepted for publication in ApJ

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations., Comment: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

We compare the sizes and luminosities of 307 faint z=6-8 sources revealed by the Hubble Frontier Fields (HFF) program with sources in the nearby universe. Making use of the latest lensing models and data from the first four HFF clusters with an extensive suite of public lens models, we measure both the sizes and luminosities for 153 z~6, 101 z~7, and 53 z~8 galaxies. The sizes range over more than a decade from ~500 to <50 pc. Extremely small sizes are inferred for many of our lowest luminosity sources, reaching individual sizes as small as 10-30 pc (the smallest is 11(-6)(+28) pc). The uncertainty in these measures ranges from 80 pc for the largest sources to typically about 20 pc for the smallest. Such sizes are smaller than extrapolations of the size-luminosity relation, and expectations for the completeness of our faint samples, suggesting a likely break in the size-luminosity relation at ~-17 mag with size proportional to L**(0.50(-0.11)(+0.10)). The sizes and luminosities of the lowest-luminosity sources are similar to those of single star cluster complexes like 30 Doradus in the lower-redshift universe and -- in a few cases -- super star clusters. Remarkably, our identification of these compact, faint star-forming sources in the z~6-8 universe also allow us to set upper limits on the proto-globular cluster LF at z~6. Comparisons with recent models allow us to rule out (with some caveats) some scenarios for proto-globular cluster formation and set useful upper limits on other less extreme ones. Our results suggest we may be very close to discovering a bona-fide population of forming globular clusters at high redshift., Comment: The authors have decided to withdraw this manuscript from the Astrophysical Journal given the unprofessional footnotes in the original astro-ph posting, offering our sincere apologies

We present an analysis of all prime HST legacy fields spanning >800 arcmin^2 for the search of z~10 galaxy candidates and the study of their UV luminosity function (LF). In particular, we present new z~10 candidates selected from the full Hubble Frontier Field (HFF) dataset. Despite the addition of these new fields, we find a low abundance of z~10 candidates with only 9 reliable sources identified in all prime HST datasets that include the HUDF09/12, the HUDF/XDF, all the CANDELS fields, and now the HFF survey. Based on this comprehensive search, we find that the UV luminosity function decreases by one order of magnitude from z~8 to z~10 at all luminosities over a four magnitude range. This also implies a decrease of the cosmic star-formation rate density by an order of magnitude within 170 Myr from z~8 to z~10. We show that this accelerated evolution compared to lower redshift can entirely be explained by the fast build-up of the dark matter halo mass function at z>8. Consequently, the predicted UV LFs from several models of galaxy formation are in good agreement with this observed trend, even though the measured UV LF lies at the low end of model predictions. In particular, the number of only 9 observed candidate galaxies is lower, by ~50%, than predicted by galaxy evolution models. The difference is generally still consistent within the Poisson and cosmic variance uncertainties. However, essentially all models predict larger numbers than observed. We discuss the implications of these results in light of the upcoming James Webb Space Telescope mission, which is poised to find much larger samples of z~10 galaxies as well as their progenitors at less than 400 Myr after the Big Bang., Comment: 13 pages, 6 figures, minor updates to match accepted version

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope., Comment: 11 pages, 6 figures, submitted to ApJ

We use the largest sample of z~6 galaxies to date from the first four Hubble Frontier Fields clusters to set constraints on the shape of the z~6 luminosity functions (LFs) to fainter than Muv=-14 mag. We quantify, for the first time, the impact of magnification uncertainties on LF results and thus provide more realistic constraints than other recent work. Our simulations reveal that for the highly-magnified sources the systematic uncertainties can become extremely large fainter than -14 mag, reaching several orders of magnitude at 95% confidence at ~-12 mag. Our new forward-modeling formalism incorporates the impact of magnification uncertainties into the LF results by exploiting the availability of many independent magnification models for the same cluster. One public magnification model is used to construct a mock high-redshift galaxy sample that is then analyzed using the other magnification models to construct a LF. Large systematic errors occur at high magnifications (mu>30) because of differences between the models. The volume densities we derive for faint (>-17 mag) sources are ~3-4x lower than one recent report and give a faint-end slope alpha=-1.92+/-0.04, which is 3.0-3.5sigma shallower (including or not including the size uncertainties, respectively). We introduce a new curvature parameter delta to model the faint end of the LF and demonstrate that the observations permit (at 68% confidence) a turn-over at z~6 in the range -15.3 to -14.2 mag, depending on the assumed lensing model. The present consideration of magnification errors and new size determinations raise doubts about previous reports regarding the form of the LF at >-14 mag. We discuss the implications of our turn-over constraints in the context of recent theoretical predictions., Comment: 33 pages, 19 figures, 5 tables, accepted for publication in the Astrophysical Journal, updated to match the version in press

We provide the first observational constraints on the sizes of the faintest galaxies lensed by the Hubble Frontier Fields (HFF) clusters. Ionizing radiation from faint galaxies likely drives cosmic reionization, and the HFF initiative provides a key opportunity to find such galaxies. Yet, we cannot really assess their ionizing emissivity without a robust measurement of their sizes, since this is key to quantifying both their prevalence and the faint-end slope to the UV luminosity function. Here we provide the first such size constraints with 2 new techniques. The first utilizes the fact that the detectability of highly-magnified galaxies as a function of shear is very dependent on a galaxy's size. Only the most compact galaxies will remain detectable in regions of high shear (vs. a larger detectable size range for low shear), a phenomenon we carefully quantify using simulations. Remarkably, however, no correlation is found between the surface density of faint galaxies and the predicted shear, using 87 faint high-magnification mu>10 z~2-8 galaxies seen behind the first 4 HFF clusters. This can only be the case if such faint (~-15 mag) galaxies have significantly smaller sizes than luminous galaxies. We constrain their half-light radii to be <~30 mas (<160-240 pc). As a 2nd size probe, we rotate and stack 26 faint high-magnification sources along the major shear axis. Less elongation is found than even for objects with an intrinsic half-light radius of 10 mas. Together these results indicate that extremely faint z~2-8 galaxies have near point-source profiles in the HFF dataset (half-light radii conservatively <30 mas and likely 5-10 mas). These results suggest smaller completeness corrections and hence much lower volume densities for faint z~2-8 galaxies and shallower faint-end slopes than have been derived in many recent studies (by factors of ~2-3 and by dalpha>~0.1-0.3)., Comment: 19 pages, 15 figures, 3 tables, accepted for publication in ApJ

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion., Comment: Accepted for publication in A&A. New version addressing last referee comments, 17 pages and 12 figures

The James Webb Space Telescope (JWST) will be an exquisite new near-infrared observatory with imaging and multi-object spectroscopy through ESA's NIRspec instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be positioned on astronomical targets by opening specific 0.002"-wide micro shutter doors. To ensure proper target acquisition, the on-sky position of the MSA needs to be verified before spectroscopic observations start. An onboard centroiding program registers the position of pre-identified guide stars in a Target Acquisition (TA) image, a short pre-spectroscopy exposure without dispersion (image mode) through the MSA with all shutters open. The outstanding issue is the availability of Galactic stars in the right luminosity range for TA relative to typical high redshift targets. We explore this here using the stars and $z\sim8$ candidate galaxies identified in the source extractor catalogs of Brightest of Reionizing Galaxies survey (BoRG[z8]), a pure-parallel program with Hubble Space Telescope Wide-Field Camera 3. We find that (a) a single WFC3 field contains enough Galactic stars to satisfy the NIRspec astrometry requirement (20 milli-arcseconds), provided its and the NIRspec TA's are $m_{lim}>24.5$ AB in WFC3 F125W, (b) a single WFC3 image can therefore serve as the pre-image if need be, (c) a WFC3 mosaic and accompanying TA image satisfy the astrometry requirement at $\sim23$ AB mag in WFC3 F125W, (d) no specific Galactic latitude requires deeper TA imaging due to a lack of Galactic stars, and (e) a depth of $\sim24$ AB mag in WFC3 F125W is needed if a guide star in the same MSA quadrant as a target is required. We take the example of a BoRG identified $z\sim8$ candidate galaxy and require a Galactic star within 20" of it. In this case, a depth of 25.5 AB in F125W is required (with $\sim$97% confidence)., Comment: 17 pages, 15 figures, to appear in the Journal of Astronomical Instrumentation

The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2{\sigma} with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective., Comment: Accepted by ApJ. 10 pages, 5 figures

We have identified 274 M-type Brown Dwarfs in the Hubble Space Telescope's Wide Field Camera 3 (WFC3) pure parallel fields from the Brightest of Reionizing Galaxies (BoRG) survey for high redshift galaxies. These are near-infrared observations with multiple lines-of-sight out of our Milky Way. Using these observed M-type Brown Dwarfs we fitted a Galactic disk and halo model with a Markov chain Monte Carlo (MCMC) analysis. This model worked best with the scale length of the disk fixed at $h$ = 2.6 kpc. For the scale height of the disk, we found $z_0 = 0.29^{+0.02}_{-0.019}$ kpc and for the central number density $\rho_0 = 0.29^{+0.20}_{-0.13}$ \#/pc$^3$. For the halo we derived a flattening parameter $\kappa$ = 0.45$\pm{0.04}$ and a power-law index $p$ = 2.4$\pm{0.07}$. We found the fraction of M-type brown dwarfs in the local density that belong to the halo to be $f_{h}$ = 0.0075$^{+0.0025}_{-0.0019}$. We found no correlation between subtype of M-dwarf and any model parameters. The total number of M-type Brown Dwarfs in the disk and halo was determined to be $58.2^{+9.81}_{-6.70} \times10^{9}$. We found an upper limit for the fraction of M-type Brown Dwarfs in the halo of 7$^{+5}_{-4}$\%. The upper limit for the total Galactic Disk mass in M-dwarfs is $4.34^{+0.73}_{-0.5}\times10^{9}$ $M_{\odot}$, assuming all M-type Brown Dwarfs have a mass of $80 M_J$., Comment: 9 Figures, 8 Tables, MNRAS

We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $\lambda=1.47\pm0.01~\mu$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that this continuum break is the Ly$\alpha$ break redshifted to ${z_\mathrm{grism}=11.09^{+0.08}_{-0.12}}$, just $\sim$400 Myr after the Big Bang. This observation extends the current spectroscopic frontier by 150 Myr to well before the Planck (instantaneous) cosmic reionization peak at z~8.8, demonstrating that galaxy build-up was well underway early in the reionization epoch at z>10. GN-z11 is remarkably and unexpectedly luminous for a galaxy at such an early time: its UV luminosity is 3x larger than L* measured at z~6-8. The Spitzer IRAC detections up to 4.5 $\mu$m of this galaxy are consistent with a stellar mass of ${\sim10^{9}~M_\odot}$. This spectroscopic redshift measurement suggests that the James Webb Space Telescope (JWST) will be able to similarly and easily confirm such sources at z>10 and characterize their physical properties through detailed spectroscopy. Furthermore, WFIRST, with its wide-field near-IR imaging, would find large numbers of similar galaxies and contribute greatly to JWST's spectroscopy, if it is launched early enough to overlap with JWST., Comment: 12 pages, 11 figures, accepted for publication in ApJ

We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from {\lambda}=0.35{\mu}m to {\lambda}=1.7{\mu}m. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization., Comment: Accepted for publication on ApJ. 21 pages, 11 figures, 4 tables

Galaxies represent one of the preferred candidate sources to drive the reionization of the universe. Even as gains are made in mapping the galaxy UV luminosity density to z>6, significant uncertainties remain regarding the conversion to the implied ionizing emissivity. The relevant unknowns are the Lyman-continuum (LyC) photon production efficiency xi_{ion} and the escape fraction f_{esc}. As we show here, the first of these unknowns is directly measureable in z=4-5 galaxies, based on the impact the Halpha line has on the observed IRAC fluxes. By computing a LyC photon production rate from the implied Halpha luminosities for a broad selection of z=4-5 galaxies and comparing this against the dust-corrected UV-continuum luminosities, we provide the first-ever direct estimates of the LyC photon production efficiency xi_{ion} for the z>~4 galaxy population. We find log_{10} xi_{ion}/[Hz/ergs] to have a mean value of 25.27_{-0.03}^{+0.03} and 25.34_{-0.02}^{+0.02} for sub-L* z=4-5 galaxies adopting Calzetti and SMC dust laws, respectively. Reassuringly, both values are consistent with standardly assumed xi_{ion}'s in reionization models, with a slight preference for higher xi_{ion}'s (by ~0.1 dex) adopting the SMC dust law. A modest ~0.03-dex increase in these estimates would result if the escape fraction for ionizing photons is non-zero and galaxies dominate the ionizing emissivity at z~4.4. High values of xi_{ion} (~25.5-25.8 dex) are derived for the bluest galaxies (beta<-2.3) in our samples, independent of dust law and consistent with results for a z=7.045 galaxy. Such elevated values of xi_{ion} would have important consequences, indicating that f_{esc} cannot be in excess of 13% unless the galaxy UV luminosity function does not extend down to -13 mag or the clumping factor is greater than 3. A low escape fraction would fit well with the low rate of LyC leakage observed at z~3., Comment: 12 pages, 6 figures, 3 tables, updated to match the version in press, both statistical and systematic errors given for xi_{ion}

We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,\sigma$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered., Comment: 16 pages, 13 figures, accepted for publication in MNRAS

The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6{\mu}m and 4.5{\mu}m with the Spitzer Infrared Array Camera (IRAC). The primary aim is to directly detect the infrared light of reionization epoch galaxies at z > 7 and to constrain their stellar populations. The observations cover the Hubble Ultra Deep Field (HUDF), including the two HUDF parallel fields, and the CANDELS/GOODS-South, and are combined with archival data from all previous deep programs into one ultradeep dataset. The resulting imaging reaches unprecedented coverage in IRAC 3.6{\mu}m and 4.5{\mu}m ranging from > 50 hour over 150 arcmin^2, > 100 hour over 60 sq arcmin2, to 200 hour over 5 - 10 arcmin$^2$. This paper presents the survey description, data reduction, and public release of reduced mosaics on the same astrometric system as the CANDELS/GOODS-South WFC3 data. To facilitate prior-based WFC3+IRAC photometry, we introduce a new method to create high signal-to-noise PSFs from the IRAC data and reconstruct the complex spatial variation due to survey geometry. The PSF maps are included in the release, as are registered maps of subsets of the data to enable reliability and variability studies. Simulations show that the noise in the ultradeep IRAC images decreases approximately as the square root of integration time over the range 20 - 200 hours, well below the classical confusion limit, reaching 1{\sigma} point source sensitivities as faint as of 15 nJy (28.5 AB) at 3.6{\mu}m and 18 nJy (28.3 AB) at 4.5{\mu}m. The value of such ultradeep IRAC data is illustrated by direct detections of z = 7 - 8 galaxies as faint as HAB = 28., Comment: Accepted for publication in ApJS

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models., Comment: 23 pages, 15 figures, 7 tables, updated to match the version in press, including some minor textual corrections identified at the proof stage

We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We identify z ~ 7.1 - 9.1 galaxies by selecting z >~ 6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating strong [OIII]+Hbeta lines in the 4.5 micron band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage, and subsequently used to select the brightest z > 7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the CANDELS survey revealed four unusually bright z ~ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z ~ 7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1500A, in the [4.5] band. Keck/MOSFIRE spectroscopy has been independently reported for two of our selected sources, showing Ly-alpha at redshifts of 7.7302+/-0.0006 and 8.683^+0.001_-0.004, respectively. We present similar Keck/MOSFIRE spectroscopy for a third selected galaxy with a probable 4.7sigma Ly-alpha line at z_spec=7.4770+/-0.0008. All three have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV) ~ -22.0) than any previously discovered z ~ 8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z > 7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance., Comment: 19 pages, 13 figures, 5 tables, accepted for publication in the Astrophysical Journal. Updated to match the version in press and reflect spectroscopic confirmation of the highest redshift candidate in our sample by Zitrin+2015

Thomson optical depth tau measurements from Planck provide new insights into the reionization of the universe. In pursuit of model-independent constraints on the properties of the ionising sources, we determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter model to map out the evolution in the emissivity at z>~6 from the new Planck optical depth tau measurements, from the constraints provided by quasar absorption spectra and from the prevalence of Ly-alpha emission in z~7-8 galaxies. We find the redshift evolution in the emissivity dot{N}_{ion}(z) required by the observations to be d(log Nion)/dz=-0.15(-0.11)(+0.08), largely independent of the assumed clumping factor C_{HII} and entirely independent of the nature of the ionising sources. The trend in dot{N}_{ion}(z) is well-matched by the evolution of the galaxy UV-luminosity density (dlog_{10} rho_UV/dz=-0.11+/-0.04) to a magnitude limit >~-13 mag, suggesting that galaxies are the sources that drive the reionization of the universe. The role of galaxies is further strengthened by the conversion from the UV luminosity density rho_UV to dot(N)_{ion}(z) being possible for physically-plausible values of the escape fraction f_{esc}, the Lyman-continuum photon production efficiency xi_{ion}, and faint-end cut-off $M_{lim}$ to the luminosity function. Quasars/AGN appear to match neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred evolution in the ionizing emissivity, we estimate that the z~10 UV-luminosity density is 8(-4)(+15)x lower than at $z~6, consistent with the observations. The present approach of contrasting the inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the growing observational evidence that faint, star-forming galaxies drive the reionization of the universe., Comment: 20 pages, 12 figures, 5 tables, Astrophysical Journal, updated to match version in press, Figure 6 shows the main result of the paper

We present a spectroscopic redshift measurement of a very bright Lyman break galaxy at z=7.7302+-0.0006 using Keck/MOSFIRE. The source was pre-selected photometrically in the EGS field as a robust z~8 candidate with H=25.0 mag based on optical non-detections and a very red Spitzer/IRAC [3.6]-[4.5] broad-band color driven by high equivalent width [OIII]+Hbeta line emission. The Lyalpha line is reliably detected at 6.1 sigma and shows an asymmetric profile as expected for a galaxy embedded in a relatively neutral inter-galactic medium near the Planck peak of cosmic reionization. The line has a rest-frame equivalent width of EW0=21+-4 A and is extended with V_FWHM=360+90-70 km/s. The source is perhaps the brightest and most massive z~8 Lyman break galaxy in the full CANDELS and BoRG/HIPPIES surveys, having assembled already 10^(9.9+-0.2) M_sol of stars at only 650 Myr after the Big Bang. The spectroscopic redshift measurement sets a new redshift record for galaxies. This enables reliable constraints on the stellar mass, star-formation rate, formation epoch, as well as combined [OIII]+Hbeta line equivalent widths. The redshift confirms that the IRAC [4.5] photometry is very likely dominated by line emission with EW0(OIII+Hbeta)= 720-150+180 A. This detection thus adds to the evidence that extreme rest-frame optical emission lines are a ubiquitous feature of early galaxies promising very efficient spectroscopic follow-up in the future with infrared spectroscopy using JWST and, later, ELTs., Comment: 6 pages, 4 figures, small updates to match ApJL accepted version

We search the complete Hubble Frontier Field dataset of Abell 2744 and its parallel field for z~10 sources to further refine the evolution of the cosmic star-formation rate density (SFRD) at z>8. We independently confirm two images of the recently discovered triply-imaged z~9.8 source by Zitrin et al. (2014) and set an upper limit for similar z~10 galaxies with red colors of J_125-H_160>1.2 in the parallel field of Abell 2744. We utilize extensive simulations to derive the effective selection volume of Lyman-break galaxies at z~10, both in the lensed cluster field and in the adjacent parallel field. Particular care is taken to include position-dependent lensing shear to accurately account for the expected sizes and morphologies of highly-magnified sources. We show that both source blending and shear reduce the completeness at a given observed magnitude in the cluster, particularly near the critical curves. These effects have a significant, but largely overlooked, impact on the detectability of high-redshift sources behind clusters, and substantially reduce the expected number of highly-magnified sources. The detections and limits from both pointings result in a SFRD which is higher by 0.4+-0.4 dex than previous estimates at z~10 from blank fields. Nevertheless, the combination of these new results with all other estimates remain consistent with a rapidly declining SFRD in the 170 Myr from z~8 to z~10 as predicted by cosmological simulations and dark-matter halo evolution in LambdaCDM. Once biases introduced by magnification-dependent completeness are accounted for, the full six cluster and parallel Frontier Field program will be an extremely powerful new dataset to probe the evolution of the galaxy population at z>8 before the advent of the JWST., Comment: 10 pages, 7 figures, changed to match accepted version to appear in ApJ

Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $\sim4\sigma$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{\odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}\sim-17.7$) are expected in lower-mass halos ($M\sim10^{8}-10^{10.5}$M$_{\odot}$). We compute the ACF of LBGs at $z\sim3.8-z\sim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($b\sim3.0$) to $z=7.2$ ($b\sim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4

The remarkable HST datasets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map out the evolution of the UV LF from z~10 to z~4. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z~4, z~5, z~6, z~7, z~8, and z~10, respectively from the ~1000 arcmin**2 area probed. The selection of z~4-8 galaxies over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are apparent at z>=7. Our new LF determinations at z~4 and z~5 span a 6-mag baseline (-22.5 to -16 AB mag). These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L* galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to 3.4 sigma significance for a steeper faint-end slope to the UV LF at z>4, with alpha evolving from alpha=-1.64+/-0.04 at z~4 to alpha=-2.06+/-0.13 at z~7 (and alpha = -2.02+/-0.23 at z~8), consistent with that expected from the evolution of the halo mass function. With our improved constraints at the bright end, we find less evolution in the characteristic luminosity M* over the redshift range z~4 to z~7; the observed evolution in the LF is now largely represented by changes in phi*. No evidence for a non-Schechter-like form to the z~4-8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio of halos ((1+z)**-1.5) provides a good representation of the observed evolution., Comment: 53 pages, 28 figures, 11 tables, accepted for publication in the Astrophysical Journal

[abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10., Comment: 20 pages, 12 figures, changed to match resubmitted version to ApJ