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

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA Band-7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 16. per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (7 show a counterpart only in the mid-IR and no HST or near-IR identification, while 2 are detected as [CII] emitters at z=5). The 7 HST and near-IR dark galaxies with mid-IR counterpart contribute for about 15 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.

We present the first [CII] 158 $\mu$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018).

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

The infrared spectral energy distributions (SEDs) of main-sequence galaxies in the early universe (z > 4) is currently unconstrained as infrared continuum observations are time consuming and not feasible for large samples. We present Atacama Large Millimetre Array (ALMA) Band 8 observations of four main-sequence galaxies at z ~ 5.5 to study their infrared SED shape in detail. Our continuum data (rest-frame 110$\rm \mu m$, close to the peak of infrared emission) allows us to constrain luminosity weighted dust temperatures and total infrared luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame 100$\rm \mu m$. The derived dust peak temperatures at z ~ 5.5 (38$\pm$8K) are elevated compared to average local galaxies, however, 5-10K below what would be predicted from an extrapolation of the trend at $z<4$. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the infrared SED at the peak to be more optically thin, making hot dust more visible to the external observer. From the 850$\rm \mu m$ dust continuum, we derive molecular gas masses between $10^{10}$ and $10^{11}\,{\rm M_{\odot}}$ and gas fractions (gas over total mass) of 30-80% (gas depletion times of 100-220Myrs). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending out to ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$\alpha$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$\pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($\Sigma$), we propose a model relating $f_{esc}$ to $\Sigma$ and find $f_{esc}\propto\Sigma^{0.4\pm0.1}$. Since $\Sigma$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs') account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($\alpha$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$\sigma$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($\beta$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($\Delta\beta$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$\lambda$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

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.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

We present an analysis of the dust attenuation of star forming galaxies at $z=2.5-4.0$ through the relationship between the UV spectral slope ($\beta$), stellar mass ($M_{\ast}$) and the infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$) based on far-infrared continuum observations from the Atacama Large Millimeter/sub-millimeter Array (ALMA). Our study exploits the full ALMA archive over the COSMOS field processed by the A$^3$COSMOS team, which includes an unprecedented sample of $\sim1500$ galaxies at $z\sim3$ as primary or secondary targets in ALMA band 6 or 7 observations with a median continuum sensitivity of 126 $\rm{\mu Jy/beam}$ (1$\sigma$). The detection rate is highly mass dependent, decreasing drastically below $\log (M_{\ast}/M_{\odot})=10.5$. The detected galaxies show that the IRX-$\beta$ relationship of massive ($\log M_{\ast}/M_{\odot} > 10$) main sequence galaxies at $z=2.5-4.0$ is consistent with that of local galaxies, while starbursts are generally offset by $\sim0.5\,{\rm dex}$ to larger IRX values. At the low mass end, we derive upper limits on the infrared luminosities through stacking of the ALMA data. The combined IRX-$M_{\ast}$ relation at $\rm{log\,(M_{\ast}/M_{\odot})>9}$ exhibits a significantly steeper slope than reported in previous studies at similar redshifts, implying little dust obscuration at $\log M_{\ast}/M_{\odot}<10$. However, our results are consistent with early measurements at $z\sim5.5$, indicating a potential redshift evolution between $z\sim2$ and $z\sim6$. Deeper observations targeting low mass galaxies will be required to confirm this finding.

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.

We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate redshift galaxies in the COSMOS field. One is near a foreground galaxy at 1.3", but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable $z_{CO}=3.329$, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-mm spectral energy distribution (SED) favors a grey $\lambda^{-0.4}$ attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/sub-mm/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical ($>29.6$ AB) and near-infrared ($>27.9$ AB) images. Using the ALMA position as a prior reveals faint $SNR\sim3$ measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850$\mu$m, and VLA 3GHz, indicating the source is real. The SED is robustly reproduced by a massive $M^*=10^{10.8}$M$_\odot$ and $M_{gas}=10^{11}$M$_\odot$, highly obscured $A_V\sim4$, star forming $SFR\sim300$ M$_{\odot}$yr$^{-1}$ galaxy at redshift $z=5.5\pm$1.1. The ultrasmall 8 arcmin$^{2}$ survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z=5) $\sim0.9\times10^{-2}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at $z>4$, below the typical detection limit of bright galaxies found in single-dish sub-mm surveys, but with larger space densities $\sim3 \times 10^{-5}$ Mpc$^{-3}$, higher duty cycles $50-100\%$, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies.

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

This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 exposures taken in 2635 orbits totaling 6.3 Msec with the HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) and the Wide Field Camera 3 (WFC3) UVIS/IR Channels in the greater GOODS-S extragalactic field, covering all major observational efforts (e.g., GOODS, GEMS, CANDELS, ERS, UVUDF and many other programs; see Illingworth et al 2019, in prep). The HLF GOODS-S catalogs include photometry in 13 bandpasses from the UV (WFC3/UVIS F225W, F275W and F336W filters), optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters), to near-infrared (WFC3/IR F098M, F105W, F125W, F140W and F160W filters). Such a data set makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range from high resolution mosaics that are largely contiguous. Here, we describe a photometric analysis of 186,474 objects in the HST imaging at wavelengths 0.2--1.6$\mu$m. We detect objects from an ultra-deep image combining the PSF-homogenized and noise-equalized F850LP, F125W, F140W and F160W images, including Gaia astrometric corrections. SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. All of the data presented herein are available through the HLF website (https://archive.stsci.edu/prepds/hlf/).

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1<z<8.0$) galaxy candidates found in Brightest of Reionizing Galaxies (BoRG) Survey fields. The original sample includes eight galaxies that are $Y$-band dropout objects with $H$-band magnitudes of $m_H<25.5$. Four of these objects were originally reported in Calvi et al. 2016. Combining new Hubble Space Telescope (HST) WFC3/F814W imaging and $Spitzer$ IRAC data with archival imaging from BoRG and other surveys, we explore the properties of these galaxies. Photometric redshift fitting places six of these galaxies in the redshift range of $7.1<z<8.0$, resulting in three new high-redshift galaxies and confirming three of the four high-redshift galaxy candidates from Calvi et al. 2016. We calculate the half-light radii of the Super Eight galaxies using the HST F160W filter and find that the Super Eight sizes are in line with typical evolution of size with redshift. The Super Eights have a mean mass of log(M$_*$/M$_\odot$) $\sim10$, which is typical for sources in this luminosity range. Finally, we place our sample on the UV $z\sim8$ luminosity function and find that the Super Eight number density is consistent with other surveys in this magnitude and redshift range.

We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2$\leq z\leq$2.0. The novelty of our work lies in the use of a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to $\geq$70% above the MS, implying that -- already at $z{\gtrsim}1$ -- starburst (SB) events ($\Delta_{\rm MS}\geq$0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS-offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of $\sim$5-10% at all masses ${\gtrsim} 10^{10}M_{\odot}$. The MS major merger fraction roughly doubles between $z=0.2$ and $z=2$, with morphological mergers driving the overall increase at $z{\gtrsim}1$. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability timescales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with $\Delta_{\rm MS}$ follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different SFR enhancements. This points to a starburst-merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.

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.

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs.

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

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.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

High-resolution multi-wavelength photometry is crucial to explore the spatial distribution of star formation in galaxies and understand how these evolve. To this aim, in this paper we exploit the deep, multi-wavelength Hubble Space Telescope (HST) data available in the central parts of the GOODS fields and study the distribution of star formation activity and mass in galaxies located at different positions with respect to the Main Sequence (MS) of star-forming galaxies. Our sample consists of galaxies with stellar mass $\geq 10^{9.5} M_{\odot}$ in the redshift range 0.2 $ \leq z \leq 1.2$. Exploiting 10-band photometry from the UV to the near-infrared at HST resolution, we derive spatially resolved maps of galaxies properties, such as stellar mass and star formation rate and specific star formation rate, with a resolution of $\sim 0.16$ arcsec. We find that the star formation activity is centrally enhanced in galaxies above the MS and centrally suppressed below the MS, with quiescent galaxies (1 dex below the MS) characterised by the highest suppression. The sSFR in the outer region does not show systematic trends of enhancement or suppression above or below the MS. The distribution of mass in MS galaxies indicates that bulges are growing when galaxies are still on the MS relation. Galaxies below the MS are more bulge-dominated with respect to MS counterparts at fixed stellar mass, while galaxies in the upper envelope are more extended and have S\'ersic indexes that are always smaller than or comparable to their MS counterparts. The suppression of star formation activity in the central region of galaxies below the MS hints at \textit{inside-out} quenching, as star formation is still ongoing in the outer regions.

We present a full data analysis of the pure-parallel Hubble Space Telescope (HST) imaging observations in the Brightest of Reionizing Galaxies Survey (BoRG[z9]) in Cycle 22. The medium-deep exposures with five HST/WFC3IR+UVIS filter bands from 79 independent sightlines (~370 arcmin^2) provide the least biased determination of number density for z>9 bright galaxies against cosmic variance. After a strict two-step selection for candidate galaxies, including dropout color and photometric redshift analyses, and revision of previous BoRG candidates, we identify one source at z~10 and two sources at z~9. The z~10 candidate shows evidence of line-of-sight lens magnification (mu~1.5), yet it appears surprisingly luminous (MUV ~ -22.6\pm0.3 mag), making it one of the brightest candidates at z > 8 known (~ 0.3 mag brighter than the z = 8.68 galaxy EGSY8p7, spectroscopically confirmed by Zitrin and collaborators). For z ~ 9 candidates, we include previous data points at fainter magnitudes and find that the data are well fitted by a Schechter luminosity function with alpha ~ -2.1, MUV ~ -21.5 mag, and log phi ~ -4.5 Mpc^-3mag^-1, for the first time without fixing any parameters. The inferred cosmic star formation rate density is consistent with unaccelerated evolution from lower redshift.

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering.

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.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

We present Keck-MOSFIRE H and K spectra for a sample of 24 candidate quiescent galaxies (QGs) at 3<z<4, identified from UVJ colors and photometric redshifts in the ZFOURGE and 3DHST surveys. We obtain spectroscopic redshifts for half of the sample, using absorption or emission lines, and confirm the high accuracy of the photometric redshifts with a median error of 1.2%. Two galaxies turn out to be dusty objects at lower redshifts (z<2.5), and are the only two detected in the sub-mm with ALMA. High equivalent-width [OIII] was observed in two galaxies, contributing up to 30% of the K-band flux and mimicking the colors of an old stellar population. This implies a failure rate of only 20% for the UVJ selection at these redshifts. Balmer absorption was identified in 4 of the brighest galaxies, confirming the absence of OB stars. Modeling all QGs with a wide range of star-formation histories, we find sSFR a factor of 10 below the main sequence (MS) for all but one galaxy, and less than 0.01 Gyr$^{-1}$ for half of the sample. This is consistent with the H$\beta$ and [OII] luminosities, and the ALMA non-detections. We then find that these QGs have quenched on average 300 Myr before observation, between z=3.5 and 5, and that they formed at z~5.5 with a mean SFR~300 Msun/yr. Considering an alternative selection of QGs based solely on the sSFR from SED modeling, we find that galaxies a factor 10 below the MS are 40% more numerous than UVJ-quiescent galaxies, implying that the UVJ selection is pure but incomplete. Current models fail at reproducing our observations and underestimate either the number density of QGs by more than an order of magnitude or the duration of their quiescence by a factor two. Overall, these results confirm the existence of an unexpected population of QGs at z>3, and offer the first insights on their formation history. [abridged]

Recent studies have revealed intense UV metal emission lines in a modest sample of z>7 Lyman-alpha emitters, indicating a hard ionizing spectrum is present. If such high ionization features are shown to be common, it may indicate that extreme radiation fields play a role in regulating the visibility of Lyman-alpha in the reionization era. Here we present deep near-infrared spectra of seven galaxies with Lyman-alpha emission at 5.4<z<8.7 (including a newly-confirmed lensed galaxy at z=6.031) and three bright z>7 photometric targets. In nine sources we do not detect UV metal lines. However in the z=8.683 galaxy EGSY8p7, we detect a 4.6 sigma emission line in the narrow spectral window expected for NV 1243. The feature is unresolved (FWHM<90 km/s) and is likely nebular in origin. A deep H-band spectrum of EGSY8p7 reveals non-detections of CIV, He II, and OIII]. The presence of NV requires a substantial flux of photons above 77 eV, pointing to a hard ionizing spectrum powered by an AGN or fast radiative shocks. Regardless of its origin, the intense radiation field of EGSY8p7 may aid the transmission of Lyman-alpha through what is likely a partially neutral IGM. With this new detection, five of thirteen known Lyman-alpha emitters at z>7 have now been shown to have intense UV line emission, suggesting that extreme radiation fields are commonplace among the Lyman-alpha population. Future observations with JWST will eventually clarify the origin of these features and explain their role in the visibility of Lyman-alpha in the reionization era.

We use deep F275W imaging from the $Hubble$ Deep UV Legacy Survey (HDUV) and G280 grism spectroscopy from $HST$/WFC3, along with new and archival optical spectra from Keck/DEIMOS, to search for candidate ionizing sources in the GOODS-N field at z ~ 2.5 - 3. Spectroscopic identification of our UV-selected sources are 99% complete to F275W = 25.5 in the region of the UV imaging, and we identify 6 potential ionizing galaxies or AGNs at z ~ 3. By far the brightest of these is a z = 2.583 AGN that totally dominates the ionizing flux in the region, with a specific ionizing volume emissivity at $912~ \AA$ of $\epsilon_{912}$ = $8.3^{27}_{1.4} \times 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$. Based on our spectroscopic data, we find four candidates are contaminated by foreground galaxies at z ~ 0.5 - 0.7. At $\epsilon_{912}$ = $2.2^{7.2}_{0.4} \times 10^{23}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$, the remaining candidate galaxy's contribution to the ionizing background lies well below the flux required to ionize the intergalactic medium at z ~ 2.5 - 3, consistent with previous observations that show AGNs provide the bulk of the ionizing background at these redshifts.

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).

The luminosity function is a fundamental observable for characterizing how galaxies form and evolve throughout the cosmic history. One key ingredient to derive this measurement from the number counts in a survey is the characterization of the completeness and redshift selection functions for the observations. In this paper we present GLACiAR, an open python tool available on GitHub to estimate the completeness and selection functions in galaxy surveys. The code is tailored for multiband imaging surveys aimed at searching for high-redshift galaxies through the Lyman Break technique, but it can be applied broadly. The code generates artificial galaxies that follow Sersic profiles with different indexes and with customizable size, redshift and spectral energy distribution properties, adds them to input images, and measures the recovery rate. To illustrate this new software tool, we apply it to quantify the completeness and redshift selection functions for J-dropouts sources (redshift z~10 galaxies) in the Hubble Space Telescope Brightest of Reionizing Galaxies Survey (BoRG). Our comparison with a previous completeness analysis on the same dataset shows overall agreement, but also highlights how different modelling assumptions for artificial sources can impact completeness estimates.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Recent work has suggested extreme [OIII] emitting star-forming galaxies are important to reionization. Relatedly, [OIII]/[OII] has been put forward as an indirect estimator of the Lyman Continuum (LyC) escape fraction ($f_{esc}$) at $z\gtrsim4.5$ when the opaque IGM renders LyC photons unobservable. Using deep archival U-band (VLT/VIMOS) imaging of a recently confirmed overdensity at $z\sim3.5$ we calculate tight constraints on $f_{esc}$ for a sample (N=73) dominated by extreme [OIII] emitters. We find no Lyman Continuum signal ($f_{esc}^{rel} < 6.3^{+0.7}_{-0.7} \%$ at $1\sigma$) in a deep U-band stack of our sample (31.98 mag at 1$\sigma$). This constraint is in agreement with recent studies of star-forming galaxies spanning $z\sim1-4$ that have found very low average $f_{esc}$. Despite the galaxies in our study having an estimated average rest-frame EW([OIII]$\lambda5007$)$\sim400\AA$ and [OIII]/[OII]$\sim 4$ from composite SED-fitting, we find no LyC detection, which brings into question the potential of [OIII]/[OII] as an effective probe of the LyC--a majority of LyC emitters have [OIII]/[OII]$>3$, but we establish here that [OIII]/[OII]$>3$ does not guarantee significant LyC leakage for a population. Since even extreme star-forming galaxies are unable to produce the $f_{esc}\sim10-15\%$ required by most theoretical calculations for star-forming galaxies to drive reionization, there must either be a rapid evolution of $f_{esc}$ between $z\sim3.5$ and the Epoch of Reionization, or hitherto observationally unstudied sources (e.g. ultra-faint low-mass galaxies with $\log(M/M_\odot)\sim7-8.5$) must make an outsized contribution to reionization.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

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.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

We present Hubble multi-wavelength photometric catalogs, including (up to) 17 filters with the Advanced Camera for Surveys and Wide Field Camera 3 from the ultra-violet to near-infrared for the Hubble Frontier Fields and associated parallels. We have constructed homogeneous photometric catalogs for all six clusters and their parallels. To further expand these data catalogs, we have added ultra-deep $K_{S}$-band imaging at 2.2~\micron\ from the Very Large Telescope HAWK-I and Keck-I MOSFIRE instruments. We also add post-cryogenic \spitzer\ imaging at 3.6~\micron\ and 4.5~\micron\ with the Infrared Array Camera (IRAC), as well as archival IRAC 5.8~\micron\ and 8.0~\micron\ imaging when available. We introduce the public release of the multi-wavelength (0.2--8~\micron) photometric catalogs, and we describe the unique steps applied for the construction of these catalogs. Particular emphasis is given to the source detection band, the contamination of light from the bright cluster galaxies and intra-cluster light. In addition to the photometric catalogs, we provide catalogs of photometric redshifts and stellar population properties. Furthermore, this includes all the images used in the construction of the catalogs, including the combined models of bright cluster galaxies and intra-cluster light, the residual images, segmentation maps and more. These catalogs are a robust data set of the Hubble Frontier Fields and will be an important aide in designing future surveys, as well as planning follow-up programs with current and future observatories to answer key questions remaining about first light, reionization, the assembly of galaxies and many more topics, most notably, by identifying high-redshift sources to target.

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

