Radial age gradients hold the cumulative record of the multitude of physical processes driving the build-up of stellar populations and the ensuing star formation (SF) quenching process in galaxy bulges, therefore potentially sensitive discriminators between competing theoretical concepts on bulge formation and evolution. Based on spectral modeling of integral field spectroscopy data from the CALIFA survey, we derive mass- and light-weighted stellar age gradients ($\nabla$(t,B)L,M) within the photometrically determined bulge radius (RB) of a representative sample of local face-on late-type galaxies that span 2.6 dex in stellar mass. Our analysis documents a trend for decreasing $\nabla$(t,B)L,M with increasing M,T, with high-mass bulges predominantly showing negative age gradients and vice versa. The inversion from positive to negative $\nabla$(t,B)L,M occurs at logM,T ~ 10, which roughly coincides with the transition from lower-mass bulges whose gas excitation is powered by SF to bulges classified as Composite, LINER or Seyfert. We discuss two limiting cases for the origin of radial age gradients in massive LTG bulges. The first assumes that the stellar age in the bulge is initially spatially uniform, thus the observed age gradients arise from an inside-out SF quenching (ioSFQ) front that is radially expanding with a mean velocity vq. In this case, the age gradients translate into a slow ioSFQ that lasts until z~2, suggesting mild negative feedback by SF or an AGN. If negative age gradients in massive bulges are not due to ioSFQ but primarily due to their inside-out formation process, then the standard hypothesis of quasi-monolithic bulge formation has to be discarded in favor of a scenario that involves gradual buildup of stellar mass over 2-3 Gyr through, e.g., inside-out SF and inward migration of SF clumps from the disk. In this case, rapid AGN-driven ioSFQ cannot be ruled out.

Taking advantage of HST CANDELS data, we analyze the lowest redshift (z<0.5) massive galaxies in order to disentangle their structural constituents and study possible faint non-axis-symmetric features. Due to the excellent HST spatial resolution for intermediate-z objects, they are hard to model by purely automatic parametric fitting algorithms. We performed careful single and double S\'ersic fits to their galaxy surface brightness profiles. We also compare the model color profiles with the observed ones and also derive multi-component global effective radii attempting to obtain a better interpretation of the mass-size relation. Additionally, we test the robustness of our measured structural parameters via simulations. We find that the S\'ersic index does not offer a good proxy for the visual morphological type for our sample of massive galaxies. Our derived multi-component effective radii give a better description of the size of our sample galaxies than those inferred from single S\'ersic models with GALFIT. Our galaxy population lays on the scatter of the local mass-size relation, indicating that these massive galaxies do not experience a significant growth in size since z~0.5. Interestingly the few outliers are late-type galaxies, indicating that spheroids must reach the local mass-size relation earlier. For most of our sample galaxies, both single and multi-component S\'ersic models with GALFIT show substantial systematic deviations from the observed SBPs in the outskirts. These residuals may be partly due to several factors, namely a non-optimal data reduction for low surface brightness features, the existence of prominent stellar haloes for massive galaxies and could also arise from conceptual shortcomings of parametric 2D image decomposition tools. They consequently propagate into galaxy color profiles.

We study the impact of local density and stellar mass on the morphology of ~500 quiescent (Q) and SFGs from the VIS3COS survey. We perform B/D decomposition of the SBPs and find ~41% of 10^10 Msun galaxies to be best fitted with 2 components. We complement our analysis with non-parametric measurements and qualitative visual classifications. We find that galaxy morphology depends on stellar mass and environment for our sample as a whole. We only find an impact of the environment on galaxy sizes for galaxies more massive than 10^11 Msun. We find higher n and B/T in high-density regions when compared to low-density counterparts at similar stellar masses. We also find that higher stellar mass galaxies have steeper light profiles compared to the lower ones. Using visual classifications, we find a morphology-density relation at z~0.84 for galaxies more massive than 10^10 Msun, with elliptical galaxies being dominant at high-density regions and disks more common in low-density regions. However, when splitting the sample into colour-colour selected SF and Q sub-populations, there are no statistically significant differences between low- and high-density regions. We find that Q galaxies are smaller, have higher n, and higher B/T when compared to SF counterparts. We confirm these trends with non-parametric quantities, finding Q galaxies to be smoother and having most of their light over smaller areas than SFGs. Overall, we find that the dependence of morphology on stellar mass is stronger than with local density and these relations are strongly correlated with the quenching fraction. The change in average morphology corresponds to a change in the relative fractions of blue disk-like galaxies and red elliptical galaxies with stellar mass and environment. We hypothesize that the processes responsible for the quenching of SF must also affect the galaxy morphology on similar timescales.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

