The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose formation, nature and role in galaxy evolution are yet to be understood. In this work we first present a new automated method based on deep neural networks to detect clumps in galaxy images and show that it is more sensitive and faster than previous proposed methods. We then use this method to systematically detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA forward modeled zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair and direct comparisons between observations and simulations. After correcting for these effects with a simple mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. We find that between $20-40\%$ of $\sim10^{10}$ solar mass star-forming galaxies contain at least one massive ($M_c>10^{7}M_\odot$) off-centered clump but only $\sim2-5\%$ of the total galaxy stellar mass is in those clumps. We also show indications that the contribution of clumps to the stellar mass is more important in extended and low mass galaxies. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function when confronted under the same conditions although tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ through violent disk instabilities.

We report the bivariate HI- and H$_2$-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for HI, H$_2$, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the HI and H$_2$ conditional distributions and the galaxy stellar mass function, GSMF. For the conditional distributions we use the compilation presented in Calette et al. 2018. For determining the GSMF from $M_{\ast}\sim3\times10^{7}$ to $3\times10^{12}$ $M_{\odot}$, we combine two spectroscopic samples from the SDSS at the redshift range $0.0033<z<0.2$. We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $\alpha\approx-1.4$, consistent with previous determinations. The obtained HI MFs agree with radio blind surveys. Similarly, the H$_2$ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $\sim96$% of the HI and H$_2$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of H$_2$ and total neutral H in late-type galaxies of $\sim 1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting H$_2$ into stars and in transforming HI gas into H$_2$. Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas--stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods enable us to look for subtle morphological deviations which are typically missed by more traditional moment-based approaches. We use a generative model to learn a representation of expected data defined by the training set and then look for deviations from the learned representation by looking for the best reconstruction of a given object. In this first proof-of-concept work, we apply our method to two different test cases. We first show that from a set of simulated galaxies, we are able to detect $\sim90\%$ of merging galaxies if we train our network only with a sample of isolated ones. We then explore how the presented approach can be used to compare observations and hydrodynamic simulations by identifying observed galaxies not well represented in the models.

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings.

The new generation of deep photometric surveys requires unprecedentedly precise shape and photometry measurements of billions of galaxies to achieve their main science goals. At such depths, one major limiting factor is the blending of galaxies due to line-of-sight projection, with an expected fraction of blended galaxies of up to 50%. Current deblending approaches are in most cases either too slow or not accurate enough to reach the level of requirements. This work explores the use of deep neural networks to estimate the photometry of blended pairs of galaxies in monochrome space images, similar to the ones that will be delivered by the Euclid space telescope. Using a clean sample of isolated galaxies from the CANDELS survey, we artificially blend them and train two different network models to recover the photometry of the two galaxies. We show that our approach can recover the original photometry of the galaxies before being blended with $\sim$7% accuracy without any human intervention and without any assumption on the galaxy shape. This represents an improvement of at least a factor of 4 compared to the classical SExtractor approach. We also show that forcing the network to simultaneously estimate a binary segmentation map results in a slightly improved photometry. All data products and codes will be made public to ease the comparison with other approaches on a common data set.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

We analyze the optical morphologies of galaxies in the IllustrisTNG simulation at $z\sim0$ with a Convolutional Neural Network trained on visual morphologies in the Sloan Digital Sky Survey. We generate mock SDSS images of a mass complete sample of $\sim12,000$ galaxies in the simulation using the radiative transfer code SKIRT and include PSF and noise to match the SDSS r-band properties. The images are then processed through the exact same neural network used to estimate SDSS morphologies to classify simulated galaxies in four morphological classes (E, S0/a, Sab, Scd). The CNN model finds that $\sim95\%$ of the simulated galaxies fall in one the four main classes with high confidence. The mass-size relations of the simulated galaxies divided by morphological type also reproduce well the slope and the normalization of observed relations which confirms the realism of optical morphologies in the TNG suite. However, the Stellar Mass Functions decomposed into different morphologies still show significant discrepancies with observations both at the low and high mass end. We find that the high mass end of the SMF is dominated in TNG by massive disk galaxies while early-type galaxies dominate in the observations according to the CNN classifications. The present work highlights the importance of detailed comparisons between observations and simulations in comparable conditions.

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.

Deep Learning (DL) algorithms for morphological classification of galaxies have proven very successful, mimicking (or even improving) visual classifications. However, these algorithms rely on large training samples of labelled galaxies (typically thousands of them). A key question for using DL classifications in future Big Data surveys is how much of the knowledge acquired from an existing survey can be exported to a new dataset, i.e. if the features learned by the machines are meaningful for different data. We test the performance of DL models, trained with Sloan Digital Sky Survey (SDSS) data, on Dark Energy survey (DES) using images for a sample of $\sim$5000 galaxies with a similar redshift distribution to SDSS. Applying the models directly to DES data provides a reasonable global accuracy ($\sim$ 90%), but small completeness and purity values. A fast domain adaptation step, consisting in a further training with a small DES sample of galaxies ($\sim$500-300), is enough for obtaining an accuracy > 95% and a significant improvement in the completeness and purity values. This demonstrates that, once trained with a particular dataset, machines can quickly adapt to new instrument characteristics (e.g., PSF, seeing, depth), reducing by almost one order of magnitude the necessary training sample for morphological classification. Redshift evolution effects or significant depth differences are not taken into account in this study.

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.

We present visual-like morphologies over 16 photometric bands, from ultra-violet to near infrared, for 8,412 galaxies in the Cluster Lensing And Supernova survey with Hubble (CLASH) obtained by a convolutional neural network (CNN) model. Our model follows the CANDELS main morphological classification scheme, obtaining the probability for each galaxy at each CLASH band of being spheroid, disk, irregular, point source, or unclassifiable. Our catalog contains morphologies for each galaxy with Hmag < 24.5 in every filter where the galaxy is observed. We trained an initial CNN model using approximately 7,500 expert eyeball labels from The Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). We created eyeball labels for 100 randomly selected galaxies per each of the 16-filters set of CLASH (1,600 galaxy images in total), where each image was classified by at least five of us. We use these labels to fine-tune the network in order to accurately predict labels for the CLASH data and to evaluate the performance of our model. We achieve a root-mean-square error of 0.0991 on the test set. We show that our proposed fine-tuning technique reduces the number of labeled images needed for training, as compared to directly training over the CLASH data, and achieves a better performance. This approach is very useful to minimize eyeball labeling efforts when classifying unlabeled data from new surveys. This will become particularly useful for massive datasets such as the ones coming from near future surveys such as EUCLID or the LSST. Our catalog consists of prediction of probabilities for each galaxy by morphology in their different bands and is made publicly available at http://www.inf.udec.cl/~guille/data/Deep-CLASH.csv.

We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

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.

Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $\Upsilon_*$ gradients than if the IMF is held fixed. If $\Upsilon_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $\Upsilon_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $\Upsilon_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $\Upsilon_*$ gradients changes the high-mass slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $\Upsilon_*$ gradients in larger samples.

We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs). We use two visual classification catalogues, GZ2 and Nair & Abraham (2010), for training CNNs with colour images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training. The catalogue is publicly released with the paper.