A randomly chosen star in today's Universe is most likely to live in a galaxy with a stellar mass between that of the Milky Way and Andromeda. Yet it remains uncertain how the structural evolution of these bulge-disk systems proceeded. Most of the unobscured star formation we observe building Andromdeda progenitors at 0.7<z<1.5 occurs in disks, but >90% of their star formation is reprocessed by dust and remains unaccounted for. Here we map 500micron dust continuum emission in an Andromeda progenitor at z=1.25 to probe where it is growing through dust-obscured star formation. Combining resolved dust measurements from the NOEMA interferometer with Hubble Space Telescope Halpha maps and multicolor imaging (including new UV data from the HDUV survey), we find a bulge growing by dust-obscured star formation: while the unobscured star formation is centrally suppressed, the dust continuum is centrally concentrated, filling in the ring-like structures evident in the Halpha and UV emission. Reflecting this, the dust emission is more compact than the optical/UV tracers of star formation with r_e(dust)=3.4kpc, r_e(Halpha)/r_e(dust)=1.4, and r_e(UV)/r_e(dust)=1.8. Crucially, however, the bulge and disk of this galaxy are building simultaneously; although the dust emission is more compact than the rest-optical emission (r_e(optical)/r_e(dust)=1.4), it is somewhat less compact than the stellar mass (r_e(M_*)/r_e(dust)=0.9). Taking the 500micron emission as a tracer of star formation, the expected structural evolution of this galaxy can be accounted for by star formation: it will grow in size by Delta(r_e)/Delta(M_*)~0.3 and central surface density by Delta(Sigma_cen)/Delta(M_*)~0.9. Finally, our observations are consistent with a picture in which merging and disk instabilities drive gas to the center of galaxies, boosting global star formation rates above the main sequence and building bulges.

We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged]

In the early Universe finding massive galaxies that have stopped forming stars present an observational challenge as their rest-frame ultraviolet emission is negligible and they can only be reliably identified by extremely deep near-infrared surveys. These have revealed the presence of massive, quiescent early-type galaxies appearing in the universe as early as z$\sim$2, an epoch 3 Gyr after the Big Bang. Their age and formation processes have now been explained by an improved generation of galaxy formation models where they form rapidly at z$\sim$3-4, consistent with the typical masses and ages derived from their observations. Deeper surveys have now reported evidence for populations of massive, quiescent galaxies at even higher redshifts and earlier times, however the evidence for their existence, and redshift, has relied entirely on coarsely sampled photometry. These early massive, quiescent galaxies are not predicted by the latest generation of theoretical models. Here, we report the spectroscopic confirmation of one of these galaxies at redshift z=3.717 with a stellar mass of 1.7$\times$10$^{11}$ M$_\odot$ whose absorption line spectrum shows no current star-formation and which has a derived age of nearly half the age of the Universe at this redshift. The observations demonstrates that the galaxy must have quickly formed the majority of its stars within the first billion years of cosmic history in an extreme and short starburst. This ancestral event is similar to those starting to be found by sub-mm wavelength surveys pointing to a possible connection between these two populations. Early formation of such massive systems is likely to require significant revisions to our picture of early galaxy assembly.

We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form.

We present the first determination of the galaxy luminosity function (LF) at z~4, 5, 6 and 7 in the rest-frame optical at lambda_rest~900 nm (z' band). The rest-frame optical light traces the content in low-mass evolved stars (~stellar mass - Mstar), minimizing potential measurement biases for Mstar. Moreover it is less affected by nebular line emission contamination and dust attenuation, is independent of stellar population models, and can be probed up to z~8 through Spitzer/IRAC. Our analysis leverages the unique full depth Spitzer/IRAC 3.6um-to-8.0um data over the CANDELS/GOODS-N, CANDELS/GOODS-S and COSMOS/UltraVISTA fields. We find that at absolute magnitudes M_z' fainter than >-23 mag, M_z' linearly correlates with M_UV,1600. At brighter M_z', M_UV,1600 presents a turnover, suggesting that the stellar mass-to-light ratio Mstar/L_UV,1600 could be characterised by a very broad range of values at high stellar masses. Median-stacking analysis recovers a Mstar/L_z' roughly independent on M_z' for M_z'>-23 mag, but exponentially increasing at brighter magnitudes. We find that the evolution of the LF marginally prefers a pure evolution in luminosity over a pure evolution in density, with the characteristic luminosity decreasing by a factor ~5x between z~4 and z~7. Direct application of the recovered Mstar/L_z' generates stellar mass functions consistent with average measurements from the literature. Measurements of the stellar-to-halo mass ratio at fixed cumulative number density show that it is roughly constant with redshift for Mh>10^12Msun. This is also supported by the fact that the evolution of the LF at 4<z<7 can be accounted for by a rigid displacement in luminosity corresponding to the evolution of the halo mass from abundance matching.

We present deep ALMA CO(5-4) observations of a main sequence, clumpy galaxy at z=1.5 in the HUDF. Thanks to the ~0.5" resolution of the ALMA data, we can link stellar population properties to the CO(5-4) emission on scales of a few kpc. We detect strong CO(5-4) emission from the nuclear region of the galaxy, consistent with the observed $L_{\rm IR}$-$L^{\prime}_{\rm CO(5-4)}$ correlation and indicating on-going nuclear star formation. The CO(5-4) gas component appears more concentrated than other star formation tracers or the dust distribution in this galaxy. We discuss possible implications of this difference in terms of star formation efficiency and mass build-up at the galaxy centre. Conversely, we do not detect any CO(5-4) emission from the UV-bright clumps. This might imply that clumps have a high star formation efficiency (although they do not display unusually high specific star formation rates) and are not entirely gas dominated, with gas fractions no larger than that of their host galaxy (~50%). Stellar feedback and disk instability torques funnelling gas towards the galaxy centre could contribute to the relatively low gas content. Alternatively, clumps could fall in a more standard star formation efficiency regime if their actual star-formation rates are lower than generally assumed. We find that clump star-formation rates derived with several different, plausible methods can vary by up to an order of magnitude. The lowest estimates would be compatible with a CO(5-4) non-detection even for main-sequence like values of star formation efficiency and gas content.

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).

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.

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]

We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.

We have submitted to MAST the 1.5 version data release of the Hubble Legacy Fields (HLF) project covering a 25 x 25 arcmin area over the GOODS-S (ECDF-S) region from the HST archival program AR-13252. The release combines exposures from Hubble's two main cameras, the Advanced Camera for Surveys (ACS/WFC) and the Wide Field Camera 3 (WFC3/IR), taken over more than a decade between mid-2002 to the end of 2016. The HLF includes essentially all optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters) and infrared (WFC3/ IR F098M, F105W, F125W, F140W and F160W filters) data taken by Hubble over the original CDF-S region including the GOODS-S, ERS, CANDELS and many other programs (31 in total). The data has been released at https://archive.stsci.edu/prepds/hlf/ as images with a common astrometric reference frame, with corresponding inverse variance weight maps. We provide one image per filter of WFC3/IR images at 60 mas per pixel resolution and two ACS/WFC images per filter, at both 30 and 60 mas per pixel. Since this comprehensive dataset combines data from 31 programs on the GOODS-S/CDF-S, the AR proposal identified the MAST products by the global name "Hubble Legacy Field", with this region being identified by "HLF-GOODS-S". This dataset complements that of the Frontier Fields program. The total incorporated in the HLF-GOODS-S is 5.8 Msec in 7211 exposures from 2442 orbits. This is ~70% of a HST full cycle!

In this study we focus on the $z\sim6.6$ Lyman-$\alpha$ CR7 consisting of clump A that is host to a potential direct collapse black hole (DCBH), and two metal enriched star forming clumps B and C. In contrast to claims that signatures of metals rule out the existence of DCBHs, we show that metal pollution of A from star forming clumps clumps B and C is inevitable, and that A can form a DCBH well before its metallicity exceeds the critical threshold of $10^{-5}-10^{-6}\ \rm Z_{\odot}$. Assuming metal mixing happens instantaneously, we derive the metallicity of A based on the star formation history of B and C. We find that treating a final accreting black hole of $10^6-10^7\ \rm M_{\odot}$ in A for nebular emission already pushes its $H_{160}$ - [3.6] and [3.6]-[4.5] colours into the 3$\sigma$ limit of observations. Hence, we show that the presence of metals in DCBH hosts is inevitable, and that it is the coevolution of the LW radiation field and metals originating from neighbouring galaxies that governs DCBH formation in a neighbouring {initially} pristine atomic cooling haloes.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

The selection of high redshift sources from broad-band photometry using the Lyman-break galaxy (LBG) technique is a well established methodology, but the characterization of its contamination for the faintest sources is still incomplete. We use the optical and near-IR data from four (ultra)deep Hubble Space Telescope legacy fields to investigate the contamination fraction of LBG samples at z~5-8 selected using a colour-colour method. Our approach is based on characterizing the number count distribution of interloper sources, that is galaxies with colors similar to those of LBGs, but showing detection at wavelengths shorter than the spectral break. Without sufficient sensitivity at bluer wavelengths, a subset of interlopers may not be properly classified, and contaminate the LBG selection. The surface density of interlopers in the sky gets steeper with increasing redshift of LBG selections. Since the intrinsic number of dropouts decreases significantly with increasing redshift, this implies increasing contamination from misclassified interlopers with increasing redshift, primarily by intermediate redshift sources with unremarkable properties (intermediate ages, lack of ongoing star formation and low/moderate dust content). Using Monte Carlo simulations, we estimate that the CANDELS deep data have contamination induced by photometric scatter increasing from ~2% at z~5 to ~6% at z~8 for a typical dropout color >1 mag, with contamination naturally decreasing for a more stringent dropout selection. Contaminants are expected to be located preferentially near the detection limit of surveys, ranging from 0.1 to 0.4 contaminants per arcmin2 at J=30, depending on the field considered. This analysis suggests that the impact of contamination in future studies of z>10 galaxies needs to be carefully considered.

Using Hubble data, including new grism spectra, Oesch et al. recently identified GN-z11, an $M_\textrm{UV}$=-21.1 galaxy at $z$=11.1 (just 400Myr after the big bang). With an estimated stellar mass of $\sim$10$^9$M$_{\odot}$, this galaxy is surprisingly bright and massive, raising questions as to how such an extreme object could form so early in the Universe. Using \Meraxes{}, a semi-analytic galaxy-formation model developed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) programme, we investigate the potential formation mechanisms and eventual fate of GN-z11. The volume of our simulation is comparable to that of the discovery observations and possesses two analogue galaxies of similar luminosity to this remarkably bright system. Existing in the two most massive subhaloes at $z$=11.1 ($M_\textrm{vir}$=1.4$\times 10^{11}$M$_{\odot}$ and 6.7$\times 10^{10}$M$_{\odot}$), our model analogues show excellent agreement with all available observationally derived properties of GN-z11. Although they are relatively rare outliers from the full galaxy population at high-$z$, they are no longer the most massive or brightest systems by $z$=5. Furthermore, we find that both objects possess relatively smooth, but extremely rapid mass growth histories with consistently high star formation rates and UV luminosities at $z{>}11$, indicating that their brightness is not a transient, merger-driven feature. Our model results suggest that future wide-field surveys with the \textit{James Webb Space Telescope} may be able to detect the progenitors of GN-z11 analogues out to $z{\sim}$13--14, pushing the frontiers of galaxy-formation observations to the early phases of cosmic reionization and providing a valuable glimpse of the first galaxies to reionize the Universe on large scales.

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.

We derive H{\alpha} fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The H{\alpha} flux is inferred based on the offset between the IRAC 3.6 {\mu}m flux and that predicted from the best-fit SED. We demonstrate that the H{\alpha} flux correlates well with dust- corrected UV star-formation rate (SFR) and therefore can serve as an independent SFR indicator. However, we also find a systematic offset in the SFR_H{\alpha}/SFR_UV ratios for z ~ 4-5 galaxies relative to local relations (assuming the same dust corrections for nebular regions and stellar light). We show that we can resolve the modest tension in the inferred SFRs by assuming bluer intrinsic UV slopes (increasing the dust correction), a rising star-formation history or assuming a low metallicity stellar population with a hard ionizing spectrum (increasing the L_H{\alpha}/SFR ratio). Using H{\alpha} as a SFR indicator, we find a normalization of the star formation main sequence in good agreement with recent SED-based determinations and also derive the SFR functions at z ~ 4-8. In addition, we assess for the first time the burstiness of star formation in z ~ 4 galaxies on <100 Myr time scales by comparing UV and H{\alpha}-based sSFRs; their one-to-one relationship argues against significantly bursty star-formation histories. Further progress will be made on these results, by incorporating new results from ALMA to constrain the dust-obscured star formation in high-redshift UV-selected samples.

The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$\mu$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm \mu_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1\sigma$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm \mu_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

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.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

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.

Recent surveys have identified a seemingly ubiquitous population of galaxies with elevated [OIII]/H$\beta$ emission line ratios at $z > 1$, though the nature of this phenomenon continues to be debated. The [OIII]/H$\beta$ line ratio is of interest because it is a main component of the standard diagnostic tools used to differentiate between active galactic nuclei (AGN) and star-forming galaxies, as well as the gas-phase metallicity indicators $O_{23}$ and $R_{23}$. Here, we investigate the primary driver of increased [OIII]/H$\beta$ ratios by median-stacking rest-frame optical spectra for a sample of star-forming galaxies in the 3D-HST survey in the redshift range $z\sim1.4-2.2$. Using $N = 4220$ star-forming galaxies, we stack the data in bins of mass and specific star formation rates (sSFR) respectively. After accounting for stellar Balmer absorption, we measure [OIII]$\lambda5007$\AA/H$\beta$ down to $\mathrm{M} \sim 10^{9.2} \ \mathrm{M_\odot}$ and sSFR $\sim 10^{-9.6} \ \mathrm{yr}^{-1}$, more than an order of magnitude lower than previous work at similar redshifts. We find an offset of $0.59\pm0.05$ dex between the median ratios at $z\sim2$ and $z\sim0$ at fixed stellar mass, in agreement with existing studies. However, with respect to sSFR, the $z \sim 2$ stacks all lie within 1$\sigma$ of the median SDSS ratios, with an average offset of only $-0.06\pm 0.05$. We find that the excitation properties of galaxies are tightly correlated with their sSFR at both $z\sim2$ and $z\sim0$, with a relation that appears to be roughly constant over the last 10 Gyr of cosmic time.

We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 7<z<9. To investigate this result in more detail, we have initiated a campaign to target UV metal emission in the four Lyman-alpha emitters as a probe of both the radiation field and the velocity offset of Lyman-alpha. Here we present the detection of intense CIII] emission in EGS-zs8-1, a galaxy from this sample previously shown to have Lyman-alpha at z=7.73. Photoionization models indicate that an intense radiation field and low metallicity are required to reproduce the intense CIII] and optical line emission. We argue that this extreme radiation field is likely to affect the local environment, increasing the transmission of Lyman-alpha through the galaxy. Moreover, the centroid of CIII] indicates that Lyman-alpha is redshifted from the systemic value by 340 km/s. This velocity offset is larger than that seen in less luminous systems, providing an additional explanation for the transmission of Lyman-alpha emission through the IGM. Since the transmission is further enhanced by the likelihood that such systems are also situated in the densest regions with the largest ionized bubbles, the visibility of Lyman-alpha at z>7 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation.

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.

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.

We present K-band spectra of rest-frame optical emission lines for 24 star-forming galaxies at z~3.2-3.7 using MOSFIRE on the Keck 1 telescope. Strong rest-frame optical [O III] and Hbeta emission lines were detected in 18 LBGs. The median flux ratio of [O III]5007 to Hbeta is 5.1+/-0.5, a factor of 5-10x higher than in local galaxies with similar stellar masses. The observed Hbeta luminosities are in good agreement with expectations from the estimated star-formation rates, and none of our sources are detected in deep X-ray stacks, ruling out significant contamination by active galactic nuclei. Combining our sample with a variety of LBGs from the literature, including 49 galaxies selected in a very similar manner, we find a high median ratio of [OIII]/Hbeta = 4.8+0.8-1.7. This high ratio seems to be an ubiquitous feature of z~3-4 LBGs, very different from typical local star-forming galaxies at similar stellar masses. The only comparable systems at z~0 are those with similarly high specific star-formation rates, though ~5x lower stellar masses. High specific star-formation rates either result in a much higher ionization parameter or other unusual conditions for the interstellar medium, which result in a much higher [OIII]/Hbeta line ratio. This implies a strong relation between a global property of a galaxy, the specific star-formation rate, and the local conditions of ISM in star-forming regions.

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.

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.

We present reduced data and data products from the 3D-HST survey, a 248-orbit HST Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper (Skelton et al. 2014) we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N. The data analysis is complicated by the fact that no slits are used: all objects in the WFC3 field are dispersed, and many spectra overlap. We developed software to automatically and optimally extract interlaced 2D and 1D spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has 98,663 measured redshifts and line strengths down to $JH_{IR}\leq 26$ and 22,548 with $JH_{IR}\leq 24$, where we comfortably detect continuum emission. Of this sample 5,459 galaxies are at $z>1.5$ and 9,621 are at $0.7<z<1.5$, where H$\alpha$ falls in the G141 wavelength coverage. Based on comparisons with ground-based spectroscopic redshifts, and on analyses of paired galaxies and repeat observations, the typical redshift error for $JH_{IR}\leq 24$ galaxies in our catalog is $\sigma_z \approx 0.003 \times (1+z)$, i.e., one native WFC3 pixel. The $3\sigma$ limit for emission line fluxes of point sources is $1.5\times10^{-17}$ ergs s$^{-1}$ cm$^{-2}$. We show various representations of the full dataset, as well as individual examples that highlight the range of spectra that we find in the survey.