We study the morphological and structural properties of the host galaxies associated with 57 optically-selected luminous type 2 AGN at $z\sim$0.3-0.4: 16 high-luminosity Seyfert 2 (HLSy2, 8.0$\le$log($L_{\rm [OIII]}/L_{\odot})<$8.3) and 41 obscured quasars (QSO2, log($L_{\rm [OIII]}/L_{\odot})\ge$8.3). With this work, the total number of QSO2 at $z<1$ with parametrized galaxies increases from $\sim$35 to 76. Our analysis is based on HST WFPC2 and ACS images that we fit with {\sc GALFIT}. HLSy2 and QSO2 show a wide diversity of galaxy hosts. The main difference lies in the higher incidence of highly-disturbed systems among QSO2. This is consistent with a scenario in which galaxy interactions are the dominant mechanism triggering nuclear activity at the highest AGN power. There is a strong dependence of galaxy properties with AGN power (assuming $L_ {\rm [OIII]}$ is an adequate proxy). The relative contribution of the spheroidal component to the total galaxy light (B/T) increases with $L_ {\rm [OIII]}$. While systems dominated by the spheoridal component spread across the total range of $L_ {\rm [OIII]}$, most disk-dominated galaxies concentrate at log($L_{\rm [OIII]}/L_{\odot})<$8.6. This is expected if more powerful AGN are powered by more massive black holes which are hosted by more massive bulges or spheroids. The average galaxy sizes ($\langle r_{\rm e} \rangle$) are 5.0$\pm$1.5 kpc for HLSy2 and 3.9$\pm$0.6 kpc for HLSy2 and QSO2 respectively. These are significantly smaller than those found for QSO1 and narrow line radio galaxies at similar $z$. We put the results of our work in context of related studies of AGN with quasar-like luminosities.

Massive Ultracompact Galaxies (MUGs) are common at z=2-3, but very rare in the nearby Universe. Simulations predict that the few surviving MUGs should reside in galaxy clusters, whose large relative velocities prevent them from merging, thus maintaining their original properties (namely stellar populations, masses, sizes and dynamical state). We take advantage of the high-completeness, large-area spectroscopic GAMA survey, complementing it with deeper imaging from the KiDS and VIKING surveys. We find a set of 22 bona-fide MUGs, defined as having high stellar mass (>8x10^10 M_Sun) and compact size (R_e<2 Kpc) at 0.02 < z < 0.3. An additional set of 7 lower-mass objects (6x10^10 < M_star/M_Sun < 8x10^10) are also potential candidates according to typical mass uncertainties. The comoving number density of MUGs at low redshift (z < 0.3) is constrained at $(1.0\pm 0.4)x 10^-6 Mpc^-3, consistent with galaxy evolution models. However, we find a mixed distribution of old and young galaxies, with a quarter of the sample representing (old) relics. MUGs have a predominantly early/swollen disk morphology (Sersic index 1<n<2.5) with high stellar surface densities (<Sigma_e> ~ 10^10 M_Sun Kpc^-2). Interestingly, a large fraction feature close companions -- at least in projection -- suggesting that many (but not all) reside in the central regions of groups. Halo masses show these galaxies inhabit average-mass groups. As MUGs are found to be almost equally distributed among environments of different masses, their relative fraction is higher in more massive overdensities, matching the expectations that some of these galaxies fell in these regions at early times. However, there must be another channel leading some of these galaxies to an abnormally low merger history because our sample shows a number of objects that do not inhabit particularly dense environments. (abridged)

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

Using the exquisite depth of the Hubble Ultra Deep Field (HUDF12 programme) dataset, we explore the ongoing assembly of the outermost regions of the most massive galaxies ($\rm M_{\rm stellar}\geq$ 5$\times$10$^{10}$ M$_{\odot}$) at $z \leq$ 1. The outskirts of massive objects, particularly Early-Types Galaxies (ETGs), are expected to suffer a dramatic transformation across cosmic time due to continuous accretion of small galaxies. HUDF imaging allows us to study this process at intermediate redshifts in 6 massive galaxies, exploring the individual surface brightness profiles out to $\sim$25 effective radii. We find that 5-20\% of the total stellar mass for the galaxies in our sample is contained within 10 $< R <$ 50 kpc. These values are in close agreement with numerical simulations, and higher than those reported for local late-type galaxies ($\lesssim$5\%). The fraction of stellar mass stored in the outer envelopes/haloes of Massive Early-Type Galaxies increases with decreasing redshift, being 28.7\% at $< z > =$ 0.1, 15.1\% at $< z > =$ 0.65 and 3.5\% at $< z > =$ 2. The fraction of mass in diffuse features linked with ongoing minor merger events is $>$ 1-2\%, very similar to predictions based on observed close pair counts. Therefore, the results for our small albeit meaningful sample suggest that the size and mass growth of the most massive galaxies have been solely driven by minor and major merging from $z =$ 1 to today.