Numerous ongoing and future large area surveys (e.g. DES, EUCLID, LSST, WFIRST), will increase by several orders of magnitude the volume of data that can be exploited for galaxy morphology studies. The full potential of these surveys can only be unlocked with the development of automated, fast and reliable analysis methods. In this paper we present DeepLeGATo, a new method for two-dimensional photometric galaxy profile modeling, based on convolutional neural networks. Our code is trained and validated on analytic profiles (HST/CANDELS F160W filter) and it is able to retrieve the full set of parameters of one- component S\'ersic models: total magnitude, effective radius, S\'ersic index, axis ratio. We show detailed comparisons between our code and GALFIT. On simulated data, our method is more accurate than GALFIT and 3000 time faster on GPU (50 times when run on the same CPU). On real data, DeepLeGATo trained on simulations behaves similarly to GALFIT on isolated galaxies. With a fast domain adaptation step made with the 0.1 - 0.8 per cent the size of the training set, our code is easily capable to reproduce the results obtained with GALFIT even on crowded regions. DeepLeGATo does not require any human intervention beyond the training step, rendering it much automated than traditional profiling methods. The development of this method for more complex models (two-component galaxies, variable PSF, dense sky regions) could constitute a fundamental tool in the era of big data in astronomy.

Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $\sigma_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm \alpha_{JAM}}$, from $M_{*}$ and $\sigma_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,\sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(\sigma_R/\sigma_e) = (R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on $M_*^{\rm \alpha_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm \alpha_{JAM}}/M_* - \sigma_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $\phi(M_*^{\rm \alpha_{JAM}})$ is an estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form.

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and S\'{e}rsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be S\'{e}rsic and the input S\'{e}rsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing S\'{e}rsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution.

We analyze the properties of a multiply-imaged Lyman-alpha (Lya) emitter at z=5.75 identified through SHARDS Frontier Fields intermediate-band imaging of the Hubble Frontier Fields (HFF) cluster Abell 370. The source, A370-L57, has low intrinsic luminosity (M_UV~-16.5), steep UV spectral index (\beta=-2.4+/-0.1), and extreme rest-frame equivalent width of Lya (EW(Lya)=420+180-120 \AA). Two different gravitational lens models predict high magnification (\mu~10--16) for the two detected counter-images, separated by 7", while a predicted third counter-image (\mu~3--4) is undetected. We find differences of ~50% in magnification between the two lens models, quantifying our current systematic uncertainties. Integral field spectroscopy of A370-L57 with MUSE shows a narrow (FWHM=204+/-10 km/s) and asymmetric Lya profile with an integrated luminosity L(Lya)~10^42 erg/s. The morphology in the HST bands comprises a compact clump (r_e<100 pc) that dominates the Lya and continuum emission and several fainter clumps at projected distances <1 kpc that coincide with an extension of the Lya emission in the SHARDS F823W17 and MUSE observations. The latter could be part of the same galaxy or an interacting companion. We find no evidence of contribution from AGN to the Lya emission. Fitting of the spectral energy distribution with stellar population models favors a very young (t<10 Myr), low mass (M*~10^6.5 Msun), and metal poor (Z<4x10^-3) stellar population. Its modest star formation rate (SFR~1.0 Msun/yr) implies high specific SFR (sSFR~2.5x10^-7 yr^-1) and SFR density (Sigma_SFR ~ 7-35 Msun/yr/kpc^2). The properties of A370-L57 make it a good representative of the population of galaxies responsible for cosmic reionization.

The study of galaxy morphology is fundamental to understand the physical processes driving the structural evolution of galaxies. It has long been known that dense environments host high fractions of early-type galaxies and low fractions of late-type galaxies, indicating that the environment affects the structural evolution of galaxies. In this paper we present an analysis of the morphological composition of red sequence galaxies in a sample of 9 galaxy clusters at $0.8<z<1.5$ drawn from the HAWK-I cluster survey (HCS), with the aim of investigating the evolutionary paths of galaxies with different morphologies. We classify galaxies according to their apparent bulge-to-total light ratio and compare with red sequence galaxies from the lower-redshift WINGS and EDisCS surveys. We find that, while the HCS red sequence is dominated by elliptical galaxies at all luminosities and stellar masses, the WINGS red sequence is dominated by elliptical galaxies only at its bright end ($M_V<-21.0$ mag), while S0s become the most frequent class at fainter luminosities. Disc-dominated galaxies comprise 10-14\% of the red sequence population in the low (WINGS) and high (HCS) redshift samples, although their fraction increases up to 40\% at $0.4 < z < 0.8$ (EDisCS). We find a 20\% increase in the fraction of S0 galaxies from $z \sim 1.5$ to $z \sim 0.05$ on the red sequence. These results suggest that elliptical and S0 galaxies follow different evolutionary histories and, in particular, that S0 galaxies result, at least at intermediate luminosities ($-22.0 < M_V < -20.0$), from the morphological transformation of quiescent spiral galaxies.

We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.

We set out to quantify the number density of quiescent massive compact galaxies at intermediate redshifts. We determine structural parameters based on i-band imaging using the CFHT equatorial SDSS Stripe 82 (CS82) survey (~170 sq. degrees) taking advantage of an exquisite median seeing of ~0.6''. We select compact massive (M > 5x10^10 M_sun) galaxies within the redshift range of 0.2<z<0.6. The large volume sampled allows to decrease the effect of cosmic variance that has hampered the calculation of the number density for this enigmatic population in many previous studies. We undertake an exhaustive analysis in an effort to untangle the various findings inherent to the diverse definition of compactness present in the literature. We find that the absolute number of compact galaxies is very dependent on the adopted definition and can change up to a factor of >10. We systematically measure a factor of ~5 more compacts at the same redshift than what was previously reported on smaller fields with HST imaging, which are more affected by cosmic variance. This means that the decrease in number density from z ~ 1.5 to z ~ 0.2 might be only of a factor of ~2-5, significantly smaller than what previously reported. This supports progenitor bias as the main contributor to the size evolution. This milder decrease is roughly compatible with the predictions from recent numerical simulations. Only the most extreme compact galaxies, with Reff < 1.5x( M/10^11 M_sun)^0.75 and M > 10^10.7 M_sun, appear to drop in number by a factor of ~20 and hence likely experience a noticeable size evolution.

The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared.

Establishing accurate morphological measurements of galaxies in a reasonable amount of time for future big-data surveys such as EUCLID, the Large Synoptic Survey Telescope or the Wide Field Infrared Survey Telescope is a challenge. Because of its high level of abstraction with little human intervention, deep learning appears to be a promising approach. Deep learning is a rapidly growing discipline that models high-level patterns in data as complex multilayered networks. In this work we test the ability of deep convolutional networks to provide parametric properties of Hubble Space Telescope like galaxies (half-light radii, Sersic indices, total flux etc..). We simulate a set of galaxies including point spread function and realistic noise from the CANDELS survey and try to recover the main galaxy parameters using deep-learning. We com- pare the results with the ones obtained with the commonly used profile fitting based software GALFIT. This way showing that with our method we obtain results at least equally good as the ones obtained with GALFIT but, once trained, with a factor 5 hundred time faster.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

We analyse the evolution of the red sequence in a sample of galaxy clusters at redshifts $0.8<z<1.5$ taken from the HAWK-I Cluster Survey (HCS). The comparison with the low-redshift ($0.04<z<0.08$) sample of the WIde-field Nearby Galaxy-cluster Survey (WINGS) and other literature results shows that the slope and intrinsic scatter of the cluster red sequence have undergone little evolution since $z=1.5$. We find that the luminous-to-faint ratio and the slope of the faint end of the luminosity distribution of the HCS red sequence are consistent with those measured in WINGS, implying that there is no deficit of red galaxies at magnitudes fainter than $M_V^*$ at high redshifts. We find that the most massive HCS clusters host a population of bright red sequence galaxies at $M_V < -22.0$ mag, which are not observed in low-mass clusters. Interestingly, we also note the presence of a population of very bright ($M_V < -23.0$ mag) and massive (${\log(M_*/M_\odot)} > 11.5$) red sequence galaxies in the WINGS clusters, which do not include only the brightest cluster galaxies and which are not present in the HCS clusters, suggesting that they formed at epochs later than $z=0.8$. The comparison with the luminosity distribution of a sample of passive red sequence galaxies drawn from the COSMOS/UltraVISTA field in the photometric redshift range $0.8<z_{phot}<1.5$ shows that the red sequence in clusters is more developed at the faint end, suggesting that halo mass plays an important role in setting the time-scales for the build-up of the red sequence.

We describe the luminosity function, based on Sersic fits to the light profiles, of CMASS galaxies at z ~ 0.55. Compared to previous estimates, our Sersic-based reductions imply more luminous, massive galaxies, consistent with the effects of Sersic- rather than Petrosian or de Vaucouleur-based photometry on the Sloan Digital Sky Survey (SDSS) main galaxy sample at z ~ 0.1. This implies a significant revision of the high mass end of the correlation between stellar and halo mass. Inferences about the evolution of the luminosity and stellar mass functions depend strongly on the assumed, and uncertain, k+e corrections. In turn, these depend on the assumed age of the population. Applying k+e corrections taken from fitting the models of Maraston et al. (2009) to the colors of both SDSS and CMASS galaxies, the evolution of the luminosity and stellar mass functions appears impressively passive, provided that the fits are required to return old ages. However, when matched in comoving number- or luminosity-density, the SDSS galaxies are less strongly clustered compared to their counterparts in CMASS. This rules out the passive evolution scenario, and, indeed, any minor merger scenarios which preserve the rank ordering in stellar mass of the population. Potential incompletenesses in the CMASS sample would further enhance this mismatch. Our analysis highlights the virtue of combining clustering measurements with number counts.

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.

We present a catalog of visual like H-band morphologies of $\sim50.000$ galaxies ($H_{f160w}<24.5$) in the 5 CANDELS fields (GOODS-N, GOODS-S, UDS, EGS and COSMOS). Morphologies are estimated with Convolutional Neural Networks (ConvNets). The median redshift of the sample is $<z>\sim1.25$. The algorithm is trained on GOODS-S for which visual classifications are publicly available and then applied to the other 4 fields. Following the CANDELS main morphology classification scheme, our model retrieves the probabilities for each galaxy of having a spheroid, a disk, presenting an irregularity, being compact or point source and being unclassifiable. ConvNets are able to predict the fractions of votes given a galaxy image with zero bias and $\sim10\%$ scatter. The fraction of miss-classifications is less than $1\%$. Our classification scheme represents a major improvement with respect to CAS (Concentration-Asymmetry-Smoothness)-based methods, which hit a $20-30\%$ contamination limit at high z. The catalog is released with the present paper via the $\href{http://rainbowx.fis.ucm.es/Rainbow_navigator_public}{Rainbow\,database}$

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We present a catalogue of 348 galaxy clusters and groups with $0.2<z<1.2$ selected in the 2.78 $deg^2$ ALHAMBRA Survey. The high precision of our photometric redshifts, close to $1\%$, and the wide spread of the seven ALHAMBRA pointings ensure that this catalogue has better mass sensitivity and is less affected by cosmic variance than comparable samples. The detection has been carried out with the Bayesian Cluster Finder (BCF), whose performance has been checked in ALHAMBRA-like light-cone mock catalogues. Great care has been taken to ensure that the observable properties of the mocks photometry accurately correspond to those of real catalogues. From our simulations, we expect to detect galaxy clusters and groups with both $70\%$ completeness and purity down to dark matter halo masses of $M_h\sim3\times10^{13}\rm M_{\odot}$ for $z<0.85$. Cluster redshifts are expected to be recovered with $\sim0.6\%$ precision for $z<1$. We also expect to measure cluster masses with $\sigma_{M_h|M^*_{CL}}\sim0.25-0.35\, dex$ precision down to $\sim3\times10^{13}\rm M_{\odot}$, masses which are $50\%$ smaller than those reached by similar work. We have compared these detections with previous optical, spectroscopic and X-rays work, finding an excellent agreement with the rates reported from the simulations. We have also explored the overall properties of these detections such as the presence of a colour-magnitude relation, the evolution of the photometric blue fraction and the clustering of these sources in the different ALHAMBRA fields. Despite the small numbers, we observe tentative evidence that, for a fixed stellar mass, the environment is playing a crucial role at lower redshifts (z$<$0.5).

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

The Next Generation Virgo Cluster Survey is an optical imaging survey covering 104 deg^2 centered on the Virgo cluster. Currently, the complete survey area has been observed in the u*giz-bands and one third in the r-band. We present the photometric redshift estimation for the NGVS background sources. After a dedicated data reduction, we perform accurate photometry, with special attention to precise color measurements through point spread function-homogenization. We then estimate the photometric redshifts with the Le Phare and BPZ codes. We add a new prior which extends to iAB = 12.5 mag. When using the u*griz-bands, our photometric redshifts for 15.5 \le i \lesssim 23 mag or zphot \lesssim 1 galaxies have a bias |\Delta z| < 0.02, less than 5% outliers, and a scatter \sigma_{outl.rej.} and an individual error on zphot that increase with magnitude (from 0.02 to 0.05 and from 0.03 to 0.10, respectively). When using the u*giz-bands over the same magnitude and redshift range, the lack of the r-band increases the uncertainties in the 0.3 \lesssim zphot \lesssim 0.8 range (-0.05 < \Delta z < -0.02, \sigma_{outl.rej} ~ 0.06, 10-15% outliers, and zphot.err. ~ 0.15). We also present a joint analysis of the photometric redshift accuracy as a function of redshift and magnitude. We assess the quality of our photometric redshifts by comparison to spectroscopic samples and by verifying that the angular auto- and cross-correlation function w(\theta) of the entire NGVS photometric redshift sample across redshift bins is in agreement with the expectations.

We study of the role of major mergers (mass ratios >1:4) in driving size growth in high-redshift (1<z<2) spheroidal galaxies (SGs) with stellar masses between 10^9.5 MSun and 10^10.7 MSun. This is a largely unexplored mass range at this epoch, containing the progenitors of more massive SGs on which the bulk of the size-evolution literature is based. We visually split our SGs into systems that are relaxed and those that exhibit tidal features indicative of a recent merger. Numerical simulations indicate that, given the depth of our images, only tidal features due to major mergers will be detectable at these epochs (features from minor mergers being too faint), making the disturbed SGs a useful route for estimating major-merger-driven size growth. The disturbed SGs are offset in size from their relaxed counterparts, lying close to the upper envelope of the local size -- mass relation. The mean size ratio of the disturbed SGs to their relaxed counterparts is ~2. Combining this observed size growth with empirical major-merger histories from the literature suggests that the size evolution of a significant fraction (around two-thirds) of SGs in this mass range could be driven by major mergers. If, as is likely, our galaxies are progenitors of more massive M > 10^10.7 MSun SGs at z<1, then major mergers are also likely to play an important role in the size growth of at least some massive SGs in this mass range.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

We quantify the systematics in the size-luminosity relation of galaxies in the SDSS main sample which arise from fitting different 1- and 2-component model profiles to the images. In objects brighter than L*, fitting a single Sersic profile to what is really a two-component SerExp system leads to biases: the half-light radius is increasingly overestimated as n of the fitted single component increases; it is also overestimated at B/T ~ 0.6. However, the net effect on the R-L relation is small, except for the most luminous tail, where it curves upwards towards larger sizes. We also study how this relation depends on morphological type. Our analysis is one of the first to use Bayesian-classifier derived weights, rather than hard cuts, to define morphology. Crudely, there appear to be only two relations: one for early-types (Es, S0s and Sa's) and another for late-types (Sbs and Scds). However, closer inspection shows that within the early-type sample S0s tend to be 15% smaller than Es of the same luminosity, and, among faint late types, Sbs are more than 25% smaller than Scds. Neither the early- nor the late-type relations are pure power-laws: both show significant curvature, which we quantify. However, the R-L relations of the bulges of early-types are almost pure power laws. Our analysis confirms that two mass scales are special for early-type galaxies: M* = 3e10 and 2e11 Msun. These same mass scales are also special for late types: there is almost no correlation between R and M* below the former, and almost no late-types above the latter. We also show that the intrinsic scatter around the relation decreases at large luminosity and/or stellar mass; this should provide additional constraints on models of how the most massive galaxies formed.

[abridged] New near-infrared surveys, using the HST, offer an unprecedented opportunity to study rest-frame optical galaxy morphologies at z>1 and to calibrate automated morphological parameters that will play a key role in classifying future massive datasets like EUCLID or LSST. We study automated parameters (e.g. CAS, Gini, M20) of massive galaxies at 1<z<3, measure their dependence on wavelength and evolution with redshift and quantify the reliability of these parameters in discriminating between visually-determined morphologies, using machine learning algorithms. We find that the relative trends between morphological types observed in the low-redshift literature are preserved at z>1: bulge-dominated systems have systematically higher concentration and Gini coefficients and are less asymmetric and rounder than disk-dominated galaxies. However, at z>1, galaxies are, on average, 50% more asymmetric and have Gini and M20 values that are 10% higher and 20% lower respectively. In bulge-dominated galaxies, morphological parameters derived from the rest-frame UV and optical wavelengths are well correlated; however late-type galaxies exhibit higher asymmetry and clumpiness when measured in the rest-frame UV. We find that broad morphological classes (e.g. bulge vs. disk dominated) can be distinguished using parameters with high (80%) purity and completeness of 80%. In a similar vein, irregular disks and mergers can also be distinguished from bulges and regular disks with a contamination lower than 20%. However, mergers cannot be differentiated from the irregular morphological class using these parameters, due to increasingly asymmetry of non-interacting late-type galaxies at z>1. Our automated procedure is applied to the CANDELS GOODS-S field and compared with the visual classification recently released on the same area getting similar results.

The ALHAMBRA (Advance Large Homogeneous Area Medium Band Redshift Astronomical) survey has observed 8 different regions of the sky, including sections of the COSMOS, DEEP2, ELAIS, GOODS-N, SDSS and Groth fields using a new photometric system with 20 contiguous ~ $300\AA$ filters covering the optical range, combining them with deep $JHKs$ imaging. The observations, carried out with the Calar Alto 3.5m telescope using the wide field (0.25 sq. deg FOV) optical camera LAICA and the NIR instrument Omega-2000, correspond to ~700hrs on-target science images. The photometric system was designed to maximize the effective depth of the survey in terms of accurate spectral-type and photo-zs estimation along with the capability of identification of relatively faint emission lines. Here we present multicolor photometry and photo-zs for ~438k galaxies, detected in synthetic F814W images, complete down to I~24.5 AB, taking into account realistic noise estimates, and correcting by PSF and aperture effects with the ColorPro software. The photometric ZP have been calibrated using stellar transformation equations and refined internally, using a new technique based on the highly robust photometric redshifts measured for emission line galaxies. We calculate photometric redshifts with the BPZ2 code, which includes new empirically calibrated templates and priors. Our photo-zs have a precision of $dz/(1+z_s)=1%$ for I<22.5 and 1.4% for 22.5<I<24.5. Precisions of less than 0.5% are reached for the brighter spectroscopic sample, showing the potential of medium-band photometric surveys. The global $P(z)$ shows a mean redshift <z>=0.56 for I<22.5 AB and <z>=0.86 for I<24.5 AB. The data presented here covers an effective area of 2.79 sq. deg, split into 14 strips of 58.5'x15.5' and represents ~32 hrs of on-target.

We show that recently documented trends in galaxy sizes with mass and redshift can be understood in terms of the influence of underlying cosmic evolution; a holistic view which is complimentary to interpretations involving the accumulation of discreet evolutionary processes acting on individual objects. Using standard cosmology theory, supported with results from the Millennium simulations, we derive expected size trends for collapsed cosmic structures, emphasising the important distinction between these trends and the assembly paths of individual regions. We then argue that the observed variation in the stellar mass content of these structures can be understood to first order in terms of natural limitations of cooling and feedback. But whilst these relative masses vary by orders of magnitude, galaxy and host radii have been found to correlate linearly. We explain how these two aspects will lead to galaxy sizes that closely follow observed trends and their evolution, comparing directly with the COSMOS and SDSS surveys. Thus we conclude that the observed minimum radius for galaxies, the evolving trend in size as a function of mass for intermediate systems, and the observed increase in the sizes of massive galaxies, may all be considered an emergent consequence of the cosmic expansion.

Our goal is to estimate empirically, for the first time, the cosmic variance that affects merger fraction studies based on close pairs. We compute the merger fraction from photometric redshift close pairs with 10h^-1 kpc <= rp <= 50h^-1 kpc and Dv <= 500 km/s, and measure it in the 48 sub-fields of the ALHAMBRA survey. We study the distribution of the measured merger fractions, that follow a log-normal function, and estimate the cosmic variance sigma_v as the intrinsic dispersion of the observed distribution. We develop a maximum likelihood estimator to measure a reliable sigma_v and avoid the dispersion due to the observational errors (including the Poisson shot noise term). The cosmic variance of the merger fraction depends mainly on (i) the number density of the populations under study, both for the principal (n_1) and the companion (n_2) galaxy in the close pair, and (ii) the probed cosmic volume V_c. We find a significant dependence on neither the search radius used to define close companions, the redshift, nor the physical selection (luminosity or stellar mass) of the samples. We provide a parametrisation of the cosmic variance with n_1, n_2, and V_c, sigma_v = 0.48 n_1^{-0.54} V_c^{-0.48} (n_2/n_1)^{-0.37}. Thanks to this prescription, future merger fraction studies based on close pairs could account properly for the cosmic variance on their results.

We present the results of a detailed analysis of galaxy properties along the red sequence in XMMU J1229+0151, an X-ray selected cluster at $z=0.98$ drawn from the HAWK-I Cluster Survey (HCS). Taking advantage of the broad photometric coverage and the availability of 77 spectra in the cluster field, we fit synthetic spectral energy distributions, and estimate stellar masses and photometric redshifts, which we use to determine the cluster membership. We investigate morphological and structural properties of red sequence galaxies and find that elliptical galaxies populate the bright end, while S0 galaxies represent the predominant population at intermediate luminosities, with their fraction decreasing at fainter magnitudes. A comparison with the low-redshift sample of the WINGS cluster survey reveals that at $z\sim1$ the bright end of the red sequence of XMMU J1229+0151 is richer in S0 galaxies. The faint end of the red sequence in XMMUJ1229+0151 appears rich in disc-dominated galaxies, which are rarer in the low redshift comparison sample at the same luminosities. Despite these differences between the morphological composition of the red sequence in XMMUJ1229+0151 and in low redshift samples, we find that to within the uncertainties, no such difference exists in the ratio of luminous to faint galaxies along the red sequence.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5\sigma limiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with \mu_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.

The early--type galaxy (ETG) mass--size relation has been largely studied to understand how these galaxies have assembled their mass. One key observational result of the last years is that massive galaxies increased their size by a factor of a few at fixed stellar mass from $z\sim2$. Hierarchical models favor minor mergers as a plausible driver of this size growth. Some of these models, predict a significant environmental dependence in the sense that galaxies residing in more massive halos tend to be larger than galaxies in lower mass halos, at fixed stellar mass and redshift. At present, observational results of this environmental dependence have been contradictory. In this paper we revisit this issue in the local Universe, by investigating how the sizes of massive ETGs depend on large-scale environment using an updated and accurate sample of ETGs in different environments - field, group, clusters - from the Sloan Digital Sky Survey DR7. Our analysis does not show any significant environmental dependence of the sizes of central and satellites ETGs at fixed stellar mass at $z\sim0$. The size-mass relation of early-type galaxies at $z\sim0$ seems to be universal, i.e., independent of the mass of the host halo and of the position of the galaxy in that halo (central or satellite). The result is robust to different galaxy selections based on star formation, morphology or central density. In fact, considering our observational errors and the size of the sample, any size ratio larger than $30-40%$ between massive galaxies ($log(M_*/M_\odot)>11$) living in clusters and in the field can be ruled out at $3\sigma$ level.

In addition to the large systematic differences arising from assumptions about the stellar mass-to-light ratio, the massive end of the stellar mass function is rather sensitive to how one fits the light profiles of the most luminous galaxies. We quantify this by comparing the luminosity and stellar mass functions based on SDSS cmodel magnitudes, and PyMorph single-Sersic and Sersic-Exponential fits to the surface brightness profiles of galaxies in the SDSS. The PyMorph fits return more light, so that the predicted masses are larger than when cmodel magnitudes are used. As a result, the total stellar mass density at z~0.1 is about 1.2x larger than in our previous analysis of the SDSS. The differences are most pronounced at the massive end, where the measured number density of objects having M* > 6 x 10^{11} Msun is ~5x larger. Alternatively, at number densities of 10^{-6} Mpc^{-3}, the limiting stellar mass is 2x larger. The differences with respect to fits by other authors, typically based on Petrosian-like magnitudes, are even more dramatic, although some of these differences are due to sky-subtraction problems, and are sometimes masked by large differences in the assumed $M_*/L$ (even after scaling to the same IMF). Our results impact studies of the growth and assembly of stellar mass in galaxies, and of the relation between stellar and halo mass, so we provide simple analytic fits to these new luminosity and stellar mass functions and quantify how they depend on morphology, as well as the binned counts in electronic format.

ALHAMBRA is a photometric survey designed to trace the cosmic evolution and cosmic variance. It covers a large area of ~ 4 sq. deg in 8 fields, where 7 fields overlap with other surveys, allowing to have complementary data in other wavelengths. All observations were carried out in 20 continuous, medium band (30 nm width) optical and 3 near-infrared (JHK) bands, providing the precise measurements of photometric redshifts. In addition, morphological classification of galaxies is crucial for any kind of galaxy formation and cosmic evolution studies, providing the information about star formation histories, their environment and interactions, internal perturbations, etc. We present a morphological classification of > 40,000 galaxies in the ALHAMBRA survey. We associate to every galaxy a probability to be early-type using the automated Bayesian code galSVM. Despite of the spatial resolution of the ALHAMBRA images (~ 1 arcsec), for 22,051 galaxies we obtained the contamination by other type of less than 10%. Of those, 1,640 and 10,322 galaxies are classified as early- (down to redshifts ~ 0.5) and late-type (down to redshifts ~ 1.0), respectively, with magnitudes F613W < 22.0. In addition, for magnitude range 22.0 < F613W < 23.0 we classified other 10,089 late-type galaxies with redshifts < 1.3. We show that the classified objects populate the expected regions in the colour-mass and colour-magnitude planes. The presented dataset is especially attractive given the homogeneous multi-wavelength coverage available in the ALHAMBRA fields, and is intended to be used in a variety of scientific applications. The low-contamination catalogue (< 10%) is made publicly available with the present paper.

[abridged] The mass-size relation of early-type galaxies (ETGs) has been largely studied in the last years to probe the mass assembly of the most massive objects in the Universe. In this paper, we focus on the mass-size relation of quiescent massive ETGs (Mstar/Msol > 3*10^10) living in massive clusters (M200 ~ 10^14 Mstar) at 0.8< z <1.5, as compared to those living in the field at the same epoch. Our sample contains ~ 400 ETGs in clusters and the same number in the field. Therefore, our sample is approximately an order of magnitude larger than previous studies in the same redshift range for galaxy clusters. We find that ETGs living in clusters are between ~30-50% larger than galaxies with the same stellar mass residing in the field. We parametrize the size using the mass-normalized size, gamma=Re/Mstar^0.57. The gamma distributions in both environments peak at the same position but the distributions in clusters are more skewed towards larger sizes. Since this size difference is not observed in the local Universe, the size evolution at fixed stellar mass from z~1.5 of cluster galaxies is less steep ((1+z)-0.53pm0.04) than the evolution of field galaxies ((1+z)-0.92pm0.04). The size difference seems to be essentially driven by the galaxies residing in the clusters cores (R<0.5*R200). If part of the size evolution is due to mergers, the difference we see between cluster and field galaxies could be due to higher merger rates in clusters at higher redshift, probably during the formation phase of the clusters when velocity dispersions are lower. We cannot exclude however that the difference is driven by newly quenched galaxies which are quenched more efficiently in clusters. The implications of these results for the hierarchical growth of ETGs will be discussed in a companion paper.

Recent independent results from numerical simulations and observations have shown that brightest cluster galaxies (BCGs) have increased their stellar mass by a factor of almost two between z~0.9 and z~0.2. The numerical simulations further suggest that more than half this mass is accreted through major mergers. Using a sample of 18 distant galaxy clusters with over 600 spectroscopically confirmed cluster members between them, we search for observational evidence that major mergers do play a significant role. We find a major merger rate of 0.38 +/- 0.14 mergers per Gyr at z~1. While the uncertainties, which stem from the small size of our sample, are relatively large, our rate is consistent with the results that are derived from numerical simulations. If we assume that this rate continues to the present day and that half of the mass of the companion is accreted onto the BCG during these mergers, then we find that this rate can explain the growth in the stellar mass of the BCGs that is observed and predicted by simulations. Major mergers therefore appear to be playing an important role, perhaps even the dominant one, in the build up of stellar mass in these extraordinary galaxies.

(Abridged) AVOCADO aims at establishing firm conclusions on the formation and evolution of dwarf galaxies by constructing and analysing a homogeneous, multiwavelength dataset for a statistically significant sample of approximately 6500 nearby dwarfs (M_i > -18 mag). We present the sample selection criteria and describe the suite of analysis tools, some of them developed in the framework of the Virtual Observatory. We make use of optical spectra and UV-to-NIR imaging of the dwarf sample to derive SFRs, stellar masses, ages and metallicities - which are further supplemented with structural parameters that are used to classify them morphologically. This dataset, coupled with a detailed characterisation of each dwarf's environment, allows for a fully comprehensive investigation of their origins and to track the potential evolutionary paths between the different dwarf types. We characterise the local environment of all dwarfs in our sample, paying special attention to trends with current star formation activity. We find that virtually all quiescent dwarfs are located in the vicinity (projected distances < 1.5 Mpc) of L* companions, consistent with recent results. While star-forming dwarfs are preferentially found at separations of the order of 1 Mpc, there appears to be a tail towards low separations (<100 kpc) in the distribution of projected distances. We speculate that, modulo projection effects, this probably represents a genuine population of late-type dwarfs caught upon first infall about their host and before environmental quenching has fully operated. In this context, these results suggest that internal mechanisms -such as gas exhaustion via star formation or feedback effects- are not sufficient to completely cease the star formation activity in dwarf galaxies, and that becoming the satellite of a massive central galaxy appears to be a necessary condition to create a quiescent dwarf.

Dwarf galaxies are by far the most abundant of all galaxy types, yet their properties are still poorly understood -especially due to the observational challenge that their intrinsic faintness represents. AVOCADO aims at establishing firm conclusions on their formation and evolution by constructing a homogeneous, multiwavelength dataset for a statistically significant sample of several thousand nearby dwarfs (-18 < Mi < -14). Using public data and Virtual Observatory tools, we have built GALEX+SDSS+2MASS spectral energy distributions that are fitted by a library of single stellar population models. Star formation rates, stellar masses, ages and metallicities are further complemented with structural parameters that can be used to classify them morphologically. This unique dataset, coupled with a detailed characterization of each dwar's environment, allows for a fully comprehensive investigation of their origins and to track the (potential) evolutionary paths between the different dwarf types.

[abridged] We study the dependence of the galaxy size evolution on morphology, stellar mass and large scale environment for a sample of 298 group and 384 field quiescent early-type galaxies from the COSMOS survey, selected from z~1 to the present, and with masses $log(M/M_\odot)>10.5$. The galaxy size growth depends on galaxy mass and early-type galaxy morphology, e.g., elliptical galaxies evolve differently than lenticular galaxies. At the low mass end -$10.5<Log(M/M_\odot)<11$, ellipticals do not show strong size growth from $z\sim1$ to the present (10% to 30% depending on the morphological classification). On the other end, massive ellipticals -log(M/M_\odot)>11.2$- approximately doubled their size. Interestingly, lenticular galaxies display different behavior: they appear more compact on average and they do show a size growth of \sim60% since z=1 independent of stellar mass. We compare our results with state-of-the art semi-analytic models. While major and minor mergers can account for most of the galaxy size growth, we find that with present data and the theoretical uncertainties in the modeling we cannot state clear evidence favoring either merger or mass loss via quasar and/or stellar winds as the primary mechanism driving the evolution. The galaxy mass--size relation and the size growth do not depend on environment in the halo mass range explored in this work (field to group mass $log(M_h/M_\odot)<14$), i.e., group and field galaxies follow the same trends, which is at variance with predictions from current hierarchical models that show a clear dependence of size growth on halo mass for massive ellipticals -$log(M_*/M_\odot)>11.2$.

We confirm the detection of 3 groups in the Lynx supercluster, at z~1.3, and give their redshifts and masses. We study the properties of the group galaxies as compared to the central clusters, RXJ0849+4452 and RXJ0848+4453, selecting 89 galaxies in the clusters and 74 galaxies in the groups. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30% -40% level by simple automated classification methods (e.g. based on Sersic index). In luminosity selected samples, both clusters and groups show high fractions of Sa galaxies. The ETG fractions never rise above ~50% in the clusters, which is low compared to the fractions observed in clusters at z~1. However, ETG plus Sa fractions are similar to those observed for ETGs in clusters at z~1. Bulge-dominated galaxies visually classified as Sas might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples, the ETG fraction show no significant evolution with respect to local clusters, suggesting that morphological transformations occur at lower masses and densities. The ETG mass-size relation shows evolution towards smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. The group ETG red sequence shows lower zero points and larger scatters than in clusters, both expected to be an indication of a younger galaxy population. The estimated age difference is small when compared to the difference in age at different galaxy masses.

The connection between AGN and their host galaxies has been widely studied over recent years, showing it to be of great importance for providing answers to some fundamental questions related with AGN fueling mechanisms, their formation and evolution. Using X-ray and one of the deepest broad-band optical data sets, we studied morphology and colours in relationship with X-ray properties for sources at redshifts z < 2.0, using a sample of 262 AGN in the Subaru/XMM-Newton Deep Survey (SXDS). Morphological classification was obtained using the galSVM code, one of the new methods useful especially when dealing with high-redshift sources and low-resolution data. Colour-magnitude diagrams were studied in relationship with redshift, morphology, X-ray obscuration, and X-ray-to-optical flux ratio. Finally, the significance of different regions was analysed on colour-magnitude diagrams, relating the observed properties of AGN populations with some models of their formation and evolution.

We revisit the scaling relations and star-forming histories of local elliptical galaxies using a novel selection method applied to the Sloan Digital Sky Survey DR7. We combine two probability-based automated spectroscopic and morphological classifications of about 600000 galaxies with z<0.25 to isolate true elliptical galaxies. Our sample selection method does not introduce artificial cuts in the parameters describing the galaxy but instead it associates to every object a weight measuring the probability of being in a given spectro-morphological class. Thus the sample minimizes the selection biases. We show that morphologically defined ellipticals are basically distributed in 3 spectral classes, which dominate at different stellar masses. The bulk of the population (about 50%) is formed by a well defined class of galaxies with old stellar populations that formed their stars at very early epochs in a short episode of star formation. They dominate the scaling relations of elliptical galaxies known from previous works and represent the canonical elliptical class. At the low mass end, we find a population of slightly larger ellipticals, with smaller velocity dispersions at fixed stellar mass, which seem to have experienced a more recent episode of star formation probably triggered by gas-rich minor mergers. The high mass end tends to be dominated by a third spectral class, slightly more metal rich and with more efficient stellar formation than the reference class. This third class contributes to the curvature of the mass-size relation at high masses reported in previous works. Our method is therefore able to isolate typical spectra of elliptical galaxies following different evolutive pathways.

We present the Kormendy and mass-size relations for early-type galaxies (ETGs) as a function of environment at z~1.3. Our sample includes 76 visually classified ETGs with masses 10^10 < M/Msun < 10^11.5, selected in the Lynx supercluster and in the GOODS/CDF-S field, 31 ETGs in clusters, 18 in groups and 27 in the field, all with multi-wavelength photometry and HST/ACS observations. The Kormendy relation, in place at z~1.3, does not depend on the environment. The mass-size relation reveals that ETGs overall appear to be more compact in denser environments: cluster ETGs have sizes on average around 30-50% smaller than those of the local universe, and a distribution with a smaller scatter, whereas field ETGs show a mass-size relation with a similar distribution than the local one. Our results imply that (1) the mass-size relation in the field did not evolve overall from z ~ 1.3 to present; this is interesting and in contrast to the trend found at higher masses from previous works; (2) in denser environments, either ETGs have increased their size by 30-50%, on average, and spread their distributions, or more ETGs have been formed within the dense environment from not ETG progenitors or larger galaxies have been accreted to a pristine compact population to reproduce the mass-size relation observed in the local Universe. Our results are driven by galaxies with masses M<2*10^11Msun and those with masses M~10^11Msun follow the same trends that the entire sample. Following Valentinuzzi et al. definition of superdense ETGs, around 35-45% of our cluster sample is made of superdense ETGs.

We compare the Hubble type and the spectroscopic class of the galaxies with spectra in SDSS/DR7. As it is long known, elliptical galaxies tend to be red whereas spiral galaxies tend to be blue, however, this relationship presents a large scatter, which we measure and quantify in detail. We compare the Automatic Spectroscopic K-means based classification (ASK) with most of the commonly used morphological classifications. All of them provide consistent results. Given a spectral class, the morphological type wavers with a standard deviation between 2 and 3 T types, and the same large dispersion characterizes the variability of spectral classes fixed the morphological type. The distributions of Hubble types given an ASK class are very skewed -- they present long tails that go to the late morphological types for the red galaxies, and to the early morphological types for the blue spectroscopic classes. The scatter is not produced by problems in the classification, and it remains when particular subsets are considered. A considerable fraction of the red galaxies are spirals (40--60 %), but they never present very late Hubble types (Sd or later). Even though red spectra are not associated with ellipticals, most ellipticals do have red spectra: 97 % of the ellipticals in the morphological catalog by Nair & Abraham, used here for reference, belong to ASK 0, 2 or 3. It contains only a 3 % of blue ellipticals. The galaxies in the green valley class (ASK~5) are mostly spirals, and the AGN class (ASK 6) presents a large scatter of Hubble types from E to Sd. From redshift 0.25 to now the galaxies redden from ASK 2 to ASK 0, as expected from the passive evolution of their stellar populations. Two of the ASK classes (1 and 4) gather edge-on spirals, a property of interest in studies requiring knowing the intrinsic shape of a galaxy (e.g., weak lensing calibration).

We have derived masses and ages for 79 early-type galaxies (ETGs) in different environments at z~1.3 in the Lynx supercluster and in the GOODS/CDF-S field using multiwavelength (0.6-4.5 $\mu$m; KPNO, Palomar, Keck, HST, Spitzer) datasets. At this redshift the contribution of the TP-AGB phase is important for ETGs, and the mass and age estimates depend on the choice of the stellar population model used in the spectral energy distribution fits. We describe in detail the differences among model predictions for a large range of galaxy ages, showing the dependence of these differences on age. Current models still yield large uncertainties. While recent models from Maraston and Charlot & Bruzual offer better modeling of the TP-AGB phase with respect to less recent Bruzual & Charlot models, their predictions do not often match. The modeling of this TP-AGB phase has a significant impact on the derived parameters for galaxies observed at high-redshift. Some of our results do not depend on the choice of the model: for all models, the most massive galaxies are the oldest ones, independent of the environment. When using Maraston and Charlot & Bruzual models, the mass distribution is similar in the clusters and in the groups, whereas in our field sample there is a deficit of massive (M $\gtrsim$ 10^11 Msun) ETGs. According to those last models, ETGs belonging to the cluster environment host on average older stars with respect to group and field populations. This difference is less significant than the age difference in galaxies of different masses.

We study the environmental dependence of stellar population properties at z ~ 1.3. We derive galaxy properties (stellar masses, ages and star formation histories) for samples of massive, red, passive early-type galaxies in two high-redshift clusters, RXJ0849+4452 and RXJ0848+4453 (with redshifts of z = 1.26 and 1.27, respectively), and compare them with those measured for the RDCS1252.9-2927 cluster at z=1.24 and with those measured for a similarly mass-selected sample of field contemporaries drawn from the GOODS-South Field. Robust estimates of the aforementioned parameters have been obtained by comparing a large grid of composite stellar population models with extensive 8-10 band photometric coverage, from the rest-frame far-ultraviolet to the infrared. We find no variations of the overall stellar population properties among the different samples of cluster early-type galaxies. However, when comparing cluster versus field stellar population properties we find that, even if the (star formation weighted) ages are similar and depend only on galaxy mass, the ones in the field do employ longer timescales to assemble their final mass. We find that, approximately 1 Gyr after the onset of star formation, the majority (75%) of cluster galaxies have already assembled most (> 80%) of their final mass, while, by the same time, fewer (35%) field ETGs have. Thus we conclude that while galaxy mass regulates the timing of galaxy formation, the environment regulates the timescale of their star formation histories.

The GALEX Ultraviolet Virgo Cluster Survey (GUViCS) is a complete blind survey of the Virgo cluster covering about 40 sq. deg. in the far UV (FUV, lambda_eff=1539A, Delta-lambda=442A) and about 120 sq. deg. in the near UV (NUV, lambda_eff=2316A, Delta-lambda=1060A). The goal of the survey is to study the ultraviolet (UV) properties of galaxies in a rich cluster environment, spanning a wide luminosity range from giants to dwarfs, and regardless of prior knowledge of their star formation activity. The UV data will be combined with those in other bands (optical: NGVS; far-infrared - submm: HeViCS; HI: ALFALFA) and with our multizone chemo-spectrophotometric models of galaxy evolution to make a complete and exhaustive study of the effects of the environment on the evolution of galaxies in high density regions. We present here the scientific objectives of the survey, describing the observing strategy and briefly discussing different data reduction techniques. Using UV data already in-hand for the central 12 sq. deg. we determine the FUV and NUV luminosity functions of the Virgo cluster core for all cluster members and separately for early- and late-type galaxies and compare it to the one obtained in the field and other nearby clusters (Coma, A1367). This analysis shows that the FUV and NUV luminosity functions of the core of the Virgo clusters are flatter (alpha about -1.1) than those determined in Coma and A1367. We discuss the possible origin of this difference

New deep z'-J data readly show a narrow red sequence co-centered with, and similary concentrated to, the extended X-ray emission of the cluster of galaxies JKCS041. The JKCS041 red sequence is 0.32+/-0.06 mag redder in z'-J than the red sequence of the zspec=1.62 IRC0218A cluster, putting JKCS041 at z<<1.62. The colour difference of the two red sequences gives a red-sequence based redshift of z=2.20+/-0.11 for JKCS041, where the uncertainty accounts for uncertainties in stellar synthesis population models, in photometric calibration and in the red sequence colour of both JKCS041 and IRC0218A clusters.

This paper aims at robustly determining the redshift of the cluster of galaxies JKCS041 and at putting constraints on the formation epoch of the color-magnitude sequence in two very high redshift clusters. New deep z'-J data show a clear narrow red sequence that is co-centered with, and similarly concentrated on, the extended X-ray emission of the cluster of galaxies JKCS041. The JKCS041 red sequence is 0.32+/-0.06 mag redder in z'-J than the red sequence of the zspec=1.62 IRC0218A cluster, putting JKCS041 at z>>1.62 and ruling out z<~1.49 the latter claimed by a recent paper. The color difference of the two red sequences gives a red-sequence-based redshift of z=2.20+/-0.11 for JKCS041, where the uncertainty accounts for uncertainties in stellar synthesis population models, in photometric calibration, and in the red sequence color of both JKCS041 and IRC0218A clusters. We do not observe any sign of truncation of the red sequence for both clusters down to J=23 mag (1.0e+11 solar masses), which suggests that it is already in place in clusters rich and massive enough to heat and retain hot gas at these high redshifts.

We present an automated morphological classification in 4 types (E,S0,Sab,Scd) of ~700.000 galaxies from the SDSS DR7 spectroscopic sample based on support vector machines. The main new property of the classification is that we associate to each galaxy a probability of being in the four morphological classes instead of assigning a single class. The classification is therefore better adapted to nature where we expect a continuos transition between different morphological types. The algorithm is trained with a visual classification and then compared to several independent visual classifications including the Galaxy Zoo first release catalog. We find a very good correlation between the automated classification and classical visual ones. The compiled catalog is intended for use in different applications and can be downloaded at http://gepicom04.obspm.fr/sdss_morphology/Morphology_2010.html and soon from the CasJobs database.

Studying outliers from the bimodal distribution of galaxies in the color-mass space, such as morphological early-type galaxies residing in the blue cloud, can help to better understand the physical mechanisms that lead galaxy migrations in this space. In this paper we study the evolution of the properties of 210 M*/Msol>10^10 blue E/S0s between z~1.4 and z~0.2 in the COSMOS field with confirmed spectroscopic redshifts from the zCOSMOS 10k release. We first observe that the threshold mass, defined at z=0 in previous studies as the mass below which the population of blue early-type galaxies starts to be abundant relative to passive E/S0s, evolves from log(M*/Msol)~10.1 at z~0.3 to log(M*/Msol)~10.9 at z~1. Second, there seems to be a turn-over mass in the nature of blue E/S0 galaxies. Above log(M*/Msol)~10.8 blue E/S0 resemble to merger remnants probably migrating to the red-sequence in a time-scale of ~3 Gyr. Below this mass, they seem to be closer to normal late-type galaxies as if they were the result of minor mergers which triggered the central star-formation and built a central bulge component or were (re)building a disk from the surrounding gas, suggesting that they are moving back or staying in the blue-cloud. This turn-over mass does not seem to evolve significantly from z~1 in contrast with the threshold mass and therefore does not seem to be linked with the relative abundance of blue E/S0s.

This paper has been withdrawn A new version is in preparation and will be soon available.

[abridged] We analyze a sample of 9 massive clusters at 0.4<z<0.6 observed with MegaCam in 4 photometric bands (g,r,i,z) from the core to a radius of 5 Mpc (~4000 galaxies). Galaxy cluster candidates are selected using photometric redshifts computed with HyperZ. Morphologies are estimated with galSVM in two broad morphological types (early-type and late-type). We examine the morphological composition of the red-sequence and the blue-cloud and study the relations between galaxies and their environment through the morphology-density relations (T-Sigma) and the morphology-radius relation (T-R) in a mass limited sample (log(M/Msol)>9.5). We find that the red sequence is already in place at z~0.5 and it is mainly composed of very massive (log(M/Msol)>11.3) early-type galaxies. These massive galaxies seem to be already formed when they enter the cluster, probably in infalling groups, since the fraction remains constant with the cluster radius. Their presence in the cluster center could be explained by a segregation effect reflecting an early assembly history. Any evolution that takes place in the galaxy cluster population occurs therefore at lower masses (10.3<log(M/Msol)<11.3). For these galaxies, the evolution, is mainly driven by galaxy-galaxy interactions in the outskirts as revealed by the T-Sigma relation. Finally, the majority of less massive galaxies (9.5<log(M/Msol)<10.3) are late-type galaxies at all locations, suggesting that they have not started the morphological transformation yet even if this low mass bin might be affected by incompleteness.

We quantify the effects of \emph{morphological k-correction} at $1<z<2$ by comparing morphologies measured in the K and I-bands in the COSMOS area. Ks-band data have indeed the advantage of probing old stellar populations for $z<2$, enabling a determination of galaxy morphological types unaffected by recent star formation. In paper I we presented a new non-parametric method to quantify morphologies of galaxies on seeing limited images based on support vector machines. Here we use this method to classify $\sim$$50 000$ $Ks$ selected galaxies in the COSMOS area observed with WIRCam at CFHT. The obtained classification is used to investigate the redshift distributions and number counts per morphological type up to $z\sim2$ and to compare to the results obtained with HST/ACS in the I-band on the same objects from other works. We associate to every galaxy with $Ks<21.5$ and $z<2$ a probability between 0 and 1 of being late-type or early-type. The classification is found to be reliable up to $z\sim2$. The mean probability is $p\sim0.8$. It decreases with redshift and with size, especially for the early-type population but remains above $p\sim0.7$. The classification is globally in good agreement with the one obtained using HST/ACS for $z<1$. Above $z\sim1$, the I-band classification tends to find less early-type galaxies than the Ks-band one by a factor $\sim$1.5 which might be a consequence of morphological k-correction effects. We argue therefore that studies based on I-band HST/ACS classifications at $z>1$ could be underestimating the elliptical population. [abridged]

We present a new non-parametric method to quantify morphologies of galaxies based on a particular family of learning machines called support vector machines. The method, that can be seen as a generalization of the classical CAS classification but with an unlimited number of dimensions and non-linear boundaries between decision regions, is fully automated and thus particularly well adapted to large cosmological surveys. The source code is available for download at http://www.lesia.obspm.fr/~huertas/galsvm.html To test the method, we use a seeing limited near-infrared ($K_s$ band, $2,16\mu m$) sample observed with WIRCam at CFHT at a median redshift of $z\sim0.8$. The machine is trained with a simulated sample built from a local visually classified sample from the SDSS chosen in the high-redshift sample's rest-frame (i band, $0.77\mu m$) and artificially redshifted to match the observing conditions. We use a 12-dimensional volume, including 5 morphological parameters and other caracteristics of galaxies such as luminosity and redshift. We show that a qualitative separation in two main morphological types (late type and early type) can be obtained with an error lower than 20% up to the completeness limit of the sample ($KAB\sim 22$) which is more than 2 times better that what would be obtained with a classical C/A classification on the same sample and indeed comparable to space data. The method is optimized to solve a specific problem, offering an objective and automated estimate of errors that enables a straightforward comparison with other surveys.

We present the results of an imaging programme of distant galaxies (z~0.8) at high spatial resolution (~0.1").We observed 7 fields of 1'*1' with the NACO Adaptive Optics system (VLT) in Ks (2.16um) band with typical V ~ 14 guide stars and 3h integration time per field. Observed fields are selected within the COSMOS survey area. High angular resolution K-band data have the advantage to probe old stellar populations in the rest-frame, enabling to determine galaxy morphological types unaffected by recent star formation, better linked to the underlying mass than classical optical morphology studies (HST). Adaptive optics on ground based telescopes is the only method today to obtain such high resolution in the K-band. In this paper we show that reliable results can be obtained and establish a first basis for larger observing programmes. We analyze the morphologies by means of B/D (Bulge/Disk) decomposition with GIM2D and CAS (Concentration-Asymmetry) estimators for 79 galaxies with magnitudes between Ks = 17-23 and classify them in three main morphological types (Late Type, Early Type and Irregulars). We obtain for the first time an estimate of the distribution of galaxy types at redshift z ~ 1 as measured from the near infrared at high spatial resolution. We show that galactic parameters (disk scale length, bulge effective radius and bulge fraction) can be estimated with a random error lower than 20% for the bulge fraction up to Ks = 19 (AB = 21) and that classification into the three main morphological types can be done up to Ks = 20 (AB = 22) with at least 70% of correct identifications. We used the known photometric redshifts to obtain a redshift distribution over 2 redshift bins (z < 0.8, 0.8 < z < 1.5) for each morphological type.