Using the most recent measurements of the ultraviolet (UV) luminosity functions (LFs) and dust estimates of early galaxies, we derive updated dust-corrected star-formation rate functions (SFRFs) at z~4-8, which we model to predict the evolution to higher redshifts, z>8. We employ abundance matching techniques to calibrate a relation between galaxy star formation rate (SFR) and host halo mass M$_{h}$ by mapping the shape of the observed SFRFs at z~4-8 to that of the halo mass function. The resulting scaling law remains roughly constant over this redshift range. We apply the average SFR-M$_{h}$ relation to reproduce the observed SFR functions at 4 <= z <= 8 and also derive the expected UV LFs at higher redshifts. At z~9 and z~10 these model LFs are in excellent agreement with current observed estimates. Our predicted number densities and UV LFs at z>10 indicate that JWST will be able to detect galaxies out to z~15 with an extensive treasury sized program. We also derive the redshift evolution of the star formation rate density and associated reionization history by galaxies. Models which integrate down to the current HUDF12/XDF detection limit (M$_{UV}$ ~ -17.7 mag) result in a SFRD that declines as (1+z)$^{-10.4\pm0.3}$ at high redshift and fail to reproduce the observed CMB electron scattering optical depth, $\tau$ ~ 0.066, to within 1$\sigma$. On the other hand, we find that the inclusion of galaxies with SFRs well below the current detection limit (M$_{UV}$ < -5.7 mag) leads to a fully reionized universe by z ~ 6.5 and an optical depth of $\tau$ ~ 0.054, consistent with the recently derived Planck value at the 1$\sigma$ level.

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes.

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0<z<2.5. They have a median NII/Halpha ratio of 0.6, are highly obscured with SFR(tot)/SFR(Halpha)~10, and have a large range of observed line widths. We infer from the kinematics and spatial distribution of Halpha that the galaxies have rotating disks of ionized gas that are a factor of ~2 more extended than the stellar distribution. By combining measurements of individual galaxies, we find that the kinematics are consistent with a nearly Keplerian fall-off from V_rot~500 km/s at 1 kpc to V_rot~250 km/s at 7 kpc, and that the total mass out to this radius is dominated by the dense stellar component. Next, we study the size and mass evolution of the progenitors of compact massive galaxies. Even though individual galaxies may have had complex histories with periods of compaction and mergers, we show that the population of progenitors likely followed a simple inside-out growth track in the size-mass plane of d(log r_e) ~ 0.3 d(log M_stars). This mode of growth gradually increases the stellar mass within a fixed physical radius, and galaxies quench when they reach a stellar density or velocity dispersion threshold. As shown in other studies, the mode of growth changes after quenching, as dry mergers take the galaxies on a relatively steep track in the size-mass plane.

Recent Spitzer/InfraRed Array Camera (IRAC) photometric observations have revealed that rest-frame optical emission lines contribute signficantly to the broadband fluxes of high-redshift galaxies. Specifically, in the narrow redshift range z~5.1-5.4 the [3.6]-[4.5] color is expected to be very red, due to contamination of the 4.5-micron band by the dominant Halpha line, while the 3.6-micron filter is free of nebular emission lines. We take advantage of new reductions of deep Spitzer/IRAC imaging over the GOODS-North+South fields (Labbe+2015) to obtain a clean measurement of the mean Halpha equivalent width from the [3.6]-[4.5] color in the redshift range z=5.1-5.4. The selected sources either have measured spectroscopic redshifts (13 sources) or lie very confidently in the redshift range z=5.1-5.4 based on the photometric redshift likelihood intervals (11 sources). Our z_{phot}=5.1-5.4 sample and z_{spec}=5.10-5.40 spectroscopic sample have a mean [3.6]-[4.5] color of 0.31+/-0.05 mag and 0.35+/-0.07 mag, implying a rest-frame equivalent width EW(Halpha+[NII]+[SII]) of 665+/-53 Angstroms and 707+/-74 Angstroms, respectively, for sources in these samples. These values are consistent albeit slightly higher than derived by Stark+2013 at z~4, suggesting an evolution to higher values of the Halpha+[NII]+[SII] EW at z>2. Using the 3.6micron band, which is free of emission line contamination, we perform robust SED fitting and find a median specific star formation rate of sSFR = 17_{-5}^{+2} Gyr^{-1}, 7_{-2}^{+1}x higher than at z~2. We find no strong correlation (<2sigma) between the Halpha+[NII]+[SII] EW and the stellar mass of sources. Before the advent of JWST, improvements in these results will come through an expansion of current spectroscopic samples and deeper Spitzer/IRAC measurements.

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.

We present a detailed analysis of an individual case of gravitational lensing of a $z\sim8$ Lyman-Break galaxy (LBG) in a blank field, identified in Hubble Space Telescope imaging obtained as part of the Brightest of Reionizing Galaxies survey. To investigate the close proximity of the bright ($m_{AB}=25.8$) $Y_{098}$-dropout to a small group of foreground galaxies, we obtained deep spectroscopy of the dropout and two foreground galaxies using VLT/X-Shooter. We detect H-$\alpha$, H-$\beta$, [OIII] and [OII] emission in the brightest two foreground galaxies (unresolved at the natural seeing of $0.8$ arcsec), placing the pair at $z=1.327$. We can rule out emission lines contributing all of the observed broadband flux in $H_{160}$ band at $70\sigma$, allowing us to exclude the $z\sim8$ candidate as a low redshift interloper with broadband photometry dominated by strong emission lines. The foreground galaxy pair lies at the peak of the luminosity, redshift and separation distributions for deflectors of strongly lensed $z\sim8$ objects, and we make a marginal detection of a demagnified secondary image in the deepest ($J_{125}$) filter. We show that the configuration can be accurately modelled by a singular isothermal ellipsoidal deflector and a S\'{e}rsic source magnified by a factor of $\mu=4.3\pm0.2$. The reconstructed source in the best-fitting model is consistent with luminosities and morphologies of $z\sim8$ LBGs in the literature. The lens model yields a group mass of $9.62\pm0.31\times10^{11} M_{\odot}$ and a stellar mass-to-light ratio for the brightest deflector galaxy of $M_{\star}/L_{B}=2.3^{+0.8}_{-0.6} M_{\odot}/L_{\odot}$ within its effective radius. The foreground galaxies' redshifts would make this one of the few strong lensing deflectors discovered at $z>1$.