We present new results on the evolution of rest-frame blue/UV sizes and Sersic indices of H$\alpha$-selected star-forming galaxies over the last 11 Gyrs. We investigate how the perceived evolution can be affected by a range of biases and systematics such as cosmological dimming and resolution effects. We use GALFIT and an artificial redshifting technique, which includes the luminosity evolution of H$\alpha$-selected galaxies, to quantify the change on the measured structural parameters with redshift. We find typical sizes of 2 to 3 kpc and Sersic indices of n~1.2, close to pure exponential disks all the way from z=2.23 to z=0.4. At z=0 we find typical sizes of 4-5 kpc. Our results show that, when using GALFIT, cosmological dimming has a negligible impact on the derived effective radius for galaxies with <10 kpc, but we find a ~20% bias on the estimate of the median Sersic indices, rendering galaxies more disk-like. Star-forming galaxies have grown on average by a factor of 2-3 in the last 11 Gyrs with $r_e\propto(1+z)^{-0.75}$. By exploring the evolution of the stellar mass-size relation we find evidence for a stronger size evolution of the most massive star-forming galaxies since z~2, as they grow faster towards z~0 when compared to the lower stellar mass counterparts. As we are tracing the rest-frame blue/UV, we are likely witnessing the growth of disks where star formation is ongoing in galaxies while their profiles remain close to exponential disks, n<1.5, across the same period.

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$\sigma$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

We present multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been obtained using our PSF-matching deconfusion code T-PHOT. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. The procedure allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of a complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K band from Hawk-I, and 3.6 - 4.5 {\mu}m from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in a companion Paper II. [abridged]

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

We present a new method to search for candidate z~>2 Herschel 500{\mu}m sources in the GOODS-North field, using a S500{\mu}m/S24{\mu}m "color deconfusion" technique. Potential high-z sources are selected against low-redshift ones from their large 500{\mu}m to 24{\mu}m flux density ratios. By effectively reducing the contribution from low-redshift populations to the observed 500{\mu}m emission, we are able to identify counterparts to high-z 500{\mu}m sources whose 24{\mu}m fluxes are relatively faint. The recovery of known z~4 starbursts confirms the efficiency of this approach in selecting high-z Herschel sources. The resulting sample consists of 34 dusty star-forming galaxies at z~>2. The inferred infrared luminosities are in the range 1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates (SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA 850{\mu}m-selected galaxy samples shows that our method is more efficient at selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10 of the sources at z~>4. We find that at a fixed luminosity, the dust temperature is ~5K cooler than that expected from the Td-LIR relation at z<1, though different temperature selection effects should be taken into account. The radio-detected subsample (excluding three strong AGN) follows the far-infrared/radio correlation at lower redshifts, and no evolution with redshift is observed out to z~5, suggesting that the far-infrared emission is star formation dominated. The contribution of the high-z Herschel 500{\mu}m sources to the cosmic SFR density is comparable to that of SMG populations at z~2.5 and at least 40% of the extinction-corrected UV samples at z~4 (abridged).

We present T-PHOT, a publicly available software aimed at extracting accurate photometry from low-resolution images of deep extragalactic fields, where the blending of sources can be a serious problem for the accurate and unbiased measurement of fluxes and colours. T-PHOT has been developed within the ASTRODEEP project and it can be considered as the next generation to TFIT, providing significant improvements above it and other similar codes. T-PHOT gathers data from a high-resolution image of a region of the sky, and uses it to obtain priors for the photometric analysis of a lower resolution image of the same field. It can handle different types of datasets as input priors: i) a list of objects that will be used to obtain cutouts from the real high-resolution image; ii) a set of analytical models; iii) a list of unresolved, point-like sources, useful e.g. for far-infrared wavelength domains. We show that T-PHOT yields accurate estimations of fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be done thanks to a flagging code given in the output). T-PHOT is many times faster than similar codes like TFIT and CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more robust and more versatile. This makes it an optimal choice for the analysis of large datasets. In addition we show how the use of different settings and methods significantly enhances the performance. Given its versatility and robustness, T-PHOT can be considered the preferred choice for combined photometric analysis of current and forthcoming extragalactic optical to far-infrared imaging surveys. [abridged]

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 have constructed a mass-selected sample of Mstar>10^11Msolar galaxies at 1<z<3 in the CANDELS UDS and COSMOS fields and have decomposed these systems into their separate bulge and disk components according to their H(160)-band morphologies. By extending this analysis to multiple bands we have been able to conduct individual bulge and disk component SED fitting which has provided us with stellar-mass and star-formation rate estimates for the separate bulge and disk components. These have been combined with size measurements to explore the evolution of these massive high-redshift galaxies. By utilising the new decomposed stellar-mass estimates, we confirm that the bulge components display a stronger size evolution than the disks. This can be seen from both the fraction of bulge components which lie below the local relation and the median sizes of the bulge components, where the bulges are a median factor of 2.93+/-0.32 times smaller than similarly massive local galaxies at 1<z<2 and 3.41+/-0.58 smaller at 2<z<3; for the disks the corresponding factors are 1.65+/-0.14 and 1.99+/-0.25. Moreover, by splitting our sample into the passive and star-forming bulge and disk sub-populations and examining their sizes as a fraction of their present-day counter-parts, we find that the star-forming and passive bulges are equally compact, star-forming disks are larger, while the passive disks have intermediate sizes. This trend is not evident when classifying galaxy morphology on the basis of single-Sersic fits and adopting the overall star-formation rates. Finally, by evolving the star-formation histories of the passive disks back to the redshifts when the passive disks were last active, we show that the passive and star-forming disks have consistent sizes at the relevant epoch. These trends need to be reproduced by any mechanisms which attempt to explain the morphological evolution of galaxies.

