We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (10^9.5 Msun for star-forming galaxies). At this early epoch, the fraction of quiescent galaxies is already highly elevated in the clusters compared to the field at the same redshift. The quenched fraction excess (QFE) represents the fraction of galaxies that would be star-forming in the field, but are quenched due to their environment. The QFE is strongly mass dependent, and increases from ~30% at Mstar=10^9.7 Msun, to ~80% at Mstar=10^11.0 Msun. Nonetheless, the shapes of the SMFs of the two individual galaxy types, star-forming and quiescent galaxies, are identical between the clusters and the field - to high statistical precision. Yet, along with the different quiescent fractions is the total galaxy SMF environmentally dependent, with a relative deficit of low-mass galaxies in the clusters. These results are in stark contrast with findings in the local Universe, and thus require a substantially different quenching mode to operate at early times. We discuss these results in the light of several popular quenching models.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of $3.1^{+0.1}_{-0.2} \times 10^{11} ~\rm{M}_\odot$ at $z = 3.493$, this galaxy is nearly three times more massive than the highest redshift spectroscopically confirmed absorption-line identified galaxy known. The star-formation history of this quiescent galaxy implies that it formed $>1000 ~\rm{M}_\odot$/yr for almost 0.5 Gyr beginning at $z\sim7.2$, strongly suggestive that it is the descendant of massive dusty star-forming galaxies at $5<z<7$ recently observed with ALMA. While galaxies with similarly extreme stellar masses are reproduced in some simulations at early times, such a lack of ongoing star formation is not seen there. This suggests the need for a more rapid quenching process than is currently prescribed, challenging our current understanding of how ultra-massive galaxies form and evolve in the early Universe.

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

We investigate the stellar and dust properties of massive (log$(M_*/M_\odot) \ge 10.5$) and dusty ($A_V \ge 1$) galaxies at $1 \le z \le 4$ by modeling their spectral energy distributions (SEDs) obtained from the combination of UltraVISTA DR3 photometry and \textit{Herschel} PACS-SPIRE data using MAGPHYS. Although the rest-frame U-V vs V-J (UVJ) diagram traces well the star-formation rates (SFR) and dust obscuration (A$_V$) out to $z \sim 3$, $\sim$15-20\% of the sample surprisingly resides in the quiescent region of the UVJ diagram, while $\sim50$\% at $3<z<4$ fall in the unobscured star-forming region. The median SED of massive dusty galaxies exhibits weaker MIR and UV emission, and redder UV slopes with increasing cosmic time. The IR emission for our sample has a significant contribution ($>20\%$) from dust heated by evolved stellar populations rather than star formation, demonstrating the need for panchromatic SED modeling. The local relation between dust mass and SFR is followed only by a sub-sample with cooler dust temperatures, while warmer objects have reduced dust masses at a given SFR. Most star-forming galaxies in our sample do not follow local IRX-$\beta$ relations, though IRX does strongly correlate with A$_V$. Our sample follows local relations, albeit with large scatter, between ISM diagnostics and sSFR. We show that FIR-detected sources represent the extreme of a continuous population of dusty galaxies rather than a fundamentally different population. Finally, using commonly adopted relations to derive SFRs from the combination of the rest-frame UV and the observed 24$\mu$m is found to overestimate the SFR by a factor of 3-5 for the galaxies in our sample.

We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 \mu$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ \alpha \sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15.

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

Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $\Gamma \approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation.

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.

The IRAC mapping of the NMBS-II fields program is an imaging survey at 3.6 and 4.5$\mu$m with the Spitzer Infrared Array Camera (IRAC). The observations cover three Canada-France-Hawaii Telescope Legacy Survey Deep (CFHTLS-D) fields, including one also imaged by AEGIS, and two MUSYC fields. These are then combined with archival data from all previous programs into deep mosaics. The resulting imaging covers a combined area of about 3 $deg^2$, with at least $\sim$2 hr integration time for each field. In this work, we present our data reduction techniques and document the resulting coverage maps at 3.6 and 4.5$\mu$m. All of the images are W-registered to the reference image, which is either the z-band stack image of the 25\% best seeing images from the CFHTLS-D for CFHTLS-D1, CFHTLS-D3, and CFHTLS-D4, or the K-band images obtained at the Blanco 4-m telescope at CTIO for MUSYC1030 and MUSYC1255. We make all images and coverage maps described herein publicly available via the Spitzer Science Center.

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.