We report the discovery of Lyman-alpha emission (Ly$\alpha$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 {\AA} which we argue is likely Ly$\alpha$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5\sigma$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 {\AA}. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$\alpha$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.

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.

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.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

Recently, a small sample of six $z\sim9-10$ candidates was discovered in CANDELS that are $\sim10-20\times$ more luminous than any of the previous $z\sim9-10$ galaxies identified over the HUDF/XDF and CLASH fields. We measure the sizes of these candidates to map out the size evolution of galaxies from the earliest observable times. Their sizes are also used to provide a valuable constraint on whether these unusual galaxy candidates are at high redshift. Using galfit to derive sizes from the CANDELS F160W images of these candidates, we find a mean size of 0.13$\pm$0.02" (or 0.5$\pm$0.1 kpc at $z\sim9-10$). This handsomely matches the 0.6 kpc size expected extrapolating lower redshift measurements to $z\sim9-10$, while being much smaller than the 0.59" mean size for lower-redshift interlopers to $z\sim9-10$ photometric selections lacking the blue IRAC color criterion. This suggests that source size may be an effective constraint on contaminants from $z\sim9-10$ selections lacking IRAC data. Assuming on the basis of the strong photometric evidence that the Oesch et al. 2014 sample is entirely at $z\sim9-10$, we can use this sample to extend current constraints on the size-luminosity, size-mass relation, and size evolution of galaxies to $z\sim10$. We find that the $z\sim9-10$ candidate galaxies have broadly similar sizes and luminosities as $z\sim6$-8 counterparts with star-formation-rate surface densities in the range of $\rm \Sigma_{SFR}=1-20\, M_\odot~ yr^{-1}\, kpc^{-2}$. The stellar mass-size relation is uncertain, but shallower than those inferred for lower-redshift galaxies. In combination with previous size measurements at z=4-7, we find a size evolution of $(1+z)^{-m}$ with $m=1.0\pm0.1$ for $>0.3L^*_{z=3}$ galaxies, consistent with the evolution previously derived from $2 < z < 8$ galaxies.

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.

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by gravitational lensing induced by foreground deflectors not only in galaxy clusters, but also in blank fields. We quantify the impact of strong magnification in the samples of $B$, $V$, $i$, $z$ $\&$ $Y$ LBGs ($4\lesssim z \lesssim8$) observed in the XDF and GOODS/CANDELS fields, by investigating the proximity of dropouts to foreground objects. We find that $\sim6\%$ of bright LBGs ($m_{H_{160}}<26$) at $z\sim7$ have been strongly lensed ($\mu>2$) by foreground objects. This fraction decreases from $\sim 3.5\%$ at $z\sim6$ to $\sim1.5\%$ at $z\sim4$. Since the observed fraction of strongly lensed galaxies is a function of the shape of the luminosity function (LF), it can be used to derive Schechter parameters, $\alpha$ and $M_{\star}$, independently from galaxy number counts. Our magnification bias analysis yields Schechter-function parameters in close agreement with those determined from galaxy counts albeit with larger uncertainties. Extrapolation of our analysis to $z\gtrsim 8$ suggests that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, we highlight how the magnification bias measurement near the detection limit can be used as probe of the population of galaxies too faint to be detected. Preliminary results using this novel idea suggest that the magnification bias at $M_{UV}\sim -18$ is not as strong as expected if $\alpha\lesssim -1.7$ extends well below the current detection limits in the XDF. At face value this implies a flattening of the LF at $M_{UV}\gtrsim-16.5$. However, selection effects and completeness estimates are difficult to quantify precisely. Thus, we do not rule out a steep LF extending to $M_{UV}\gtrsim -15$.

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.

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.

We measure the UV-continuum slope beta for over 4000 high-redshift galaxies over a wide range of redshifts z~4-8 and luminosities from the HST HUDF/XDF, HUDF09-1, HUDF09-2, ERS, CANDELS-N, and CANDELS-S data sets. Our new beta results reach very faint levels at z~4 (-15.5 mag: 0.006 L*(z=3)), z~5 (-16.5 mag: 0.014L*(z=3)), and z~6 and z~7 (-17 mag: 0.025 L*(z=3)). Inconsistencies between previous studies led us to conduct a comprehensive review of systematic errors and develop a new technique for measuring beta that is robust against biases that arise from the impact of noise. We demonstrate, by object-by-object comparisons, that all previous studies, including our own and those done on the latest HUDF12 dataset, suffer from small systematic errors in beta. We find that after correcting for the systematic errors (typically d(beta) ~0.1-0.2) all beta results at z~7 from different groups are in excellent agreement. The mean beta we measure for faint (-18 mag: 0.1L*(z=3)) z~4, z~5, z~6, and z~7 galaxies is -2.03+/-0.03+/-0.06 (random and systematic errors), -2.14+/-0.06+/-0.06, -2.24+/-0.11+/-0.08, and -2.30+/-0.18+/-0.13, respectively. Our new beta values are redder than we have reported in the past, but bluer than other recent results. Our previously reported trend of bluer beta's at lower luminosities is confirmed, as is the evolution to bluer beta's at high redshifts. beta appears to show only a mild luminosity dependence faintward of M(UV,AB) ~ -19 mag, suggesting that the mean beta asymptotes to ~ -2.2 to -2.4 for faint z>~4 galaxies. At z~7, the observed beta's suggest non-zero, but low dust extinction, and they agree well with values predicted in cosmological hydrodynamical simulations.

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<z<6$, followed by a slight rise in halo masses at $z\sim7$ ($\sim2\sigma$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7\pm1.2$) implying a constant mass-to-light ratio $\frac{dlogM}{dlogL}\sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $\epsilon_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.

The dense interiors of massive galaxies are among the most intriguing environments in the Universe. In this paper we ask when these dense cores were formed and determine how galaxies gradually assembled around them. We select galaxies that have a stellar mass >3x10^10 Msun inside r=1 kpc out to z=2.5, using the 3D-HST survey and data at low redshift. Remarkably, the number density of galaxies with dense cores appears to have decreased from z=2.5 to the present. This decrease is probably mostly due to stellar mass loss and the resulting adiabatic expansion, with some contribution from merging. We infer that dense cores were mostly formed at z>2.5, consistent with their largely quiescent stellar populations. While the cores appear to form early, the galaxies in which they reside show strong evolution: their total masses increase by a factor of 2-3 from z=2.5 to z=0 and their effective radii increase by a factor of 5-6. As a result, the contribution of dense cores to the total mass of the galaxies in which they reside decreases from ~50% at z=2.5 to ~15% at z=0. Because of their early formation, the contribution of dense cores to the total stellar mass budget of the Universe is a strong function of redshift. The stars in cores with M_1kpc>3x10^10 Msun make up ~0.1% of the stellar mass density of the Universe today but 10%-20% at z~2, depending on their IMF. The formation of these cores required the conversion of ~10^11 Msun of gas into stars within ~1 kpc, while preventing significant star formation at larger radii.

The 3D-HST and CANDELS programs have provided WFC3 and ACS spectroscopy and photometry over ~900 square arcminutes in five fields: AEGIS, COSMOS, GOODS-North, GOODS-South, and the UKIDSS UDS field. All these fields have a wealth of publicly available imaging datasets in addition to the HST data, which makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range. In this paper we describe a photometric analysis of the CANDELS and 3D-HST HST imaging and the ancillary imaging data at wavelengths 0.3um to 8um. Objects were selected in the WFC3 near-IR bands, and their SEDs were determined by carefully taking the effects of the point spread function in each observation into account. A total of 147 distinct imaging datasets were used in the analysis. The photometry is made available in the form of six catalogs: one for each field, as well as a master catalog containing all objects in the entire survey. We also provide derived data products: photometric redshifts, determined with the EAZY code, and stellar population parameters determined with the FAST code. We make all the imaging data that were used in the analysis available, including our reductions of the WFC3 imaging in all five fields. 3D-HST is a spectroscopic survey with the WFC3 and ACS grisms, and the photometric catalogs presented here constitute a necessary first step in the analysis of these grism data. All the data presented in this paper are available through the 3D-HST website.

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $\alpha = -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} \phi^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} \epsilon [\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$\sigma$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

We present direct observational evidence for star formation quenching in galaxy groups in the redshift range 0<z<2.5. We utilize a large sample of nearly 6000 groups, selected by fixed cumulative number density from three photometric catalogs, to follow the evolving quiescent fractions of central and satellite galaxies over roughly 11 Gyr. At z~0, central galaxies in our sample range in stellar mass from Milky Way/M31 analogs (M=6.5x10^10 M\solar) to nearby massive ellipticals (M=1.5x10^11 M\solar). Satellite galaxies in the same groups reach masses as low as twice that of the Large Magellanic Cloud (M=6.5x10^9 M\solar). Using statistical background subtraction, we measure the average rest-frame colors of galaxies in our groups and calculate the evolving quiescent fractions of centrals and satellites over seven redshift bins. Our analysis shows clear evidence for star formation quenching in group halos, with a different quenching onset for centrals and their satellite galaxies. Using halo mass estimates for our central galaxies, we find that star formation shuts off in centrals when typical halo masses reach between 10^12 and 10^13 M\solar, consistent with predictions from the halo quenching model. In contrast, satellite galaxies in the same groups most likely undergo quenching by environmental processes, whose onset is delayed with respect to their central galaxy. Although star formation is suppressed in all galaxies over time, the processes that govern quenching are different for centrals and satellites. While mass plays an important role in determining the star formation activity of central galaxies, quenching in satellite galaxies is dominated by the environment in which they reside.

[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.

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IRAC data increase the coverage at 3.6 micron and 4.5 micron to ~120h over the HUDF reaching depths of ~28 (AB,1 sigma). The improved depth and inclusion of brighter candidates result in direct >3 sigma IRAC detections of 20/63 sources, of which 11/63 are detected at > 5 sigma. The average [3.6]-[4.5] colors of IRAC detected galaxies at z~8 are markedly redder than those at z~7, observed only 130Myr later. The simplest explanation is that we witness strong rest-frame optical emission lines (in particular [OIII]4959,5007+Hbeta) moving through the IRAC bandpasses with redshift. Assuming that the average rest-frame spectrum is the same at both z~7 and z~8 we estimate a rest-frame equivalent width of W([OIII]4959,5007+Hbeta) = 670 (+260,-170) Angstrom contributing 0.56 (+0.16,-0.11) mag to the 4.5 micron filter at z~8. The corresponding W(Halpha) = 430 (+160,-110) Angstrom implies an average specific star formation rate of sSFR = 11 (+11,-5) Gyr^-1 and a stellar population age of 100 (+100,-50) Myr. Correcting the spectral energy distribution for the contribution of emission lines lowers the average best-fit stellar masses and mass-to-light ratios by x3, decreasing the integrated stellar mass density to rho*(z=8,MUV<-18)=0.6 (+0.4,-0.3) x 10^6 Msun Mpc^-3.

We measure the evolution of the specific star formation rate (sSFR = SFR / Mstellar) between redshift 4 and 6 to investigate the previous reports of "constant" sSFR at z>2. We obtain photometry on a large sample of galaxies at z~4-6 located in the GOODS-S field that have high quality imaging from HST and Spitzer. We have derived stellar masses and star formation rates (SFRs) through stellar population modeling of their spectral energy distributions (SEDs). We estimate the dust extinction from the observed UV colors. In the SED fitting process we have studied the effects of assuming a star formation history (SFH) both with constant SFR and one where the SFR rises exponentially with time. The latter SFH is chosen to match the observed evolution of the UV luminosity function. We find that neither the mean SFRs nor the mean stellar masses change significantly when the rising SFR (RSF) model is assumed instead of the constant SFR model. When focusing on galaxies with Mstar ~ 5x10^9 Msun, we find that the sSFR evolves weakly with redshift (sSFR(z) \propto (1+z)^(0.6+/-0.1) Gyr^-1), consistent with previous results and with recent estimates of the sSFR at z~2-3 using similar assumptions. We have also investigated the impact of optical emission lines on our results. We estimate that the contribution of emission lines to the rest-frame optical fluxes is only modest at z~4 and 5 but it could reach ~50% at z~6. When emission lines of this strength are taken into account, the sSFR shows somewhat higher values at high redshifts, according to the relation sSFR(z) \propto (1+z)^(1.0+/-0.1) Gyr^-1, i.e., ~2.3x higher at z~6 than at z~2. However, the observed evolution is substantially weaker than that found at z<2 or that expected from current models (which corresponds to sSFR(z) \propto (1+z)^(2.5) Gyr^-1). -abridged-

We present the analysis of photometric, spectroscopic, and morphological properties for differently selected samples of passive galaxies up to z=1 extracted from the zCOSMOS-20k spectroscopic survey. This analysis intends to explore the dependence of galaxy properties on the selection criterion adopted, study the degree of contamination due to star-forming outliers, and provide a comparison between different commonly used selection criteria. We extracted from the zCOSMOS-20k catalog six different samples of passive galaxies, based on morphology, optical colors, specific star-formation rate, a best fit to the observed spectral energy distribution, and a criterion that combines morphological, spectroscopic, and photometric information. The morphological sample has the higher percentage of contamination in colors, specific star formation rate and presence of emission lines, while the red & passive ETGs sample is the purest, with properties mostly compatible with no star formation activity; however, it is also the less economic criterion in terms of information used. The best performing among the other criteria are the red SED and the quiescent ones, providing a percentage of contamination only slightly higher than the red & passive ETGs criterion (on average of a factor of ~2) but with absolute values of the properties of contaminants still compatible with a red, passively evolving population. We also provided two revised definitions of early type galaxies based on restframe color-color and color-mass criteria, that better reproduce the observed bimodalities. The analysis of the number densities shows evidences of mass-assembly downsizing, with galaxies at 10.25<log(M/Msun)<10.75 increasing their number by a factor ~2-4 from z=0.6 to z=0.2, by a factor ~2-3 from z=1 to z=0.2 at 10.75<log(M/Msun)<11, and by only ~10-50% from z=1 to z=0.2 at 11<log(M/Msun)<11.5.

The eXtreme Deep Field (XDF) combines data from ten years of observations with the HST Advanced Camera for Surveys (ACS) and the Wide-Field Camera 3 Infra-Red (WFC3/IR) into the deepest image of the sky ever in the optical/near-IR. Since the initial observations on the Hubble Ultra-Deep Field (HUDF) in 2003, numerous surveys and programs, including supernova followup, HUDF09, CANDELS, and HUDF12 have contributed additional imaging data across the HUDF region. Yet these have never been combined and made available as one complete ultra-deep optical and near-infrared image dataset. We do so now for the eXtreme Deep Field (XDF) program. Our new and improved processing techniques provide higher quality reductions of the total dataset. All WFC3 near-IR and optical ACS data sets have been fully combined and accurately matched, resulting in the deepest imaging ever taken at these wavelengths ranging from 29.1 to 30.3 AB mag (5sigma in a 0.35" diameter aperture) in 9 filters. The gains in the optical for the four filters done in the original ACS HUDF correspond to a typical improvement of 0.15 mag, with gains of 0.25 mag in the deepest areas. Such gains are equivalent to adding ~130 to ~240 orbits of ACS data to the HUDF. Improved processing alone results in a typical gain of ~0.1 mag. Our 5sigma (optical+near-IR) SExtractor catalogs reveal about 14140 sources in the full field and about 7121 galaxies in the deepest part of the XDF (the HUDF09 region). The XDF is the deepest image of the universe ever taken, reaching, in the combined image for a flat f_nu source, to 31.2 AB mag 5sigma (32.9 at 1sigma) in a 0.35" diameter aperture.

We show that a bursty model of high redshift star formation explains several puzzling observations of the high redshift galaxy population. We begin by pointing out that the observed specific star formation rate requires a duty-cycle of ~10%, which is much lower than found in many hydro-dynamical simulations. This value follows directly from the fact that the observed star formation rate in galaxies integrated over a Hubble time exceeds the observed stellar mass by an order of magnitude. We use the large observed specific star formation rate to calibrate the efficiency of feedback in a model for the high redshift star formation rate which includes merger driven star formation regulated by SNe feedback. This model reproduces the star formation rate density function and the stellar mass function of galaxies at 4<z<7. A prediction of the model is that the specific star formation rate does not evolve very rapidly with either mass or redshift, in agreement with observation. This is in contrast to results from hydrodynamical simulations where the star formation closely follows the accretion rate, and so increases strongly towards high redshift. The bursty star formation model naturally explains the observation that there is not enough stellar mass at z~2-4 to account for all of the star-formation observed, without invoking properties like an evolving initial mass function of stars. The finding of a duty cycle that is ~10% implies that there should be ten times the number of known galaxies at fixed stellar mass that have not yet been detected through standard UV selection at high redshift. We therefore predict the existence of a large undetected population of UV-faint galaxies that accounts for most of the stellar mass density at z~4-8.

We present an analysis of cool outflowing gas around galaxies, traced by MgII absorption lines in the co-added spectra of a sample of 486 zCOSMOS galaxies at 1 < z < 1.5. These galaxies span a range of stellar masses (9.45< log[M*/Msun]<10.7) and star formation rates (0.14 < log [SFR/Msun/yr] < 2.35). We identify the cool outflowing component in the MgII absorption and find that the equivalent width of the outflowing component increases with stellar mass. The outflow equivalent width also increases steadily with the increasing star formation rate of the galaxies. At similar stellar masses the blue galaxies exhibit a significantly higher outflow equivalent width as compared to red galaxies. The outflow equivalent width shows strong effect with star formation surface density ({\Sigma}SFR) of the sample. For the disk galaxies, the outflow equivalent width is higher for the face-on systems as compared to the edge-on ones, indicating that for the disk galaxies, the outflowing gas is primarily bipolar in geometry. Galaxies typically exhibit outflow velocities ranging from -200 km/s to -300 km/s and on average the face-on galaxies exhibit higher outflow velocity as compared to the edge-on ones. Galaxies with irregular morphologies exhibit outflow equivalent width as well as outflow velocities comparable to face on disk galaxies. These galaxies exhibit minimum mass outflow rates > 5-7 Msun/yr and a mass loading factor ({\eta} = dMout/dt /SFR) comparable to the star formation rates of the galaxies.

We explore the role of environment in the evolution of galaxies over 0.1<z<0.7 using the final zCOSMOS-bright data set. Using the red fraction of galaxies as a proxy for the quenched population, we find that the fraction of red galaxies increases with the environmental overdensity and with the stellar mass, consistent with previous works. As at lower redshift, the red fraction appears to be separable in mass and environment, suggesting the action of two processes: mass and environmental quenching. The parameters describing these appear to be essentially the same at z~0.7 as locally. We explore the relation between red fraction, mass and environment also for the central and satellite galaxies separately, paying close attention to the effects of impurities in the central-satellite classification and using carefully constructed samples matched in stellar mass. There is little evidence for a dependence of the red fraction of centrals on overdensity. Satellites are consistently redder at all overdensities, and the satellite quenching efficiency increases with overdensity at 0.1<z<0.4. This is less marked at higher redshift, but both are nevertheless consistent with the equivalent local measurements. At a given stellar mass, the fraction of galaxies that are satellites also increases with the overdensity. At a given overdensity and mass, the obtained relation between the environmental quenching and the satellite fraction agrees well with the satellite quenching efficiency, demonstrating that the environmental quenching in the overall population is consistent with being entirely produced through the satellite quenching process at least up to z=0.7. However, despite the unprecedented size of our high redshift samples, the associated statistical uncertainties are still significant and our statements should be understood as approximations to physical reality, rather than physically exact formulae.

We investigate the link between AGN activity, star-formation and stellar mass of the host galaxy at 0<z<1, looking for spectroscopic traces of AGN and aging of the host. This work provides an extension of the existing studies at z<0.1 and contributes to shed light on galaxy evolution at intermediate redshifts. We used the zCOSMOS 20k data to create a sample of galaxies at z<1. We divided the sample in several mass-redshift bins to obtain stacked galaxy spectra with an improved S/N. We exploited emission-line diagnostic diagrams to separate AGN from star-forming galaxies. We found indication of a role for the total galaxy stellar mass in leading galaxy classification. Stacked spectra show AGN signatures above the log M_*/M_sun>10.2 threshold. Moreover, the stellar populations of AGN hosts are found to be older with respect to star-forming and composites galaxies. This could be due to the the tendency of AGN to reside in massive hosts. The dependence of the AGN classification on the stellar mass is in agreement with what has been already found in previous studies. It is consistent with, together with the evidence of older stellar populations inhabiting the AGN-like galaxies, the downsizing scenario. In particular, our evidence points to an evolutionary scenario where the AGN-feedback is capable of quenching the star formation in the most massive galaxies. Therefore, the AGN-feedback is the best candidate for initiating the passive evolutionary phase of galaxies.

Ultra-deep WFC3/IR observations on the HUDF from the HUDF09 program revealed just one plausible z~10 candidate UDFj-39546284. UDFj-39546284 had all the properties expected of a galaxy at z~10 showing (1) no detection in the deep ACS+WFC3 imaging data blueward of the F160W band, exhibiting (2) a blue spectral slope redward of the break, and showing (3) no prominent detection in deep IRAC observations. The new, similarly deep WFC3/IR HUDF12 F160W observations over the HUDF09/XDF allow us to further assess this candidate. These observations show that this candidate, previously only detected at ~5.9 sigma in a single band, clearly corresponds to a real source. It is detected at ~5.3 sigma in the new H-band data and at ~7.8 sigma in the full 85-orbit H-band stack. Interestingly, the non-detection of the source (<1 sigma) in the new F140W observations suggests a higher redshift. Formally, the best-fit redshift of the source utilizing all the WFC3+ACS (and IRAC+K-band) observations is 11.8+/-0.3. However, we consider the z~12 interpretation somewhat unlikely, since the source would either need to be ~20x more luminous than expected or show very high-EW Ly-alpha emission (which seems improbable given the extensive neutral gas prevalent early in the reionization epoch). Lower-redshift solutions fail if only continuum models are allowed. Plausible lower-redshift solutions require that the H-band flux be dominated by line emission such as Halpha or [OIII] with extreme EWs. The tentative detection of line emission at 1.6 microns in UDFj-39546284 in a companion paper suggests that such emission may have already been found.

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 {\AA} break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a \Lambda CDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2.

We present a comprehensive analysis of z>8 galaxies based on ultra-deep WFC3/IR data. We constrain the evolution of the UV luminosity function (LF) and luminosity densities from z~11 to z~8 by exploiting all the WFC3/IR data over the Hubble Ultra-Deep Field from the HUDF09 and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area WFC3/IR imaging over GOODS-South. Galaxies are selected based on the Lyman Break Technique in three samples centered around z~9, z~10 and z~11, with seven z~9 galaxy candidates, and one each at z~10 and z~11. We confirm a new z~10 candidate (with z=9.8+-0.6) that was not convincingly identified in our first z~10 sample. The deeper data over the HUDF confirms all our previous z>~7.5 candidates as genuine high-redshift candidates, and extends our samples to higher redshift and fainter limits (H_160~29.8 mag). We perform one of the first estimates of the z~9 UV LF and improve our previous constraints at z~10. Extrapolating the lower redshift UV LF evolution should have revealed 17 z~9 and 9 z~10 sources, i.e., a factor ~3x and 9x larger than observed. The inferred star-formation rate density (SFRD) in galaxies above 0.7 M_sun/yr decreases by 0.6+-0.2 dex from z~8 to z~9, in good agreement with previous estimates. The low number of sources found at z>8 is consistent with a very rapid build-up of galaxies across z~10 to z~8. From a combination of all current measurements, we find a best estimate of a factor 10x decrease in the SFRD from z~8 to z~10, following (1+z)^(-11.4+-3.1). Our measurements thus confirm our previous finding of an accelerated evolution beyond z~8, and signify a rapid build-up of galaxies with M_UV<-17.7 within only ~200 Myr from z~10 to z~8, in the heart of cosmic reionization.

We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < {\lambda} < 1.65{\mu}m. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is unique to slitless spectroscopy at HST resolution. We find a 2.7{\sigma} detection of an emission line at 1.599 {\mu}m---just redward of the JH140 filter---with flux 3.5 +/- 1.3x10-18 erg/s/cm2. Assuming the line is real, it contributes 110 +/- 40% of the observed H160 flux and has an observed equivalent width > 7300 {\AA}. If the line is confirmed, it could be Ly-{\alpha} at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift feature such as [O III]{\lambda}4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km/s of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery of a bright "[O III] blob" with a secure G141 grism redshift of z = 1.605. This object has a strikingly large observed equivalent width of nearly 9000 {\AA} that results in similar "dropout" colors as UDFj-39546284.

We present stellar mass surface density profiles of a mass-selected sample of 177 galaxies at 0.5 < z < 2.5, obtained using very deep HST optical and near-infrared data over the GOODS-South field, including recent CANDELS data. Accurate stellar mass surface density profiles have been measured for the first time for a complete sample of high-redshift galaxies more massive than 10^10.7 M_sun. The key advantage of this study compared to previous work is that the surface brightness profiles are deconvolved for PSF smoothing, allowing accurate measurements of the structure of the galaxies. The surface brightness profiles account for contributions from complex galaxy structures such as rings and faint outer disks. Mass profiles are derived using radial rest-frame u-g color profiles and a well-established empirical relation between these colors and the stellar mass-to-light ratio. We derive stellar half-mass radii from the mass profiles, and find that these are on average ~25% smaller than rest-frame g band half-light radii. This average size difference of 25% is the same at all redshifts, and does not correlate with stellar mass, specific star formation rate, effective surface density, Sersic index, or galaxy size. Although on average the difference between half-mass size and half-light size is modest, for approximately 10% of massive galaxies this difference is more than a factor two. These extreme galaxies are mostly extended, disk-like systems with large central bulges. These results are robust, but could be impacted if the central dust extinction becomes high. ALMA observations can be used to explore this possibility. These results provide added support for galaxy growth scenarios wherein massive galaxies at these epochs grow by accretion onto their outer regions.

In this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M_star >= 10^11 M_Sun) in the COSMOS field, and study their role in mass and size evolution. We use the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close pairs with a separation 10h^-1 kpc <= r_p <= 30h-1 kpc and a relative velocity Delta v <= 500 km s^-1. We measure both major (stellar mass ratio mu = M_star,2/M_star,1 >= 1/4) and minor (1/10 <= mu < 1/4) merger fractions of massive galaxies, and study their dependence on redshift and on morphology. The merger fraction and rate of massive galaxies evolves as a power-law (1+z)^n, with major mergers increasing with redshift, n_MM = 1.4, and minor mergers showing little evolution, n_mm ~ 0. When split by their morphology, the minor merger fraction for early types is higher by a factor of three than that for spirals, and both are nearly constant with redshift. Our results show that massive early-type galaxies have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a mass growth of ~30%. We find that mu >= 1/10 mergers can explain ~55% of the observed size evolution of these galaxies since z ~ 1. Another ~20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (mu < 1/10) could contribute with an extra ~20%. The remaining ~5% should come from other processes (e.g., adiabatic expansion or observational effects). This picture also reproduces the mass growth and velocity dispersion evolution of these galaxies. We conclude from these results that merging is the main contributor to the size evolution of massive ETGs at z <= 1, accounting for ~50-75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (mu < 1/4) events.

We present a study of rest-frame UV-to-optical color distributions for z~4 galaxies based on the combination of deep HST/ACS+WFC3/IR data with Spitzer/IRAC imaging. In particular, we use new, ultra-deep data from the IRAC Ultradeep Field program (IUDF10). Our sample contains a total of ~2600 galaxies selected as B-dropout Lyman Break Galaxies (LBGs) in the HUDF and one of its deep parallel fields, the HUDF09-2, as well as GOODS-North and South. This sample is used to investigate the UV continuum slopes beta and Balmer break colors (J_125-[4.5]) as a function of rest-frame optical luminosity. The [4.5] filter is chosen to avoid potential contamination by strong rest-frame optical emission lines. We find that galaxies at M_z<-21.5 (roughly corresponding to L*[z~4]) are significantly redder than their lower luminosity counterparts. The UV continuum slopes and the J_125-[4.5] colors are well correlated. The most simple explanation for this correlation is that the dust reddening at these redshifts is better described by an SMC-like extinction curve, rather than the typically assumed Calzetti reddening. After correcting for dust, we find that the galaxy population shows mean stellar population ages in the range 10^8.5 to 10^9 yr, with a dispersion of ~0.5 dex, and only weak trends as a function of luminosity. In contrast to some results from the literature, we find that only a small fraction of galaxies shows Balmer break colors which are consistent with extremely young ages, younger than 100 Myr. Under the assumption of smooth star-formation histories, this fraction is only 12-19% for galaxies at M_z<-19.75. Our results are consistent with a gradual build-up of stars and dust in galaxies at z>4, with only a small fraction of stars being formed in short, intense bursts of star-formation.

We present new z~8 galaxy candidates from a search over ~95 arcmin^2 of WFC3/IR data, tripling the previous search area for bright z~8 galaxies. Our analysis uses newly acquired WFC3/IR imaging data from the CANDELS Multi-Cycle Treasury program over the GOODS South field. These new data are combined with existing deep optical ACS imaging to search for relatively bright (M_UV < -19.5 mag) z~8 galaxy candidates using the Lyman Break technique. These new candidates are used to determine the bright end of the UV luminosity function (LF) of star-forming galaxies at z~8. To minimize contamination from lower redshift galaxies, we make full use of all optical ACS data and impose strict non-detection criteria based on an optical chi^2_opt flux measurement. In the whole search area we identify 16 candidate z~8 galaxies, spanning a magnitude range H_160 = 25.7-27.9 mag. The new data show that the UV LF is a factor ~1.7x lower at M_UV < -19.5 mag than determined from the HUDF09 and ERS data alone. Combining this new sample with the previous candidates from the HUDF09 and ERS data allows us to perform the most accurate measurement of the z~8 UV LF yet. Schechter function fits to the combined data result in a best-fit characteristic magnitude of M*(z=8) = -20.04+-0.46 mag. The faint-end slope is very steep, though quite uncertain, with alpha = -2.06+-0.32. A combination of wide area data with additional ultra-deep imaging will be required to significantly reduce the uncertainties on these parameters in the future.

We present an optical group catalog between 0.1 < z < 1 based on 16,500 high-quality spectroscopic redshifts in the completed zCOSMOS-bright survey. The catalog published herein contains 1498 groups in total and 192 groups with more than five observed members. The catalog includes both group properties and the identification of the member galaxies. Based on mock catalogs, the completeness and purity of groups with three and more members should be both about 83% with respect to all groups that should have been detectable within the survey, and more than 75% of the groups should exhibit a one-to-one correspondence to the "real" groups. Particularly at high redshift, there are apparently more galaxies in groups in the COSMOS field than expected from mock catalogs. We detect clear evidence for the growth of cosmic structure over the last seven billion years in the sense that the fraction of galaxies that are found in groups (in volume-limited samples) increases significantly with cosmic time. In the second part of the paper, we develop a method for associating galaxies that only have photo-z to our spectroscopically identified groups. We show that this leads to improved definition of group centers, improved identification of the most massive galaxies in the groups, and improved identification of central and satellite galaxies, where we define the former to be galaxies at the minimum of the gravitational potential wells. Subsamples of centrals and satellites in the groups can be defined with purities up to 80%, while a straight binary classification of all group and non-group galaxies into centrals and satellites achieves purities of 85% and 75%, respectively, for the spectroscopic sample.

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. \phi(L) = \phi_* (L/L_*)^{\alpha}\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive \phi_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope \alpha = -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between \phi_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe.

For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 \AA \ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer/IRAC) and UV fluxes. We show that the stellar mass-size relation of Lyman-break galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass increases towards lower redshifts. For galaxies with stellar masses of 9.5<Log(M*/Msun)<10.4 sizes evolve as $(1+z)^{-1.20\pm0.11}$. This evolution is very similar for galaxies with lower stellar masses of 8.6<Log(M*/Msun)<9.5 which is $r_{e} \propto (1+z)^{-1.18\pm0.10}$, in agreement with simple theoretical galaxy formation models at high z. Our results are consistent with previous measurements of the LBGs mass-size relation at lower redshifts (z ~ 1-3).

We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the H{\alpha} and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, \Delta\alpha\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time.

Ultra-deep ACS and WFC3/IR HUDF+HUDF09 data, along with the wide-area GOODS+ERS+CANDELS data over the CDF-S GOODS field, are used to measure UV colors, expressed as the UV-continuum slope beta, of star-forming galaxies over a wide range in luminosity (0.1L*(z=3) to 2L*(z=3)) at high redshift (z~7 to z~4). Beta is measured using all ACS and WFC3/IR passbands uncontaminated by Ly_alpha and spectral breaks. Extensive tests show that our beta measurements are only subject to minimal biases. Using a different selection procedure, Dunlop et al. recently found large biases in their beta measurements. To reconcile these different results, we simulated both approaches and found that beta measurements for faint sources are subject to large biases if the same passbands are used both to select the sources and to measure beta. High-redshift galaxies show a well-defined rest-frame UV color-magnitude (CM) relationship that becomes systematically bluer towards fainter UV luminosities. No evolution is seen in the slope of the UV CM relationship in the first 1.5 Gyr, though there is a small evolution in the zero-point to redder colors from z~7 to z~4. This suggests that galaxies are evolving along a well-defined sequence in the L(UV)-color (beta) plane (a "star-forming sequence"?). Dust appears to be the principal factor driving changes in the UV color (beta) with luminosity. These new larger beta samples lead to improved dust extinction estimates at z~4-7 and confirm that the extinction is still essentially zero at low luminosities and high redshifts. Inclusion of the new dust extinction results leads to (i) excellent agreement between the SFR density at z~4-8 and that inferred from the stellar mass density, and (ii) to higher SSFRs at z>~4, suggesting the SSFR may evolve modestly (by factors of ~2) from z~4-7 to z~2.

The HUDF09 data are the deepest near-IR observations ever, reaching to 29.5 mag. Luminosity functions (LF) from these new HUDF09 data for 132 z\sim7 and z\sim8 galaxies are combined with new LFs for z\sim5-6 galaxies and the earlier z\sim4 LF to reach to very faint limits (<0.05 L*(z=3)). The faint-end slopes alpha are steep: -1.79+/-0.12 (z\sim5), -1.73+/-0.20 (z\sim6), -2.01+/-0.21 (z\sim7), and -1.91+/-0.32 (z\sim8). Slopes alpha\lesssim-2 lead to formally divergent UV fluxes, though galaxies are not expected to form below \sim-10 AB mag. These results have important implications for reionization. The weighted mean slope at z\sim6-8 is -1.87+/-0.13. For such steep slopes, and a faint-end limit of -10 AB mag, galaxies provide a very large UV ionizing photon flux. While current results show that galaxies can reionize the universe by z\sim6, matching the Thomson optical depths is more challenging. Extrapolating the current LF evolution to z>8, taking alpha to be -1.87+/-0.13 (the mean value at z\sim6-8), and adopting typical parameters, we derive Thomson optical depths of 0.061_{-0.006}^{+0.009}. However, this result will change if the faint-end slope alpha is not constant with redshift. We test this hypothesis and find a weak, though uncertain, trend to steeper slopes at earlier times (dalpha/dz\sim-0.05+/-0.04), that would increase the Thomson optical depths to 0.079_{-0.017}^{+0.063}, consistent with recent WMAP estimates (tau=0.088+/-0.015). It may thus not be necessary to resort to extreme assumptions about the escape fraction or clumping factor. Nevertheless, the uncertainties remain large. Deeper WFC3/IR+ACS observations can further constrain the ionizing flux from galaxies.

Using the group catalog obtained from zCOSMOS spectroscopic data and the complementary photometric data from the COSMOS survey, we explore segregation effects occurring in groups of galaxies at intermediate/high redshifts. We built two composite groups at intermediate (0.2 <= z <= 0.45) and high (0.45 < z <= 0.8) redshifts, and we divided the corresponding composite group galaxies into three samples according to their distance from the group center. We explored how galaxy stellar masses and colors - working in narrow bins of stellar masses - vary as a function of the galaxy distance from the group center. We found that the most massive galaxies in our sample (Log(M_gal/M_sun) >= 10.6) do not display any strong group-centric dependence of the fractions of red/blue objects. For galaxies of lower masses (9.8 <= Log(M_gal/M_sun) <= 10.6) there is a radial dependence in the changing mix of red and blue galaxies. This dependence is most evident in poor groups, whereas richer groups do not display any obvious trend of the blue fraction. Interestingly, mass segregation shows the opposite behavior: it is visible only in rich groups, while poorer groups have a a constant mix of galaxy stellar masses as a function of radius. We suggest a simple scenario where color- and mass-segregation originate from different physical processes. While dynamical friction is the obvious cause for establishing mass segregation, both starvation and galaxy-galaxy collisions are plausible mechanisms to quench star formation in groups at a faster rate than in the field. In poorer groups the environmental effects are caught in action superimposed to secular galaxy evolution. Their member galaxies display increasing blue fractions when moving from the group center to more external regions, presumably reflecting the recent accretion history of these groups.

The success of future large scale weak lensing surveys will critically depend on the accurate estimation of photometric redshifts of very large samples of galaxies. This in turn depends on both the quality of the photometric data and the photo-z estimators. In a previous study, (Bordoloi et al. 2010) we focussed primarily on the impact of photometric quality on photo-z estimates and on the development of novel techniques to construct the N(z) of tomographic bins at the high level of precision required for precision cosmology, as well as the correction of issues such as imprecise corrections for Galactic reddening. We used the same set of templates to generate the simulated photometry as were then used in the photo-z code, thereby removing any effects of "template error". In this work we now include the effects of "template error" by generating simulated photometric data set from actual COSMOS photometry. We use the trick of simulating redder photometry of galaxies at higher redshifts by using a bluer set of passbands on low z galaxies with known redshifts. We find that "template error" is a rather small factor in photo-z performance, at the photometric precision and filter complement expected for all-sky surveys. With only a small sub-set of training galaxies with spectroscopic redshifts, it is in principle possible to construct tomographic redshift bins whose mean redshift is known, from photo-z alone, to the required accuracy of 0.002(1+z).

(Abridged) We analyze the environments and galactic properties (morphologies and star-formation histories) of a sample of 153 close kinematic pairs in the redshift range 0.2 < z < 1 identified in the zCOSMOS-bright 10k spectroscopic sample of galaxies. Correcting for projection effects, the fraction of close kinematic pairs is three times higher in the top density quartile than in the lowest one. This translates to a three times higher merger rate because the merger timescales are shown, from mock catalogues based on the Millennium simulation, to be largely independent of environment once the same corrections for projection is applied. We then examine the morphologies and stellar populations of galaxies in the pairs, comparing them to control samples that are carefully matched in environment so as to remove as much as possible the well-known effects of environment on the properties of the parent population of galaxies. Once the environment is properly taken into account in this way, we find that the early-late morphology mix is the same as for the parent population, but that the fraction of irregular galaxies is boosted by 50-75%, with a disproportionate increase in the number of irregular-irregular pairs (factor of 4-8 times), due to the disturbance of disk galaxies. Future dry-mergers, involving elliptical galaxies comprise less than 5% of all close kinematic pairs. In the closest pairs, there is a boost in the specific star-formation rates of star-forming galaxies of a factor of 2-4, and there is also evidence for an increased incidence of post star-burst galaxies. Although significant for the galaxies involved, the "excess" star-formation associated with pairs represents only about 5% of the integrated star-formation activity in the parent sample. Although most pair galaxies are in dense environments, the effects of interaction appear to be largest in the lower density environments.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

We use the ultra-deep HUDF09 and the deep ERS data from the HST WFC3/IR camera, along with the wide area Spitzer/IRAC data from GOODS-S to derive SEDs of star-forming galaxies from the rest-frame UV to the optical over a wide luminosity range (M_1500 ~ -21 to M_1500 ~ -18) from z ~ 7 to z ~ 4. The sample contains ~ 400 z ~ 4, ~ 120 z ~ 5, ~ 60 z ~ 6, and 36 prior z ~ 7 galaxies. Median stacking enables the first comprehensive study of very faint high-z galaxies at multiple redshifts (e.g., [3.6] = 27.4 +/- 0.1 AB mag for the M_1500 ~ -18 sources at z ~ 4). At z ~ 4 our faint median-stacked SEDs reach to ~ 0.06 L*(z=4) and are combined with recently published results at high luminosity L > L* that extend to M_1500 ~ -23. We use the observed SEDs and template fits to derive rest frame UV-to-optical colors (U - V) at all redshifts and luminosities. We find that this color does not vary significantly with redshift at a fixed luminosity. The UV-to-optical color does show a weak trend with luminosity, becoming redder at higher luminosities. This is most likely due to dust. At z >~ 5 we find blue colors [3.6]-[4.5] ~ -0.3 mag that are most likely due to rest-frame optical emission lines contributing to the flux in the IRAC filter bandpasses. The scatter across our derived SEDs remains substantial, but the results are most consistent with a lack of any evolution in the SEDs with redshift at a given luminosity. The similarity of the SEDs suggests a self-similar mode of evolution over a timespan from 0.7 Gyr to 1.5 Gyr that encompasses very substantial growth in the stellar mass density in the universe (from ~ 4x10^6 to ~ 2x10^7 Msun Mpc^-3).

We present a photometric and spectroscopic study of galaxies at 0.5<z<1 as a function of environment based on data from the zCOSMOS survey. There is a fair amount of evidence that galaxy properties depend on mass of groups and clusters, in the sense that quiescent galaxies prefer more massive systems. We base our analysis on a mass-selected environment using X-ray groups of galaxies and define the group membership using a large number of spectroscopic redshifts from zCOSMOS. We show that the fraction of red galaxies is higher in groups than in the field at all redshifts probed in our study. Interestingly, the fraction of [OII] emitters on the red sequence increases at higher redshifts in groups, while the fraction does not strongly evolve in the field. This is due to increased dusty star formation activities and/or increased activities of active galactic nuclei (AGNs) in high redshift groups. We study these possibilities using the 30-band photometry and X-ray data. We find that the stellar population of the red [OII] emitters in groups is old and there is no clear hint of dusty star formation activities in those galaxies. The observed increase of red [OII] emitters in groups is likely due to increased AGN activities. However, our overall statistics is poor and any firm conclusions need to be drawn from a larger statistical sample of z~1 groups.

The zCOSMOS SINFONI project is aimed at studying the physical and kinematical properties of a sample of massive z~1.4-2.5 star-forming galaxies, through SINFONI near-IR integral field spectroscopy (IFS), combined with the multi-wavelength information from the zCOSMOS (COSMOS) survey. The project is based on 1 hour of natural-seeing observations per target, and Adaptive Optics (AO) follow-up for a major part of the sample, which includes 30 galaxies selected from the zCOSMOS/VIMOS spectroscopic survey. This first paper presents the sample selection, and the global physical characterization of the target galaxies from multicolor photometry, i.e., star formation rate (SFR), stellar mass, age, etc. The Halpha integrated properties such as, flux, velocity dispersion, and size, are derived from the natural-seeing observations, while the follow up AO observations will be presented in the next paper of this series. Our sample appears to be well representative of star-forming galaxies at z~2, covering a wide range in mass and SFR. The Halpha integrated properties of the 25 Halpha detected galaxies are similar to those of other IFS samples at the same redshifts. Good agreement is found among the SFRs derived from Halpha luminosity and other diagnostic methods, provided the extinction affecting the Halpha luminosity is about twice that affecting the continuum. A preliminary kinematic analysis, based on the maximum observed velocity difference across the source, and on the integrated velocity dispersion, indicates that the sample splits nearly 50-50 into rotation-dominated and velocity dispersion-dominated galaxies, in good agreement with previous surveys.

Our goal is to develop a new and reliable statistical method to classify galaxies from large surveys. We probe the reliability of the method by comparing it with a three-dimensional classification cube, using the same set of spectral, photometric and morphological parameters.We applied two different methods of classification to a sample of galaxies extracted from the zCOSMOS redshift survey, in the redshift range 0.5 < z < 1.3. The first method is the combination of three independent classification schemes, while the second method exploits an entirely new approach based on statistical analyses like Principal Component Analysis (PCA) and Unsupervised Fuzzy Partition (UFP) clustering method. The PCA+UFP method has been applied also to a lower redshift sample (z < 0.5), exploiting the same set of data but the spectral ones, replaced by the equivalent width of H$\alpha$. The comparison between the two methods shows fairly good agreement on the definition on the two main clusters, the early-type and the late-type galaxies ones. Our PCA-UFP method of classification is robust, flexible and capable of identifying the two main populations of galaxies as well as the intermediate population. The intermediate galaxy population shows many of the properties of the green valley galaxies, and constitutes a more coherent and homogeneous population. The fairly large redshift range of the studied sample allows us to behold the downsizing effect: galaxies with masses of the order of $3\cdot 10^{10}$ Msun mainly are found in transition from the late type to the early type group at $z>0.5$, while galaxies with lower masses - of the order of $10^{10}$ Msun - are in transition at later epochs; galaxies with $M <10^{10}$ Msun did not begin their transition yet, while galaxies with very large masses ($M > 5\cdot 10^{10}$ Msun) mostly completed their transition before $z\sim 1$.

We map the radial and azimuthal distribution of Mg II gas within 200 kpc (physical) of 4000 galaxies at redshifts 0.5 < z < 0.9 using co-added spectra of more than 5000 background galaxies at z > 1. We investigate the variation of Mg II rest frame equivalent width as a function of the radial impact parameter for different subsets of foreground galaxies selected in terms of their rest-frame colors and masses. Blue galaxies have a significantly higher average Mg II equivalent width at close galactocentric radii as compared to the red galaxies. Amongst the blue galaxies, there is a correlation between Mg II equivalent width and galactic stellar mass of the host galaxy. We also find that the distribution of Mg II absorption around group galaxies is more extended than that for non-group galaxies, and that groups as a whole have more extended radial profiles than individual galaxies. Interestingly, these effects can be satisfactorily modeled by a simple superposition of the absorption profiles of individual member galaxies, assuming that these are the same as those of non-group galaxies, suggesting that the group environment may not significantly enhance or diminish the Mg II absorption of individual galaxies. We show that there is a strong azimuthal dependence of the Mg II absorption within 50 kpc of inclined disk-dominated galaxies, indicating the presence of a strongly bipolar outflow aligned along the disk rotation axis. There is no significant dependence of Mg II absorption on the apparent inclination angle of disk-dominated galaxies.

We search for z~10 galaxies over ~160 arcmin^2 of WFC3/IR data in the Chandra Deep Field South, using the public HUDF09, ERS, and CANDELS surveys, that reach to 5sigma depths ranging from 26.9 to 29.4 in H_160 AB mag. z>~9.5 galaxy candidates are identified via J_125-H_160>1.2 colors and non-detections in any band blueward of J_125. Spitzer IRAC photometry is key for separating the genuine high-z candidates from intermediate redshift (z~2-4) galaxies with evolved or heavily dust obscured stellar populations. After removing 16 sources of intermediate brightness (H_160~24-26 mag) with strong IRAC detections, we only find one plausible z~10 galaxy candidate in the whole data set, previously reported in Bouwens et al. (2011). The newer data cover a 3x larger area and provide much stronger constraints on the evolution of the UV luminosity function (LF). If the evolution of the z~4-8 LFs is extrapolated to z~10, six z~10 galaxies are expected in our data. The detection of only one source suggests that the UV LF evolves at an accelerated rate before z~8. The luminosity density is found to increase by more than an order of magnitude in only 170 Myr from z~10 to z~8. This increase is >=4x larger than expected from the lower redshift extrapolation of the UV LF. We are thus likely witnessing the first rapid build-up of galaxies in the heart of cosmic reionization. Future deep HST WFC3/IR data, reaching to well beyond 29 mag, can enable a more robust quantification of the accelerated evolution around z~10.

We use the high angular resolution in the near-infrared of the WFC3 on HST to determine YHVz color-color selection criteria to identify and characterize 1.5<z<3.5 galaxies in the HUDF09 and ERS (GOODS-South) fields. The WFC3 NIR images reveal galaxies at these redshifts that were undetected in the rest-frame UV HUDF/GOODS images, as well as true centers and regular disks in galaxies classified as highly irregular in rest-frame UV light. Across the 1.5<z<2.15 redshift range, regular disks are unveiled in the WFC3 images of ~25% of both intermediate and high mass galaxies, i.e., above 10^10 Msun. Meanwhile, galaxies maintaining diffuse and/or irregular morphologies in the rest-frame optical light---i.e., not yet dynamically settled---at these epochs are almost entirely restricted to masses below 10^11 Msun. In contrast at 2.25 < z < 3.5 these diffuse and/or irregular structures overwhelmingly dominate the morphological mix in both the intermediate and high mass regimes, while no regular disks, and only a small fraction (25%) of smooth spheroids, are evident above 10^11 Msun. Strikingly, by 1.5 < z < 2.25 roughly 2 out of every 3 galaxies at the highest masses are spheroids. In our small sample, the fraction of star-forming galaxies at these mass scales decreases concurrently from ~60% to ~5%. If confirmed, this indicates that z~2 is the epoch of both the morphological transformation and quenching of star-formation which assemble the first substantial population of massive ellipticals.

We identify 73 z~7 and 59 z~8 candidate galaxies in the reionization epoch, and use this large 26-29.4 AB mag sample of galaxies to derive very deep luminosity functions to <-18 AB mag and the star formation rate density at z~7 and z~8. The galaxy sample is derived using a sophisticated Lyman-Break technique on the full two-year WFC3/IR and ACS data available over the HUDF09 (~29.4 AB mag, 5 sigma), two nearby HUDF09 fields (~29 AB mag, 14 arcmin) and the wider area ERS (~27.5 AB mag) ~40 arcmin**2). The application of strict optical non-detection criteria ensures the contamination fraction is kept low (just ~7% in the HUDF). This very low value includes a full assessment of the contamination from lower redshift sources, photometric scatter, AGN, spurious sources, low mass stars, and transients (e.g., SNe). From careful modelling of the selection volumes for each of our search fields we derive luminosity functions for galaxies at z~7 and z~8 to <-18 AB mag. The faint-end slopes alpha at z~7 and z~8 are uncertain but very steep at alpha = -2.01+/-0.21 and alpha=-1.91+/-0.32, respectively. Such steep slopes contrast to the local alpha<~-1.4 and may even be steeper than that at z~4 where alpha=-1.73+/-0.05. With such steep slopes (alpha<~-1.7) lower luminosity galaxies dominate the galaxy luminosity density during the epoch of reionization. The star formation rate densities derived from these new z~7 and z~8 luminosity functions are consistent with the trends found at later times (lower redshifts). We find reasonable consistency, with the SFR densities implied from reported stellar mass densities, being only ~40% higher at z<7. This suggests that (1) the stellar mass densities inferred from the Spitzer IRAC photometry are reasonably accurate and (2) that the IMF at very high redshift may not be very different from that at later times.