We present the results of a new and improved study of the morphological and spectral evolution of massive galaxies over the redshift range 1<z<3. Our analysis is based on a bulge-disk decomposition of 396 galaxies with Mstar>10^11 Msolar from the CANDELS WFC3/IR imaging within the COSMOS and UKIDSS UDS survey fields. We find that, by modelling the H(160) image of each galaxy with a combination of a de Vaucouleurs bulge (Sersic index n=4) and an exponential disk (n=1), we can then lock all derived morphological parameters for the bulge and disk components, and successfully reproduce the shorter-wavelength J(125), i(814), v(606) HST images simply by floating the magnitudes of the two components. This then yields sub-divided 4-band HST photometry for the bulge and disk components which, with no additional priors, is well described by spectrophotometric models of galaxy evolution. Armed with this information we are able to properly determine the masses and star-formation rates for the bulge and disk components, and find that: i) from z=3 to z=1 the galaxies move from disk-dominated to increasingly bulge-dominated, but very few galaxies are pure bulges/ellipticals by z=1; ii) while most passive galaxies are bulge-dominated, and most star-forming galaxies disk-dominated, 18+/-5% of passive galaxies are disk-dominated, and 11+/-3% of star-forming galaxies are bulge-dominated, a result which needs to be explained by any model purporting to connect star-formation quenching with morphological transformations; iii) there exists a small but significant population of pure passive disks, which are generally flatter than their star-forming counterparts (whose axial ratio distribution peaks at b/a~0.7); iv) flatter/larger disks re-emerge at the highest star-formation rates, consistent with recent studies of sub-mm galaxies, and with the concept of a maximum surface-density for star-formation activity.

There is cumulative evidence showing that, for the most massive galaxies, the fraction of disk-like objects compared to those with spheroidal properties increases with redshift. However, this evolution is thus far based on the surface brightness study of these objects. To explore the consistency of this scenario, it is necessary to measure the dynamical status of these galaxies. With this aim we have obtained seeing-limited near-infrared integral field spectra in the H-band for 10 massive galaxies (M_{stellar} >10^{11} h_{70}^-2 M_{Sun}) at z~1.4 with SINFONI at the VLT. Our sample is selected by their stellar mass and EW[OII] > 15 \AA, to secure their kinematic measurements, but without accounting for any morphological or flux criteria a priori. Through this 3D kinematic spectroscopy analysis we find that half (i.e. 50+/-7%) of our galaxies are compatible with being rotationally supported disks, in agreement with previous morphological expectations. This is a factor of approximately two higher than what is observed in the present Universe for objects of the same stellar mass. Strikingly, the majority of our sample of massive galaxies show extended and fairly high rotational velocity maps, implying that massive galaxies acquire rapidly rotational support and hence gravitational equilibrium. Our sample also show evidence for ongoing interactions and mergers. Summarizing, massive galaxies at high-z show a significant diversity and must have continued evolution beyond the fading of stellar populations, to become their present day counterparts.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

In this paper we describe the first data release of the UltraVISTA near-infrared imaging survey of the COSMOS field. We summarise the key goals and design of the survey and provide a detailed description of our data reduction techniques . We provide stacked, sky-subtracted images in $YJHK_{\rm s}$ and narrow-band filters constructed from data collected during the first year of UltraVISTA observations. Our stacked images reach $5\sigma$ $AB$ depths in an aperture of $2\arcsec$ diameter of $\sim 25$ in $Y$ and $\sim 24$ in $JHK_{\rm s}$ bands and all have sub-arcsecond seeing. To this $5\sigma$ limit, our $K_{\rm s}$ catalogue contains 216,268 sources. We carry out a series of quality assessment tests on our images and catalogues, comparing our stacks with existing catalogues. The $1\sigma$ astrometric RMS in both directions for stars selected with $17.0<K_{\rm s}\rm {(AB)} <19.5$ is $\sim 0.08\arcsec$ in comparison to the publicly-available COSMOS ACS catalogues. Our images are resampled to the same pixel scale and tangent point as the publicly available COSMOS data and so may be easily used to generate multi-colour catalogues using this data. All images and catalogues presented in this paper are publicly available through ESO's "phase 3" archiving and distribution system and from the UltraVISTA web site.

We investigate the resolved star formation properties of a sample of 45 massive galaxies (M_*>10^11M_solar) within a redshift range of 1.5 < z < 3 detected in the GOODS NICMOS Survey (Conselice et al. 2011), a HST H-band imaging program. We derive the star formation rate as a function of radius using rest frame UV data from deep z_{850} ACS imaging. The star formation present at high redshift is then extrapolated to z=0, and we examine the stellar mass produced in individual regions within each galaxy. We also construct new stellar mass profiles of the in-situ stellar mass at high redshift from Sersic fits to rest-frame optical, H_{160}-band, data. We combine the two stellar mass profiles to produce a modelled evolved stellar mass profile. We then fit a new Sersic profile to the evolved profile, from which we examine what effect the resulting stellar mass distribution added via star formation has on the structure and size of each individual galaxy. We conclude that due to the lack of sufficient size growth and Sersic evolution by star formation other mechanisms such as merging must contribute a large proportion to account for the observed structural evolution from z>1 to the present day.