Distant luminous Lyman-alpha emitters (LAEs) are excellent targets for spectroscopic observations of galaxies in the epoch of reionisation (EoR). We present deep high-resolution (R=5000) VLT/X-SHOOTER observations, along with an extensive collection of photometric data of `COLA1', a proposed double peaked LAE at z=6.6 (Hu et al. 2016). We rule out that COLA1's emission line is an [OII] doublet at z=1.475 on the basis of i) the asymmetric red line-profile and flux ratio of the peaks (blue/red=$0.31\pm0.03$) and ii) an unphysical [OII]/Ha ratio ([OII]/Ha > 22). We show that COLA1's observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lya and [OIII] at z=2.142. We thus conclude that COLA1 is a real double-peaked LAE at z=6.593, the first discovered at z>6, confirming the result from Hu et al. (2016). COLA1 is UV luminous (M$_{1500}=-21.6\pm0.3$), has a high equivalent width (EW$_{0}$~120 \AA) and very compact Lya emission (r$_{50} = 0.3$ kpc). Relatively weak inferred Hb+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z<1/20 Z$_{\odot}$ or reduced strength of nebular lines due to high escape of ionising photons. The small Lya peak separation of $220\pm20$ km/s implies a low HI column density and an ionising photon escape fraction of ~15-30 %, providing the first direct evidence that such galaxies contribute actively to the reionisation of the Universe at z>6. Based on simple estimates, we find that COLA1 could have provided just enough photons to reionise its own ~0.3 pMpc bubble, allowing the blue Lya line to be observed. However, we also discuss alternative scenarios explaining the detected double peaked nature of COLA1. Our results show that future high-resolution observations of statistical samples of double peaked LAEs at z>5 are a promising probe of the occurrence of ionised regions around galaxies in the EoR.

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 targeted follow-up Hubble Space Telescope WFC3 F160W imaging study of very massive galaxies $(\log(M_{\rm{star}}/M_{\odot})> 11.2)$ selected from a combination of ground-based near-infrared galaxy surveys (UltraVISTA, NMBS-II, UKIDSS UDS) at $1.5<z<3$. We find that these galaxies are diverse in their structures, with $\sim1/3$ of the targets being composed of close pairs, and span a wide range in sizes. At $1.5<z<2.5$, the sizes of both star-forming and quiescent galaxies are consistent with the extrapolation of the stellar mass-size relations determined at lower stellar masses. At $2.5<z<3.0$, however, we find evidence that quiescent galaxies are systematically larger than expected based on the extrapolation of the relation derived using lower stellar mass galaxies. We used the observed light profiles of the blended systems to decompose their stellar masses and investigate the effect of the close pairs on the measured number densities of very massive galaxies in the early universe. We estimate correction factors to account for close-pair blends and apply them to the observed stellar mass functions measured using ground-based surveys. Given the large uncertainties associated with this extreme population of galaxies, there is currently little tension between the (blending-corrected) number density estimates and predictions from theoretical models. Although we currently lack the statistics to robustly correct for close-pair blends, we show that this is a systematic effect which can reduce the observed number density of very massive galaxies by up to a factor of $\sim1.5$, and should be accounted for in future studies of stellar mass functions.

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.

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 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 rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history.

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 examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

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 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 present the analysis and results of a spectroscopic follow-up program of a mass-selected sample of six galaxies at 3 < z < 4 using data from Keck-NIRSPEC and VLT-Xshooter. We confirm the z > 3 redshifts for half of the sample through the detection of strong nebular emission lines, and improve the zphot accuracy for the remainder of the sample through the combination of photometry and spectra. The modeling of the emission-line-corrected spectral energy distributions (SEDs) adopting improved redshifts confirms the very large stellar masses of the sample (M_* ~ 1.5-4 x 10^11 Msun) in the first 2 Gyrs of cosmic history, with a diverse range in stellar ages, star formation rates and dust content. From the analysis of emission line luminosities and widths, and far-infrared (FIR) fluxes we confirm that >80% of the sample are hosts to luminous hidden active galactic nuclei (AGNs), with bolometric luminosities of ~10^(44-46) erg/s. We find that the MIPS 24um photometry is largely contaminated by AGN continuum, rendering the SFRs derived using only 24um photometry to be severely overestimated. By including the emission from the AGN in the modeling of the UV-to-FIR SEDs, we confirm that the presence of the AGN does not bias considerably the stellar masses (< 0.3 dex at 1\sigma). We show evidence for a rapid increase of the AGN fraction from ~30% to ~60-100% over the 1 Gyr between z~2 and z~3. Although we cannot exclude some enhancement of the AGN fraction for our sample due to selection effects, the small measured [OIII] contamination to the observed K-band fluxes suggests that our sample is not significantly biased toward massive galaxies hosting AGNs.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

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.

In this paper, we use stacking analysis to trace the mass-growth, colour evolution, and structural evolution of present-day massive galaxies ($\log(M_{*}/M_{\odot})=11.5$) out to $z=5$. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by $\approx1.5~\mathrm{dex}$ since $z=5$. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by $z\sim2$, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build up at $r<3~\mathrm{kpc}$ beyond $z>4$ probing an era of significant mass assembly in the interiors of present day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with $z$ at $r<3~\mathrm{kpc}$. We identify $z\sim1.5$ as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in-situ star formation in disks to a period of efficient growth in $r_{e}$ consistent with the minor merger scenario.

We present the results of a pilot near-infrared (NIR) spectroscopic campaign of five very massive galaxies ($\log(\text{M}_\star/\text{M}_\odot)>11.45$) in the range of $1.7<z<2.7$. We measure an absorption feature redshift for one galaxy at $z_\text{spec}=2.000\pm0.006$. For the remaining galaxies, we combine the photometry with the continuum from the spectra to estimate continuum redshifts and stellar population properties. We define a continuum redshift ($z_{\rm cont}$ ) as one in which the redshift is estimated probabilistically using EAZY from the combination of catalog photometry and the observed spectrum. We derive the uncertainties on the stellar population synthesis properties using a Monte Carlo simulation and examine the correlations between the parameters with and without the use of the spectrum in the modeling of the spectral energy distributions (SEDs). The spectroscopic constraints confirm the extreme stellar masses of the galaxies in our sample. We find that three out of five galaxies are quiescent (star formation rate of $\lesssim 1 M_\odot~yr^{-1}$) with low levels of dust obscuration ($A_{\rm V} < 1$) , that one galaxy displays both high levels of star formation and dust obscuration (${\rm SFR} \approx 300 M_\odot~{\rm yr}^{-1}$, $A_{\rm V} \approx 1.7$~mag), and that the remaining galaxy has properties that are intermediate between the quiescent and star-forming populations.

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

Our goal is to study the evolution of the $B-$band luminosity function (LF) since $z=1$ using ALHAMBRA data. We used the photometric redshift and the $I-$band selection magnitude probability distribution functions (PDFs) of those ALHAMBRA galaxies with $I\leq24$ mag to compute the posterior LF. We statistically studied quiescent and star-forming galaxies using the template information encoded in the PDFs. The LF covariance matrix in redshift-magnitude-galaxy type space was computed, including the cosmic variance. That was estimated from the intrinsic dispersion of the LF measurements in the 48 ALHAMBRA sub-fields. The uncertainty due to the photometric redshift prior is also included in our analysis. We modelled the LF with a redshift-dependent Schechter function affected by the same selection effects than the data. The measured ALHAMBRA LF at $0.2\leq z<1$ and the evolving Schechter parameters both for quiescent and star-forming galaxies agree with previous results in the literature. The estimated redshift evolution of $M_B^* \propto Qz$ is $Q_{\rm SF}=-1.03\pm0.08$ and $Q_{\rm Q}=-0.80\pm0.08$, and of $\log \phi^* \propto Pz$ is $P_{\rm SF}=-0.01\pm0.03$ and $P_{\rm Q}=-0.41\pm0.05$. The measured faint-end slopes are $\alpha_{\rm SF}=-1.29\pm0.02$ and $\alpha_{\rm Q}=-0.53\pm0.04$. We find a significant population of faint quiescent galaxies, modelled by a second Schechter function with slope $\beta=-1.31\pm0.11$. We find a factor $2.55\pm0.14$ decrease in the luminosity density $j_B$ of star-forming galaxies, and a factor $1.25\pm0.16$ increase in the $j_B$ of quiescent ones since $z=1$, confirming the continuous build-up of the quiescent population with cosmic time. The contribution of the faint quiescent population to $j_B$ increases from 3% at $z=1$ to 6% at $z=0$. The developed methodology will be applied to future multi-filter surveys such as J-PAS.

The original ALHAMBRA catalogue contained over 400,000 galaxies selected using a synthetic F814W image, to the magnitude limit AB(F814W)$\approx$24.5. Given the photometric redshift depth of the ALHAMBRA multiband data (<z>=0.86) and the approximately $I$-band selection, there is a noticeable bias against red objects at moderate redshift. We avoid this bias by creating a new catalogue selected in the $K_s$ band. This newly obtained catalogue is certainly shallower in terms of apparent magnitude, but deeper in terms of redshift, with a significant population of red objects at $z>1$. We select objects using the $K_s$ band images, which reach an approximate AB magnitude limit $K_s \approx 22$. We generate masks and derive completeness functions to characterize the sample. We have tested the quality of the photometry and photometric redshifts using both internal and external checks. Our final catalogue includes $\approx 95,000$ sources down to $K_s \approx 22$, with a significant tail towards high redshift. We have checked that there is a large sample of objects with spectral energy distributions that correspond to that of massive, passively evolving galaxies at $z > 1$, reaching as far as $z \approx 2.5$. We have tested the possibility of combining our data with deep infrared observations at longer wavelengths, particularly Spitzer IRAC data.

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.

Using the UltraVISTA DR1 and 3D-HST catalogs, we construct a stellar-mass-complete sample, unique for its combination of surveyed volume and depth, to study the evolution of the fractions of quiescent galaxies, moderately unobscured star-forming galaxies, and dusty star-forming galaxies as a function of stellar mass over the redshift interval $0.2 \le z \le 3.0$. We show that the role of dusty star-forming galaxies within the overall galaxy population becomes more important with increasing stellar mass, and grows rapidly with increasing redshift. Specifically, dusty star-forming galaxies dominate the galaxy population with $\log{(M_{\rm star}/M_{\odot})} \gtrsim 10.3$ at $z\gtrsim2$. The ratio of dusty and non-dusty star-forming galaxies as a function of stellar mass changes little with redshift. Dusty star-forming galaxies dominate the star-forming population at $\log{(M_{\rm star}/M_{\odot})} \gtrsim 10.0-10.5$, being a factor of $\sim$3-5 more common, while unobscured star-forming galaxies dominate at $\log{(M_{\rm star}/M_{\odot})} \lesssim 10$. At $\log{(M_{\rm star}/M_{\odot})} > 10.5$, red galaxies dominate the galaxy population at all redshift $z<3$, either because they are quiescent (at late times) or dusty star-forming (in the early universe).

We present an overview of the "KIFF" project, which provides ultra-deep Ks-band imaging of all six of the Hubble Frontier Fields clusters Abell 2744, MACS-0416, Abell S1063, Abell 370, MACS-0717 and MACS-1149. All of these fields have recently been observed with large allocations of Directors' Discretionary Time with the HST and Spitzer telescopes covering 0.4 < lambda < 1.6 microns and 3.6--4.5 microns, respectively. VLT/HAWK-I integrations of the first four fields reach 5-sigma limiting depths of Ks~26.0 (AB, point sources) and have excellent image quality (FWHM ~ 0."4). Shorter Keck/MOSFIRE integrations of the MACS-0717 (MACS-1149) field better observable in the north reach limiting depths Ks=25.5 (25.1) with seeing FWHM ~0."4 (0."5). In all cases the Ks-band mosaics cover the primary cluster and parallel HST/ACS+WFC3 fields. The total area of the Ks-band coverage is 490 arcmin^2. The Ks-band at 2.2 microns crucially fills the gap between the reddest HST filter (1.6 micron ~ H-band) and the IRAC 3.6 micron passband. While reaching the full depths of the space-based imaging is not currently feasible from the ground, the deep Ks-band images provide important constraints on both the redshifts and the stellar population properties of galaxies extending well below the characteristic stellar mass across most of the age of the universe, down to, and including, the redshifts of the targeted galaxy clusters (z < 0.5).

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.

Using a sample of 98 galaxy clusters recently imaged in the near infra-red with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from the ESO archive, we measure how the stellar mass of the most massive galaxies in the universe, namely Brightest Cluster Galaxies (BCG), increases with time. Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$, which has been noted in recent works to mark an epoch over which the growth in the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a subsample of 102 systems that includes only those clusters that have estimates of the cluster mass. We combine the BCGs in this subsample with BCGs from the literature, and find that the growth in stellar mass of BCGs from 10 billion years ago to the present epoch is broadly consistent with recent semi-analytic and semi-empirical models. As in other recent studies, tentative evidence indicates that the stellar mass growth rate of BCGs may be slowing in the past 3.5 billion years. Further work in collecting larger samples, and in better comparing observations with theory using mock images is required if a more detailed comparison between the models and the data is to be made.

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 study the clustering of galaxies as a function of spectral type and redshift in the range $0.35 < z < 1.1$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The data cover 2.381 deg$^2$ in 7 fields, after applying a detailed angular selection mask, with accurate photometric redshifts [$\sigma_z < 0.014(1+z)$] down to $I_{AB} < 24$. From this catalog we draw five fixed number density, redshift-limited bins. We estimate the clustering evolution for two different spectral populations selected using the ALHAMBRA-based photometric templates: quiescent and star-forming galaxies. For each sample, we measure the real-space clustering using the projected correlation function. Our calculations are performed over the range $[0.03,10.0] h^{-1}$ Mpc, allowing us to find a steeper trend for $r_p \lesssim 0.2 h^{-1}$ Mpc, which is especially clear for star-forming galaxies. Our analysis also shows a clear early differentiation in the clustering properties of both populations: star-forming galaxies show weaker clustering with evolution in the correlation length over the analysed redshift range, while quiescent galaxies show stronger clustering already at high redshifts, and no appreciable evolution. We also perform the bias calculation where similar segregation is found, but now it is among the quiescent galaxies where a growing evolution with redshift is clearer. These findings clearly corroborate the well known colour-density relation, confirming that quiescent galaxies are mainly located in dark matter halos that are more massive than those typically populated by star-forming galaxies.

The growth of galaxies is a key problem in understanding the structure and evolution of the universe. Galaxies grow their stellar mass by a combination of star formation and mergers, with a relative importance that is redshift dependent. Theoretical models predict quantitatively different contributions from the two channels; measuring these from the data is a crucial constraint. Exploiting the UltraVISTA catalog and a unique sample of progenitors of local ultra massive galaxies selected with an abundance matching approach, we quantify the role of the two mechanisms from z=2 to 0. We also compare our results to two independent incarnations of semi-analytic models. At all redshifts, progenitors are found in a variety of environments, ranging from being isolated to having 5-10 companions with mass ratio at least 1:10 within a projected radius of 500 kpc. In models, progenitors have a systematically larger number of companions, entailing a larger mass growth for mergers than in observations, at all redshifts. Generally, in both observations and models, the inferred and the expected mass growth roughly agree, within the uncertainties. Overall, our analysis confirms the model predictions, showing how the growth history of massive galaxies is dominated by in situ star formation at z~2, both star-formation and mergers at 1<z<2, and by mergers alone at z<1. Nonetheless, detailed comparisons still point out to tensions between the expected mass growth and our results, which might be due to either an incorrect progenitors-descendants selection, uncertainties on star formation rate and mass estimates, or the adopted assumptions on merger rates.

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 present the first spectroscopic confirmation of an ultra-massive galaxy at redshift z>3 using data from Keck-NIRSPEC, VLT-Xshooter, and GTC-Osiris. We detect strong [OIII] and Ly$\alpha$ emission, and weak [OII], CIV, and HeII, placing C1-23152 at a spectroscopic redshift of $z_{spec}$=3.351. The modeling of the emission-line corrected spectral energy distribution results in a best-fit stellar mass of $M_{*}=3.1^{+0.6}_{-0.7}\times10^{11} M_{\odot}$, a star-formation rate of <7 $M_{\odot} yr^{-1}$, and negligible dust extinction. The stars appear to have formed in a short intense burst ~300-500 Myr prior to the observation epoch, setting the formation redshift of this galaxy at z~4.1. From the analysis of the line ratios and widths, and the observed flux at 24$\mu$m, we confirm the presence of a luminous hidden active galactic nucleus (AGN), with bolometric luminosity of ~$10^{46}erg s^{-1}$. Potential contamination to the observed SED from the AGN continuum is constrained, placing a lower limit on the stellar mass of $2\times10^{11} M_{\odot}$. HST/WFC3 $H_{160}$ and ACS $I_{814}$ images are modeled, resulting in an effective radius of $r_{e}$~1 kpc in the $H_{160}$ band and a Sersic index n~4.4. This object may be a prototype of the progenitors of local most massive elliptical galaxies in the first 2 Gyr of cosmic history, having formed most of its stars at z>4 in a highly dissipative, intense, and short burst of star formation. C1-23152 is completing its transition to a post-starburst phase while hosting a powerful AGN, potentially responsible for the quenching of the star formation activity.

We build a Spitzer IRAC complete catalog of objects, obtained by complementing the $K_\mathrm{s}$-band selected UltraVISTA catalog with objects detected in IRAC only. With the aim of identifying massive (i.e., $\log(M_*/M_\odot)>11$) galaxies at $4<z<7$, we consider the systematic effects on the measured photometric redshifts from the introduction of an old and dusty SED template and from the introduction of a bayesian prior taking into account the brightness of the objects, as well as the systematic effects from different star formation histories (SFHs) and from nebular emission lines in the recovery of stellar population parameters. We show that our results are most affected by the bayesian luminosity prior, while nebular emission lines and SFHs only introduce a small dispersion in the measurements. Specifically, the number of $4<z<7$ galaxies ranges from 52 to 382 depending on the adopted configuration. Using these results we investigate, for the first time, the evolution of the massive end of the stellar mass functions (SMFs) at $4<z<7$. Given the rarity of very massive galaxies in the early universe, major contributions to the total error budget come from cosmic variance and poisson noise. The SMF obtained without the introduction of the bayesian luminosity prior does not show any evolution from $z\sim6.5$ to $z\sim 3.5$, implying that massive galaxies could already be present when the Universe was $\sim0.9$~Gyr old. However, the introduction of the bayesian luminosity prior reduces the number of $z>4$ galaxies with best fit masses $\log(M_*/M_\odot)>11$ by 83%, implying a rapid growth of very massive galaxies in the first 1.5 Gyr of cosmic history. From the stellar-mass complete sample, we identify one candidate of a very massive ($\log(M_*/M_\odot)\sim11.5$), quiescent galaxy at $z\sim5.4$, with MIPS $24\mu$m detection suggesting the presence of a powerful obscured AGN.

Using the UltraVISTA catalogs, we investigate the evolution in the 11.4~Gyr since $z=3$ of the progenitors of local ultra-massive galaxies ($\log{(M_{\rm star}/M_{\odot})}\approx11.8$; UMGs), providing a complete and consistent picture of how the most massive galaxies at $z=0$ have assembled. By selecting the progenitors with a semi-empirical approach using abundance matching, we infer a growth in stellar mass of 0.56$^{+0.35}_{-0.25}$ dex, 0.45$^{+0.16}_{-0.20}$~dex, and 0.27$^{+0.08}_{-0.12}$ dex from $z=3$, $z=2$, and $z=1$, respectively, to $z=0$. At $z<1$, the progenitors of UMGs constitute a homogeneous population of only quiescent galaxies with old stellar populations. At $z>1$, the contribution from star-forming galaxies progressively increases, with the progenitors at $2<z<3$ being dominated by massive ($M_{\rm star} \approx 2 \times 10^{11}$M$_{\odot}$), dusty ($A_{\rm V}\sim$1--2.2 mag), star-forming (SFR$\sim$100--400~M$_{\odot}$ yr$^{-1}$) galaxies with a large range in stellar ages. At $z=2.75$, $\sim$15\% of the progenitors are quiescent, with properties typical of post-starburst galaxies with little dust extinction and strong Balmer break, and showing a large scatter in color. Our findings indicate that at least half of the stellar content of local UMGs was assembled at $z>1$, whereas the remaining was assembled via merging from $z\sim 1$ to the present. Most of the quenching of the star-forming progenitors happened between $z=2.75$ and $z=1.25$, in good agreement with the typical formation redshift and scatter in age of $z=0$ UMGs as derived from their fossil records. The progenitors of local UMGs, including the star-forming ones, never lived on the blue cloud since $z=3$. We propose an alternative path for the formation of local UMGs that refines previously proposed pictures and that is fully consistent with our findings.

We present an observational study of the stellar mass function of satellite galaxies around central galaxies at 0.2<z<1.2. Using statistical background subtraction of contaminating sources we derive satellite stellar mass distributions in four bins of central galaxy mass in three redshift ranges. Our results show that the stellar mass function of satellite galaxies increases with central galaxy mass, and that the distribution of satellite masses at fixed central mass is at most weakly dependent on redshift. We conclude that the average mass distribution of galaxies in groups is remarkably universal even out to z=1.2 and that it can be uniquely characterized by the group central galaxy mass. This further suggests that as central galaxies grow in stellar mass, they do so in tandem with the mass growth of their satellites. Finally, we classify all galaxies as either star forming or quiescent, and derive the mass functions of each subpopulation separately. We find that the mass distribution of both star forming and quiescent satellites show minimal redshift dependence at fixed central mass. However, while the fraction of quiescent satellite galaxies increases rapidly with increasing central galaxy mass, that of star forming satellites decreases.

We study the clustering of galaxies as function of luminosity and redshift in the range $0.35 < z < 1.25$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The ALHAMBRA data used in this work cover $2.38 \mathrm{deg}^2$ in 7 independent fields, after applying a detailed angular selection mask, with accurate photometric redshifts, $\sigma_z \lesssim 0.014 (1+z)$, down to $I_{\rm AB} < 24$. Given the depth of the survey, we select samples in $B$-band luminosity down to $L^{\rm th} \simeq 0.16 L^{*}$ at $z = 0.9$. We measure the real-space clustering using the projected correlation function, accounting for photometric redshifts uncertainties. We infer the galaxy bias, and study its evolution with luminosity. We study the effect of sample variance, and confirm earlier results that the COSMOS and ELAIS-N1 fields are dominated by the presence of large structures. For the intermediate and bright samples, $L^{\rm med} \gtrsim 0.6L^{*}$, we obtain a strong dependence of bias on luminosity, in agreement with previous results at similar redshift. We are able to extend this study to fainter luminosities, where we obtain an almost flat relation, similar to that observed at low redshift. Regarding the evolution of bias with redshift, our results suggest that the different galaxy populations studied reside in haloes covering a range in mass between $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 11.5$ for samples with $L^{\rm med} \simeq 0.3 L^{*}$ and $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 13.0$ for samples with $L^{\rm med} \simeq 2 L^{*}$, with typical occupation numbers in the range of $\sim 1 - 3$ galaxies per halo.

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

Comparing galaxies across redshifts at fixed cumulative number density is a popular way to estimate the evolution of specific galaxy populations. This method ignores scatter in mass accretion histories and galaxy-galaxy mergers, which can lead to errors when comparing galaxies over large redshift ranges (Delta z > 1). We use abundance matching in the LCDM paradigm to estimate the median change in number density with redshift and provide a simple fit (+0.16 dex per unit Delta z) for progenitors of z = 0 galaxies. We find that galaxy descendants do not evolve in the same way as galaxy progenitors, largely due to scatter in mass accretion histories. We also provide estimates for the 1-sigma range of number densities corresponding to galaxy progenitors and descendants. Finally, we discuss some limits on number density comparisons, which arise due to difficulties measuring physical quantities (e.g., stellar mass) consistently across redshifts. A public tool to calculate number density evolution for galaxies, as well as approximate halo masses, is available online.

We present measurements of the stellar mass functions (SMFs) of star-forming and quiescent galaxies to z = 4 using a sample of 95 675 galaxies in the COSMOS/UltraVISTA field. Sources have been selected from the DR1 UltraVISTA K_{s}-band imaging which covers a unique combination of a wide area (1.62 deg^2), to a significant depth (K_{s,tot} = 23.4). The SMFs of the combined population are in good agreement with previous measurements and show that the stellar mass density of the universe was only 50%, 10% and 1% of its current value at z ~ 0.75, 2.0, and 3.5, respectively. The quiescent population drives most of the overall growth, with the stellar mass density of these galaxies increasing by 2.71^{+0.93}_{-0.22} dex since z = 3.5. At z > 2.5, star-forming galaxies dominate the total SMF at all stellar masses, although a nonzero population of quiescent galaxies persists to z = 4. Comparisons of the K_{s}-selected star-forming galaxy SMFs to UV-selected SMFs at 2.5 < z < 4 show reasonable agreement and suggests UV-selected samples are representative of the majority of the stellar mass density at z > 3.5. We estimate the average mass growth of individual galaxies by selecting galaxies at fixed cumulative number density. The average galaxy with Log(M_{*}/M_{sun}) = 11.5 at z = 0.3 has grown in mass by only 0.2 dex (0.3 dex) since z = 2.0(3.5), whereas those with Log(M_{*}/M_{sun}) = 10.5 have grown by > 1.0 dex since z = 2. At z < 2, the time derivatives of the mass growth are always larger for lower-mass galaxies, which demonstrates that the mass growth in galaxies since that redshift is mass-dependent and primarily bottom-up. Lastly, we examine potential sources of systematic uncertainties on the SMFs and find that those from photo-z templates, SPS modeling, and the definition of quiescent galaxies dominate the total error budget in the SMFs.

We present the stellar mass functions (SMFs) of star-forming and quiescent galaxies from observations of 10 rich clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS) in the redshift range 0.86<z<1.34. We compare our results with field measurements at similar redshifts using data from a Ks-band selected catalogue of the COSMOS/UltraVISTA field. We construct a Ks-band selected multi-colour catalogue for the clusters in 11 photometric bands covering u-8um, and estimate photometric redshifts and stellar masses using SED fitting techniques. To correct for interlopers in our cluster sample, we use the deep spectroscopic component of GCLASS, which contains spectra for 1282 identified cluster and field galaxies taken with Gemini/GMOS. Both the photometric and spectroscopic samples are sufficiently deep that we can probe the SMF down to masses of 10^10 Msun. We distinguish between star-forming and quiescent galaxies using the rest-frame U-V versus V-J diagram, and find that the best-fitting Schechter parameters alpha and M* are similar within the uncertainties for these galaxy types within the different environments. However, there is a significant difference in the shape and normalisation of the total SMF between the clusters and the field sample. This difference in the total SMF is primarily a reflection of the increased fraction of quiescent galaxies in high-density environments. We apply a simple quenching model that includes components of mass- and environment-driven quenching, and find that in this picture 45% of the star-forming galaxies, which normally would be forming stars in the field, are quenched by the cluster. If galaxies in clusters and the field quench their star formation via different mechanisms, these processes have to conspire in such a way that the shapes of the quiescent and star-forming SMF remain similar in these different environments.

We present a catalog covering 1.62 deg^2 of the COSMOS/UltraVISTA field with PSF-matched photometry in 30 photometric bands. The catalog covers the wavelength range 0.15um - 24um including the available GALEX, Subaru, CFHT, VISTA and Spitzer data. Catalog sources have been selected from the DR1 UltraVISTA Ks band imaging that reaches a depth of K_{s,tot} = 23.4 AB (90% completeness). The PSF-matched catalog is generated using position-dependent PSFs ensuring accurate colors across the entire field. Also included is a catalog of photometric redshifts (z_phot) for all galaxies computed with the EAZY code. Comparison with spectroscopy from the zCOSMOS 10k bright sample shows that up to z ~ 1.5 the z_phot are accurate to dz/(1 + z) = 0.013, with a catastrophic outlier fraction of only 1.6%. The z_phot also show good agreement with the z_phot from the NEWFIRM Medium Band Survey (NMBS) out to z ~ 3. A catalog of stellar masses and stellar population parameters for galaxies determined using the FAST spectral energy distribution fitting code is provided for all galaxies. Also included are rest-frame U-V and V-J colors, L_2800 and L_IR. The UVJ color-color diagram confirms that the galaxy bi-modality is well-established out to z ~ 2. Star-forming galaxies also obey a star forming "main sequence" out to z ~ 2.5, and this sequence evolves in a manner consistent with previous measurements. The COSMOS/UltraVISTA Ks-selected catalog covers a unique parameter space in both depth, area, and multi-wavelength coverage and promises to be a useful tool for studying the growth of the galaxy population out to z ~ 3 - 4.

We study the structural evolution of massive galaxies by linking progenitors and descendants at a constant cumulative number density of n_c=1.4x10^{-4} Mpc^{-3} to z~3. Structural parameters were measured by fitting Sersic profiles to high resolution CANDELS HST WFC3 J_{125} and H_{160} imaging in the UKIDSS-UDS at 1<z<3 and ACS I_{814} imaging in COSMOS at 0.25<z<1. At a given redshift, we selected the HST band that most closely samples a common rest-frame wavelength so as to minimize systematics from color gradients in galaxies. At fixed n_c, galaxies grow in stellar mass by a factor of ~3 from z~3 to z~0. The size evolution is complex: galaxies appear roughly constant in size from z~3 to z~2 and then grow rapidly to lower redshifts. The evolution in the surface mass density profiles indicates that most of the mass at r<2 kpc was in place by z~2, and that most of the new mass growth occurred at larger radii. This inside-out mass growth is therefore responsible for the larger sizes and higher Sersic indices of the descendants toward low redshift. At z<2, the effective radius evolves with the stellar mass as r_e M^{2.0}, consistent with scenarios that find dissipationless minor mergers to be a key driver of size evolution. The progenitors at z~3 were likely star-forming disks with r_e~2 kpc, based on their low Sersic index of n~1, low median axis ratio of b/a~0.52, and typical location in the star-forming region of the U-V versus V-J diagram. By z~1.5, many of these star-forming disks disappeared, giving rise to compact quiescent galaxies. Toward lower redshifts, these galaxies continued to assemble mass at larger radii and became the local ellipticals that dominate the high mass end of the mass function at the present epoch.

[Abridged] Using public data from the NMBS and CANDELS surveys, we study the population of massive galaxies at z>3 to identify the potential progenitors of z~2 compact, massive, quiescent (CMQ) galaxies, furthering our understanding of the evolution of massive galaxies. Our work is enabled by high-resolution CANDELS images and accurate photometric redshifts, stellar masses and star formation rates (SFRs) from 37-band NMBS photometry. The total number of z>3 massive galaxies is consistent with the number of massive quiescent (MQ) galaxies at z~2, implying that the SFRs for all of these galaxies must be much lower by z~2. We discover 4 CMQ galaxies at z>3, pushing back the time for which such galaxies have been observed. However, the volume density for these galaxies is significantly less than that of galaxies at z<2 with similar masses, SFRs, and sizes, implying that additional CMQ galaxies must be created in the ~1 Gyr between z=3 and z=2. We find 5 star-forming galaxies at z~3 that are compact (Re<1.4 kpc) and have stellar mass M*>10^(10.6)Msun, likely to become members of the CMQ galaxy population at z~2. We evolve the stellar masses and SFRs of each individual z>3 galaxy adopting 5 different star formation histories (SFHs) and studying the resulting population of massive galaxies at z=2.3. We find that declining or truncated SFHs are necessary to match the observed number density of MQ galaxies at z~2, whereas a constant SFH results in a number density significantly smaller than observed. All of our assumed SFHs imply number densities of CMQ galaxies at z~2 that are consistent with the observed number density. Better agreement with the observed number density of CMQ galaxies at z~2 is obtained if merging is included in the analysis and better still if star formation quenching is assumed to shortly follow the merging event, as implied by recent models of formation of MQ galaxies.

We present the rest-frame J- and H-band luminosity function (LF) of field galaxies, based on a deep multi-wavelength composite sample from the MUSYC, FIRES and FIREWORKS survey public catalogues, covering a total area of 450 arcmin^2. The availability of flux measurements in the Spitzer IRAC 3.6, 4.5, 5.8, and 8 um channels allows us to compute absolute magnitudes in the rest-frame J and H bands up to z=3.5 minimizing the dependence on the stellar evolution models. We compute the LF in the four redshift bins 1.5<z<2.0, 2.0<z<2.5, 2.5<z<3.0 and 3.0<z<3.5. Combining our results with those already available at lower redshifts, we find that (1) the faint end slope is consistent with being constant up to z=3.5, with alpha=-1.05+/-0.03 for the rest-frame J band and alpha=-1.15+/-0.02 for the rest-frame H band; (2) the normalization phi* decreases by a factor of 6 between z=0 and z~1.75 and by a factor 3 between z~1.75 and z=3.25; (3) the characteristic magnitude M* shows a brightening from z=0 to z~2 followed by a slow dimming to z=3.25. We finally compute the luminosity density (LD) in both rest-frame J and H bands. The analysis of our results together with those available in the literature shows that the LD is approximately constant up to z~1, and it then decreases by a factor of 6 up to z=3.5.

We report the discovery of a massive ultra-compact quiescent galaxy that has been strongly-lensed into multiple images by a foreground galaxy at z = 0.960. This system was serendipitously discovered as a set of extremely Ks-bright high-redshift galaxies with red J - Ks colors using new data from the UltraVISTA YJHKs near-infrared survey. The system was also previously identified as an optically-faint lens/source system using the COSMOS ACS imaging by Faure et al. (2008, 2011). Photometric redshifts for the three brightest images of the source galaxy determined from twenty-seven band photometry place the source at z = 2.4 +/- 0.1. We provide an updated lens model for the system which is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model implies that the magnification of the three brightest images is a factor of 4 - 5. We use the lens model, combined with the Ks-band image to constrain the size and Sersic profile of the galaxy. The best-fit model is an ultra-compact galaxy (Re = 0.64^{+0.08}_{-0.18} kpc, lensing-corrected), with a Sersic profile that is intermediate between a disk and bulge profile (n = 2.2^{+2.3}_{-0.9}). We present aperture photometry for the source galaxy images which have been corrected for flux contamination from the central lens. The best-fit stellar population model is a massive galaxy (Log(M_{star}/M_{sol}) = 10.8^{+0.1}_{-0.1}, lensing-corrected) with an age of 1.0^{+1.0}_{-0.4} Gyr, moderate dust extinction (Av = 0.8^{+0.5}_{-0.6}), and a low specific star formation rate (Log(SSFR) < -11.0 yr^{-1}). This is typical of massive "red-and-dead" galaxies at this redshift and confirms that this source is the first bona fide strongly-lensed massive ultra-compact quiescent galaxy to be discovered. We conclude with a discussion of the prospects of finding a larger sample of these galaxies.

We present the rest-frame V-band luminosity function (LF) of galaxies at 0.4<z<4.0, measured from a near-infrared selected sample constructed from the NMBS, the FIRES, the FIREWORKS, and the ultra-deep NICMOS and WFC3 observations in the HDFN, HUDF, and GOODS-CDFS, all having high-quality optical to mid-infrared data. This unique sample combines data from surveys with a large range of depths and areas in a self-consistent way, allowing us to (1) minimize the uncertainties due to cosmic variance; and (2) simultaneously constrain the bright and faint ends with unprecedented accuracy over the targeted redshift range, probing the LF down to 0.1 L* at z~3.9. We find that (1) the faint end is fairly flat and with a constant slope from z=4, with alpha= -1.27 +/- 0.05; (2) the characteristic magnitude has dimmed by 1.3 mag from z~3.7 to z=0.1; (3) the characteristic density has increased by a factor of ~8 from z~3.7 to z=0.1, with 50% of this increase from z~4 to z~1.8; and (4) the luminosity density peaks at z~1-1.5, increasing by a factor of ~4 from z=4.0 to z~1-1.5, and subsequently decreasing by a factor of ~1.5 by z=0.1. We find no evidence for a steepening of the faint-end slope with redshift out to z=4, in contrast with previous observational claims and theoretical predictions. The constant faint-end slope suggests that the efficiency of stellar feedback may evolve with redshift. Alternative interpretations are discussed, such as different masses of the halos hosting faint galaxies at low and high redshifts and/or environmental effects.

We have developed a new method, close in philosophy to the photometric redshift technique, which can be applied to spectral data of very low signal-to-noise ratio. Using it we intend to measure redshifts while minimising the dangers posed by the usual extraction techniques. GRB afterglows have generally very simple optical spectra over which the separate effects of absorption and reddening in the GRB host, the intergalactic medium, and our own Galaxy are superimposed. We model all these effects over a series of template afterglow spectra to produce a set of clean spectra that reproduce what would reach our telescope. We also model carefully the effects of the telescope-spectrograph combination and the properties of noise in the data, which are then applied on the template spectra. The final templates are compared to the two-dimensional spectral data, and the basic parameters (redshift, spectral index, Hydrogen absorption column) are estimated using statistical tools. We show how our method works by applying it to our data of the NIR afterglow of GRB090423. At z ~ 8.2, this was the most distant object ever observed. We use the spectrum taken by our team with the Telescopio Nazionale Galileo to derive the GRB redshift and its intrinsic neutral Hydrogen column density. Our best fit yields z=8.4^+0.05/-0.03 and N(HI)<5x10^20 cm^-2, but with a highly non-Gaussian uncertainty including the redshift range z [6.7, 8.5] at the 2-sigma confidence level. Our method will be useful to maximise the recovered information from low-quality spectra, particularly when the set of possible spectra is limited or easily parameterisable while at the same time ensuring an adequate confidence analysis.

Context: Multiwavelength observations of gamma-ray burst prompt and afterglow emission are a key tool to disentangle the various possible emission processes and scenarios proposed to interpret the complex gamma-ray burst phenomenology. Aims: We collected a large dataset on GRB060908 in order to carry out a comprehensive analysis of the prompt emission as well as the early and late afterglow. Methods: Data from Swift-BAT, -XRT and -UVOT together with data from a number of different ground-based optical/NIR and millimeter telescopes allowed us to follow the afterglow evolution from about a minute from the high-energy event down to the host galaxy limit. We discuss the physical parameters required to model these emissions. Results: The prompt emission of GRB060908 was characterized by two main periods of activity, spaced by a few seconds of low intensity, with a tight correlation between activity and spectral hardness. Observations of the afterglow began less than one minute after the high-energy event, when it was already in a decaying phase, and it was characterized by a rather flat optical/NIR spectrum which can be interpreted as due to a hard energy-distribution of the emitting electrons. On the other hand, the X-ray spectrum of the afterglow could be fit by a rather soft electron distribution. Conclusions: GRB060908 is a good example of a gamma-ray burst with a rich multi-wavelength set of observations. The availability of this dataset, built thanks to the joint efforts of many different teams, allowed us to carry out stringent tests for various interpretative scenarios showing that a satisfactorily modeling of this event is challenging. In the future, similar efforts will enable us to obtain optical/NIR coverage comparable in quality and quantity to the X-ray data for more events, therefore opening new avenues to progress gamma-ray burst research.

The approach of Observational Astronomy is mainly aimed at the construction of larger aperture telescopes, more sensitive detectors and broader wavelength coverage. Certainly fruitful, this approach turns out to be not completely fulfilling the needs when phenomena related to the formation of black holes (BH), neutron stars (NS) and relativistic stars in general are concerned. Recently, mainly through the Vela, Beppo-SAX and Swift satellites, we reached a reasonable knowledge of the most violent events in the Universe and of some of the processes we believe are leading to the formation of black holes (BH). We plan to open a new window of opportunity to study the variegated physics of very fast astronomical transients, particularly the one related to extreme compact objects. The innovative approach is based on three cornerstones: 1) the design (the conceptual design has been already completed) of a 3m robotic telescope and related focal plane instrumentation characterized by the unique features: "No telescope points faster"; 2) simultaneous multi-wavelengths observations (photometry, spectroscopy o\& polarimetry); 3) high time resolution observations. The conceptual design of the telescope and related instrumentation is optimized to address the following topics: High frequency a-periodic variability, Polarization, High z GRBs, Short GRBs, GRB-Supernovae association, Multi-wavelengths simultaneous photometry and rapid low dispersion spectroscopy. This experiment will turn the "exception" (like the optical observations of GRB 080319B) to "routine".

The Telescopi ROBotic de ARas (TROBAR) is a new robotic facility built at Aras de Los Olmos (Valencia-Spain). This is a 60cm telescope equipped with a 4kx4k optical camera, corresponding to 30x30 arcmin2 FoV, and it will be primarily used for a systematic search of Ha emitting stars in the Galactic Plane to a depth of ~14mag. Both data acquisition and reduction will be performed automatically. The robotization of data acquisition is now entering its final phase while the development of the data reduction pipeline has just started.

Gamma-ray burst are thought to be produced by highly relativistic outflows. Although upper and lower limits for the outflow initial Lorentz factor $\Gamma_0$ are available, observational efforts to derive a direct determination of $\Gamma_0$ have so far failed or provided ambiguous results. As a matter of fact, the shape of the early-time afterglow light curve is strongly sensitive on $\Gamma_0$ which determines the time of the afterglow peak, i.e. when the outflow and the shocked circumburst material share a comparable amount of energy. We now comment early-time observations of the near-infrared afterglows of GRB 060418 and GRB 060607A performed by the REM robotic telescope. For both events, the afterglow peak was singled out and allowed us to determine the initial fireball Lorentz, $\Gamma_0\sim 400$.

The Variable Star One-shot Project (VSOP) is aimed at (1) providing the variability type and spectral type of all unstudied variable stars, (2) process, publish, and make the data available as automatically as possible, and (3) generate serendipitous discoveries. This first paper describes the project itself, the acquisition of the data, the dataflow, the spectroscopic analysis and the on-line availability of the fully calibrated and reduced data. We also present the results on the 221 stars observed during the first semester of the project. We used the high-resolution echelle spectrographs HARPS and FEROS in the ESO La Silla Observatory (Chile) to survey known variable stars. Once reduced by the dedicated pipelines, the radial velocities are determined from cross correlation with synthetic template spectra, and the spectral types are determined by an automatic minimum distance matching to synthetic spectra, with traditional manual spectral typing cross-checks. The variability types are determined by manually evaluating the available light curves and the spectroscopy. In the future, a new automatic classifier, currently being developed by members of the VSOP team, based on these spectroscopic data and on the photometric classifier developed for the COROT and Gaia space missions, will be used. We confirm or revise spectral types of 221 variable stars from the GCVS. We identify 26 previously unknown multiple systems, among them several visual binaries with spectroscopic binary individual components. We present new individual results for the multiple systems V349 Vel and BC Gru, for the composite spectrum star V4385 Sgr, for the T-Tauri star V1045 Sco, and for DM Boo which we re-classify as a BY Draconis variable. The complete data release can be accessed via the VSOP web site.

The Variable Star One-shot Project (VSOP) aimed at providing to the world-wide stellar community the necessary one-shot spectrum of unstudied variable stars, too often classified as such by an analysis of photometric data only. The VSOP has established an new kind of observational model, where all steps from observations to spectral analysis, are automatized (or are underway to be fully automatized). The project is centralized on a collaborative wiki website. The VSOP operational model is very successful, data is continously flowing and being analyszed, and VSOP is now a worldwide open collaboration of people with very different and complementary skills and expertise. The idea of a central wiki website has been extended by one of us to propose a new service to the whole astronomical community, called Wikimbad. Wikimbad is an open wiki website aimed at collecting, organizing and making publicly available all kind of reduced and published astronomical data. Its strengths and a comparison with the Virtual Observatory are discussed. See: http://vsop.sc.eso.org and http://wikimbad.org

Gamma-ray burst (GRB) emission is believed to originate in highly relativistic fireballs. Currently, only lower limits were securely set to the initial fireball Lorentz factor Gamma_0. We aim to provide a direct measure of Gamma_0. The early-time afterglow light curve carries information about Gamma_0, which determines the time of the afterglow peak. We have obtained early observations of the near-infrared afterglows of GRB 060418 and GRB 060607A with the REM robotic telescope. For both events, the afterglow peak could be clearly singled out, allowing a firm determination of the fireball Lorentz of Gamma_0 ~ 400, fully confirming the highly relativistic nature of GRB fireballs. The deceleration radius was inferred to be R_dec ~ 10^17 cm. This is much larger than the internal shocks radius (believed to power the prompt emission), thus providing further evidence for a different origin of the prompt and afterglow stages of the GRB.

We present ground-based and Hubble Space Telescope optical observations of the X-ray flash (XRF) 020903, covering 300 days. The afterglow showed a very rapid rise in the first day, followed by a relatively slow decay in the next few days. There was a clear bump in the light curve after ~25 days, accompanied by a drastic change in the spectral energy distribution. The light curve and the spectral energy distribution are naturally interpreted as the emergence -- and subsequent decay -- of a supernova (SN), similar to SN 1998bw. At peak luminosity, the SN is estimated to be 0.8 +/- 0.1 mag fainter than SN1998bw. This argues in favor of the existence of a supernova associated with this X-ray flash. A spectrum obtained 35 days after the burst shows emission lines from the host galaxy. We use this spectrum to put an upper limit on the oxygen abundance of the host at [O/H] < -0.6 dex. We also discuss a possible trend between the softness of several bursts and the early behavior of the optical afterglow, in the sense that XRFs and X-ray rich GRBs seem to have a plateau phase or even a rising light curve. This can be naturally explained in models where XRFs are similar to GRBs but seen off the jet axis.

Multiband VRIJHK photometry of the Blazar PKS 0537-441 obtained with the REM telescope from December 2004 to March 2005 is presented. A major period of activity is found with more than four magnitudes variability in the V filter in 50 days and of 2.5 in 10 days. In intensity and duration the activity is similar to that of 1972 reported by Eggen (1973), but it is much better documented. No clear evidence of variability on time-scale of minutes is found. The spectral energy distribution is roughly described by a power-law, with the weaker state being the softer.

GRB 021004 is one of the best sampled gamma-ray bursts (GRB) to date, although the nature of its light curve is still being debated. Here we present a large amount (107) of new optical, near-infrared (NIR) and millimetre observations, ranging from 2 hours to more than a year after the burst. Fitting the multiband data to a model based on multiple energy injections suggests that at least 7 refreshed shocks took place during the evolution of the afterglow, implying a total energy release (collimated within an angle of 1.8 deg) of ~ 8x10^51 erg. Analysis of the late photometry reveals that the GRB 021004 host is a low extinction (Av ~ 0.1) starburst galaxy with M_B ~ -22.0.

We present optical and near-infrared (NIR) photometry of the bright afterglow of GRB 020813. Our data span from 3 hours to 4 days after the GRB event. A rather sharp achromatic break is present in the light curve, 14 hours after the trigger. In the framework of jetted fireballs, this break corresponds to a jet half opening angle of (1.9 +- 0.2) deg, the smallest value ever inferred for a GRB. We discuss our results in the framework of currently available models, and find that they have problems in explaining the joint temporal and spectral properties, and in particular the slow decay before the break.

Polarization measurements for the optical counterpart to GRB 020405 are presented and discussed. Our observations were performed with the VLT-UT3 (Melipal) during the second and third night after the gamma-ray burst discovery. The polarization degree (and the position angle) appears to be constant between our two observations at a level around (1.5 - 2)%. The polarization can be intrinsic but it is not possible to unambiguously exclude that a substantial fraction of it is induced by dust in the host galaxy.

The results of a CCD survey for variability of stars in the nearby galaxy NGC 6822 are presented. The goal of the survey was to obtain good light curves of Cepheids for Fourier decomposition and to detect shorter period Cepheids. Since the program was carried out with a relatively small telescope, the Dutch 0.9 m at ESO-La Silla, the observations were unfiltered (white light, or Wh-band). The analysis revealed the presence of more than 130 variable stars. 21 population I Cepheids are detected; 6 of them were already known from previous works (Kayser, 1967). For at least three Cepheids, however, the previous identification or period was wrong. Some probable population II (W Vir) stars are also identified. The dispersion of the fundamental mode Cepheid PL relation appears to be small.