In this paper, we present a derivation of the rest-frame 1400A luminosity function (LF) at redshift six from a new application of the maximum likelihood method by exploring the five deepest HST/ACS fields, i.e., the HUDF, two UDF05 fields, and two GOODS fields. We work on the latest improved data products, which makes our results more robust than those of previous studies. We use un-binned data and thereby make optimal use of the information contained in the dataset. We focus on the analysis to a magnitude limit where the completeness is larger than 50% to avoid possibly large errors in the faint end slope that are difficult to quantify. We also take into account scattering in and out of the dropout sample due to photometric errors by defining for each object a probability that it belongs to the dropout sample. We find the best fit Schechter parameters to the z~6 LF are: alpha = 1.87 +/- 0.14, M* = -20.25 +/- 0.23, and phi*=1.77^{+0.62}_{-0.49} * 10^{-3} Mpc^{-3}. Such a steep slope suggests that galaxies, especially the faint ones, are possibly the main sources of ionizing photons in the universe at redshift six. We also combine results from all studies at z~6 to reach an agreement in 95% confidence level that -20.45<M*<-20.05 and -1.90<alpha<-1.55. The luminosity density has been found not to evolve significantly between z~6 and z~5, but considerable evolution is detected from z~6 to z~3.

The contribution of major mergers to galaxy mass assembly along cosmic time is an important ingredient to the galaxy evolution scenario. We aim to measure the evolution of the merger rate for both luminosity/mass selected galaxy samples and investigate its dependence with the local environment. We use a sample of 10644 spectroscopically observed galaxies from the zCOSMOS redshift survey to identify pairs of galaxies destined to merge, using only pairs for which the velocity difference and projected separation of both components with a confirmed spectroscopic redshift indicate a high probability of merging. We have identified 263 spectroscopically confirmed pairs with r_p^{max} = 100 h^{-1} kpc. We find that the density of mergers depends on luminosity/mass, being higher for fainter/less massive galaxies, while the number of mergers a galaxy will experience does not depends significantly on its intrinsic luminosity but rather on its stellar mass. We find that the pair fraction and merger rate increase with local galaxy density, a property observed up to redshift z=1. We find that the dependence of the merger rate on the luminosity or mass of galaxies is already present up to redshifts z=1, and that the evolution of the volumetric merger rate of bright (massive) galaxies is relatively flat with redshift with a mean value of 3*10^{-4} (8*10^{-5} respectively) mergers h^3 Mpc^{-3} Gyr^{-1}. The dependence of the merger rate with environment indicates that dense environments favors major merger events as can be expected from the hierarchical scenario. The environment therefore has a direct impact in shapping-up the mass function and its evolution therefore plays an important role on the mass growth of galaxies along cosmic time.