Present-day massive galaxies are composed mostly of early-type objects. It is unknown whether this was also the case at higher redshifts. In a hierarchical assembling scenario the morphological content of the massive population is expected to change with time from disk-like objects in the early Universe to spheroid-like galaxies at present. In this paper we have probed this theoretical expectation by compiling a large sample of massive (M_{stellar}>10^{11} h_{70}^{-2} M_{Sun}$) galaxies in the redshift interval 0 < z < 3. Our sample of 1082 objects comprises 207 local galaxies selected from SDSS plus 875 objects observed with the HST belonging to the POWIR/DEEP2 and GNS surveys. 639 of our objects have spectroscopic redshifts. Our morphological classification is performed as close as possible to the optical restframe according to the photometric bands available in our observations both quantitatively (using the Sersic index as a morphological proxy) and qualitative (by visual inspection). Using both techniques we find an enormous change on the dominant morphological class with cosmic time. The fraction of early-type galaxies among the massive galaxy population has changed from ~20-30% at z~3 to ~70% at z=0. Early type galaxies have been the predominant morphological class for massive galaxies since only z~1.

We report the results of a comprehensive study of the relationship between galaxy size, stellar mass and specific star-formation rate (sSFR) at redshifts 1.3<z<1.5. Based on a mass complete (M_star >= 6x10^10 Msun), spectroscopic sample from the UKIDSS Ultra-deep Survey (UDS), with accurate stellar-mass measurements derived from spectro photometric fitting, we find that at z~1.4 the location of massive galaxies on the size-mass plane is determined primarily by their sSFR. At this epoch we find that massive galaxies which are passive (sSFR <= 0.1 Gyr^-1) follow a tight size-mass relation, with half-light radii a factor f=2.4+/-0.2 smaller than their local counterparts. Moreover, amongst the passive sub-sample we find no evidence that the off-set from the local size-mass relation is a function of stellar population age. Based on a sub-sample with dynamical mass estimates we also derive an independent estimate of f=2.3+/-0.3 for the typical growth in half-light radius between z~1.4 and the present day. Focusing on the passive sub-sample, we conclude that to produce the necessary evolution predominantly via major mergers would require an unfeasible number of merger events and over populate the high-mass end of the local stellar mass function. In contrast, we find that a scenario in which mass accretion is dominated by minor mergers can produce the necessary evolution, whereby an increase in stellar mass by a factor of ~2, accompanied by an increase in size by a factor of ~3.5, is sufficient to reconcile the size-mass relation at z~1.4 with that observed locally. Finally, we note that a significant fraction (44+/-12%) of the passive galaxies in our sample have a disk-like morphology, providing additional evidence that separate physical processes are responsible for the quenching of star-formation and the morphological transformation of massive galaxies (abridged).

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 a study of the resolved star-forming properties of a sample of distant massive M_*>10^11M_solar galaxies in the GOODS NICMOS Survey (GNS). We derive dust corrected UV star formation rates (SFRs) as a function of radius for 45 massive galaxies within the redshift range 1.5<z<3 in order to measure the spatial location of ongoing star formation. We find that the star formation rates present in different regions of a galaxy reflect the already existent stellar mass density, i.e. high density regions have higher star formation rates than lower density regions, on average. This observed star formation is extrapolated in several ways to the present day, and we measure the amount of new stellar mass that is created in individual portions of each galaxy to determine how the stellar mass added via star formation changes the observed stellar mass profile, the Sersic index (n) and effective radius (R_e) over time. We find that these massive galaxies fall into three broad classifications of star formation distribution. These different star formation distributions increase the effective radii over time, which are on average a factor of ~16pm5% larger, with little change in n (average Delta n=-0.9pm0.9) after evolution. We also implement a range of simple stellar migration models into the simulated evolutionary path of these galaxies in order to gauge its effect on the properties of our sample. This yields a larger increase in the evolved R_e than the pure static star formation model, with a maximum average increase of Delta R_e~54pm19%, but with little change in n, Delta n ~-1.1pm1.3. These results are not in agreement with the observed change in the R_e and n between z~2.5 and 0 obtained via various observational studies. We conclude that star formation and stellar migration alone cannot account for the observed change in structural parameters for this galaxy population (abridged).

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