Searches for very-high-redshift galaxies over the past decade have yielded a large sample of more than 6,000 galaxies existing just 900-2,000 million years (Myr) after the Big Bang (redshifts 6 > z > 3; ref. 1). The Hubble Ultra Deep Field (HUDF09) data have yielded the first reliable detections of z ~ 8 galaxies that, together with reports of a gamma-ray burst at z ~ 8.2 (refs 10, 11), constitute the earliest objects reliably reported to date. Observations of z ~ 7-8 galaxies suggest substantial star formation at z > 9-10. Here we use the full two-year HUDF09 data to conduct an ultra-deep search for z ~ 10 galaxies in the heart of the reionization epoch, only 500 Myr after the Big Bang. Not only do we find one possible z ~ 10 galaxy candidate, but we show that, regardless of source detections, the star formation rate density is much smaller (~10%) at this time than it is just ~200 Myr later at z ~ 8. This demonstrates how rapid galaxy build-up was at z ~ 10, as galaxies increased in both luminosity density and volume density from z ~ 8 to z ~ 10. The 100-200 Myr before z ~ 10 is clearly a crucial phase in the assembly of the earliest galaxies.

We present the first results on the search for very bright (M_AB -21) galaxies at redshift z~8 from the Brightest of Reionizing Galaxies (BoRG) survey. BoRG is a Hubble Space Telescope Wide Field Camera 3 pure-parallel survey that is obtaining images on random lines of sight at high Galactic latitudes in four filters (F606W, F098M, F125W, F160W), with integration times optimized to identify galaxies at z>7.5 as F098M-dropouts. We discuss here results from a search area of approximately 130 arcmin^2 over 23 BoRG fields, complemented by six other pure-parallel WFC3 fields with similar filters. This new search area is more than two times wider than previous WFC3 observations at z~8. We identify four F098M-dropout candidates with high statistical confidence (detected at greater than 8sigma confidence in F125W). These sources are among the brightest candidates currently known at z~8 and approximately ten times brighter than the z=8.56 galaxy UDFy-38135539. They thus represent ideal targets for spectroscopic followup observations and could potentially lead to a redshift record, as our color selection includes objects up to z~9. However, the expected contamination rate of our sample is about 30% higher than typical searches for dropout galaxies in legacy fields, such as the GOODS and HUDF, where deeper data and additional optical filters are available to reject contaminants.

We present UV luminosity functions (LFs) at 1500 Angstrom derived from the HST Early Release Science WFC3/UVIS data acquired over ~50 arcmin^2 of the GOODS-South field. The LFs are determined over the entire redshift range z=0.75-2.5 using two methods, similar to those used at higher redshifts for Lyman Break Galaxies (LBGs): (1) 13-band UV+optical+NIR photometric redshifts to study galaxies in the range z=0.5-2 in three bins of dz=0.5, and (2) dropout samples in three redshift windows centered at z~1.5, z~1.9, and z~2.5. The characteristic luminosity dims by ~1.5 mag from z=2.5 to z=0.75, consistent with earlier work. However, the other Schechter function parameters, the faint-end slope and the number density, are found to be remarkably constant over the range z=0.75-2.5. Using these LF determinations we find the UV luminosity density to increase by ~1.4 dex according to (1+z)^{2.58+-0.15} from z~0 to its peak at z~2.5. Strikingly, the inferred faint-end slopes for our LFs are all steeper than alpha=-1.5, in agreement with higher-redshift LBG studies. Since the faint-end slope in the local universe is found to be much flatter with alpha~=-1.2, this poses the question as to when and how the expected flattening occurs. Despite relatively large uncertainties, our data suggest alpha~=-1.7 at least down to z~1. These new results from such a shallow early dataset demonstrate very clearly the remarkable potential of WFC3/UVIS for the thorough characterization of galaxy evolution over the full redshift range z~0.5 to z~3.

We measure the spatial clustering of galaxies as a function of their morphological type at z~0.8, for the first time in a deep redshift survey with full morphological information. This is obtained by combining high-resolution HST imaging and VLT spectroscopy for about 8,500 galaxies to I_AB=22.5 with accurate spectroscopic redshifts from the zCOSMOS-Bright redshift survey. At this epoch, early-type galaxies already show a significantly stronger clustering than late-type galaxies on all probed scales. A comparison to the SDSS at z~0.1, shows that the relative clustering strength between early and late morphological classes tends to increase with cosmic time at small separations, while on large scales it shows no significant evolution since z~0.8. This suggests that most early-type galaxies had already formed in intermediate and dense environments at this epoch. Our results are consistent with a picture in which the relative clustering of different morphological types between z~1 and z~0, reflects the evolving role of environment in the morphological transformation of galaxies, on top of the global mass-driven evolution.

We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1<z<1.0 over an area of ~1.5 deg2. The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated by SED fitting of the multiwavelength optical to NIR photometry. The evolution of the galaxy stellar mass function (GSMF) in high and low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from different galaxy types, as defined by their SEDs and their morphologies, can be quantified. At redshift z~1, the GSMF is only slightly dependent on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high density regions exhibiting a marked bimodality. As a result, we infer that galaxy evolution depends on both the stellar mass and the environment, the latter setting the probability of a galaxy to have a given mass: all the galaxy properties related to the stellar mass show a dependence on environment, reflecting the difference observed in the mass functions. The shapes of the GSMFs of early- and late-type galaxies are almost identical for the extremes of the density contrast we consider. The evolution toward z=0 of the mass at which the early- and late-type GSMFs match each other is more rapid in high density environments. The comparison of the observed GSMFs to the same quantities derived from a set of mock catalogues shows that blue galaxies in sparse environments are overproduced in the semi-analytical models at intermediate and high masses, because of a deficit of star formation suppression, while at z<0.5 an excess of red galaxies is present in dense environments at intermediate and low masses, because of the overquenching of satellites. ABRIDGED

We explore the inter-relationships between mass, star-formation rate and environment in the SDSS, zCOSMOS and other surveys. The differential effects of mass and environment are completely separable to z ~ 1, indicating that two distinct processes are operating, "mass-quenching" and "environment-quenching". Environment-quenching, at fixed over-density, evidently does not change with epoch to z ~ 1, suggesting that it occurs as large-scale structure develops in the Universe. The observed constancy of the mass-function shape for star-forming galaxies, demands that the mass-quenching of galaxies around and above M*, must be proportional to their star-formation rates at all z < 2. We postulate that this simple mass-quenching law also holds over a much broader range of stellar mass and epoch. These two simple quenching processes, plus some additional quenching due to merging, then naturally produce (a) a quasi-static Schechter mass function for star-forming galaxies with a value of M* that is set by the proportionality between the star-formation and mass-quenching rates, (b) a double Schechter function for passive galaxies with two components: the dominant one is produced by mass-quenching and has exactly the same M* as the star-forming galaxies but an alpha shallower by +1, while the other is produced by environment effects and has the same M* and alpha as the star-forming galaxies, and is larger in high density environments. Subsequent merging of quenched galaxies modifies these predictions somewhat in the denser environments, slightly increasing M* and making alpha more negative. All of these detailed quantitative relationships between the Schechter parameters are indeed seen in the SDSS, lending strong support to our simple empirically-based model. The model naturally produces for passive galaxies the "anti-hierarchical" run of mean ages and alpha-element abundances with mass.

We present the analysis of the U-V rest-frame color distribution and some spectral features as a function of mass and environment for two sample of early-type galaxies up to z=1 extracted from the zCOSMOS spectroscopic survey. The first sample ("red galaxies") is defined with a photometric classification, while the second ("ETGs") by combining morphological, photometric, and spectroscopic properties to obtain a more reliable sample. We find that the color distribution of red galaxies is not strongly dependent on environment for all mass bins, with galaxies in overdense regions redder than galaxies in underdense regions with a difference of 0.027\pm0.008 mag. The dependence on mass is far more significant, with average colors of massive galaxies redder by 0.093\pm0.007 mag than low-mass galaxies throughout the entire redshift range. We study the color-mass relation, finding a mean slope 0.12\pm0.005, while the color-environment relation is flatter, with a slope always smaller than 0.04. The spectral analysis that we perform on our ETGs sample is in good agreement with our photometric results: we find for D4000 a dependence on mass between high and low-mass galaxies, and a much weaker dependence on environment (respectively a difference of of 0.11\pm0.02 and of 0.05\pm0.02); for the equivalent width of H{\delta}we measure a difference of 0.28\pm0.08 {\AA}across the same mass range and no significant dependence on environment.By analyzing the lookback time of early-type galaxies, we support the possibility of a downsizing scenario, in which massive galaxies with a stronger D4000 and an almost constant equivalent width of $H\delta$ formed their mass at higher redshift than lower mass ones. We also conclude that the main driver of galaxy evolution is the galaxy mass, the environment playing a subdominant role.