We investigate the total major (> 1:4 by stellar mass) and minor (> 1:100 by stellar mass) merger history of a population of 80 massive (M_* > 10^11 M_sol) galaxies at high redshifts (z = 1.7 - 3). We utilize extremely deep and high resolution HST H-band imaging from the GOODS NICMOS Survey (GNS), which corresponds to rest-frame optical wavelengths at the redshifts probed. We find that massive galaxies at high redshifts are often morphologically disturbed, with a CAS deduced merger fraction f_m = 0.23 +/- 0.05 at z = 1.7 - 3. We find close accord between close pair methods (within 30 kpc apertures) and CAS methods for deducing major merger fractions at all redshifts. We deduce the total (minor + major) merger history of massive galaxies with M_* > 10^9 M_sol galaxies, and find that this scales roughly linearly with log-stellar-mass and magnitude range. We test our close pair methods by utilizing mock galaxy catalogs from the Millennium Simulation. We compute the total number of mergers to be (4.5 +/- 2.9) / <\tau_m> from z = 3 to the present, to a stellar mass sensitivity threshold of ~ 1:100 (where \tau_m is the merger timescale in Gyr which varies as a function of mass). This corresponds to an average mass increase of (3.4 +/- 2.2) x 10^11 M_sol over the past 11.5 Gyrs due to merging. We show that the size evolution observed for these galaxies may be mostly explained by this merging.

We present a near-infrared spectroscopic study of a stellar mass selected sample of galaxies at z~1 utilising the LIRIS multi-object spectrograph on the WHT. We detect continuum, and the H alpha line for our sample, which is one of the better direct tracers of star formation in external galaxies. We spectroscopically measure the H alpha emission from 41 massive (M_{*}>10^{10.5} Msol) galaxies taken from the POWIR Survey with spectroscopic redshifts 0.4<z_{spec}<1.4. We correct our H alpha fluxes for dust extinction by using multi-wavelength data, and investigate SFR trends with mass and colour. We find a drop in the fraction of massive galaxies with M_{*}>10^{11} Msol which are detected in H alpha emission at z<0.9. We furthermore find that the fraction of galaxies with H alpha emission drops steadily and significantly with redder (U-B) colours at z~1, and that the SSFR drops with increasing (U-B) colour for galaxies at all masses. By investigating the SFR-mass relation we find that the SFR is roughly constant with mass, in possible contrast to previous work, and that the specific star formation rate (SSFR) is lower in the most massive galaxies. The scatter in the SFR vs. mass relationship is very small for those systems with ongoing star formation which suggests that star formation in the most massive galaxies at z~1 shuts off rather abruptly over <1 Gyr, without an obvious gradual decline. We furthermore investigate the SFR as a function of (U-B) colour divided into different mass bins, revealing a tracer of the epoch of transition from star forming to passive, as a form of star formation "downsizing". This suggests that the shut off of star formation occurs before the change in a galaxy's colour. We find that galaxy stellar mass is the primary driving mechanisms behind the star formation history for these galaxies and discuss several possible mechanisms for regulating this process.

We present our results on the structure and activity of massive galaxies at z=1-3 using one of the largest (166 with M_star>=5e10 M_sun) and most diverse samples of massive galaxies derived from the GOODS-NICMOS survey: (1) Sersic fits to deep NIC3/F160W images indicate that the rest-frame optical structures of massive galaxies are very different at z=2-3 compared to z~0. Approximately 40% of massive galaxies are ultra-compact (r_e<=2 kpc), compared to less than 1% at z~0. Furthermore, most (~65%) systems at z=2-3 have a low Sersic index n<=2, compared to ~13% at z~0. We present evidence that the n<=2 systems at z=2-3 likely contain prominent disks, unlike most massive z~0 systems. (2) There is a correlation between structure and star formation rates (SFR). The majority (~85%) of non-AGN massive galaxies at z=2-3, with SFR high enough to yield a 5 sigma (30 micro Jy) 24 micron Spitzer detection have low n<=2. Such n<=2 systems host the highest SFR. (3) The frequency of AGN is ~40% at z=2-3. Most (~65%) AGN hosts have disky (n<=2) morphologies. Ultra-compact galaxies appear quiescent in terms of both AGN activity and star formation. (4) Large stellar surface densities imply massive galaxies at z=2-3 formed via rapid, highly dissipative events at z>2. The large fraction of n<=2 disky systems suggests cold mode accretion complements gas-rich major mergers at z>2. In order for massive galaxies at z=2-3 to evolve into present-day massive E/S0s, they need to significantly increase (n, r_e). Dry minor and major mergers may play an important role in this process.

We investigate the relation between star formation rates and local galaxy environment for a stellar mass selected galaxy sample in the redshift range 1.5 < z < 3. We use near-infra-red imaging from an extremely deep Hubble Space Telescope survey, the GOODS-NICMOS Survey (GNS) to measure local galaxy densities based on the nearest neighbour approach, while star-formation rates are estimated from rest-frame UV-fluxes. Due to our imaging depth we can examine galaxies down to a colour-independent stellar mass completeness limit of log M\ast = 9.5 M\odot at z ~ 3. We find a strong dependence of star formation activity on galaxy stellar mass over the whole redshift range, which does not depend on local environment. The average star formation rates are largely independent of local environment apart from in the highest relative over-densities. Galaxies in over-densities of a factor of > 5 have on average lower star formation rates by a factor of 2 - 3, but only up to redshifts of z ~ 2. We do not see any evidence for AGN activity influencing these relations. We also investigate the influence of the very local environment on star-formation activity by counting neighbours within 30 kpc radius. This shows that galaxies with two or more close neighbours have on average significantly lower star formation rates as well as lower specific star formation rates up to z ~ 2.5. We suggest that this might be due to star formation quenching induced by galaxy merging processes.