[Abridged] With the first 10000 spectra of the flux limited zCOSMOS sample (I<=22.5) we study the evolution of environmental effects on galaxy properties since z=1.0, and disentangle the dependence among galaxy colour, stellar mass and local density (3D local density contrast `delta', computed with the 5th nearest neighbour approach). We confirm that within a luminosity-limited sample (M_B<=-20.5-z) the fraction of red (U-B>=1) galaxies 'f_red' depends on delta at least up to z=1, with red galaxies residing mainly in high densities. This trend weakens for increasing z, and it is mirrored by the behaviour of the fraction of galaxies with D4000A break >=1.4. We also find that up to z=1 the fraction of galaxies with log(EW[OII]) >=1.15 is higher for lower delta, and also this dependence weakens for increasing z. Given the triple dependence among galaxy colours, stellar mass and delta, the colour-delta relation found in the luminosity-selected sample can be due to the broad range of stellar masses. Thus, we fix the stellar mass and we find that in this case the colour-delta relation is flat up to z=1 for galaxies with log(M/M_sun)>=10.7. This means that for these masses the colour-delta relation found in a luminosity-selected sample is the result of the combined colour-mass and mass-delta relations. In contrast, we find that for 0.1<=z<=0.5 and log(M/M_sun)<=10.7 'f_red' depends on delta even at fixed mass. In these mass and z ranges, environment affects directly also galaxy colours. We suggest a scenario in which the colour depends primarily on stellar mass, but for relatively low mass galaxies the local density modulates this dependence. These galaxies formed more recently, in an epoch when evolved structures were already in place, and their longer SFH allowed environment-driven physical processes to operate during longer periods of time.

The increasing abundance of passive "red-sequence" galaxies since z=1-2 is mirrored by a coincident rise in the number of galaxies with spheroidal morphologies. In this paper, however, we show that in detail the correspondence between galaxy morphology and color is not perfect, providing insight into the physical origin of this evolution. Using the COSMOS survey, we study a significant population of red sequence galaxies with disk-like morphologies. These passive disks typically have Sa-Sb morphological types with large bulges, but they are not confined to dense environments. They represent nearly one-half of all red-sequence galaxies and dominate at lower masses (log Mstar < 10) where they are increasingly disk-dominated. As a function of time, the abundance of passive disks with log Mstar < 11 increases, but not as fast as red-sequence spheroidals in the same mass range. At higher mass, the passive disk population has declined since z~1, likely because they transform into spheroidals. We estimate that as much as 60% of galaxies transitioning onto the red sequence evolve through a passive disk phase. The origin of passive disks therefore has broad implications for understanding how star formation shuts down. Because passive disks tend to be more bulge-dominated than their star-forming counterparts, a simple fading of blue disks does not fully explain their origin. We explore several more sophisticated explanations, including environmental effects, internal stabilization, and disk regrowth during gas-rich mergers. While previous work has sought to explain color and morphological transformations with a single process, these observations open the way to new insight by highlighting the fact that galaxy evolution may actually proceed through several separate stages.

We present a galaxy group catalogue spanning the redshift range 0.1 <~ z <~ 1 in the ~1.7 deg^2 COSMOS field, based on the first ~10,000 zCOSMOS spectra. The performance of both the Friends-of-Friends (FOF) and Voronoi-Delaunay-Method (VDM) approaches to group identification has been extensively explored and compared using realistic mock catalogues. We find that the performance improves substantially if groups are found by progressively optimizing the group-finding parameters for successively smaller groups, and that the highest fidelity catalogue, in terms of completeness and purity, is obtained by combining the independently created FOF and VDM catalogues. The final completeness and purity of this catalogue, both in terms of the groups and of individual members, compares favorably with recent results in the literature. The current group catalogue contains 102 groups with N >= 5 spectroscopically confirmed members, with a further ~700 groups with 2 <= N <= 4. Most of the groups can be assigned a velocity dispersion and a dark-matter mass derived from the mock catalogues, with quantifiable uncertainties. The fraction of zCOSMOS galaxies in groups is about 25% at low redshift and decreases toward ~15% at z ~ 0.8. The zCOSMOS group catalogue is broadly consistent with that expected from the semi-analytic evolution model underlying the mock catalogues. Not least, we show that the number density of groups with a given intrinsic richness increases from redshift z ~ 0.8 to the present, consistent with the hierarchical growth of structure.

We investigate how the shape of the galaxy two-point correlation function as measured in the zCOSMOS survey depends on local environment, quantified in terms of the density contrast on scales of 5 Mpc/h. We show that the flat shape previously observed at redshifts between z=0.6 and z=1 can be explained by this volume being simply 10% over-abundant in high-density environments, with respect to a Universal density probability distribution function. When galaxies corresponding to the top 10% tail of the distribution are excluded, the measured w_p(r_p) steepens and becomes indistinguishable from LCDM predictions on all scales. This is the same effect recognised by Abbas & Sheth in the SDSS data at z~0 and explained as a natural consequence of halo-environment correlations in a hierarchical scenario. Galaxies living in high-density regions trace dark matter halos with typically higher masses, which are more correlated. If the density probability distribution function of the sample is particularly rich in high-density regions because of the variance introduced by its finite size, this produces a distorted two-point correlation function. We argue that this is the dominant effect responsible for the observed "peculiar" clustering in the COSMOS field.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We present the Galaxy Stellar Mass Function (MF) up to z~1 from the zCOSMOS-bright 10k spectroscopic sample. We investigate the total MF and the contribution of ETGs and LTGs, defined by different criteria (SED, morphology or star formation). We unveil a galaxy bimodality in the global MF, better represented by 2 Schechter functions dominated by ETGs and LTGs, respectively. For the global population we confirm that low-mass galaxies number density increases later and faster than for massive galaxies. We find that the MF evolution at intermediate-low values of Mstar (logM<10.6) is mostly explained by the growth in stellar mass driven by smoothly decreasing star formation activities. The low residual evolution is consistent with ~0.16 merger per galaxy per Gyr (of which fewer than 0.1 are major). We find that ETGs increase in number density with cosmic time faster for decreasing Mstar, with a median "building redshift" increasing with mass, in contrast with hierarchical models. For LTGs we find that the number density of blue or spiral galaxies remains almost constant from z~1. Instead, the most extreme population of active star forming galaxies is rapidly decreasing in number density. We suggest a transformation from blue active spirals of intermediate mass into blue quiescent and successively (1-2 Gyr after) into red passive types. The complete morphological transformation into red spheroidals, required longer time-scales or follows after 1-2 Gyr. A continuous replacement of blue galaxies is expected by low-mass active spirals growing in stellar mass. We estimate that on average ~25% of blue galaxies is transforming into red per Gyr for logM<11. We conclude that the build-up of galaxies and ETGs follows the same downsizing trend with mass as the formation of their stars, converse to the trend predicted by current SAMs. We expect a negligible evolution of the global Galaxy Baryonic MF.

We investigate the Spitzer/IRAC properties of 36 z~7 z-dropout galaxies and 3 z~8 Y-dropout galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the SEDs to derive mean redshifts, stellar masses, and ages. The z~7 galaxies are best characterized by substantial ages (>100 Myr) and M/L_V~0.2. The main trend with decreasing luminosity is that of bluing of the far-UV slope from beta~-2.0 to beta~-3.0. This can be explained by decreasing metallicity, except for the lowest luminosity galaxies (0.1 L*{z=3}), where low metallicity and smooth SFHs fail to match the blue far-UV and moderately red H-[3.6] color. Such colors may require episodic SFHs with short periods of activity and quiescence ("on-off" cycles) and/or a contribution from emission lines. The stellar mass of our sample of z~7 star forming galaxies correlates with SFR according to log M* = 8.70 (+-0.09) + 1.06(+-0.10) log SFR, implying star formation may have commenced at z>10. No galaxies are found with SFRs much higher or lower than the past averaged SFR suggesting that the typical star formation timescales are probably a substantial fraction of the Hubble time. We report the first IRAC detection of Y-dropout galaxies at z~8. The average rest-frame U-V=0.3 (AB) of the 3 galaxies are similar to faint z~7 galaxies, implying similar M/L. The stellar mass density to M_{UV,AB}<-18 is rho*(z=8) = 1.8^{+0.7}_{-1.0} x 10^6 M_sun Mpc^-3, following log rho*(z)= 10.6(+-0.6) - 4.4(+-0.7) log(1+z) [M_sun Mpc^-3] over 3<z< 8.

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopic redshift z=1.91 and has been identified as an extremely compact galaxy by Daddi et al. 2005. We use new H-F160W imaging data obtained with Hubble Space Telescope/WFC3 to measure the deconvolved surface brightness profile to H = 28 mag arcsec**-2. We find that the surface brightness profile is well approximated by an n=3.7 Sersic profile. Our deconvolved profile is constructed by a new technique which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: r_e=0.42 +- 0.14 kpc in the observed H-F160W-band. We show that this result is robust to deviations from the Sersic model used in the fit. We test the sensitivity of our analysis to faint "wings" in the profile using simulated galaxy images consisting of a bright compact component and a faint extended component. We find that due to the combination of the WFC3 imaging depth and our method's sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, and that we can therefore confidently rule out the existence of a faint extended envelope around the observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a factor of 10 off from the local mass-size relation.

The new Wide Field Camera 3/IR observations on the Hubble Ultra-Deep Field started investigating the properties of galaxies during the reionization epoch. To interpret these observations, we present a novel approach inspired by the conditional luminosity function method. We calibrate our model to observations at z=6 and assume a non-evolving galaxy luminosity versus halo mass relation. We first compare model predictions against the luminosity function measured at z=5 and z=4. We then predict the luminosity function at z>=7 under the sole assumption of evolution in the underlying dark-matter halo mass function. Our model is consistent with the observed z>6.5 galaxy number counts in the HUDF survey and suggests a possible steepening of the faint-end slope of the luminosity function: alpha(z>8)< -1.9 compared to alpha=-1.74 at z=6. Although we currently see only the brightest galaxies, a hidden population of lower luminosity objects (L/L_{*}> 10^{-4}) might provide >75% of the total reionizing flux. Assuming escape fraction f_{esc}~0.2, clumping factor C~5, top heavy-IMF and low metallicity, galaxies below the detection limit produce complete reionization at z>8. For solar metallicity and normal stellar IMF, reionization finishes at z>6, but a smaller C/f_{esc} is required for an optical depth consistent with the WMAP measurement. Our model highlights that the star formation rate in sub-L_* galaxies has a quasi-linear relation to dark-matter halo mass, suggesting that radiative and mechanical feedback were less effective at z>6 than today.

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.

The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, appears to be represented by post-starburst galaxies. To investigate their impact on galaxy evolution, we initiated a multiwavelength study of galaxies with k+a spectral features in the COSMOS field. We examine a mass-selected sample of k+a galaxies at z=0.48-1.2 using the spectroscopic zCOSMOS sample. K+a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (<8%) of these galaxies have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~8Msun/yr). Over the entire redshift range explored, the class of k+a galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution does not vary with the strength of the Hdelta absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology relation. Although k+a galaxies are also found in underdense regions, they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~2-8Mpc, and in groups they show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical significant evolution in either the number/mass density of k+a galaxies at intermediate redshift with respect to the local values, or the spectral properties. Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of the red-sequence by up to a non-negligible level of ~10%.

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.

We use the ultra-deep WFC3/IR data over the HUDF and the Early Release Science WFC3/IR data over the CDF-South GOODS field to quantify the broadband spectral properties of candidate star-forming galaxies at z~7. We determine the UV-continuum slope beta in these galaxies, and compare the slopes with galaxies at later times to measure the evolution in beta. For luminous L*(z=3) galaxies, we measure a mean UV-continuum slope beta of -2.0+/-0.2, which is comparable to the beta~-2 derived at similar luminosities at z~5-6. However, for the lower luminosity 0.1L*(z=3) galaxies, we measure a mean beta of -3.0+/-0.2. This is substantially bluer than is found for similar luminosity galaxies at z~4, just 800 Myr later, and even at z~5-6. In principle, the observed beta of -3.0 can be matched by a very young, dust-free stellar population, but when nebular emission is included the expected beta becomes >~-2.7. To produce these very blue beta's (i.e., beta~-3), extremely low metallicities and mechanisms to reduce the red nebular emission are likely required. For example, a large escape fraction (i.e., f_{esc}>~0.3) could minimize the contribution from this red nebular emission. If this is correct and the escape fraction in faint z~7 galaxies is >~0.3, it may help to explain how galaxies reionize the universe.

We utilize the newly-acquired, ultra-deep WFC3/IR observations over the HUDF to search for star-forming galaxies at z~8-8.5, only 600 million years from recombination, using a Y_{105}-dropout selection. The new 4.7 arcmin**2 WFC3/IR observations reach to ~28.8 AB mag (5 sigma) in the Y_{105}J_{125}H_{160} bands. These remarkable data reach ~1.5 AB mag deeper than the previous data over the HUDF, and now are an excellent match to the HUDF optical ACS data. For our search criteria, we use a two-color Lyman-Break selection technique to identify z~8-8.5 Y_{105}-dropouts. We find 5 likely z~8-8.5 candidates. The sources have H_{160}-band magnitudes of ~28.3 AB mag and very blue UV-continuum slopes, with a median estimated beta of <~-2.5 (where f_{\lambda}\propto \lambda^{\beta}). This suggests that z~8 galaxies are not only essentially dust free but also may have very young ages or low metallicities. The observed number of Y_{105}-dropout candidates is smaller than the 20+/-6 sources expected assuming no evolution from z~6, but is consistent with the 5 expected extrapolating the Bouwens et al. 2008 LF results to z~8. These results provide evidence that the evolution in the LF seen from z~7 to z~3 continues to z~8. The remarkable improvement in the sensitivity of WFC3/IR has enabled HST to cross a threshold, revealing star-forming galaxies at z~8-9.

We present a sample of 16 robust z~7 z-drop galaxies detected by the newly installed WFC3/IR camera on the Hubble Space Telescope. Our analysis is based on the first epoch data of the HUDF09 program covering the Hubble Ultra Deep Field with 60 orbits of Y, J, and H observations. These remarkable data cover 4.7 arcmin^2 and are the deepest NIR images ever taken, reaching to ~29 mag AB (5 sigma). The 16 z~6.5-7.5 galaxies have been identified based on the Lyman Break technique utilizing (z-Y) vs. (Y-J) colors. They have magnitudes J = 26.0-29.0 (AB), an average apparent half-light radius of ~0.16 arcsec (<~1 kpc), and show very blue colors (some even beta<~-2.5), in particular at low luminosities. The WFC3/IR data confirms previous NICMOS detections indicating that the dropout selection at z~7 is very reliable. Our data allow a first determination of the faint end slope of the z~7 luminosity function, reaching down to M_UV ~ -18, a full magnitude fainter than previous measurements. When fixing phi*=1.4e-3/Mpc^3/mag to the value previously measured at z~6, we find a best-fit value of alpha=-1.77+-0.20, with a characteristic luminosity of M_*=-19.91+-0.09. This steep slope is similar to what is seen at z~2-6 and indicates that low luminosity galaxies could potentially provide adequate flux to reionize the universe. The remarkable depth and resolution of these new images provide insights into the coming power of JWST.

We study the Spitzer Infrared Array Camera (IRAC) mid-infrared (rest-frame optical) fluxes of 14 newly WFC3/IR-detected z=7 z_{850}-dropout galaxies and 5 z=8 Y_{105}-dropout galaxies. The WFC3/IR depth and spatial resolution allow accurate removal of contaminating foreground light, enabling reliable flux measurements at 3.6 micron and 4.5 micron. None of the galaxies are detected to [3.6]=26.9 (AB, 2 sigma), but a stacking analysis reveals a robust detection for the z_{850}-dropouts and an upper limit for the Y_{105}-dropouts. We construct average broadband SEDs using the stacked ACS, WFC3, and IRAC fluxes and fit stellar population synthesis models to derive mean redshifts, stellar masses, and ages. For the z_{850}-dropouts, we find z=6.9^{+0.1}_{-0.1}, (U-V)_{rest}=0.4, reddening A_V=0, stellar mass M*=1.2^{+0.3}_{-0.6} x 10^9 M_sun (Salpeter IMF). The best-fit ages ~300Myr, M/L_V=0.2, and SSFR=1.7Gyr^{-1} are similar to values reported for luminous z=7 galaxies, indicating the galaxies are smaller but not younger. The sub-L* galaxies observed here contribute significantly to the stellar mass density and under favorable conditions may have provided enough photons for sustained reionization at 7<z<11. In contrast, the z=8.3^{+0.1}_{-0.2} Y_{105}-dropouts have stellar masses that are uncertain by 1.5 dex due to the near-complete reliance on far-UV data. Adopting the 2 sigma upper limit on the M/L(z=8), the stellar mass density to M_{UV,AB} < -18 declines from rho*(z=7)=3.7^{+1.0}_{-1.8} x 10^6 M_sun Mpc^{-3} to rho*(z=8) < 8 x 10^5 M_sun Mpc^{-3}, following (1+z)^{-6} over 3<z<8. Lower masses at z=8 would signify more dramatic evolution, which can be established with deeper IRAC observations, long before the arrival of the James Webb Space Telescope.

We present a first morphological study of z~7-8 Lyman Break galaxies (LBGs) from Oesch et al. 2009 and Bouwens et al. 2009 detected in ultra-deep near-infrared imaging of the Hubble Ultra Deep field (HUDF) by the HUDF09 program. With an average intrinsic size of 0.7+-0.3 kpc these galaxies are found to be extremely compact having an average observed surface brightness of mu_J ~= 26 mag arcsec^(-2), and only two out of the full sample of 16 z~7 galaxies show extended features with resolved double cores. By comparison to lower redshift LBGs it is found that only little size evolution takes place from z~7 to z~6, while galaxies between z~4-5 show more extended wings in their apparent profiles. The average size scales as (1+z)^(-m) with m=1.12+-0.17 for galaxies with luminosities in the range (0.3-1)L*_{z=3} and with m=1.32+-0.52 for (0.12-0.3)L*_{z=3}, consistent with galaxies having constant comoving sizes. The peak of the size distribution changes only slowly from z~7 to z~4. However, a tail of larger galaxies (>~ 1.2 kpc) is gradually built up towards later cosmic times, possibly via hierarchical build-up or via enhanced accretion of cold gas. Additionally, the average star-formation surface density of LBGs with luminosities (0.3-1)L*_{z=3} is nearly constant at Sigma_{SFR}=1.9 Msun/yr/kpc^2 over the entire redshift range z~4-7 suggesting similar star-formation efficiencies at these early epochs. The above evolutionary trends seem to hold out to z~8 though the sample is still small and possibly incomplete.

We present a catalog of 213 type-2 AGN selected from the zCOSMOS survey. The selected sample covers a wide redshift range (0.15<z<0.92) and is deeper than any other previous study, encompassing the luminosity range 10^{5.5} < Lsun< L[OIII] < 10^{9.1} Lsun. We explore the intrinsic properties of these AGN and the relation to their X-ray emission (derived from the XMM-COSMOS observations). We study their evolution by computing the [OIII]5007A line luminosity function (LF) and we constrain the fraction of obscured AGN as a function of luminosity and redshift. The sample was selected on the basis of the optical emission line ratios, after applying a cut to the signal-to-noise ratio (S/N) of the relevant lines. We used the standard diagnostic diagrams [OIII]/Hbeta versus [NII]/Halpha and ([OIII]/Hbeta versus [SII]/Halpha) to isolate AGN in the redshift range 0.15<z<0.45 and the diagnostic diagram [OIII]/Hbeta versus [OII]/Hbeta to extend the selection to higher redshift (0.5<z<0.92). Combining our sample with one drawn from SDSS, we found that the best description of the evolution of type-2 AGN is a luminosity-dependent density evolution model. Moreover, using the type-1 AGN LF we were able to constrain the fraction of type-2 AGN to the total (type-1 + type-2) AGN population. We found that the type-2 fraction decreases with luminosity, in agreement with the most recent results, and shows signs of a slight increase with redshift. However, the trend with luminosity is visible only after combining the SDSS+zCOSMOS samples. From the COSMOS data points alone, the type-2 fraction seems to be quite constant with luminosity.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $\mu$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24-micron-selected galaxies with 0.06< S(24um)< 0.50 mJy and I(AB)<22.5, over 1.5 deg^2 of the COSMOS field, and studied different spectral properties, depending on redshift. At 0.2<z<0.3, we found that different reddening laws of common use in the literature explain the dust extinction properties of around 80% of our infrared (IR) sources, within the error bars. For up to 16% of objects, instead, the Halpha/Hbeta ratios are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust distributions. In only a few of our galaxies at 0.2<z<0.3 the IR emission could be mainly produced by dust heated by old rather than young stars. Besides, the line ratios of ~22% of our galaxies suggest that they might be star-formation/nuclear-activity composite systems. At 0.5<z<0.7, we estimated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined analysis of the Hdelta equivalent-width versus Dn(4000) diagram for 1722 faint and bright 24um galaxies at 0.6<z<1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, quiescent, exponentially-declining star formation histories can well reproduce the spectral properties of ~40% of the less luminous sources. Our results suggest a transition in the possible modes of star formation at total IR luminosities L(TIR)=(3 +/-2)x10^11 Lsun.

We use the overdensity field reconstructed in the volume of the COSMOS area to study the nonlinear biasing of the zCOSMOS galaxies. The galaxy overdensity field is reconstructed using the current sample of ~8500 accurate zCOSMOS redshifts at I(AB)<22.5 out to z~1 on scales R from 8 to 12 Mpc/h. By comparing the probability distribution function (PDF) of galaxy density contrast delta_g to the lognormal approximation of the PDF of the mass density contrast delta, we obtain the mean biasing function b(delta,z,R) between the galaxy and matter overdensity field and its second moments b(hat) and b(tilde) up to z~1. Over the redshift interval 0.4<z<1 the conditional mean function <delta_g|delta> = b(delta,z,R) delta is of the following characteristic shape. The function vanishes in the most underdense regions and then sharply rises in a nonlinear way towards the mean densities. <delta_g|delta> is almost a linear tracer of the matter in the overdense regions, up to the most overdense regions in which it is nonlinear again and the local effective slope of <delta_g|delta> vs. delta is smaller than unity. The <delta_g|delta> function is evolving only slightly over the redshift interval 0.4<z<1. The linear biasing parameter increases from b(hat)=1.24+/-0.11 at z=0.4 to b(hat)=1.64+/-0.15 at z=1 for the M_B<-20-z sample of galaxies. b(hat) does not show any dependence on the smoothing scale from 8 to 12 Mpc/h, but increases with luminosity. The measured nonlinearity parameter b(tilde)/b(hat) is of the order of a few percent (but it can be consistent with 0) and it does not change with redshift, the smoothing scale or the luminosity. By matching the linear bias of galaxies to the halo bias, we infer that the M_B<-20-z galaxies reside in dark matter haloes with a characteristic mass of about 3-6 x 10^12 Msol, depending on the halo bias fit.

(Abridged) We studied the evolution in the B band luminosity function to z~1 in the zCOSMOS 10k sample, for which both accurate galaxy classifications and a detailed description of the local density field are available. The global LF exhibits a brightening of ~0.7 mag in M* from z~0.2 to z~0.9. At low z, late types dominate at faint magnitudes, while the bright end is populated mainly by early types. At higher z, late-type galaxies evolve significantly and, at z~1, the contribution from the various types to the bright end of the LF is comparable. The evolution for early types is in both luminosity and normalization. A similar behaviour is exhibited by late types, but with an opposite trend for the normalization. Studying the role of the environment, we find that the global LF of galaxies in overdense regions has always a brighter M* and a flatter slope. In low density environments, the main contribution to the LF is from blue galaxies, while for high density environments there is an important contribution from red galaxies to the bright end. The differences between the global LF in the two environments are not due to only a difference in the relative numbers of red and blue galaxies, but also to their relative luminosity distributions: the value of M* for both types in underdense regions is always fainter than in overdense environments. The "specular" evolution of late- and early-type galaxies is consistent with a scenario where a part of blue galaxies is transformed in red galaxies with increasing cosmic time, without significant changes in the fraction of intermediate-type galaxies. The bulk of this tranformation in overdense regions probably happened before z~1, while it is still ongoing at lower z in underdense environments.

We study the evolution of galaxies inside and outside of the group environment since z=1 using a large well defined set of groups and galaxies from the zCOSMOS-bright redshift survey in the COSMOS field. The fraction of galaxies with early-type morphologies increases monotonically with M_B luminosity and stellar mass and with cosmic epoch. It is higher in the groups than elsewhere, especially at later epochs. The emerging environmental effect is superposed on a strong global mass-driven evolution, and at z~0.5 and log(M*/Msol)~10.2, the "effect" of group environment is equivalent to (only) about 0.2 dex in stellar mass or 2 Gyr in time. The stellar mass function of galaxies in groups is enriched in massive galaxies. We directly determine the transformation rates from late to early morphologies, and for transformations involving colour and star formation indicators. The transformation rates are systematically about twice as high in the groups as outside, or up to 3-4 times higher correcting for infall and the appearance of new groups. The rates reach values, for masses around the crossing mass 10^10.5 Msol, as high as (0.3-0.7)/Gyr in the groups, implying transformation timescales of 1.4-3 Gyr, compared with less than 0.2/Gyr, i.e. timescales >5 Gyr, outside of groups. All three transformation rates decrease at higher stellar masses, and must decrease also at the lower masses below 10^10 Msol which we cannot well probe. The rates involving colour and star formation are consistently higher than those for morphology, by a factor of about 50%. Our conclusion is that the transformations which drive the evolution of the overall galaxy population since z~1 must occur at a rate 2-4 times higher in groups than outside of them.

We took advantage of the wealth of information provided by the first ~10000 galaxies of the zCOSMOS-bright survey and its group catalogue to study the complex interplay between group environment and galaxy properties. The classical indicator F_blue (fraction of blue galaxies) proved to be a simple but powerful diagnostic tool. We studied its variation for different luminosity and mass selected galaxy samples. Using rest-frame B-band selected samples, the groups galaxy population exhibits significant blueing as redshift increases, but maintains a lower F_blue with respect both to the global and the isolated galaxy population. However moving to mass selected samples it becomes apparent that such differences are largely due to the biased view imposed by the B-band luminosity selection, being driven by the population of lower mass, bright blue galaxies for which we miss the redder, equally low mass, counterparts. By focusing the analysis on narrow mass bins such that mass segregation becomes negligible we find that only for the lowest mass bin explored (logMass <= 10.6) does a significant residual difference in color remain as a function of environment, while this difference becomes negligible toward higher masses. Our results indicate that red galaxies of logMass >= 10.8 are already in place at z ~ 1 and do not exhibit any strong environmental dependence, possibly originating from so-called 'nature'/internal mechanisms. In contrast, for lower galaxy masses and redshifts lower than z ~ 1, we observe the emergence in groups of a population of 'nurture' red galaxies: slightly deviating from the trend of the downsizing scenario followed by the global galaxy population, and more so with cosmic time. These galaxies exhibit signatures of group-related secular physical mechanisms directly influencing galaxy evolution.

We present in these proceedings some preliminary results we have obtained studying the evolution of the specific star formation rate as a function of surface mass density and Sersic indices at z<0.7. These results are based on the consistent comparison of the properties of ~ 650 massive zCOSMOS galaxies in a mass-complete sample at 0.5<z<0.7 with a mass-complete sample of ~ 21500 SDSS local galaxies.

One of the key unanswered questions in the study of galaxy evolution is what physical processes inside galaxies drive the changes in the SFRs in individual galaxies that, taken together, produce the large decline in the global star-formation rate density (SFRD) to redshifts since z~2. Many studies of the SFR at intermediate redshifts have been made as a function of the integrated stellar mass of galaxies but these did not use information on the internal structural properties of the galaxies. In this paper we present a comparative study of the dependence of SFRs on the average surface mass densities (SigmaM) of galaxies of different morphological types up to z~1 using the zCOSMOS and SDSS surveys. The main findings about the evolution of these relatively massive galaxies are: 1) There is evidence that, for both SDSS ans zCOSMOS galaxies, the mean specific SFR within a given population (either disk-dominated or bulge-dominated) is independent of SigmaM; 2) The observed SSFR - SigmaM step-function relation is due, at all investigated redshifts, to the changing mix of disk-dominated and bulge-dominated galaxies as surface density increases and the strong difference in the average SSFR between disks and bulges. We also find a modest differential evolution in the size-mass relations of disk and spheroid galaxies; 3) The shape of the median SSFR - SigmaM relation is similar, but with median SSFR values that are about 5-6 times higher in zCOSMOS galaxies than for SDSS, across the whole range of SigmaM, and in both spheroid and disk galaxies. This increase matches that of the global SFRD of the Universe as a whole, emphasizing that galaxies of all types are contributing, proportionally, to the global increase in SFRD in the Universe back to these redshifts (abridged).

For more than two decades we have known that galaxy morphological segregation is present in the Local Universe. It is important to see how this relation evolves with cosmic time. To investigate how galaxy assembly took place with cosmic time, we explore the evolution of the morphology-density relation up to redshift z~1 using about 10000 galaxies drawn from the zCOSMOS Galaxy Redshift Survey. Taking advantage of accurate HST/ACS morphologies from the COSMOS survey, of the well-characterised zCOSMOS 3D environment, and of a large sample of galaxies with spectroscopic redshift, we want to study here the evolution of the morphology-density relation up to z~1 and its dependence on galaxy luminosity and stellar mass. The multi-wavelength coverage of the field also allows a first study of the galaxy morphological segregation dependence on colour. We further attempt to disentangle between processes that occurred early in the history of the Universe or late in the life of galaxies. The zCOSMOS field benefits of high-resolution imaging in the F814W filter from the Advanced Camera for Survey (ACS). We use standard morphology classifiers, optimised for being robust against band-shifting and surface brightness dimming, and a new, objective, and automated method to convert morphological parameters into early, spiral, and irregular types. We use about 10000 galaxies down to I_AB=22.5 with a spectroscopic sampling rate of 33% to characterise the environment of galaxies up to z~1 from the 100 kpc scales of galaxy groups up to the 100 Mpc scales of the cosmic web. ABRIDGED

We study the dependence of galaxy clustering on luminosity and stellar mass at redshifts z ~ [0.2-1] using the first zCOSMOS 10K sample. We measure the redshift-space correlation functions xi(rp,pi) and its projection wp(rp) for sub-samples covering different luminosity, mass and redshift ranges. We quantify in detail the observational selection biases and we check our covariance and error estimate techniques using ensembles of semi-analytic mock catalogues. We finally compare our measurements to the cosmological model predictions from the mock surveys. At odds with other measurements, we find a weak dependence of galaxy clustering on luminosity in all redshift bins explored. A mild dependence on stellar mass is instead observed. At z~0.7, wp(rp) shows strong excess power on large scales. We interpret this as produced by large-scale structure dominating the survey volume and extending preferentially in direction perpendicular to the line-of-sight. We do not see any significant evolution with redshift of the amplitude of clustering for bright and/or massive galaxies. The clustering measured in the zCOSMOS data at 0.5<z<1 for galaxies with log(M/M_\odot)>=10 is only marginally consistent with predictions from the mock surveys. On scales larger than ~2 h^-1 Mpc, the observed clustering amplitude is compatible only with ~1% of the mocks. Thus, if the power spectrum of matter is LCDM with standard normalization and the bias has no unnatural scale-dependence, this result indicates that COSMOS has picked up a particularly rare, ~2-3 sigma positive fluctuation in a volume of ~10^6 h^-1 Mpc^3. These findings underline the need for larger surveys of the z~1 Universe to appropriately characterize the level of structure at this epoch.

We use the current sample of ~10,000 zCOSMOS spectra of sources selected with I(AB) < 22.5 to define the density field out to z~1, with much greater resolution in the radial dimension than has been possible with either photometric redshifts or weak lensing. We apply new algorithms that we have developed (ZADE) to incorporate objects not yet observed spectroscopically by modifying their photometric redshift probability distributions using the spectroscopic redshifts of nearby galaxies. This strategy allows us to probe a broader range of galaxy environments and reduce the Poisson noise in the density field. The reconstructed overdensity field of the 10k zCOSMOS galaxies consists of cluster-like patterns surrounded by void-like regions, extending up to z~1. Some of these structures are very large, spanning the ~50 Mpc/h transverse direction of the COSMOS field and extending up to Delta z~0.05 in redshift. We present the three dimensional overdensity maps and compare the reconstructed overdensity field to the independently identified virialised groups of galaxies and clusters detected in the visible and in X-rays. The distribution of the overdense structures is in general well traced by these virialised structures. A comparison of the large scale structures in the zCOSMOS data and in the mock catalogues reveals an excellent agreement between the fractions of the volume enclosed in structures of all sizes above a given overdensity between the data and the mocks in 0.2<z<1.

The impact of environment on AGN activity up to z~1 is assessed by utilizing a mass-selected sample of galaxies from the 10k catalog of the zCOSMOS spectroscopic redshift survey. We identify 147 AGN by their X-ray emission as detected by XMM-Newton from a parent sample of 7234 galaxies. We measure the fraction of galaxies with stellar mass M_*>2.5x10^10 Msun that host an AGN as a function of local overdensity using the 5th, 10th and 20th nearest neighbors that cover a range of physical scales (~1-4 Mpc). Overall, we find that AGNs prefer to reside in environments equivalent to massive galaxies with substantial levels of star formation. Specifically, AGNs with host masses between 0.25-1x10^11 Msun span the full range of environments (i.e., field-to-group) exhibited by galaxies of the same mass and rest-frame color or specific star formation rate. Host galaxies having M_*>10^11 Msun clearly illustrate the association with star formation since they are predominantly bluer than the underlying galaxy population and exhibit a preference for lower density regions analogous to SDSS studies of narrow-line AGN. To probe the environment on smaller physical scales, we determine the fraction of galaxies (M_*>2.5x10^10 Msun) hosting AGNs inside optically-selected groups, and find no significant difference with field galaxies. We interpret our results as evidence that AGN activity requires a sufficient fuel supply; the probability of a massive galaxy to have retained some sufficient amount of gas, as evidence by its ongoing star formation, is higher in underdense regions where disruptive processes (i.e., galaxy harrassment, tidal stripping) are lessened.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if accretion onto supermassive black holes and star formation are explicitly linked up to z~1. We identify 152 galaxies that harbor AGN, based on XMM observations of 7543 galaxies (i<22.5). Star formation rates (SFRs), including those weighted by stellar mass, are determined using the [OII]3727 emission-line, corrected for an AGN contribution. We find that the majority of AGN hosts have significant levels of star formation with a distribution spanning ~1-100 Msun yr^-1. The close association between AGN activity and star formation is further substantiated by an increase in the AGN fraction with the youthfulness of their stars as indicated by the rest-frame color (U-V) and spectral index Dn(4000); we demonstrate that mass-selection alleviates an artifical peak falling in the transition region when using luminosity-limited samples. We also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the overall galaxy population and explains the low SFRs in AGNs (z<0.3) from the SDSS. We conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR is shown to be weak although the average ratio is constant with redshift, effectively shifting the evidence for a co-evolution scenario in a statistical manner to smaller physical scales. Our findings illustrate an intermittent scenario with an AGN lifetime substantially shorter than that of star formation and underlying complexities regarding fueling over vastly different physical scales yet to be determined [Abridged].

Aims. We investigate the relationships between three main optical galaxy observables (spectral properties, colours, and morphology), exploiting the data set provided by the COSMOS/zCOSMOS survey. The purpose of this paper is to define a simple galaxy classification cube, using a carefully selected sample of around 1000 galaxies. Methods. Using medium resolution spectra of the first 1k zCOSMOS-bright sample, optical photometry from the Subaru/COSMOS observations, and morphological measurements derived from ACS imaging, we analyze the properties of the galaxy population out to z~1. Applying three straightforward classification schemes (spectral, photometric, and morphological), we identify two main galaxy types, which appear to be linked to the bimodality of galaxy population. The three parametric classifications constitute the axes of a "classification cube". Results. A very good agreement exists between the classification from spectral data (quiescent/star-forming galaxies) and that based on colours (red/blue galaxies). The third parameter (morphology) is less well correlated with the first two: in fact a good correlation between the spectral classification and that based on morphological analysis (early-/late-type galaxies) is achieved only after partially complementing the morphological classification with additional colour information. Finally, analyzing the 3D-distribution of all galaxies in the sample, we find that about 85% of the galaxies show a fully concordant classification, being either quiescent, red, bulge-dominated galaxies (~20%) or star-forming, blue, disk-dominated galaxies (~65%). These results imply that the galaxy bimodality is a consistent behaviour both in morphology, colour and dominant stellar population, at least out to z~1.

[Abridged] We detect three (plus one less certain) z-dropout sources in two separate fields of our UDF05 HST NICMOS images. These z~7 Lyman-Break Galaxy (LBG) candidates allow us to constrain the Luminosity Function (LF) of the star forming galaxy population at those epochs. By assuming a change in only M* and adopting a linear evolution in redshift, anchored to the measured values at z~6, the best fit evolution coefficient is found to be 0.43+-0.19 mag per unit redshift (0.36+-0.18, if including all four candidates), which provides a value of M*(z=7.2)=-19.7+-0.3. This implies a steady evolution for the LBG LF out to z~7, at the same rate that is observed throughout the z~3 to 6 period. This puts a strong constraint on the star-formation histories of z~6 galaxies, whose ensemble star-formation rate density must be lower by a factor 2 at ~170 Myr before the epoch at which they are observed. In particular, a large fraction of stars in the z~6 LBG population must form at redshifts well above z~7. Extrapolating this steady evolution of the LF out to higher redshifts, we estimate that galaxies would be able to reionize the universe by z~6, provided that the faint-end slope of the z>7 LF steepens to alpha~-1.9, and that faint galaxies, with luminosities below the current detection limits, contribute a substantial fraction of the required ionizing photons. This scenario gives however an integrated optical depth to electron scattering that is ~2sigma below the WMAP-5 measurement. Therefore, altogether, our results indicate that, should galaxies be the primary contributors to reionization, either the currently detected evolution of the galaxy population slows down at z>7, or the LF evolution must be compensated by a decrease in metallicity and a corresponding increase in ionization efficiency at these early epochs.

We investigate the close environment of 203 Spitzer 24 micron-selected sources at 0.6<z<1.0 using zCOSMOS-bright redshifts and spectra of I<22.5 AB mag galaxies, over 1.5 sq. deg. of the COSMOS field. We quantify the degree of passivity of the LIRG and ULIRG environments by analysing the fraction of close neighbours with Dn(4000)>1.4. We find that LIRGs at 0.6<z<0.8 live in more passive environments than those of other optical galaxies that have the same stellar mass distribution. Instead, ULIRGs inhabit more active regions (e.g. LIRGs and ULIRGs at 0.6<z<0.8 have, respectively, (42.0 +/- 4.9)% and (24.5 +/- 5.9)% of neighbours with Dn (4000)>1.4 within 1 Mpc and +/- 500 km/s). The contrast between the activities of the close environments of LIRGs and ULIRGs appears especially enhanced in the COSMOS field density peak at z~0.67, because LIRGs on this peak have a larger fraction of passive neighbours, while ULIRGs have as active close environments as those outside the large-scale structure. The differential environmental activity is related to the differences in the distributions of stellar mass ratios between LIRGs/ULIRGs and their close neighbours, as well as in the general local density fields. At 0.8<z<1.0, instead, we find no differences in the environment densities of ULIRGs and other similarly massive galaxies, in spite of the differential activities. We discuss a possible scenario to explain these findings.

ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer, is a tool for estimating redshifts and template types of galaxies using medium- and broad-band photometric data. ZEBRA employs novel techniques within the template-fitting approach to produce high-quality Maximum-Likelihood and Bayesian redshift estimates. This manuscript serves as a user guide to ZEBRA. It explains how to use ZEBRA, specifies input and output formats, and gives a short account of the available options. ZEBRA is a free and open-source software distributed under the GNU Public License 3 and available at http://www.exp-astro.phys.ethz.ch/ZEBRA . Several upgrades are currently being implemented in ZEBRA. An updated documentation will be provided at each new release. Any problems, comments and suggestions on the code and the manual should be sent via e-mail to zebra@phys.ethz.ch.

Accurate photometric redshifts are among the key requirements for precision weak lensing measurements. Both the large size of the Sloan Digital Sky Survey (SDSS) and the existence of large spectroscopic redshift samples that are flux-limited beyond its depth have made it the optimal data source for developing methods to properly calibrate photometric redshifts for lensing. Here, we focus on galaxy-galaxy lensing in a survey with spectroscopic lens redshifts, as in the SDSS. We develop statistics that quantify the effect of source redshift errors on the lensing calibration and on the weighting scheme, and show how they can be used in the presence of redshift failure and sampling variance. We then demonstrate their use with 2838 source galaxies with spectroscopy from DEEP2 and zCOSMOS, evaluating several public photometric redshift algorithms, in two cases including a full p(z) for each object, and find lensing calibration biases as low as 1% (due to fortuitous cancellation of two types of bias) or as high as 20% for methods in active use (despite the small mean photoz bias of these algorithms). Our work demonstrates that lensing-specific statistics must be used to reliably calibrate the lensing signal, due to asymmetric effects of (frequently non-Gaussian) photoz errors. We also demonstrate that large-scale structure (LSS) can strongly impact the photoz calibration and its error estimation, due to a correlation between the LSS and the photoz errors, and argue that at least two independent degree-scale spectroscopic samples are needed to suppress its effects. Given the size of our spectroscopic sample, we can reduce the galaxy-galaxy lensing calibration error well below current SDSS statistical errors.

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP) and NICMOS (F110W, F160W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z~5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V-i, i-z selection criterion, we identify a sample of 101 and 133 candidate z~5 galaxies down to z_{850}=28.5 and 29.25 magnitudes in the NICP12 and in the HUDF fields, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 A LBG LF over the range M_{1400}=[-21.4, -17.1], i.e. down to ~0.04 L* at z~5, and complement it with data from the Subaru Deep Field (SDF) from Yoshida et al. (2006) to extend it to the brighter end (M_{1400}>-22.2). We show that: (i) Different assumptions regarding the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z~5; (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z~5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha ~ -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z~3 and z~6.