We present a study of the star-forming properties of a stellar mass-selected sample of galaxies in the GOODS NICMOS Survey (GNS), based on deep Hubble Space Telescope imaging of the GOODS North and South fields. Using a stellar mass selected sample, combined with HST/ACS and Spitzer data to measure both UV and infrared derived star formation rates (SFR), we investigate the star forming properties of a complete sample of ~1300 galaxies down to log M*=9.5 at redshifts 1.5<z<3. Eight percent of the sample is made up of massive galaxies with M*>10^11 Msun. We derive optical colours, dust extinctions, and ultraviolet and infrared SFR to determine how the star formation rate changes as a function of both stellar mass and time. Our results show that SFR increases at higher stellar mass such that massive galaxies nearly double their stellar mass from star formation alone over the redshift range studied, but the average value of SFR for a given stellar mass remains constant over this 2 Gyr period. Furthermore, we find no strong evolution in the SFR for our sample as a function of mass over our redshift range of interest, in particular we do not find a decline in the SFR among massive galaxies, as is seen at z < 1. The most massive galaxies in our sample (log M*>11) have high average SFRs with values, SFR(UV,corr) = 103+/-75 Msun/yr, yet exhibit red rest-frame (U-B) colours at all redshifts. We conclude that the majority of these red high-redshift massive galaxies are red due to dust extinction. We find that A(2800) increases with stellar mass, and show that between 45% and 85% of massive galaxies harbour dusty star formation. These results show that even just a few Gyr after the first galaxies appear, there are strong relations between the global physical properties of galaxies, driven by stellar mass or another underlying feature of galaxies strongly related to the stellar mass.

We examine in this paper a stellar mass selected sample of galaxies at 1 < z < 3 within the Hubble Ultra Deep Field, utilising WFC3 imaging to study the rest-frame optical morphological distribution of galaxies at this epoch. We measure how apparent morphologies (disk, elliptical, peculiar) correlate with physical properties, such as quantitative structure and spectral-types. One primary result is that apparent morphology does not correlate strongly with stellar populations, nor with galaxy structure at this epoch, suggesting a chaotic formation history for Hubble types at z > 1. By using a locally defined definition of disk and elliptical galaxies based on structure and spectral-type, we find no true ellipticals at z > 2, and a fraction of 3.2+/-2.3% at 1.5 < z < 2. Local counterparts of disk galaxies are at a similar level of 7-10%, much lower than the 75% fraction at lower redshifts. We further compare WFC3 images with the rest-frame UV view of galaxies from ACS imaging, showing that galaxies imaged with ACS that appear peculiar often contain an `elliptical' like morphology in WFC3. We show through several simulations that this larger fraction of elliptical-like galaxies is partially due to the courser PSF of WFC3, and that the `elliptical' class very likely includes early-type disks. We also measure the merger history for our sample using CAS parameters, finding a redshift evolution increasing with redshift, and a peak merger fraction of ~30% at z~2 for the most massive galaxies with M_*> 10^{10} M_sol, consistent with previous results from ACS and NICMOS. We compare our results to semi-analytical model results and find a relatively good agreement between our morphological break-down and the predictions. Finally, we argue that the peculiars, ellipticals and peculiar ellipticals have similar properties, suggesting similar formation modes, likely driven by major mergers.

We use a sample of 8298 galaxies observed in the HST GOODS NICMOS Survey (GNS) to construct the galaxy stellar mass function as a function of both redshift and stellar mass up to z=3.5 and down to masses of Mstar=10^8.5 Msun at z~1. We discover that a significant fraction of all massive Mstar>10^11 Msun galaxies are in place up to the highest redshifts we probe, with a decreasing fraction of lower mass galaxies present at all redshifts. This is an example of `galaxy mass downsizing', and is the result of massive galaxies forming before lower mass ones, and not just simply ending their star formation earlier as in traditional downsizing scenarios. We find that the faint end slope is significantly steeper than what is found in previous investigations. We demonstrate that this steeper mass function better matches the stellar mass added due to star formation, thereby alleviating some of the mismatch between these two measures of the evolution of galaxy mass. We furthermore examine the stellar mass function divided into blue/red systems, as well as for star forming and non-star forming galaxies. We find a similar mass downsizing present for both blue/red and star-forming/non-star forming galaxies, and that the low mass galaxies are mostly all blue, and are therefore creating the steep mass functions. We furthermore show that, although there is a downsizing such that high mass galaxies are nearer their z=0 values at high redshift, this turns over at masses Mstar~10^10 Msun, such that the lowest mass galaxies are more common than galaxies at slight higher masses, creating a `dip' in the observed galaxy mass function. We argue that the galaxy assembly process may be driven by different mechanisms at low and high masses, and that the efficiency of the galaxy formation process is lowest at masses Mstar~10^10 Msun at 1<z<3. (Abridged)

We present a study of the relationship between galaxy colour, stellar mass, and local galaxy density in a deep near-infrared imaging survey up to a redshift of z~3 using the GOODS NICMOS Survey (GNS). The GNS is a very deep, near-infrared Hubble Space Telescope survey imaging a total of 45 arcmin^2 in the GOODS fields, reaching a stellar mass completeness limit of M* = 10^9.5 M_sun at z=3. Using this data we measure galaxy local densities based on galaxy counts within a fixed aperture, as well as the distance to the 3rd, 5th and 7th nearest neighbour. We find a strong correlation between colour and stellar mass at all redshifts up to z~3. We do not find a strong correlation between colour and local density, however, the highest overdensities might be populated by a higher fraction of blue galaxies than average or underdense areas, indicating a possible reversal of the colour-density relation at high redshift. Our data suggests that the possible higher blue fraction at extreme overdensities might be due to a lack of massive red galaxies at the highest local densities.

We study the relationship between galaxy colour, stellar mass, and local galaxy density in a deep near-infrared imaging survey up to a redshift of z~3 using the GOODS NICMOS Survey (GNS). The GNS is a deep near-infrared Hubble Space Telescope survey imaging a total of 45 square arcminutes of the GOODS fields, reaching a stellar mass completeness limit of log M* = 9.5 M_sun at z=3. Using this data we measure galaxy local densities based on galaxy counts within a fixed aperture, as well as the distance to the 3rd, 5th and 7th nearest neighbour. We compare the average rest-frame (U-B) colour and fraction of blue galaxies in different local densities and at different stellar masses. We find a strong correlation between colour and stellar mass at all redshifts up to z~3. Massive red galaxies are already in place at z~3 at the expected location of the red-sequence in the colour-magnitude diagram, although they are star forming. We do not find a strong correlation between colour and local density, however, there may be evidence that the highest overdensities are populated by a higher fraction of blue galaxies than average or underdense areas. This could indicating that the colour-density relation at high redshift is reversed with respect to lower redshifts (z<1), where higher densities are found to have lower blue fractions. Our data suggests that the possible higher blue fraction at extreme overdensities might be due to a lack of massive red galaxies at the highest local densities.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

The population of compact massive galaxies observed at z > 1 are hypothesised, both observationally and in simulations, to be merger remnants of gas-rich disc galaxies. To probe such a scenario we analyse a sample of 12 gas-rich and active star forming sub-mm galaxies (SMGs) at 1.8 < z < 3. We present a structural and size measurement analysis for all of these objects using very deep ACS and NICMOS imaging in the GOODS-North field. Our analysis reveals a heterogeneous mix of morphologies and sizes. We find that four galaxies (33% \pm 17%) show clear signs of mergers or interactions, which we classify as early-stage mergers. The remaining galaxies are divided into two categories: five of them (42% \pm 18%) are diffuse and regular disc-like objects, while three (25% \pm 14%) are very compact, spheroidal systems. We argue that these three categories can be accommodated into an evolutionary sequence, showing the transformation from isolated, gas-rich discs with typical sizes of 2-3 kpc, into compact (< 1 kpc) galaxies through violent major merger events, compatible with the scenario depicted by theoretical models. Our findings that some SMGs are already dense and compact provides strong support to the idea that SMGs are the precursors of the compact, massive galaxies found at slightly lower redshift.

We measure and analyse the sizes of 82 massive (M >= 10^11 M_Sun) galaxies at 1.7<z<3 utilizing deep HST NICMOS data taken in the GOODS North and South fields. Our sample provides the first statistical study of massive galaxy sizes at z>2. We split our sample into disk-like (Sersic index n<=2) and spheroid-like (Sersic index n>2) galaxies, and find that at a given stellar mass, disk-like galaxies at z~2.3 are a factor of 2.6+/-0.3 smaller than present day equal mass systems, and spheroid-like galaxies at the same redshift are 4.3+/-0.7 times smaller than comparatively massive elliptical galaxies today. We furthermore show that the stellar mass densities of very massive galaxies at z~2.5 are similar to present-day globular clusters with values ~2x10^10 M_Sun kpc^-3

We measure the sizes of 82 massive (M>10^11 M_sun) galaxies at 1.7<z<3 utilizing deep HST NICMOS data taken in the GOODS North and South fields. Our sample is almost an order of magnitude larger than previous studies at these redshifts, providing the first statistical study of massive galaxy sizes at z>2, confirming the extreme compactness of these galaxies. We split our sample into disk-like (n<2) and spheroid-like (n>2) galaxies based on their Sersic indices, and find that at a given stellar mass disk-like galaxies at z~2.3 are a factor of 2.6+/-0.3 smaller than present day equal mass systems, and spheroid-like galaxies at the same redshifts are 4.3+/-0.7 smaller than comparatively massive elliptical galaxies today. At z>2 our results are compatible with both a leveling off, or a mild evolution in size. Furthermore, the high density (~2x10^10 M_sun kpc^-3) of massive galaxies at these redshifts, which are similar to present day globular clusters, possibly makes any further evolution in sizes beyond z=3 unlikely.