sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

De Lucia, Gabriella

Normalized to: De Lucia, G.

186 article(s) in total. 609 co-authors, from 1 to 46 common article(s). Median position in authors list is 4,0.

[1]  [pdf] - 2071928
The influence of environment on satellite galaxies in the GAEA semi-analytic model
Comments: 19 pages, 14 figures; submitted to MNRAs; comments are welcome
Submitted: 2020-03-28
Reproducing the observed quenched fraction of satellite galaxies has been a long-standing issue for galaxy formation models. We modify the treatment of environmental effects in our state-of-the-art GAlaxy Evolution and Assembly (GAEA) semi-analytic model to improve our modelling of satellite galaxies. Specifically, we implement gradual stripping of hot gas, ram-pressure stripping of cold gas, and an updated algorithm to account for angular momentum exchanges between the gaseous and stellar disc components of model galaxies. Our updated model is in good agreement with the observed quenched fractions for central and satellite galaxies, and their dependencies on stellar mass and halo mass. We also find consistency between model predictions and observational estimates of quenching time for satellite galaxies; HI, H$_2$ fractions of central galaxies; and deficiencies of HI, H$_2$, SFR of galaxies in cluster haloes. Analysing results from our updated model, we find that the dominant quenching mechanisms are hot-gas stripping for low-mass galaxies, and AGN feedback for massive galaxies. The Ram-pressure stripping of cold gas has a limited impact on the predicted quenched fractions, nonetheless, is needed to reproduce the observed HI deficiencies.
[2]  [pdf] - 2063806
The GOGREEN survey: The environmental dependence of the star-forming galaxy main sequence at $1.0<z<1.5$
Comments: Accepted for publication in MNRAS on February 25 2020
Submitted: 2020-02-26, last modified: 2020-03-13
We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3\sigma$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales.
[3]  [pdf] - 2054270
The Rise of Active Galactic Nuclei in the GAlaxy Evolution and Assembly semi-analytic model
Comments: 16 pages, 12 figures, 1 table, MNRAS submitted
Submitted: 2020-02-24
We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disk instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} {\rm M}_\odot$) without changing their overall assembly history. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 {\rm M}_\odot$ are already in a self-regulated state as inferred from observations.
[4]  [pdf] - 1994452
Modelling the quenching of star formation activity from the evolution of the colour-magnitude relation in VIPERS
Comments: 21 pages, submitted to MNRAS
Submitted: 2019-11-06
We study the evolution of the colour-magnitude relation for galaxies in the VIMOS Public Extragalactic Redshift Survey (VIPERS) by introducing the concept of the bright edge, and use this to derive constraints on the quenching of star formation activity in galaxies over the redshift range $0.5 < z < 1.1$. The bright-edge of the colour-magnitude diagram evolves with little dependence on galaxy colour, and therefore on the amount of star formation taking place in bright galaxies. We modelled this evolution with delayed exponential star formation histories (SFHs), to better understand the time-scale of the turn-off in star formation activity. We show that using SFHs without quenching, the transition from the blue cloud to the red sequence is too slow. This indicates that a scenario purely driven by the consumption of the gas inside each galaxy does not reproduce the observed evolution of the colour-magnitude bright edge. Among the quenching scenarios explored, the one that best matches the observations assumes that galaxies stop their star formation at a randomly selected time with a uniform distribution up to $2.5$ Gyr. We argue that quenching is required over a wide range of stellar masses. Qualitatively similar evolution of the bright edge is found in the predictions of a semi-analytical galaxy formation model, but quantitatively there are marked differences with the observations. This illustrates the utility of the bright edge as a test of galaxy formation models. The evolution changes and no longer matches the observed trend if feedback from heating by active galactic nuclei is turned off.
[5]  [pdf] - 1987343
Preprocessing Among the Infalling Galaxy Population of EDisCS Clusters
Comments: Accepted for publication in ApJ
Submitted: 2015-06-05, last modified: 2019-10-29
We present results from a low-resolution spectroscopic survey for 21 galaxy clusters at $0.4 < z < 0.8$ selected from the ESO Distant Cluster Survey. We measured spectra using the low-dispersion prism in IMACS on the Magellan Baade telescope and calculate redshifts with an accuracy of $\sigma_z = 0.007$. We find 1763 galaxies that are brighter than $R = 22.9$ in the large-scale cluster environs. We identify the galaxies expected to be accreted by the clusters as they evolve to $z = 0$ using spherical infall models and find that $\sim30\%$ to $\sim70\%$ of the $z = 0$ cluster population lies outside the virial radius at $z \sim 0.6$. For analogous clusters at $z = 0$, we calculate that the ratio of galaxies that have fallen into the clusters since $z \sim 0.6$ to those that were already in the core at that redshift is typically between $\sim0.3$ and $1.5$. This wide range of ratios is due to intrinsic scatter and is not a function of velocity dispersion, so a variety of infall histories is to be expected for clusters with current velocity dispersions of $300 \lesssim\sigma\lesssim 1200$ km s$^{-1}$. Within the infall regions of $z \sim 0.6$ clusters, we find a larger red fraction of galaxies than in the field and greater clustering among red galaxies than blue. We interpret these findings as evidence of "preprocessing", where galaxies in denser local environments have their star formation rates affected prior to their aggregation into massive clusters, although the possibility of backsplash galaxies complicates the interpretation.
[6]  [pdf] - 1979738
The VIMOS Public Extragalactic Redshift Survey (VIPERS) -- Environment-size relation of massive passive galaxies at 0.5 < z < 0.8
Comments: 12 pages, accepted for publication in A&A
Submitted: 2019-09-05
We use the statistics of the VIPERS survey to investigate the relation between the surface mean stellar mass density Sigma=Mstar/(2*pi*Re^2) of massive passive galaxies (MPGs, Mstar>10^11 Msun) and their environment in the redshift range 0.5<z<0.8. Passive galaxies were selected on the basis of their NUVrK colors (~900 objects), and the environment was defined as the galaxy density contrast, delta, using the fifth nearest-neighbor approach. The analysis of Sigma vs. delta was carried out in two stellar mass bins. In galaxies with Mstar<2*10^11 Msun, no correlation between Sigma and delta is observed. This implies that the accretion of satellite galaxies, which is more frequent in denser environments and efficient in reducing the galaxy Sigma, is not relevant in the formation and evolution of these systems. Conversely, in galaxies with Mstar>2*10^11 Msun, we find an excess of MPGs with low Sigma and a deficit of high-Sigma MPGs in the densest regions wrt other environments. We interpret this result as due to the migration of some high-Sigma MPGs (<1% of the total population of MPGs) into low-Sigma MPGs, probably through mergers or cannibalism of small satellites. In summary, our results imply that the accretion of satellite galaxies has a marginal role in the mass-assembly history of most MPGs. We have previously found that the number density of VIPERS massive star-forming galaxies (MSFGs) declines rapidily from z=0.8 to z=0.5, which mirrors the rapid increase in the number density of MPGs. This indicates that the MSFGs at z>0.8 migrate to the MPG population. Here, we investigate the Sigma-delta relation of MSFGs at z>0.8 and find that it is consistent within 1 sigma with that of low-Sigma MPGs at z<0.8. Thus, the results of this and our previous paper show that MSFGs at z>0.8 are consistent in terms of number and environment with being the progenitors of low-Sigma MPGs at z<0.8.
[7]  [pdf] - 2061569
The atomic Hydrogen content of the post-reionization Universe
Comments: 21 pages, 20 figures, comments welcome
Submitted: 2019-09-05
We present a comprehensive analysis of atomic hydrogen (HI) properties using a semi-analytical model of galaxy formation and N-body simulations covering a large cosmological volume at high resolution. We examine the HI mass function and the HI density, characterizing both their redshift evolution and their dependence on hosting halo mass. We analyze the HI content of dark matter haloes in the local Universe and up to redshift $z=5$, discussing the contribution of different galaxy properties. We find that different assembly history plays a crucial role in the scatter of this relation. We propose new fitting functions useful for constructing mock HI maps with HOD techniques. We investigate the HI clustering properties relevant for future $21$~cm Intensity Mapping (IM) experiments, including the HI bias and the shot noise level. The HI bias increases with redshift and it is roughly flat on the largest scales probed. The scale dependency is found at progressively larger scales with increasing redshift, apart from a dip feature at $z=0$. The shot-noise values are consistent with the ones inferred by independent studies, confirming that shot-noise will not be a limiting factor for IM experiments. We detail the contribution from various galaxy properties on the HI power spectrum and their relation to the halo bias. We find that HI poor satellite galaxies play an important role at the scales of the 1-halo term. Finally, we present the $21$~cm signal in redshift space, a fundamental prediction to be tested against data from future radio telescopes such as SKA.
[8]  [pdf] - 1928142
The Rest-frame $H$-band Luminosity Function of Red Sequence Galaxies in Clusters at $1.0 < z < 1.3$
Comments: Accepted for publication in ApJ. 24 pages, 13 figures
Submitted: 2019-06-25, last modified: 2019-07-11
We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 \mu$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ \alpha \sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15.
[9]  [pdf] - 1905883
Euclid preparation III. Galaxy cluster detection in the wide photometric survey, performance and algorithm selection
Comments: A&A, in press. Only changes in affiliations and names between v1,v2,v3
Submitted: 2019-06-11, last modified: 2019-06-19
Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]
[10]  [pdf] - 1905693
The Evolution of Sizes and Specific Angular Momenta in Hierarchical Models of Galaxy Formation and Evolution
Comments: Accepted for publication in MNRAS, 17 pages, 11 figures
Submitted: 2019-02-27, last modified: 2019-06-13
We extend our previous work focused at $z\sim0$, studying the redshift evolution of galaxy dynamical properties using the state-of-the-art semi-analytic model GAEA: we show that the predicted size-mass relation for disky/star forming and quiescent galaxies is in good agreement with observational estimates, up to $z\sim2$. Bulge dominated galaxies have sizes that are offset low with respect to observational estimates, mainly due to our implementation of disk instability at high redshift. At large masses, both quiescent and bulge dominated galaxies have sizes smaller than observed. We interpret this as a consequence of our most massive galaxies having larger gas masses than observed, and therefore being more affected by dissipation. We argue that a proper treatment of quasar driven winds is needed to alleviate this problem. Our model compact galaxies have number densities in agreement with observational estimates and they form most of their stars in small and low angular momentum high-$z$ halos. GAEA predicts that a significant fraction of compact galaxies forming at high-$z$ is bound to merge with larger structures at lower redshifts: therefore they are not the progenitors of normal-size passive galaxies at $z=0$. Our model also predicts a stellar-halo size relation that is in good agreement with observational estimates. The ratio between stellar size and halo size is proportional to the halo spin and does not depend on stellar mass but for the most massive galaxies, where AGN feedback leads to a significant decrease of the retention factor (from about 80 per cent to 20 per cent).
[11]  [pdf] - 1849164
Starburst galaxies in semi-analytic models of galaxy formation and evolution
Comments: 14 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2018-09-17, last modified: 2019-03-14
We study the shape and evolution of the star formation main sequence in three independently developed semi-analytic models of galaxy formation. We focus, in particular, on the characterization of the model galaxies that are significantly above the main sequence, and that can be identified with galaxies classified as `starburst' in recent observational work. We find that, in all three models considered, star formation triggered by merger events (both minor and major) contribute to only a very small fraction of the cosmic density of star formation. While mergers are associated to bursts of star formation in all models, galaxies that experienced recent merger events are not necessarily significantly above the main sequence. On the other hand, `starburst galaxies' are not necessarily associated with merger episodes, especially at the low-mass end. Galaxies that experienced recent mergers can have relatively low levels of star formation when/if the merger is gas-poor, and galaxies with no recent merger can experience episodes of starbursts due to a combination of large amount of cold gas available from cooling/accretion events and/or small disk radii which increases the cold gas surface density.
[12]  [pdf] - 1818881
Physical properties of SDSS satellite galaxies in projected phase-space
Comments: 24 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2019-01-14
We investigate how environment affects satellite galaxies using their location within the projected phase-space of their host haloes from the Wang et al.'s group catalogue. Using the Yonsei Zoom in Cluster Simulations, we derive zones of constant mean infall time T_inf in projected phase-space, and catalogue in which zone each observed galaxy falls. Within each zone we compute the mean observed galaxy properties including specific star formation rate, luminosity-weighted age, stellar metallicity and [alpha/Fe] abundance ratio. By comparing galaxies in different zones, we inspect how shifting the mean infall time from recent infallers (mean T_inf < 3 Gyr) to ancient infallers (mean T_inf > 5 Gyr) impacts galaxy properties at fixed stellar and halo mass. Ancient infallers are more quenched, and the impact of environmental quenching is visible down to low host masses (< group masses). Meanwhile, the quenching of recent infallers is weakly dependent on host mass, indicating they have yet to respond strongly to their current environment. [alpha/Fe] and especially metallicity are less dependent on host mass, but show a dependence on mean T_inf. We discuss these results in the context of longer exposure times for ancient infallers to environmental effects, which grow more efficient in hosts with a deeper potential well and a denser intracluster medium. We also compare our satellites with a control field sample, and find that even the most recent infallers (mean T_inf < 2 Gyr) are more quenched than field galaxies, in particular for cluster mass hosts. This supports the role of pre-processing and/or faster quenching in satellites.
[13]  [pdf] - 1799841
The VIMOS Public Extragalactic Redshift Survey (VIPERS). The complexity of galaxy populations at 0.4< z<1.3 revealed with unsupervised machine-learning algorithms
Comments: published in A&A, 27 pages, 10 figures, 4 appendixes. Please contact Malgorzata Siudek ( if you are interested in the catalogue
Submitted: 2018-05-24, last modified: 2018-12-18
Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter space of 12 rest-frame magnitudes and spectroscopic redshift. The model (DBk) and the number of classes (12) were established based on the joint analysis of standard statistical criteria and confirmed by the analysis of the galaxy distribution with respect to a number of classes and their properties. This new approach allows us to classify galaxies based just on their redshifts and UV-NIR spectral energy distributions. The FEM unsupervised algorithm has automatically distinguished 12 classes: 11 classes of VIPERS galaxies and an additional class of broad-line AGNs. After a first broad division into blue, green and red categories we obtained a further sub-division into three red, three green, and five blue galaxy classes. The FEM classes follow the galaxy sequence from the earliest to the latest types that is reflected in their colours (which are constructed from rest-frame magnitudes used in classification procedure) but also their morphological, physical, and spectroscopic properties (not included in the classification scheme). We demonstrate that the members of each class share similar physical and spectral properties. In particular, we are able to find three different classes of red passive galaxy populations. Thus, we demonstrate the potential of an unsupervised approach to galaxy classification and we retrieve the complexity of galaxy populations at z~0.7, a task that usual simpler colour-based approaches cannot fulfil.
[14]  [pdf] - 1799842
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Unsupervised classification with photometric redshifts: a method to accurately classify large galaxy samples without spectroscopic information
Comments: the letter submitted to MNRAS, 6 pages, 4 figures. Please contact Malgorzata Siudek ( if you are interested in the catalogue
Submitted: 2018-05-24, last modified: 2018-12-18
Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et al. (2018, hereafter: S1) for galaxies at z~0.7 could be reproduced if only photometric data are available. We applied the Fisher Expectation-Maximization (FEM) unsupervised clustering algorithm to 54,293 VIPERS galaxies working in a parameter space of reliable photometric redshifts and 12 corresponding rest-frame magnitudes. The FEM algorithm distinguishes four main groups: (1) red, (2) green, (3) blue, and (4) outliers. Each group is further divided into 3, 3, 4, and 2 subclasses, respectively. The accuracy of reproducing galaxy classes using spectroscopic data is high: 92%, 84%, 96% for red, green, and blue classes, respectively, except for dusty star-forming galaxies. The presented verification of the photometric classification demonstrates that large photometric samples can be used to distinguish different galaxy classes at z > 0.5 with an accuracy provided so far only by spectroscopic data except for particular galaxy classes.
[15]  [pdf] - 1790743
Nature versus nurture: what regulates star formation in satellite galaxies?
Comments: 11 pages, 6 figures, 1 appendix. Accepted for publication in MNRAS
Submitted: 2018-06-01, last modified: 2018-11-07
We use our state-of-the-art Galaxy Evolution and Assembly (GAEA) semi-analytic model to study how and on which time-scales star formation is suppressed in satellite galaxies. Our fiducial stellar feedback model, implementing strong stellar driven outflows, reproduces relatively well the variations of passive fractions as a function of galaxy stellar mass and halo mass measured in the local Universe, as well as the `quenching' time-scales inferred from the data. We show that the same level of agreement can be obtained by using an alternative stellar feedback scheme featuring lower ejection rates at high redshift, and modifying the treatment for hot gas stripping. This scheme over-predicts the number densities of low to intermediate mass galaxies. In addition, a good agreement with the observed passive fractions can be obtained only by assuming that cooling can continue on satellites, at the rate predicted considering halo properties at infall, even after their parent dark matter substructure is stripped below the resolution of the simulation. For our fiducial model, the better agreement with the observed passive fractions can be ascribed to: (i) a larger cold gas fraction of satellites at the time of accretion, and (ii) a lower rate of gas reheating by supernovae explosions and stellar winds with respect to previous versions of our model. Our results suggest that the abundance of passive galaxies with stellar mass larger than ~10^10 Msun is primarily determined by the self-regulation between star formation and stellar feedback, with environmental processes playing a more marginal role.
[16]  [pdf] - 1742077
The VIMOS Public Extragalactic Redshift Survey (VIPERS): Unbiased clustering estimate with VIPERS slit assignment
Comments: 15 pages, 13 figures - Accepted for publication on A&A
Submitted: 2018-07-16, last modified: 2018-08-30
The VIPERS galaxy survey has measured the clustering of $0.5<z<1.2$ galaxies, enabling a number of measurements of galaxy properties and cosmological redshift-space distortions (RSD). Because the measurements were made using one-pass of the VIMOS instrument on the Very Large Telescope (VLT), the galaxies observed only represent approximately 47\% of the parent target sample, with a distribution imprinted with the pattern of the VIMOS slitmask. Correcting for the effect on clustering has previously been achieved using an approximate approach developed using mock catalogues. Pairwise inverse probability (PIP) weighting has recently been proposed by Bianchi & Percival to correct for missing galaxies, and we apply it to mock VIPERS catalogues to show that it accurately corrects the clustering for the VIMOS effects, matching the clustering measured from the observed sample to that of the parent. We then apply PIP-weighting to the VIPERS data, and fit the resulting monopole and quadrupole moments of the galaxy two-point correlation function with respect to the line-of-sight, making measurements of RSD. The results are close to previous measurements, showing that the previous approximate methods used by the VIPERS team are sufficient given the errors obtained on the RSD parameter.
[17]  [pdf] - 1738887
Structural and Dynamical Properties of Galaxies in a Hierarchical Universe: Sizes and Specific Angular Momenta
Comments: 26 pages, 14 figures, 4 appendices. Accepted for publication in MNRAS
Submitted: 2018-03-21, last modified: 2018-08-27
We use a state-of-the-art semi-analytic model to study the size and the specific angular momentum of galaxies. Our model includes a specific treatment for the angular momentum exchange between different galactic components. Disk scale radii are estimated from the angular momentum of the gaseous/stellar disk, while bulge sizes are estimated assuming energy conservation. The predicted size--mass and angular momentum--mass relations are in fair agreement with observational measurements in the local Universe, provided a treatment for gas dissipation during major mergers is included. Our treatment for disk instability leads to unrealistically small radii of bulges formed through this channel, and predicts an offset between the size--mass relations of central and satellite early-type galaxies, that is not observed. The model reproduces the observed dependence of the size--mass relation on morphology, and predicts a strong correlation between specific angular momentum and cold gas content. This correlation is a natural consequence of galaxy evolution: gas-rich galaxies reside in smaller halos, and form stars gradually until present day, while gas-poor ones reside in massive halos, that formed most of their stars at early epochs, when the angular momentum of their parent halos is low. The dynamical and structural properties of galaxies can be strongly affected by a different treatment for stellar feedback, as this would modify their star formation history. A higher angular momentum for gas accreted through rapid mode does not affect significantly the properties of massive galaxies today, but has a more important effect on low-mass galaxies at higher redshift.
[18]  [pdf] - 1728094
On the influence of environment on star forming galaxies
Comments: 15 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2018-08-05
We use our state-of-the-art semi analytic model for GAlaxy Evolution and Assembly (GAEA), and observational measurements of nearby galaxies to study the influence of the environment on the gas content and gaseous/stellar disc sizes of star-forming galaxies. We analyse the origin of differences between physical properties of satellites and those of their central counterparts, identified by matching the Vmax of their host haloes at the accretion time of the satellites. Our model reproduces nicely the differences between centrals and satellites measured for the HI mass, size of the star-forming region, and stellar radii. In contrast, our model predicts larger differences with respect to data for the molecular gas mass and star formation rate. By analysing the progenitors of central and satellite model galaxies, we find that differences in the gas content arise after accretion, and can be entirely ascribed to the instantaneous stripping of the hot gas reservoir. The suppression of cold gas replenishment via cooling and star formation leads to a reduction of the cold gas and of its density. Therefore, more molecular gas is lost than lower density HI gas, and model satellites have less molecular gas and lower star formation rates than observed satellites. We argue that these disagreements could be largely resolved with the inclusion of a proper treatment for ram-pressure stripping of cold gas and a more gradual stripping of the hot gas reservoir. A more sophisticated treatment of angular momentum exchanges, accounting for the multi-phase nature of the gaseous disc is also required.
[19]  [pdf] - 1717236
On the shape and evolution of a cosmic ray regulated galaxy-wide stellar initial mass function
Comments: 8 pages, 5 figures, MNRAS accepted
Submitted: 2018-07-03
In this paper, we present a new derivation of the shape and evolution of the integrated galaxy-wide initial mass function (IGIMF), incorporating explicitly the effects of cosmic rays (CRs) as regulators of the chemical and thermal state of the gas in the dense cores of molecular clouds. We predict the shape of the IGIMF as a function of star formation rate (SFR) and CR density, and show that it can be significantly different with respect to local estimates. In particular, we focus on the physical conditions corresponding to IGIMF shapes that are simultaneously shallower at high-mass end and steeper at the low-mass end than a Kroupa IMF. These solutions can explain both the levels of $\alpha$-enrichment and the excess of low-mass stars as a function of stellar mass, observed for local spheroidal galaxies. As a preliminary test of our scenario, we use idealized star formation histories to estimate the mean IMF shape for galaxies of different $z=0$ stellar mass. We show that the fraction of low-mass stars as a function of galaxy stellar mass predicted by these mean IMFs agrees with the values derived from high-resolution spectroscopic surveys.
[20]  [pdf] - 1634308
Variations of the stellar initial mass function in semi-analytical models II: the impact of Cosmic Ray regulation
Comments: 10 Pages (+3 in Appendices), 6 figures (+4 in Appendices), 2 tables, MNRAS in press
Submitted: 2018-01-08
Recent studies proposed that cosmic rays (CR) are a key ingredient in setting the conditions for star formation, thanks to their ability to alter the thermal and chemical state of dense gas in the UV-shielded cores of molecular clouds. In this paper, we explore their role as regulators of the stellar initial mass function (IMF) variations, using the semi-analytic model for GAlaxy Evolution and Assembly (GAEA). The new model confirms our previous results obtained using the integrated galaxy-wide IMF (IGIMF) theory: both variable IMF models reproduce the observed increase of $\alpha$-enhancement as a function of stellar mass and the measured $z=0$ excess of dynamical mass-to-light ratios with respect to photometric estimates assuming a universal IMF. We focus here on the mismatch between the photometrically-derived ($M^{\rm app}_{\star}$) and intrinsic ($M_{\star}$) stellar masses, by analysing in detail the evolution of model galaxies with different values of $M_{\star}/M^{\rm app}_{\star}$. We find that galaxies with small deviations (i.e. formally consistent with a universal IMF hypothesis) are characterized by more extended star formation histories and live in less massive haloes with respect to the bulk of the galaxy population. While the IGIMF theory does not change significantly the mean evolution of model galaxies with respect to the reference model, a CR-regulated IMF implies shorter star formation histories and higher peaks of star formation for objects more massive than $10^{10.5} M_\odot$. However, we also show that it is difficult to unveil this behaviour from observations, as the key physical quantities are typically derived assuming a universal IMF.
[21]  [pdf] - 1799807
The VIMOS Public Extragalactic Redshift Survey (VIPERS). AGN feedback in [NeV] emitters
Comments: 13 pages, 9 figures, submitted to A&A
Submitted: 2017-12-21
Using an unconventional single line diagnostic that unambiguously identifies AGNs in composite galaxies we report statistical differences in the properties (stellar age, [OII] luminosity, colour) between active and inactive galaxies at 0.62<z<1.2 extracted from the VIMOS Public Extragalactic Redshift Survey (VIPERS). The nuclear activity is probed by the high-ionization [NeV] emission line and along with their parent samples, the galaxies are properly selected according to their stellar mass, redshift, and colour distributions. We report younger underlying stellar ages and higher [OII] luminosities of active galaxies in the green valley and in the blue cloud compared to control samples. We observe higher fractions of green galaxies hosting AGN activity at progressively bluer (r-K) colours. Depending on the location of the host galaxy in the NUVrK colour diagram we find higher AGN fractions in massive blue galaxies and in the least massive red galaxies, in agreement with the picture that black holes vary their properties when hosted in either star-forming or passive galaxies. Exactly where the fast quenching processes are expected to play a role, we identify a novel class of active galaxies in the blue cloud with signatures typical for a suddenly suppression of their star formation activity after a burst happening in the recent past. Their optical spectra resemble those of post-starburst galaxies, that would never be identified in a spectroscopic search using classical selection techniques. Broadly, these active galaxies selected on the [NeV] line are not commonly represented in shallow X-ray, mid-IR, or classical line diagnostics. If we consider that our results are limited by the shallow observational limits and rapid AGN variability, the impact of AGN feedback on galaxy formation and evolution may represent an important channel of fast-transiting galaxies moving to the red sequence.
[22]  [pdf] - 1641261
The VIMOS Public Extragalactic Redshift Survey (VIPERS): An unbiased estimate of the growth rate of structure at $\mathbf{\left<z\right>=0.85}$ using the clustering of luminous blue galaxies
Comments: 20 pages, 21 figures, accepted for publication in A&A
Submitted: 2017-07-31, last modified: 2017-10-23
We use the VIPERS final data release to investigate the performance of colour-selected populations of galaxies as tracers of linear large-scale motions. We empirically select volume-limited samples of blue and red galaxies as to minimise the systematic error on the estimate of the growth rate $f\sigma_8$ from the anisotropy of the two-point correlation function. To this end, rather than rigidly splitting the sample into two colour classes we define the red/blue fractional contribution of each object through a weight based on the $(U-V)$ colour distribution. Using mock surveys that are designed to reproduce the observed properties of VIPERS galaxies, we find the systematic error in recovering the fiducial value of $f\sigma_8$ to be minimized when using a volume-limited sample of luminous blue galaxies. We model non-linear corrections via the Scoccimarro extension of the Kaiser model, finding systematic errors on $f\sigma_8$ of below $1-2\%$, using scales as small as 5 $h^{-1}\mathrm{Mpc}$. We interpret this result as indicating that selection of luminous blue galaxies maximises the fraction that are central objects in their dark matter haloes; this in turn minimises the contribution to the measured $\xi(r_p,\pi)$ from the 1-halo term, which is dominated by non-linear motions. The gain is inferior if one uses the full magnitude-limited sample of blue objects, consistent with the presence of a significant fraction of blue, fainter satellites dominated by non-streaming, orbital velocities. We measure a value of $f\sigma_8=0.45 \pm 0.11$ over the single redshift range $0.6\le z\le 1.0$, corresponding to an effective redshift for the blue galaxies $\left<z\right>=0.85$. Including in the likelihood the potential extra information contained in the blue-red galaxy cross-correlation function does not lead to an appreciable improvement in the error bars, while it increases the systematic error.
[23]  [pdf] - 1602626
Determining the Halo Mass Scale where Galaxies Lose Their Gas
Comments: 30 pages, 14 figures, 3 tables. Accepted for publication in the Astrophysical Journal
Submitted: 2017-09-12
A major question in galaxy formation is how the gas supply that fuels activity in galaxies is modulated by their environment. We use spectroscopy of a set of well characterized clusters and groups at $0.4<z<0.8$ from the ESO Distant Cluster Survey (EDisCS) and compare it to identically selected field galaxies. Our spectroscopy allows us to isolate galaxies that are dominated by old stellar populations. Here we study a stellar-mass limited sample ($\log(M_*/M_\odot)>10.4$) of these old galaxies with weak [OII] emission. We use line ratios and compare to studies of local early type galaxies to conclude that this gas is likely excited by post-AGB stars and hence represents a diffuse gas component in the galaxies. For cluster and group galaxies the fraction with EW([OII])$>5$\AA\ is $f_{[OII]}=0.08^{+0.03}_{-0.02}$ and $f_{[OII]}=0.06^{+0.07}_{-0.04}$ respectively. For field galaxies we find $f_{[OII]}=0.27^{+0.07}_{-0.06}$, representing a 2.8$\sigma$ difference between the [OII] fractions for old galaxies between the different environments. We conclude that a population of old galaxies in all environments has ionized gas that likely stems from stellar mass loss. In the field galaxies also experience gas accretion from the cosmic web and in groups and clusters these galaxies have had their gas accretion shut off by their environment. Additionally, galaxies with emission preferentially avoid the virialized region of the cluster in position-velocity space. We discuss the implications of our results, among which is that gas accretion shutoff is likely effective at group halo masses (log~${\cal M}/$\msol$>12.8$) and that there are likely multiple gas removal processes happening in dense environments.
[24]  [pdf] - 1597995
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Gravity test from the combination of redshift-space distortions and galaxy-galaxy lensing at $0.5 < z < 1.2$
Comments: 21 pages, 19 figures, accepted for publication in A&A (referee report received seven months after submission)
Submitted: 2016-12-16, last modified: 2017-08-28
We carry out a joint analysis of redshift-space distortions and galaxy-galaxy lensing, with the aim of measuring the growth rate of structure; this is a key quantity for understanding the nature of gravity on cosmological scales and late-time cosmic acceleration. We make use of the final VIPERS redshift survey dataset, which maps a portion of the Universe at a redshift of $z \simeq 0.8$, and the lensing data from the CFHTLenS survey over the same area of the sky. We build a consistent theoretical model that combines non-linear galaxy biasing and redshift-space distortion models, and confront it with observations. The two probes are combined in a Bayesian maximum likelihood analysis to determine the growth rate of structure at two redshifts $z=0.6$ and $z=0.86$. We obtain measurements of $f\sigma_8(0.6) = 0.48 \pm 0.12$ and $f\sigma_8(0.86) = 0.48 \pm 0.10$. The additional galaxy-galaxylensing constraint alleviates galaxy bias and $\sigma_8$ degeneracies, providing direct measurements of $[f(0.6),\sigma_8(0.6)] = [0.93 \pm 0.22, 0.52 \pm 0.06]$ and $f(0.86),\sigma_8(0.86)] = [0.99 \pm 0.19, 0.48 \pm 0.04]$. These measurements are statistically consistent with a Universe where the gravitational interactions can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives. Finally, as a complementary test we measure the gravitational slip parameter, $E_G$ , for the first time at $z>0.6$. We find values of $\smash{\overline{E}_G}(0.6) = 0.16 \pm 0.09$ and $\smash{\overline{E}_G}(0.86) = 0.09 \pm 0.07$, when $E_G$ is averaged over scales above $3 h^{-1} \rm{Mpc}$. We find that our $E_G$ measurements exhibit slightly lower values than expected for standard relativistic gravity in a {\Lambda}CDM background, although the results are consistent within $1-2\sigma$.
[25]  [pdf] - 1580718
H$_2$-based star formation laws in hierarchical models of galaxy formation
Comments: 28 pages, 21 figures, MNRAS accepted
Submitted: 2016-11-28, last modified: 2017-07-19
We update our recently published model for GAlaxy Evolution and Assembly (GAEA), to include a self-consistent treatment of the partition of cold gas in atomic and molecular hydrogen. Our model provides significant improvements with respect to previous ones used for similar studies. In particular, GAEA (i) includes a sophisticated chemical enrichment scheme accounting for non-instantaneous recycling of gas, metals, and energy; (ii) reproduces the measured evolution of the galaxy stellar mass function; (iii) reasonably reproduces the observed correlation between galaxy stellar mass and gas metallicity at different redshifts. These are important prerequisites for models considering a metallicity dependent efficiency of molecular gas formation. We also update our model for disk sizes and show that model predictions are in nice agreement with observational estimates for the gas, stellar and star forming disks at different cosmic epochs. We analyse the influence of different star formation laws including empirical relations based on the hydro-static pressure of the disk, analytic models, and prescriptions derived from detailed hydro-dynamical simulations. We find that modifying the star formation law does not affect significantly the global properties of model galaxies, neither their distributions. The only quantity showing significant deviations in different models is the cosmic molecular-to-atomic hydrogen ratio, particularly at high redshift. Unfortunately, however, this quantity also depends strongly on the modelling adopted for additional physical processes. Useful constraints on the physical processes regulating star formation can be obtained focusing on low mass galaxies and/or at higher redshift. In this case, self-regulation has not yet washed out differences imprinted at early time.
[26]  [pdf] - 1583093
Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) I: Survey Description
Comments: Accepted for publication in MNRAS May 26, 2017
Submitted: 2017-05-03, last modified: 2017-06-27
We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.
[27]  [pdf] - 1580669
The VIMOS Public Extragalactic Redshift Survey (VIPERS). The matter density and baryon fraction from the galaxy power spectrum at redshift $0.6<z<1.1$
Comments: 16 pages, 10 figures
Submitted: 2016-11-21, last modified: 2017-06-15
We use the final catalogue of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to measure the power spectrum of the galaxy distribution at high redshift, presenting results that extend beyond $z=1$ for the first time. We apply an FFT technique to four independent sub-volumes comprising a total of $51,728$ galaxies at $0.6<z<1.1$ (out of the nearly $90,000$ included in the whole survey). We concentrate here on the shape of the direction-averaged power spectrum in redshift space, explaining the level of modelling of redshift-space anisotropies and the anisotropic survey window function that are needed to deduce this in a robust fashion. We then use covariance matrices derived from a large ensemble of mock datasets in order to fit the spectral data. The results are well matched by a standard $\Lambda$CDM model, with density parameter $\Omega_M h =\smash{0.227^{+0.063}_{-0.050}}$ and baryon fraction $\smash{f_B=\Omega_B/\Omega_M=0.220^{+0.058}_{-0.072}}$. These inferences from the high-$z$ galaxy distribution are consistent with results from local galaxy surveys, and also with the Cosmic Microwave Background. Thus the $\Lambda$CDM model gives a good match to cosmic structure at all redshifts so far accessible to observational study.
[28]  [pdf] - 1579873
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Exploring the dependence of the three-point correlation function on stellar mass and luminosity at 0.5<z<1.1
Comments: 12 pages, 11 figures, 2 tables, accepted for publication in A&A
Submitted: 2016-03-29, last modified: 2017-05-22
The three-point correlation function (3PCF) is a powerful probe to investigate the clustering of matter in the Universe in a complementary way with respect to lower-order statistics, providing additional information with respect to the two-point correlation function and allowing us to shed light on biasing, nonlinear processes, and deviations from Gaussian statistics. In this paper, we analyse the first data release of the VIMOS Public Extragalactic Redshift Survey (VIPERS), determining the dependence of the three-point correlation function on luminosity and stellar mass at $z=[0.5,1.1]$. We exploit the VIPERS Public Data Release 1, consisting of more than 50,000 galaxies with B-band magnitudes in the range $-21.6\lesssim M_{\rm B}-5\log(h)\lesssim-19.9$ and stellar masses in the range $9.8\lesssim\log(M_\star[h^{-2}\,M_\odot])\lesssim 10.7$. We measure both the connected 3PCF and the reduced 3PCF in redshift space, probing different configurations and scales, in the range $2.5<r\,$[Mpc/h]$<20$. We find a significant dependence of the reduced 3PCF on scales and triangle shapes, with stronger anisotropy at larger scales ($r\sim10$ Mpc/h) and an almost flat trend at smaller scales, $r\sim2.5$ Mpc/h. Massive and luminous galaxies present a larger connected 3PCF, while the reduced 3PCF is remarkably insensitive to magnitude and stellar masses in the range we explored. These trends, already observed at low redshifts, are confirmed for the first time to be still valid up to $z=1.1$, providing support to the hierarchical scenario for which massive and bright systems are expected to be more clustered. The possibility of using the measured 3PCF to provide independent constraints on the linear galaxy bias $b$ has also been explored, showing promising results in agreement with other probes.
[29]  [pdf] - 1581830
Strong stellar-driven outflows shape the evolution of galaxies at cosmic dawn
Comments: 6 pages, 2 figures, accepted on ApJ Letters
Submitted: 2017-03-08, last modified: 2017-05-19
We study galaxy mass assembly and cosmic star formation rate (SFR) at high-redshift (z$\gt$4), by comparing data from multiwavelength surveys with predictions from the GAlaxy Evolution and Assembly (GAEA) model. GAEA implements a stellar feedback scheme partially based on cosmological hydrodynamical simulations, that features strong stellar driven outflows and mass-dependent timescale for the re-accretion of ejected gas. In previous work, we have shown that this scheme is able to correctly reproduce the evolution of the galaxy stellar mass function (GSMF) up to $z\sim3$. We contrast model predictions with both rest-frame Ultra-Violet (UV) and optical luminosity functions (LF), which are mostly sensible to the SFR and stellar mass, respectively. We show that GAEA is able to reproduce the shape and redshift evolution of both sets of LFs. We study the impact of dust on the predicted LFs and we find that the required level of dust attenuation is in qualitative agreement with recent estimates based on the UV continuum slope. The consistency between data and model predictions holds for the redshift evolution of the physical quantities well beyond the redshift range considered for the calibration of the original model. In particular, we show that GAEA is able to recover the evolution of the GSMF up to z$\sim$7 and the cosmic SFR density up to z$\sim$10.
[30]  [pdf] - 1581441
nIFTy Cosmology: the clustering consistency of galaxy formation models
Comments: 17 pages, 7 figures
Submitted: 2017-02-08, last modified: 2017-04-13
We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.
[31]  [pdf] - 1580672
The VIMOS Public Extragalactic Redshift Survey (VIPERS). The decline of cosmic star formation: quenching, mass, and environment connections
Comments: 20 pages, 13 figures, accepted to be published in A&A
Submitted: 2016-11-21, last modified: 2017-02-08
[Abridged] We use the final data of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to investigate the effect of environment on the evolution of galaxies between $z=0.5$ and $z=0.9$. We characterise local environment in terms of the density contrast smoothed over a cylindrical kernel, the scale of which is defined by the distance to the $5^{th}$ nearest neighbour. We find that more massive galaxies tend to reside in higher-density environments over the full redshift range explored. Defining star-forming and passive galaxies through their (NUV$-r$) vs ($r-K$) colours, we then quantify the fraction of star-forming over passive galaxies, $f_{\rm ap}$, as a function of environment at fixed stellar mass. $f_{\rm ap}$ is higher in low-density regions for galaxies with masses ranging from $\log(\mathcal{M}/\mathcal{M}_\odot)=10.38$ (the lowest value explored) to at least $\log(\mathcal{M}/\mathcal{M}_\odot)\sim11.3$, although with decreasing significance going from lower to higher masses. This is the first time that environmental effects on high-mass galaxies are clearly detected at redshifts as high as $z\sim0.9$. We compared these results to VIPERS-like galaxy mock catalogues based on the galaxy formation model of De Lucia & Blaizot. The model correctly reproduces $f_{\rm ap}$ in low-density environments, but underpredicts it at high densities. The discrepancy is particularly strong for the lowest-mass bins. We find that this discrepancy is driven by an excess of low-mass passive satellite galaxies in the model. Looking at the accretion history of these model galaxies, i.e. the times when they become satellites, a better (yet not perfect) agreement with observations can be obtained in high density regions by assuming either that a not-negligible fraction of satellites is destroyed, or that their quenching time-scale is longer than $\sim 2$ Gyr.
[32]  [pdf] - 1580934
The VIMOS Public Extragalactic Redshift Survey (VIPERS): The growth of structures at $0.5<z<1.2$ from redshift-space distortions in the clustering of the PDR-2 final sample
Comments: 17 pages, 19 figures
Submitted: 2016-12-16, last modified: 2017-01-26
We present measurements of the growth rate of cosmological structure from the modelling of the anisotropic galaxy clustering measured in the final data release of the VIPERS survey. The analysis is carried out in configuration space and based on measurements of the first two even multipole moments of the anisotropic galaxy auto-correlation function, in two redshift bins spanning the range $0.5 < z < 1.2$. We provide robust and cosmology-independent corrections for the VIPERS angular selection function, allowing recovery of the underlying clustering amplitude at the percent level down to the Mpc scale. We discuss several improvements on the non-linear modelling of redshift-space distortions (RSD) and perform detailed tests of a variety of approaches against a set of realistic VIPERS-like mock realisations. This includes using novel fitting functions to describe the velocity divergence and density power spectra $P_{\theta\theta}$ and $P_{\delta\theta}$ that appear in RSD models. These tests show that we are able to measure the growth rate with negligible bias down to separations of $5h^{-1}Mpc$. Interestingly, the application to real data shows a weaker sensitivity to the details of non-linear RSD corrections compared to mock results. We obtain consistent values for the growth rate times the matter power spectrum normalisation parameter of $f\sigma_8=0.55\pm 0.12$ and $0.40\pm0.11$ at effective redshifts of $z = 0.6$ and $z=0.86$ respectively. These results are in agreement with standard cosmology predictions assuming Einstein gravity in a $\Lambda \rm{CDM}$ background.
[33]  [pdf] - 1553949
The VIMOS Public Extragalactic Redshift Survey (VIPERS): PCA-based automatic cleaning and reconstruction of survey spectra
Comments: Find the VIPERS data release at
Submitted: 2016-12-06
Identifying spurious reduction artefacts in galaxy spectra is a challenge for large surveys. We present an algorithm for identifying and repairing residual spurious features in sky-subtracted galaxy spectra with application to the VIPERS survey. The algorithm uses principal component analysis (PCA) applied to the galaxy spectra in the observed frame to identify sky line residuals imprinted at characteristic wavelengths. We further model the galaxy spectra in the rest-frame using PCA to estimate the most probable continuum in the corrupted spectral regions, which are then repaired. We apply the method to 90,000 spectra from the VIPERS survey and compare the results with a subset where careful editing was performed by hand. We find that the automatic technique does an extremely good job in reproducing the time-consuming manual cleaning and does it in a uniform and objective manner across a large data sample. The mask data products produced in this work are released together with the VIPERS second public data release (PDR-2).
[34]  [pdf] - 1532617
AGN feedback and the origin of the $\alpha$ enhancement in early type galaxies - insights from the GAEA model
Comments: 5 pages, 3 figures, accepted for publication MNRAS Letters
Submitted: 2016-11-14, last modified: 2016-11-29
We take advantage of our recently published model for Galaxy Evolution and Assembly (GAEA) to study the origin of the observed correlation between [$\alpha$/Fe] and galaxy stellar mass. In particular, we analyse the role of radio mode AGN feedback, that recent work has identified as a crucial ingredient to reproduce observations. In GAEA, this process introduces the observed trend of star formation histories extending over shorter time-scales for more massive galaxies, but does not provide a sufficient condition to reproduce the observed $\alpha$ enhancements of massive galaxies. In the framework of our model, this is possible only assuming that any residual star formation is truncated for galaxies more massive than $10^{10.5}\,{\rm M}_{\odot}$. This results, however, in even shorter star formation time-scales for the most massive galaxies, that translate in total stellar metallicities significantly lower than observed. Our results demonstrate that (i) trends of [$\alpha$/Fe] ratios cannot be simply converted in relative time-scale indicators, and (ii) AGN feedback cannot explain alone the positive correlation between [$\alpha$/Fe] and galaxy mass/velocity dispersion. Reproducing simultaneously the mass-metallicity relation and the $\alpha$ enhancements observed poses a challenge for hierarchical models, unless more exotic solutions are adopted such as metal-rich winds or a variable IMF.
[35]  [pdf] - 1580673
The VIMOS Public Extragalactic Redshift Survey (VIPERS): Downsizing of the blue cloud and the influence of galaxy size on mass quenching over the last eight billion years
Comments: 20 pages, 15 figures. Submitted to A&A. References updated. Comments welcome
Submitted: 2016-11-21, last modified: 2016-11-24
We use the full VIPERS redshift survey in combination with SDSS-DR7 to explore the relationships between star-formation history (using d4000), stellar mass and galaxy structure, and how these relationships have evolved since z~1. We trace the extents and evolutions of both the blue cloud and red sequence, by fitting double Gaussians to the d4000 distribution of galaxies in narrow stellar mass bins, for four redshift intervals over 0<z<1. This reveals downsizing in star formation, as the high-mass limit of the blue cloud retreats steadily with time from M*~10^11.2 M_sun at z~0.9 to M*~10^10.7 M_sun by the present day. The number density of massive blue-cloud galaxies (M*>10^11 M_sun, d4000<1.55) drops sharply by a factor five between z~0.8 and z~0.5. These galaxies are becoming quiescent at a rate that largely matches the increase in the numbers of massive passive galaxies seen over this period. We examine the size-mass relation of blue cloud galaxies, finding that its high-mass boundary runs along lines of constant M*/r_e or equivalently inferred velocity dispersion. Larger galaxies can continue to form stars to higher stellar masses than smaller galaxies. As blue cloud galaxies approach this high-mass limit, they start to be quenched, their d4000 values increasing to push them towards the green valley. In parallel, their structures change, showing higher Sersic indices and central stellar mass densities. For these galaxies, bulge growth is necessary for them to reach the high-mass limit of the blue cloud and be quenched by internal mechanisms. The blue cloud galaxies that are being quenched at z~0.8 lie along the same size-mass relation as present day quiescent galaxies, and seem the likely progenitors of today's S0s.
[36]  [pdf] - 1622182
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Full spectroscopic data and auxiliary information release (PDR-2)
Comments: Submitted to Astronomy and Astrophysics; 14 pages, 14 figures. High resolution version of Fig. 14 available at
Submitted: 2016-11-21
We present the full public data release (PDR-2) of the VIMOS Public Extragalactic Redshift Survey (VIPERS), performed at the ESO VLT. We release redshifts, spectra, CFHTLS magnitudes and ancillary information (as masks and weights) for a complete sample of 86,775 galaxies (plus 4,732 other objects, including stars and serendipitous galaxies); we also include their full photometrically-selected parent catalogue. The sample is magnitude limited to i_AB < 22.5, with an additional colour-colour pre-selection devised as to exclude galaxies at z<0.5. This practically doubles the effective sampling of the VIMOS spectrograph over the range 0.5<z<1.2 (reaching 47% on average), yielding a final median local galaxy density close to 5 10^{-3} h^3 Mpc}^{-3}. The total area spanned by the final data set is ~ 23.5 deg^2, corresponding to 288 VIMOS fields with marginal overlaps, split over two regions within the CFHTLS-Wide W1 and W4 equatorial fields (at R.A. ~2 and ~22 hours, respectively). Spectra were observed at a resolution R=220, covering a wavelength range 5500-9500 Angstrom. Data reduction and redshift measurements were performed through a fully automated pipeline; all redshift determinations were then visually validated and assigned a quality flag. Measurements with a quality flag >= 2 are shown to have a confidence level of 96% or larger and make up 88% of all measured galaxy redshifts (76,552 out of 86,775), constituting the VIPERS prime catalogue for statistical investigations. For this sample the rms redshift error, estimated using repeated measurements of about 3,000 galaxies, is found to be sigma_z = 0.00054(1+z). All data are available at and on the ESO Archive.
[37]  [pdf] - 1580671
VIMOS Public Extragalactic Redshift Survey (VIPERS). The distinct build-up of dense and normal massive passive galaxies
Comments: 14 pages, 13 figures, submitted to A&A
Submitted: 2016-11-21
We use the final data from the VIPERS redshift survey to extract an unparalleled sample of more than 2000 massive M > 10^11 M_sun passive galaxies (MPGs) at redshift 0.5 < z < 1.0, based on their NUVrK colours. This enables us to investigate how the population of these objects was built up over cosmic time. We find that the evolution of the number density depends on the galaxy mean surface stellar mass density, Sigma. In particular, dense (Sigma > 2000 M_sun pc^-2) MPGs show a constant comoving number density over this redshift range, whilst this increases by a factor ~ 4 for the least dense objects, defined as having Sigma < 1000 M_sun pc^-2. We estimate stellar ages for the MPG population both fitting the Spectral Energy Distribution (SED) and through the D4000_n index, obtaining results in good agreement. Our findings are consistent with passive ageing of the stellar content of dense MPGs. We show that at any redshift the less dense MPGs are younger than dense ones and that their stellar populations evolve at a slower rate than predicted by passive evolution. This points to a scenario in which the overall population of MPGs was built up over the cosmic time by continuous addition of less dense galaxies: on top of an initial population of dense objects that passively evolves, new, larger, and younger MPGs continuously join the population at later epochs. Finally, we demonstrate that the observed increase in the number density of MPGs is totally accounted for by the observed decrease in the number density of correspondingly massive star forming galaxies (i.e. all the non-passive M > 10^11 M_sun objects). Such systems observed at z ~ 1 in VIPERS, therefore, represent the most plausible progenitors of the subsequent emerging class of larger MPGs.
[38]  [pdf] - 1580670
The VIMOS Public Extragalactic Redshift Survey: Measuring the growth rate of structure around cosmic voids
Comments: 14 pages, 12 figures
Submitted: 2016-11-21
We identified voids in the completed VIMOS Public Extragalactic Redshift Survey (VIPERS), using an algorithm based on searching for empty spheres. We measured the cross-correlation between the centres of voids and the complete galaxy catalogue. The cross-correlation function exhibits a clear anisotropy in both VIPERS fields (W1 and W4), which is characteristic of linear redshift space distortions. By measuring the projected cross-correlation and then deprojecting it we are able to estimate the undistorted cross-correlation function. We propose that given a sufficiently well measured cross-correlation function one should be able to measure the linear growth rate of structure by applying a simple linear Gaussian streaming model for the redshift space distortions (RSD). Our study of voids in 306 mock galaxy catalogues mimicking the VIPERS fields would suggest that VIPERS is capable of measuring $\beta$ with an error of around $25\%$. Applying our method to the VIPERS data, we find a value for the redshift space distortion parameter, $\beta = 0.423^{+0.104}_{-0.108}$, which given the bias of the galaxy population we use gives a linear growth rate of $f\sigma_8 = 0.296^{+0.075}_{-0.078}$ at $z = 0.727$. These results are consistent with values observed in parallel VIPERS analysis using standard techniques.
[39]  [pdf] - 1521475
The VIMOS Public Extragalactic Redshift Survey (VIPERS): Galaxy segregation inside filaments at $z \simeq 0.7$
Comments: 7 pages, 3 figures, accepted for publication in MNRAS
Submitted: 2016-11-21
We present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.
[40]  [pdf] - 1530717
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Star formation history of passive galaxies
Comments: 17 pages, 12 figures, 1 table, accepted to A&A
Submitted: 2016-05-18, last modified: 2016-11-17
We trace the evolution and the star formation history of passive galaxies, using a subset of the VIMOS Public Extragalactic Redshift Survey (VIPERS). We extracted from the VIPERS survey a sample of passive galaxies in the redshift range 0.4<z<1.0 and stellar mass range 10<$log(M_{star}/M_{\odot})$<12. The sample was selected using an evolving cut in the rest-frame U-V color distribution and additional quality-ensuring cuts. We use the stacked spectra to measure the 4000$\AA$ break (D4000) and the $H\delta$ Lick index ($H\delta_{A}$) with high precision. We compare the results with a grid of synthetic spectra to constrain the star formation epochs of these galaxies. We characterize the formation redshift-stellar mass relation for intermediate-redshift passive galaxies. We find that at $z\sim1$ stellar populations in low-mass passive galaxies are younger than in high-mass passive galaxies, similarly to what is observed at the present epoch. Over the full analyzed redshift and stellar mass range, the $D4000$ index increases with redshift, while $H\delta_{A}$ gets lower. This implies that the stellar populations are getting older with increasing stellar mass. Comparison to the spectra of passive galaxies in the SDSS survey shows that the shape of the relations of $D4000$, and $H\delta_{A}$ with stellar mass has not changed significantly with redshift. Assuming a single burst formation, this implies that high-mass passive galaxies formed their stars at $z_{form}\sim2$, while low-mass galaxies formed their main stellar population more recently, at $z_{form}\sim1$. The consistency of these results, obtained using two independent estimator of the formation redshift ($D4000$ and $H\delta_{A}$), further strengthens a scenario in which star formation proceeds from higher- to lower-mass systems as time passes, i.e. what has become known as the 'downsizing' picture.
[41]  [pdf] - 1531932
HI-selected Galaxies in Hierarchical Models of Galaxy Formation and Evolution
Comments: 20 pages, 15 figures, accepted for publication in MNRAS
Submitted: 2016-10-06, last modified: 2016-11-07
In this work, we study the basic statistical properties of HI-selected galaxies extracted from six different semi-analytic models, all run on the same cosmological N-body simulation. One model includes an explicit treatment for the partition of cold gas into atomic and molecular hydrogen. All models considered agree nicely with the measured HI mass function in the local Universe, with the measured scaling relations between HI and galaxy stellar mass, and with the predicted 2-point correlation function for HI rich galaxies. One exception is given by one model that predicts very little HI associated with galaxies in haloes above 10^12 Msun: we argue this is due to a too efficient radio-mode feedback for central galaxies, and to a combination of efficient stellar feedback and instantaneous stripping of hot gas for satellites. We demonstrate that treatment of satellite galaxies introduces large uncertainties at low HI masses. While models assuming non instantaneous stripping of hot gas tend to form satellite galaxies with HI masses slightly smaller than those of centrals with the same stellar mass, instantaneous gas stripping does not translate necessarily in lower HI masses. In fact, the adopted stellar feedback and star formation affect the satellites too. We analyze the relation between HI content and spin of simulated haloes: low spin haloes tend to host HI poor galaxies while high spin haloes are populated by galaxies in a wide range of HI mass. In our simulations, this is due to a correlation between the initial gas disk size and the halo spin.
[42]  [pdf] - 1528070
Variations of the initial mass function in semi-analytical models: implications for the mass assembly and the chemical enrichment of galaxies in the GAEA model
Comments: 14 pages, 8 figures, 2 tables, submitted to MNRAS
Submitted: 2016-06-06, last modified: 2016-10-07
In this work, we investigate the implications of the Integrated Galaxy-wide stellar Initial Mass Function (IGIMF) approach in the framework of the semi-analytic model GAEA (GAlaxy Evolution and Assembly), which features a detailed treatment of chemical enrichment and stellar feedback. The IGIMF provides an analytic description of the dependence of the stellar IMF shape on the rate of star formation in galaxies. We find that our model with a universal IMF predicts a rather flat [$\alpha$/Fe]-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of $\alpha$-enhancement with stellar mass, in agreement with previous studies. This is mainly due to the fact that massive galaxies are characterized by larger star formation rates at high-redshift, leading to stronger $\alpha$-enhancement with respect to low-mass galaxies. At the same time, the IGIMF hypothesis does not affect significantly the trend for shorter star formation timescales for more massive galaxies. We argue that in the IGIMF scenario the [$\alpha$/Fe] ratios are good tracers of the highest star formation events. The final stellar masses and mass-to-light-ratio of our model massive galaxies are larger than those estimated from the synthetic photometry assuming a universal IMF, providing a self-consistent interpretation of similar recent results, based on dynamical analysis of local early type galaxies.
[43]  [pdf] - 1488839
The fundamental plane of EDisCS galaxies - The effect of size evolution
Comments: Revised according to the corrections described in Saglia et al. 2016, "ERRATA: The fundamental plane of EDisCS galaxies. The effect of size evolution", A&A, in press
Submitted: 2010-09-03, last modified: 2016-09-29
We study the evolution of spectral early-type galaxies in clusters, groups and the field up to redshift 0.9 using the EDisCS dataset. We measure Re, Ie, and sigma for 154 cluster and 68 field galaxies. We study the evolution of the zero point of the fundamental plane (FP) and confirm results in the literature, but now also for the low cluster velocity dispersion regime. The mass-to-light ratio varies as Delta log M/L_B=(-0.54+-0.01)z=(-1.61+-0.01)log(1+z) in clusters, independent of their velocity dispersion. The evolution is stronger (Delta log M/L_B=(-0.76+-0.01)z=(-2.27+-0.03)log(1+z)) for field galaxies. The FP residuals correlate with galaxy mass and become progressively negative at low masses. The effect is visible at z>=0.7 for cluster galaxies and at z>=0.5 for field galaxies. We investigate the size evolution of our galaxy sample. We find that the half-luminosity radius for a galaxy with a dynamical or stellar mass of 2x10^11 Msol varies as (1+z)^{-1.0+-0.3} for both cluster and field galaxies. At the same time, stellar velocity dispersions grow with redshift, as (1+z)^{0.59+-0.10} at constant dynamical mass, and as (1+z)^{0.34+- 0.14} at constant stellar mass. The measured size evolution reduces to Re (1+z)^{-0.5+- 0.2} and sigma (1+z)^{0.41+-0.08}, at fixed dynamical masses, and Re (1+z)^{-0.68+-0.4} and sigma (1+z)^{0.19+-0.10}, at fixed stellar masses, when the progenitor bias (galaxies that locally are of spectroscopic early-type, but not very old, disappear from the EDisCS high-redshift sample; these galaxies tend to be large in size) is taken into account. Taken together, the variations in size and velocity dispersion imply that the luminosity evolution with redshift derived from the zero point of the FP is somewhat milder than that derived without taking these variations into account.
[44]  [pdf] - 1462530
Numerical simulations challenged on the prediction of massive subhalo abundance in galaxy clusters: the case of Abell 2142
Comments: 8 pages, 2 figures. Modified to match the version published in ApJL
Submitted: 2016-07-04, last modified: 2016-08-24
In this Letter we compare the abundance of member galaxies of a rich, nearby ($z=0.09$) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12) database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant ($\gtrsim 7$ sigma confidence level) smaller amount of massive (circular velocity above $200\,{\rm km\, s^{-1}}$) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 \citep{grillo2015} and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.
[45]  [pdf] - 1498198
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Measuring nonlinear galaxy bias at z~0.8
Comments: 20 pages, 15 figures
Submitted: 2014-06-25, last modified: 2016-07-07
We use the first release of the VImos Public Extragalactic Redshift Survey of galaxies (VIPERS) of ~50,000 objects to measure the biasing relation between galaxies and mass in the redshift range z=[0.5,1.1]. We estimate the 1-point distribution function [PDF] of VIPERS galaxies from counts in cells and, assuming a model for the mass PDF, we infer their mean bias relation. The reconstruction of the bias relation from PDFs is performed through a novel method that accounts for Poisson noise, redshift distortions, inhomogeneous sky coverage and other selection effects. With this procedure we constrain galaxy bias and its deviations from linearity down to scales as small as 4 Mpc/h and out to z=1.1. We detect small (~3%) but significant deviations from linear bias. The mean biasing function is close to linear in regions above the mean density. The mean slope of the biasing relation is a proxy to the linear bias parameter. It increases both with luminosity, in agreement with results of previous analyses, and with redshift. However, we detect a strong bias evolution only for z>0.9 in agreement with some, but not all, previous studies. We also detected a significant increase of the bias with the scale, from 4 to 8 Mpc/h, now seen for the first time out to z=1. The amplitude of nonlinearity depends on redshift, luminosity and on scales but no clear trend is detected. Thanks to the large cosmic volume probed by VIPERS we find that the mismatch between the previous estimates of bias at z~1 from zCOSMOS and VVDS-Deep galaxy samples is fully accounted for by cosmic variance. The results of our work confirm the importance of going beyond the over-simplistic linear bias hypothesis showing that non-linearities can be accurately measured through the applications of the appropriate statistical tools to existing datasets like VIPERS.
[46]  [pdf] - 1470664
Clustering-based redshift estimation: application to VIPERS/CFHTLS
Submitted: 2016-05-18, last modified: 2016-06-17
We explore the accuracy of the clustering-based redshift estimation proposed by M\'enard et al. (2013) when applied to VIPERS and CFHTLS real data. This method enables us to reconstruct redshift distributions from measurement of the angular clus- tering of objects using a set of secure spectroscopic redshifts. We use state of the art spectroscopic measurements with iAB < 22.5 from the VIMOS Public Extragalactic Redshift Survey (VIPERS) as reference population to infer the redshift distribution of galaxies from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) T0007 release. VIPERS provides a nearly representative sample to the flux limit iAB < 22.5 at redshift > 0.5 which allows us to test the accuracy of the clustering-based red- shift distributions. We show that this method enables us to reproduce the true mean color-redshift relation when both populations have the same magnitude limit. We also show that this technique allows the inference of redshift distributions for a population fainter than the one of reference and we give an estimate of the color-redshift mapping in this case. This last point is of great interest for future large redshift surveys which suffer from the need of a complete faint spectroscopic sample.
[47]  [pdf] - 1579982
The Frontier Fields Lens Modeling Comparison Project
Comments: 38 pages, 25 figures, submitted to MNRAS, version with full resolution images can be found at
Submitted: 2016-06-14
Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.
[48]  [pdf] - 1422193
Galaxy assembly, stellar feedback and metal enrichment: the view from the GAEA model
Comments: 30 pages, 19 figures, accepted for publication by MNRAS; note that corresponding new galaxy catalogues (FIRE model) will soon be made publicly available at
Submitted: 2015-12-14, last modified: 2016-06-13
One major problem of current theoretical models of galaxy formation is given by their inability to reproduce the apparently `anti-hierarchical' evolution of galaxy assembly: massive galaxies appear to be in place since $z\sim 3$, while a significant increase of the number densities of low mass galaxies is measured with decreasing redshift. In this work, we perform a systematic analysis of the influence of different stellar feedback schemes, carried out in the framework of GAEA, a new semi-analytic model of galaxy formation. It includes a self-consistent treatment for the timings of gas, metal and energy recycling, and for the chemical yields. We show this to be crucial to use observational measurements of the metallicity as independent and powerful constraints for the adopted feedback schemes. The observed trends can be reproduced in the framework of either a strong ejective or preventive feedback model. In the former case, the gas ejection rate must decrease significantly with cosmic time (as suggested by parametrizations of the cosmological `FIRE' simulations). Irrespective of the feedback scheme used, our successful models always imply that up to 60-70 per cent of the baryons reside in an `ejected' reservoir and are unavailable for cooling at high redshift. The same schemes predict physical properties of model galaxies (e.g. gas content, colour, age, and metallicity) that are in much better agreement with observational data than our fiducial model. The overall fraction of passive galaxies is found to be primarily determined by internal physical processes, with environment playing a secondary role.
[49]  [pdf] - 1418947
Radiative feedback and cosmic molecular gas: the role of different radiative sources
Comments: 19 pages on MNRAS
Submitted: 2016-05-31
We present results from multifrequency radiative hydrodynamical chemistry simulations addressing primordial star formation and related stellar feedback from various populations of stars, stellar energy distributions (SEDs) and initial mass functions. Spectra for massive stars, intermediate-mass stars and regular solar-like stars are adopted over a grid of 150 frequency bins and consistently coupled with hydrodynamics, heavy-element pollution and non-equilibrium species calculations. Powerful massive population III stars are found to be able to largely ionize H and, subsequently, He and He$^+$, causing an inversion of the equation of state and a boost of the Jeans masses in the early intergalactic medium. Radiative effects on star formation rates are between a factor of a few and 1 dex, depending on the SED. Radiative processes are responsible for gas heating and photoevaporation, although emission from soft SEDs has minor impacts. These findings have implications for cosmic gas preheating, primordial direct-collapse black holes, the build-up of "cosmic fossils" such as low-mass dwarf galaxies, the role of AGNi during reionization, the early formation of extended disks and angular-momentum catastrophe.
[50]  [pdf] - 1530716
The VIMOS Public Extragalactic Redshift Survey (VIPERS). The coevolution of galaxy morphology and colour to z~1
Comments: 18 pages, 13 figures
Submitted: 2016-05-18
We explore the evolution of the statistical distribution of galaxy morphological properties and colours over the redshift range $0.5<z<1$, combining high-quality imaging data from the CFHT Legacy Survey with the large number of redshifts and extended photometry from the VIPERS survey. Galaxy structural parameters are measured by fitting S\'ersic profiles to $i$-band images and then combined with absolute magnitudes, colours and redshifts, to trace the evolution in a multi-parameter space. We analyse, using a new method, the combination of colours and structural parameters of early- and late-type galaxies in luminosity--redshift space. We found that both the rest-frame colour distributions in the (U-B) vs. (B-V) plane and the S\'ersic index distributions are well fitted by a sum of two Gaussians, with a remarkable consistency of red-spheroidal and blue-disky galaxy populations, over the explored redshift ($0.5<z<1$) and luminosity ($-1.5<B-B_*<1.0$) ranges. The combination of the UBV rest-frame colour and S\'ersic index $n$ as a function of redshift and luminosity allows us to present the structure of early- and late-type galaxies and their evolution. We found that early type galaxies display only a slow change of their concentrations since $z\sim1$; it is already established by $z\sim1$ and depends much more strongly on their luminosities. In contrast, late-type galaxies get clearly more concentrated with cosmic time since $z\sim1$, with only little evolution in colour, which remains dependent mainly on their luminosity. This flipped luminosity (mass) and redshift dependence likely reflects different evolutionary tracks of early- and late-type galaxies before and after $z\sim1$. The combination of rest-frame colours and S\'ersic index $n$ as a function of redshift and luminosity leads to a precise statistical description of the structure of galaxies and their evolution.
[51]  [pdf] - 1396820
Disc colours in field and cluster spiral galaxies at 0.5 < z < 0.8
Comments: Accepted for publication in A&A
Submitted: 2016-01-20
We present a detailed study of the colours in late-type galaxy discs for ten of the EDisCS galaxy clusters with 0.5 < z < 0.8. Our cluster sample contains 172 spiral galaxies, and our control sample is composed of 96 field disc galaxies. We deconvolve their ground-based V and I images obtained with FORS2 at the VLT with initial spatial resolutions between 0.4 and 0.8 arcsec to achieve a final resolution of 0.1 arcsec with 0.05 arcsec pixels, which is close to the resolution of the ACS at the HST. After removing the central region of each galaxy to avoid pollution by the bulges, we measured the V-I colours of the discs. We find that 50% of cluster spiral galaxies have disc V-I colours redder by more than 1 sigma of the mean colours of their field counterparts. This is well above the 16% expected for a normal distribution centred on the field disc properties. The prominence of galaxies with red discs depends neither on the mass of their parent cluster nor on the distance of the galaxies to the cluster cores. Passive spiral galaxies constitute 20% of our sample. These systems are not abnormally dusty. They are are made of old stars and are located on the cluster red sequences. Another 24% of our sample is composed of galaxies that are still active and star forming, but less so than galaxies with similar morphologies in the field. These galaxies are naturally located in the blue sequence of their parent cluster colour-magnitude diagrams. The reddest of the discs in clusters must have stopped forming stars more than ~5 Gyr ago. Some of them are found among infalling galaxies, suggesting preprocessing. Our results confirm that galaxies are able to continue forming stars for some significant period of time after being accreted into clusters, and suggest that star formation can decline on seemingly long (1 to 5 Gyr) timescales.
[52]  [pdf] - 1341300
Semi-Analytic Model Predictions of the Galaxy Population in Proto-clusters
Comments: 14 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2015-12-04
We investigate the galaxy population in simulated proto-cluster regions using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations. We select the most massive clusters at redshift $z=0$ from our set of simulations, and follow their main progenitors back in time. The analysis shows that proto-cluster regions are dominated by central galaxies and their number decreases with time as many become satellites, clustering around the central object. In agreement with observations, we find an increasing velocity dispersion with cosmic time, the increase being faster for satellites. The analysis shows that proto-clusters are very extended regions, $\gtrsim 20 \, Mpc$ at $z \gtrsim 1$. The fraction of galaxies in proto-cluster regions that are not progenitor of cluster galaxies varies with redshift, stellar mass and area considered. It is about 20-30 per cent for galaxies with stellar mass $\sim 10^9\,{\rm M}_{\sun}$, while negligible for the most massive galaxies considered. Nevertheless, these objects have properties similar to those of progenitors. We investigate the building-up of the passive-sequence in clusters, and find that their progenitors are on average always active at any redshift of interest of proto-clusters. The main mechanism which quenches their star formation is the removal of the hot gas reservoir at the time of accretion. The later galaxies are accreted (become satellite), and the more the cold gas available, the longer the time spent as active. Central galaxies are active over all redshift range considered, although a non-negligible fraction of them become passive at redshift $z<1$, due to strong feedback from Active Galactic Nuclei.
[53]  [pdf] - 1347510
Mergers and Star Formation: The environment and Stellar Mass Growth of the Progenitors of Ultra-Massive Galaxies since z = 2
Comments: 18 pages, 14 figures, accepted for publication in ApJ
Submitted: 2015-09-01, last modified: 2015-11-19
The growth of galaxies is a key problem in understanding the structure and evolution of the universe. Galaxies grow their stellar mass by a combination of star formation and mergers, with a relative importance that is redshift dependent. Theoretical models predict quantitatively different contributions from the two channels; measuring these from the data is a crucial constraint. Exploiting the UltraVISTA catalog and a unique sample of progenitors of local ultra massive galaxies selected with an abundance matching approach, we quantify the role of the two mechanisms from z=2 to 0. We also compare our results to two independent incarnations of semi-analytic models. At all redshifts, progenitors are found in a variety of environments, ranging from being isolated to having 5-10 companions with mass ratio at least 1:10 within a projected radius of 500 kpc. In models, progenitors have a systematically larger number of companions, entailing a larger mass growth for mergers than in observations, at all redshifts. Generally, in both observations and models, the inferred and the expected mass growth roughly agree, within the uncertainties. Overall, our analysis confirms the model predictions, showing how the growth history of massive galaxies is dominated by in situ star formation at z~2, both star-formation and mergers at 1<z<2, and by mergers alone at z<1. Nonetheless, detailed comparisons still point out to tensions between the expected mass growth and our results, which might be due to either an incorrect progenitors-descendants selection, uncertainties on star formation rate and mass estimates, or the adopted assumptions on merger rates.
[54]  [pdf] - 1347593
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Environmental effects shaping the galaxy stellar mass function
Comments: A&A accepted; 20 pages, 12 figures, 2 tables
Submitted: 2015-11-03
We exploit the first public data release of VIPERS to investigate environmental effects in galaxy evolution between $z\sim0.5$ and $0.9$. The large number of spectroscopic redshifts over an area of about $10\,\mathrm{deg}^2$ provides a galaxy sample with high statistical power. The accurate redshift measurements, with $\sigma_z = 0.00047(1+z_\mathrm{spec})$, allow us to robustly isolate galaxies living in the lowest- and highest-density environments, as defined in terms of spatial 3D density contrast. We estimate the stellar mass function (SMF) of galaxies residing in these two environments, and constrain its high-mass end with unprecedented precision. We find that the galaxy SMF in the densest regions has a different shape than that measured at low densities, with an enhancement of massive galaxies and a hint of a flatter (less negative) slope at $z<0.8$. We normalise each SMF to the comoving volume occupied by the corresponding environment, and relate estimates from different redshift bins. We observe an evolution of the SMF of VIPERS galaxies in high densities, while the low-density one is nearly constant. We compare these results to semi-analytical models and find consistent environmental signatures. We discuss how the halo mass function and fraction of central/satellite galaxies depend on the environments considered, making intrinsic and environmental properties of galaxies physically coupled, and therefore difficult to disentangle. The evolution of our low-density regions is well described by the formalism introduced by Peng et al.~(2010), and is consistent with the idea that galaxies become progressively passive because of internal physical processes. The same formalism could also describe the evolution of the SMF in the high density regions, but only if a significant contribution from dry mergers is considered. [Abridged]
[55]  [pdf] - 1342878
CLASH-VLT: Environment-driven evolution of galaxies in the z=0.209 cluster Abell 209
Comments: 17 pages, 20 figures, A&A in press
Submitted: 2015-10-19
The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.
[56]  [pdf] - 1249163
On the dependence of galaxy morphologies on galaxy mergers
Comments: MNRAS in press, 11 pages, 5 figures
Submitted: 2015-07-16
The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass.
[57]  [pdf] - 1125363
Star Formation in Herschel's Monsters versus Semi-Analytic Models
Comments: 9 pages; 3 figures
Submitted: 2015-06-04
We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.
[58]  [pdf] - 1546192
A high-dimensional look at VIPERS galaxies
Comments: 3 pages. To appear in Statistical Challenges in 21st Century Cosmology, Proceedings IAU Symposium No. 306, A. H. Heavens & J.-L. Starck, eds
Submitted: 2015-05-23
We investigate how galaxies in VIPERS (the VIMOS Public Extragalactic Redshift Survey) inhabit the cosmological density field by examining the correlations across the observable parameter space of galaxy properties and clustering strength. The high-dimensional analysis is made manageable by the use of group-finding and regression tools. We find that the major trends in galaxy properties can be explained by a single parameter related to stellar mass. After subtracting this trend, residual correlations remain between galaxy properties and the local environment pointing to complex formation dependencies. As a specific application of this work we build subsamples of galaxies with specific clustering properties for use in cosmological tests.
[59]  [pdf] - 1300189
The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
Comments: Submitted to Astronomy & Astrophysics
Submitted: 2015-05-23
Aims. Using the VIMOS Public Extragalactic Redshift Survey (VIPERS) we aim to jointly estimate the key parameters that describe the galaxy density field and its spatial correlations in redshift space. Methods. We use the Bayesian formalism to jointly reconstruct the redshift-space galaxy density field, power spectrum, galaxy bias and galaxy luminosity function given the observations and survey selection function. The high-dimensional posterior distribution is explored using the Wiener filter within a Gibbs sampler. We validate the analysis using simulated catalogues and apply it to VIPERS data taking into consideration the inhomogeneous selection function. Results. We present joint constraints on the anisotropic power spectrum as well as the bias and number density of red and blue galaxy classes in luminosity and redshift bins as well as the measurement covariances of these quantities. We find that the inferred galaxy bias and number density parameters are strongly correlated although these are only weakly correlated with the galaxy power spectrum. The power spectrum and redshift-space distortion parameters are in agreement with previous VIPERS results with the value of the growth rate $f\sigma_8 = 0.38$ with 18% uncertainty at redshift 0.7.
[60]  [pdf] - 1239453
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Hierarchical scaling and biasing
Comments: 20 pages, 14 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2015-05-20
We investigate the higher-order correlation properties of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to test the hierarchical scaling hypothesis at z~1 and the dependence on galaxy luminosity, stellar mass, and redshift. We also aim to assess deviations from the linearity of galaxy bias independently from a previously performed analysis of our survey (Di Porto et al. 2014). We have measured the count probability distribution function in cells of radii 3 < R < 10 Mpc/h, deriving $\sigma_{8g}$, the volume-averaged two-,three-,and four-point correlation functions and the normalized skewness $S_{3g}$ and kurtosis $S_{4g}$ for volume-limited subsamples covering the ranges $-19.5 \le M_B(z=1.1)-5log(h) \le -21.0$, $9.0 < log(M*/M_{\odot} h^{-2}) \le 11.0$, $0.5 \le z < 1.1$. We have thus performed the first measurement of high-order correlations at z~1 in a spectroscopic redshift survey. Our main results are the following. 1) The hierarchical scaling holds throughout the whole range of scale and z. 2) We do not find a significant dependence of $S_{3g}$ on luminosity (below z=0.9 $S_{3g}$ decreases with luminosity but only at 1{\sigma}-level). 3) We do not detect a significant dependence of $S_{3g}$ and $S_{4g}$ on scale, except beyond z~0.9, where the dependence can be explained as a consequence of sample variance. 4) We do not detect an evolution of $S_{3g}$ and $S_{4g}$ with z. 5) The linear bias factor $b=\sigma_{8g}/\sigma_{8m}$ increases with z, in agreement with previous results. 6) We quantify deviations from the linear bias by means of the Taylor expansion parameter $b_2$. Our results are compatible with a null non-linear bias term, but taking into account other available data we argue that there is evidence for a small non-linear bias term.
[61]  [pdf] - 1037708
nIFTy Cosmology: Comparison of Galaxy Formation Models
Comments: 35 pages, 20 figures, accepted for publication in MNRAS
Submitted: 2015-05-18
We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.
[62]  [pdf] - 1378565
The VIMOS Public Extragalactic Redshift Survey (VIPERS): On the correct recovery of the count-in-cell probability distribution function
Comments: 14 pages, 11 figures, 2 tables
Submitted: 2015-05-03
We compare three methods to measure the count-in-cell probability density function of galaxies in a spectroscopic redshift survey. From this comparison we found that when the sampling is low (the average number of object per cell is around unity) it is necessary to use a parametric method to model the galaxy distribution. We used a set of mock catalogues of VIPERS, in order to verify if we were able to reconstruct the cell-count probability distribution once the observational strategy is applied. We find that in the simulated catalogues, the probability distribution of galaxies is better represented by a Gamma expansion than a Skewed Log-Normal. Finally, we correct the cell-count probability distribution function from the angular selection effect of the VIMOS instrument and study the redshift and absolute magnitude dependency of the underlying galaxy density function in VIPERS from redshift $0.5$ to $1.1$. We found very weak evolution of the probability density distribution function and that it is well approximated, independently from the chosen tracers, by a Gamma distribution.
[63]  [pdf] - 974530
CLASH-VLT: Substructure in the galaxy cluster MACS J1206.2-0847 from kinematics of galaxy populations
Comments: A&A accepted, 19 pages, 30 figures, minor language changes
Submitted: 2015-03-18, last modified: 2015-04-02
In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.
[64]  [pdf] - 1450548
Galaxy sizes as a function of environment at intermediate redshift from the ESO Distant Cluster Survey
Comments: 11 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2015-03-27
In order to assess whether the environment has a significant effect on galaxy sizes, we compare the mass--size relations of cluster and field galaxies in the $0.4 < z < 0.8$ redshift range from the ESO Distant Cluster Survey (EDisCS) using HST images. We analyse two mass-selected samples, one defined using photometric redshifts ($10.2 \le \log M_\ast/M_{\odot} \le 12.0$), and a smaller more robust subsample using spectroscopic redshifts ($10.6 \le \log M_\ast/M_{\odot} \le 11.8$). We find no significant difference in the size distributions of cluster and field galaxies of a given morphology. Similarly, we find no significant difference in the size distributions of cluster and field galaxies of similar rest-frame $B-V$ colours. We rule out average size differences larger than $10$--$20$\% in both cases. Consistent conclusions are found with the spectroscopic and photometric samples. These results have important consequences for the physical process(es) responsible for the size evolution of galaxies, and in particular the effect of the environment. The remarkable growth in galaxy size observed from $z\sim2.5$ has been reported to depend on the environment at higher redshifts ($z>1$), with early-type/passive galaxies in higher density environments growing earlier. Such dependence disappears at lower redshifts. Therefore, if the reported difference at higher-$z$ is real, the growth of field galaxies has caught up with that of cluster galaxies by $z\sim1$. Any putative mechanism responsible for galaxy growth has to account for the existence of environmental differences at high redshift and their absence (or weakening) at lower redshifts.
[65]  [pdf] - 1223718
The early phases of galaxy clusters formation in IR: coupling hydrodynamical simulations with GRASIL3D
Comments: 14 pages, 11 figures, MNRAS in press. Minor editorial improvements
Submitted: 2014-12-18, last modified: 2015-03-25
We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.
[66]  [pdf] - 1224382
The galaxy-halo connection from a joint lensing, clustering and abundance analysis in the CFHTLenS/VIPERS field
Comments: 31 pages, 18 figures, 4 table. Accepted for publication in MNRAS. Online material available at
Submitted: 2015-02-10
We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies.
[67]  [pdf] - 1215992
CLASH-VLT: Insights on the mass substructures in the Frontier Fields Cluster MACS J0416.1-2403 through accurate strong lens modeling
Comments: 26 pages, 22 figures, 7 tables; accepted for publication in the Astrophysical Journal. A high-resolution version is available at
Submitted: 2014-07-29, last modified: 2014-12-22
We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.
[68]  [pdf] - 1222827
The definition of environment and its relation to the quenching of galaxies at z=1-2 in a hierarchical Universe
Comments: 19 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2014-10-27
A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density.
[69]  [pdf] - 1215993
Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model
Comments: 20 pages, 13 figures, minor changes to match version in press in MNRAS
Submitted: 2014-07-29, last modified: 2014-09-29
We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.
[70]  [pdf] - 1215806
The influence of the environmental history on quenching star formation in a $\Lambda$CDM universe
Comments: 24 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2014-07-21, last modified: 2014-09-12
We present a detailed analysis of the influence of the environment and of the environmental history on quenching star formation in central and satellite galaxies in the local Universe. We take advantage of publicly available galaxy catalogues obtained from applying a galaxy formation model to the Millennium simulation. In addition to halo mass, we consider the local density of galaxies within various fixed scales. Comparing our model predictions to observational data (SDSS), we demonstrate that the models are failing to reproduce the observed density dependence of the quiescent galaxy fraction in several aspects: for most of the stellar mass ranges and densities explored, models cannot reproduce the observed similar behaviour of centrals and satellites, they slightly under-estimate the quiescent fraction of centrals and significantly over-estimate that of satellites. We show that in the models, the density dependence of the quiescent central galaxies is caused by a fraction of "backsplash" centrals which have been satellites in the past (and were thus suffering from environmental processes). Turning to satellite galaxies, the density dependence of their quiescent fractions reflects a dependence on the time spent orbiting within a parent halo of a particular mass, correlating strongly with halo mass and distance from the halo centre. Comparisons with observational estimates suggest relatively long gas consumption time scales of roughly 5 Gyr in low mass satellite galaxies. The quenching time scales decrease with increasing satellite stellar mass. Overall, a change in modelling both internal processes (star formation and feedback) and environmental processes (e.g. making them dependent on dynamical friction time-scales and preventing the re-accretion of gas onto backsplash galaxies) is required for improving currently used galaxy formation models.
[71]  [pdf] - 891981
The VIMOS Public Extragalactic Redshift Survey - Searching for Cosmic Voids
Submitted: 2014-07-10, last modified: 2014-09-12
The characterisation of cosmic voids gives unique information about the large-scale distribution of galaxies, their evolution and the cosmological model. We identify and characterise cosmic voids in the VIMOS Public Extragalactic Redshift Survey (VIPERS) at redshift 0.55 < z < 0.9. A new void search method is developed based upon the identification of empty spheres that fit between galaxies. The method can be used to characterise the cosmic voids despite the presence of complex survey boundaries and internal gaps. We investigate the impact of systematic observational effects and validate the method against mock catalogues. We measure the void size distribution and the void-galaxy correlation function. We construct a catalogue of voids in VIPERS. The distribution of voids is found to agree well with the distribution of voids found in mock catalogues. The void-galaxy correlation function shows indications of outflow velocity from the voids.
[72]  [pdf] - 1216044
What Regulates Galaxy Evolution? Open Questions in Our Understanding of Galaxy Formation and Evolution
Comments: 36 pages, 6 Figures, submitted to New Astronomy Reviews
Submitted: 2014-07-31
In April 2013, a workshop entitled "What Regulates Galaxy Evolution" was held at the Lorentz Center. The aim of the workshop was to bring together the observational and theoretical community working on galaxy evolution, and to discuss in depth of the current problems in the subject, as well as to review the most recent observational constraints. A total of 42 astrophysicists attended the workshop. A significant fraction of the time was devoted to identifying the most interesting "open questions" in the field, and to discuss how progress can be made. This review discusses the four questions (one for each day of the workshop) that, in our opinion, were the focus of the most intense debate. We present each question in its context, and close with a discussion of what future directions should be pursued in order to make progress on these problems.
[73]  [pdf] - 1208933
Simulating the evolution of disc galaxies in a group environment. II. The influence of close-encounters between galaxies
Comments: 16 pages, 6 figures, 4 tables. Shortened abstract and minor additions. Accepted for publication in MNRAS
Submitted: 2014-04-10, last modified: 2014-06-25
We study the evolution of disc galaxies in group environments under the effect of both the global tidal field and close-encounters between galaxies, using controlled N-body simulations of isolated mergers. We find that close-range encounters between galaxies are less frequent and less damaging to disc galaxies than originally expected, since they mostly occur when group members have lost a significant fraction of their initial mass to tidal stripping. We also find that group members mostly affect disc galaxies indirectly by modifying their common global tidal field. Different initial orbital parameters of group members introduce a significant "scatter" in the evolution of general properties of disc galaxies around a "median" evolution that is similar to when only the effect of the global tidal field is included. Close-encounters introduce a high variability in the properties of disc galaxies, even slowing their evolution in some cases, and could wash out correlations between galaxy properties and the group total mass. The combined effect of the global tidal field and close-encounters appears to be inefficient at forming/enhancing central stellar bulges. This implies that bulges of S0 galaxies should be mostly composed by young stars, which is consistent with recent observations.
[74]  [pdf] - 1209163
What do simulations predict for the galaxy stellar mass function and its evolution in different environments?
Comments: ApJ in press
Submitted: 2014-04-24
We present a comparison between the observed galaxy stellar mass function and the one predicted from the De Lucia & Blaizot (2007) semi-analytic model applied to the Millennium Simulation, for cluster satellites and galaxies in the field (meant as a wide portion of the sky, including all environments), in the local universe (z~0.06) and at intermediate redshift (z~0.6), with the aim to shed light on the processes which regulate the mass distribution in different environments. While the mass functions in the field and in its finer environments (groups, binary and single systems) are well matched in the local universe down to the completeness limit of the observational sample, the model over-predicts the number of low mass galaxies in the field at z~0.6 and in clusters at both redshifts. Above M_*=10^10.25 M_sun, it reproduces the observed similarity of the cluster and field mass functions, but not the observed evolution. Our results point out two shortcomings of the model: an incorrect treatment of cluster-specific environmental effects and an over-efficient galaxy formation at early times (as already found by e.g. Weinmann et al. 2012). Next, we consider only simulations. Using also the Guo et al. (2011) model, we find that the high mass end of the mass functions depends on halo mass: only very massive halos host massive galaxies, with the result that their mass function is flatter. Above M_*=10^9.4 M_sun, simulations show an evolution in the number of the most massive galaxies in all the environments. Mass functions obtained from the two prescriptions are different, however results are qualitatively similar, indicating that the adopted recipes to model the evolution of central and satellite galaxies still have to be better implemented in semi-analytic models.
[75]  [pdf] - 1208783
A semi-analytic model comparison: testing cooling models against hydrodynamical simulations
Comments: 20 pages, 12 figures, accepted by MNRAS
Submitted: 2014-04-03
We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.
[76]  [pdf] - 1202962
The VIMOS Public Extragalactic Redshift Survey (VIPERS): A quiescent formation of massive red-sequence galaxies over the past 9 Gyr
Comments: 27 pages, 21 figures, accepted for publication in A&A
Submitted: 2014-01-23
We explore the evolution of the Colour-Magnitude Relation (CMR) and Luminosity Function (LF) at 0.4<z<1.3 from the VIMOS Public Extragalactic Redshift Survey (VIPERS) using ~45,000 galaxies with precise spectroscopic redshifts down to i'_AB<22.5 over ~10.32 deg^2 in two fields. From z=0.5 to z=1.3 the LF and CMR are well defined for different galaxy populations and M^*_B evolves by ~1.04(1.09)+/-0.06(0.10) mag for the total (red) galaxy sample. We compare different criteria for selecting early-type galaxies (ETGs): (1) fixed cut in rest-frame (U-V) colours, (2) evolving cut in (U-V) colours, (3) rest-frame (NUV-r')-(r'-K) colour selection, and (4) SED classification. Regardless of the method we measure a consistent evolution of the red-sequence (RS). Between 0.4<z<1.3 we find a moderate evolution of the RS intercept of Delta(U-V)=0.28+/-0.14 mag, favouring exponentially declining star formation (SF) histories with SF truncation at 1.7<=z<=2.3. Together with the rise in the ETG number density by 0.64 dex since z=1, this suggests a rapid build-up of massive galaxies (M>10^11 M_sun) and expeditious RS formation over a short period of ~1.5 Gyr starting before z=1. This is supported by the detection of ongoing SF in ETGs at 0.9<z<1.0, in contrast with the quiescent red stellar populations of ETGs at 0.5<z<0.6. There is an increase in the observed CMR scatter with redshift, two times larger than in galaxy clusters and at variance with theoretical models. We discuss possible physical mechanisms that support the observed evolution of the red galaxy population. Our findings point out that massive galaxies have experienced a sharp SF quenching at z~1 with only limited additional merging. In contrast, less-massive galaxies experience a mix of SF truncation and minor mergers which build-up the low- and intermediate-mass end of the CMR.
[77]  [pdf] - 791740
The VIMOS Public Extragalactic Redshift Survey (VIPERS): $\Omega_{\rm m_0}$ from the galaxy clustering ratio measured at $z \sim 1$
Comments: 18 pages, 17 figures, accepted for publication in A&A, references added
Submitted: 2013-10-12, last modified: 2013-12-12
We use a sample of about 22,000 galaxies at $0.65<z<1.2$ from the VIPERS PDR-1 catalogue, to constrain the cosmological model through a measurement of the galaxy {\it clustering ratio} $\eta_{g,R}$. This statistic has favourable properties, being defined as the ratio of two quantities characterizing the smoothed density field in spheres of given radius $R$: the value of its correlation function on a multiple of this scale, $\xi(nR)$, and its variance $\sigma^2(R)$. For sufficiently large values of $R$, this is a universal number, capturing 2-point clustering information independently of the linear bias and linear redshift-space distortions of the specific galaxy tracers. In this paper we discuss in detail how to extend the application of $\eta_{g,R}$ to quasi-linear scales and how to control and remove observational selection effects which are typical of redshift surveys as VIPERS. We verify the accuracy and efficiency of these procedures using mock catalogues that match the survey selection process. These results evidence the robustness of $\eta_{g,R}$ to non-linearities and observational effects, which is related to its very definition as a ratio of quantities that are similarly affected. We measure $\Omega_{m,0}=0.270_{-0.025}^{+0.029}$. In addition to the great precision achieved on our estimation of $\Omega_m$ using VIPERS PDR-1, this result is remarkable because it appears to be in good agreement with a recent estimate $z\simeq 0.3$, obtained applying the same technique to the SDSS-LRG catalogue. It, therefore, suports the robustness of the present analysis. Moreover, the combination of these two measurements at $z\sim 0.3$ and $z\sim 0.9$ provides us with a very precise estimate $\Omega_{m,0}=0.274\pm0.017$ which highlights the great consistency between our estimation and other cosmological probes such as BAOs, CMB and Supernovae.
[78]  [pdf] - 1179699
The VIMOS Public Extragalactic Survey (VIPERS): First Data Release of 57 204 spectroscopic measurements
Comments: Accepted for publication in A&A Added and/or replaced some figure, added section on DataBase interface, expaned Introduction
Submitted: 2013-10-03, last modified: 2013-12-10
We present the first Public Data Release (PDR-1) of the VIMOS Public Extragalactic Survey (VIPERS). It comprises 57 204 spectroscopic measurements together with all additional information necessary for optimal scientific exploitation of the data, in particular the associated photometric measurements and quantification of the photometric and survey completeness. VIPERS is an ESO Large Programme designed to build a spectroscopic sample of ' 100 000 galaxies with iAB < 22.5 and 0.5 < z < 1.5 with high sampling rate (~45%). The survey spectroscopic targets are selected from the CFHTLS-Wide five-band catalogues in the W1 and W4 fields. The final survey will cover a total area of nearly 24 deg2, for a total comoving volume between z = 0.5 and 1.2 of ~4x10^7 h^(-3)Mpc^3 and a median galaxy redshift of z~0.8. The release presented in this paper includes data from virtually the entire W4 field and nearly half of the W1 area, thus representing 64% of the final dataset. We provide a detailed description of sample selection, observations and data reduction procedures; we summarise the global properties of the spectroscopic catalogue and explain the associated data products and their use, and provide all the details for accessing the data through the survey database ( where all information can be queried interactively.
[79]  [pdf] - 1179731
Characterizing Diffused Stellar Light in simulated galaxy clusters
Comments: 18 Pages, 15 figures. Matches to MNRAS published version
Submitted: 2013-10-04, last modified: 2013-11-11
[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.
[80]  [pdf] - 1180600
On the formation and physical properties of the Intra-Cluster Light in hierarchical galaxy formation models
Comments: Accepted for Publication in MNRAS, 19 pages, 13 figures, 1 table
Submitted: 2013-11-08
We study the formation of the Intra-Cluster Light (ICL) using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations of groups and clusters. We assume that the ICL forms by (1) stellar stripping of satellite galaxies and (2) relaxation processes that take place during galaxy mergers. The fraction of ICL in groups and clusters predicted by our models ranges between 10 and 40 per cent, with a large halo-to-halo scatter and no halo mass dependence. We note, however, that our predicted ICL fractions depend on the resolution: for a set of simulations with particle mass one order of magnitude larger than that adopted in the high resolution runs used in our study, we find that the predicted ICL fractions are ~30-40 per cent larger than those found in the high resolution runs. On cluster scale, large part of the scatter is due to a range of dynamical histories, while on smaller scale it is driven by individual accretion events and stripping of very massive satellites, $M_{*} \gtrsim 10^{10.5} M_{\odot}$, that we find to be the major contributors to the ICL. The ICL in our models forms very late (below $z\sim 1$), and a fraction varying between 5 and 25 per cent of it has been accreted during the hierarchical growth of haloes. In agreement with recent observational measurements, we find the ICL to be made of stars covering a relatively large range of metallicity, with the bulk of them being sub-solar.
[81]  [pdf] - 1180115
Searching for large-scale structures around high-redshift radio galaxies with Herschel
Comments: 13 pages, 9 figures; Accepted for publication in MNRAS
Submitted: 2013-10-21
This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.
[82]  [pdf] - 725633
The VIMOS Public Extragalactic Redshift Survey (VIPERS). A precise measurement of the galaxy stellar mass function and the abundance of massive galaxies at redshifts 0.5<z<1.3
Comments: 19 pages, 14 figures, 2 tables. A&A in press (accepted for publication on July 12th)
Submitted: 2013-03-15, last modified: 2013-09-02
We measure the evolution of the galaxy stellar mass function from z=1.3 to z=0.5 using the first 53,608 redshifts of the ongoing VIMOS Public Extragalactic Survey (VIPERS). We estimate the galaxy stellar mass function at several epochs discussing in detail the amount of cosmic variance affecting our estimate. We find that Poisson noise and cosmic variance of the galaxy mass function in the VIPERS survey are comparable with the statistical uncertainties of large surveys in the local universe. VIPERS data allow us to determine with unprecedented accuracy the high-mass tail of the galaxy stellar mass function, which includes a significant number of galaxies that are usually too rare to detect with any of the past spectroscopic surveys. At the epochs sampled by VIPERS, massive galaxies had already assembled most of their stellar mass. We apply a photometric classification in the (U-V) rest-frame colour to compute the mass function of blue and red galaxies, finding evidence for the evolution of their contribution to the total number density budget: the transition mass above which red galaxies dominate is found to be about 10^10.4 M_sun at z=0.55 and evolves proportionally to (1+z)^3. We are able to trace separately the evolution of the number density of blue and red galaxies with masses above 10^11.4 M_sun, in a mass range barely studied in previous work. We find that for such large masses, red galaxies show a milder evolution with redshift, when compared to objects at lower masses. At the same time, we detect a population of similarly massive blue galaxies, which are no longer detectable below z=0.7. These results show the improved statistical power of VIPERS data, and give initial promising indications of mass-dependent quenching of galaxies at z~1. [Abridged]
[83]  [pdf] - 711290
The VIMOS Public Extragalactic Redshift Survey (VIPERS). A Support Vector Machine classification of galaxies, stars and AGNs
Comments: 14 pages, 11 figures, 9 tables, published in A&A
Submitted: 2013-03-11, last modified: 2013-08-27
The aim of this work is to develop a comprehensive method for classifying sources in large sky surveys and we apply the techniques to the VIMOS Public Extragalactic Redshift Survey (VIPERS). Using the optical (u*, g', r', i') and NIR data (z', Ks), we develop a classifier, based on broad-band photometry, for identifying stars, AGNs and galaxies improving the purity of the VIPERS sample. Support Vector Machine (SVM) supervised learning algorithms allow the automatic classification of objects into two or more classes based on a multidimensional parameter space. In this work, we tailored the SVM for classifying stars, AGNs and galaxies, and applied this classification to the VIPERS data. We train the SVM using spectroscopically confirmed sources from the VIPERS and VVDS surveys. We tested two SVM classifiers and concluded that including NIR data can significantly improve the efficiency of the classifier. The self-check of the best optical + NIR classifier has shown a 97% accuracy in the classification of galaxies, 97 for stars, and 95 for AGNs in the 5-dimensional colour space. In the test on VIPERS sources with 99% redshift confidence, the classifier gives an accuracy equal to 94% for galaxies, 93% for stars, and 82% for AGNs. The method was applied to sources with low quality spectra to verify their classification, and thus increasing the security of measurements for almost 4 900 objects. We conclude that the SVM algorithm trained on a carefully selected sample of galaxies, AGNs, and stars outperforms simple colour-colour selection methods, and can be regarded as a very efficient classification method particularly suitable for modern large surveys.
[84]  [pdf] - 1165170
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Luminosity and stellar mass dependence of galaxy clustering at 0.5<z<1.1
Comments: 17 pages, 10 figures, 4 tables. Accepted for publication in A&A
Submitted: 2013-03-11, last modified: 2013-07-20
We investigate the dependence of galaxy clustering on luminosity and stellar mass in the redshift range 0.5<z<1.1, using the first ~55000 redshifts from the VIMOS Public Extragalactic Redshift Survey (VIPERS). We measured the redshift-space two-point correlation functions (2PCF), and the projected correlation function, in samples covering different ranges of B-band absolute magnitudes and stellar masses. We considered both threshold and binned galaxy samples, with median B-band absolute magnitudes -21.6<MB-5log(h)<-19.5 and median stellar masses 9.8<log(M*[Msun/h^2])<10.7. We assessed the real-space clustering in the data from the projected correlation function, which we model as a power law in the range 0.2<r_p[Mpc/h]<20. Finally, we estimated the galaxy bias as a function of luminosity, stellar mass, and redshift, assuming a flat LCDM model to derive the dark matter 2PCF. We provide the best-fit parameters of the power-law model assumed for the real-space 2PCF -- the correlation length and the slope -- as well as the linear bias parameter, as a function of the B-band absolute magnitude, stellar mass, and redshift. We confirm and provide the tightest constraints on the dependence of clustering on luminosity at 0.5<z<1.1. We prove the complexity of comparing the clustering dependence on stellar mass from samples that are originally flux-limited and discuss the possible origin of the observed discrepancies. Overall, our measurements provide stronger constraints on galaxy formation models, which are now required to match, in addition to local observations, the clustering evolution measured by VIPERS galaxies between z=0.5 and z=1.1 for a broad range of luminosities and stellar masses.
[85]  [pdf] - 1166361
Detection of galaxy assembly bias
Comments: 6 pages, 4 figures, accepted for publication in MNRAS, Fig.4 updated
Submitted: 2013-05-02, last modified: 2013-07-09
Assembly bias describes the finding that the clustering of dark matter haloes depends on halo formation time at fixed halo mass. In this paper, we analyse the influence of assembly bias on galaxy clustering using both semi-analytical models (SAMs) and observational data. At fixed stellar mass, SAMs predict that the clustering of {\it central} galaxies depends on the specific star formation rate (sSFR), with more passive galaxies having a higher clustering amplitude. We find similar trends using SDSS group catalogues, and verify that these are not affected by possible biases due to the group finding algorithm. Low mass central galaxies reside in narrow bins of halo mass, so the observed trends of higher clustering amplitude for galaxies with lower sSFR is not driven by variations of the parent halo mass. We argue that the clustering dependence on sSFR represent a direct detection of assembly bias. In addition, contrary to what expected based on clustering of dark matter haloes, we find that low-mass central galaxies in SAMs with larger host halo mass have a {\it lower} clustering amplitude than their counter-parts residing in lower mass haloes. This results from the fact that, at fixed stellar mass, assembly bias has a stronger influence on clustering than the dependence on the parent halo mass.
[86]  [pdf] - 1165168
The VIMOS Public Extragalactic Redshift Survey (VIPERS). Galaxy clustering and redshift-space distortions at z=0.8 in the first data release
Comments: 19 pages, 19 figures, accepted for publication in A&A
Submitted: 2013-03-11, last modified: 2013-07-09
We present in this paper the general real- and redshift-space clustering properties of galaxies as measured in the first data release of the VIPERS survey. VIPERS is a large redshift survey designed to probe the distant Universe and its large-scale structure at 0.5 < z < 1.2. We describe in this analysis the global properties of the sample and discuss the survey completeness and associated corrections. This sample allows us to measure the galaxy clustering with an unprecedented accuracy at these redshifts. From the redshift-space distortions observed in the galaxy clustering pattern we provide a first measurement of the growth rate of structure at z = 0.8: f\sigma_8 = 0.47 +/- 0.08. This is completely consistent with the predictions of standard cosmological models based on Einstein gravity, although this measurement alone does not discriminate between different gravity models.
[87]  [pdf] - 1157886
On the scatter in the relation between stellar mass and halo mass: random or halo formation time dependent?
Comments: 10 pages, 6 figures, published in MNRAS
Submitted: 2012-11-19, last modified: 2013-07-09
The empirical HOD model of Wang et al. 2006 fits, by construction, both the stellar mass function and correlation function of galaxies in the local Universe. In contrast, the semi-analytical models of De Lucia & Blazoit 2007 (DLB07) and Guo et al. 2011 (Guo11), built on the same dark matter halo merger trees than the empirical model, still have difficulties in reproducing these observational data simultaneously. We compare the relations between the stellar mass of galaxies and their host halo mass in the three models, and find that they are different. When the relations are rescaled to have the same median values and the same scatter as in Wang et al., the rescaled DLB07 model can fit both the measured galaxy stellar mass function and the correlation function measured in different galaxy stellar mass bins. In contrast, the rescaled Guo11 model still over-predicts the clustering of low-mass galaxies. This indicates that the detail of how galaxies populate the scatter in the stellar mass -- halo mass relation does play an important role in determining the correlation functions of galaxies. While the stellar mass of galaxies in the Wang et al. model depends only on halo mass and is randomly distributed within the scatter, galaxy stellar mass depends also on the halo formation time in semi-analytical models. At fixed value of infall mass, galaxies that lie above the median stellar mass -- halo mass relation reside in haloes that formed earlier, while galaxies that lie below the median relation reside in haloes that formed later. This effect is much stronger in Guo11 than in DLB07, which explains the over-clustering of low mass galaxies in Guo11. Our results illustrate that the assumption of random scatter in the relation between stellar and halo mass as employed by current HOD and abundance matching models may be problematic in case a significant assembly bias exists in the real Universe.
[88]  [pdf] - 1171689
The Hierarchical Origins of Observed Galaxy Morphology
Comments: latex, 20 pages, 13 figures. Accepted by Monthly Notices. Source package includes full version of Table 1 from paper (file
Submitted: 2013-05-30
Galaxies grow primarily via accretion-driven star formation in discs and merger-driven growth of bulges. These processes are implicit in semi-analytical models of galaxy formation, with bulge growth in particular relating directly to the hierarchical build-up of halos and their galaxies. In this paper, we consider several implementations of two semi-analytical models. Focusing on implementations in which bulges are formed during mergers only, we examine the fractions of elliptical galaxies and both passive and star-forming disk galaxies as functions of stellar and halo mass, for central and satellite systems. This is compared to an observational cross-matched SDSS+RC3 z ~ 0 sample of galaxies with accurate visual morphological classifications and M_{stellar} > 10^10.5 M_{sol}. The models qualitatively reproduce the observed increase of elliptical fraction with stellar mass, and with halo mass for central galaxies, supporting the idea that observed ellipticals form during major mergers. However, the overall elliptical fraction produced by the models is much too high compared with the z ~ 0 data. Since the "passive" -- i.e. non-star-forming -- fractions are approximately reproduced, and since the fraction which are star-forming disc galaxies is also reproduced, the problem is that the models overproduce ellipticals at the expense of passive S0 and spiral galaxies. Bulge-growth implementations (tuned to reproduce simulations) which allow the survival of residual discs in major mergers still destroy too much of the disc. Increasing the lifetime of satellites, or allowing significant disc regrowth around merger remnants, merely increases the fraction of star-forming disc galaxies. Instead, it seems necessary to reduce the mass ratios of merging galaxies, so that most mergers produce modest bulge growth in disc-galaxy remnants instead of ellipticals. [Abridged]
[89]  [pdf] - 1164643
Isolated galaxies in hierarchical galaxy formation models - present-day properties and environmental histories
Comments: 15 pages, 13 figures, minor changes in the text, accepted for publication by MNRAS
Submitted: 2013-02-14, last modified: 2013-05-09
In this study, we have carried out a detailed, statistical analysis of isolated model galaxies, taking advantage of publicly available hierarchical galaxy formation models. To select isolated galaxies, we employ 2D methods widely used in the observational literature, as well as a more stringent 3D isolation criterion that uses the full 3D-real space information. In qualitative agreement with observational results, isolated model galaxies have larger fractions of late-type, star forming galaxies with respect to randomly selected samples of galaxies with the same mass distribution. We also find that the samples of isolated model galaxies typically contain a fraction of less than 15 per cent of satellite galaxies, that reside at the outskirts of their parent haloes where the galaxy number density is low. Projection effects cause a contamination of 2D samples of about 18 per cent, while we estimate a typical completeness of 65 per cent. Our model isolated samples also include a very small (few per cent) fraction of bulge dominated galaxies (B/T > 0.8) whose bulges have been built mainly by minor mergers. Our study demonstrates that about 65-70 per cent of 2D isolated galaxies that are classified as isolated at z = 0 have indeed been completely isolated since z = 1 and only 7 per cent have had more than 3 neighbours within a comoving radius of 1 Mpc. Irrespectively of the isolation criteria, roughly 45 per cent of isolated galaxies have experienced at least one merger event in the past (most of the mergers are minor, with mass ratios between 1:4 and 1:10). The latter point validates the approximation that isolated galaxies have been mainly influenced by internal processes.
[90]  [pdf] - 1165128
An improved prescription for merger time-scales from controlled simulations
Comments: 6 pages, 3 figures, 2 tables. Replaced to match the version accepted for publication in MNRAS
Submitted: 2013-03-08, last modified: 2013-05-04
We compare three analytical prescriptions for merger times available from the literature to simulations of isolated mergers. We probe three different redshifts, and several halo concentrations, mass ratios, orbital circularities and orbital energies of the satellite. We find that prescriptions available in the literature significantly under-predict long timescales for mergers at high redshift. We argue that these results have not been highlighted previously either because the evolution of halo concentration of satellite galaxies has been neglected (in previous isolated merger simulations), or because long merger times and mergers with high initial orbital circularities are under-represented (for prescriptions based on cosmological simulations). Motivated by the evolution of halo concentration at fixed mass, an explicit dependence on redshift added as t_merger,modified(z) = (1+z)^0.44 t_merger to the prescription based on isolated mergers gives a significant improvement in the predicted merger times up to ~20 t_dyn in the redshift range 0<z<2. When this modified prescription is used to compute galaxy stellar mass functions, we find that it leads up to a 25 per cent increase in the number of low mass galaxies surviving at z=0, and a 10 per cent increase for more massive galaxies. This worsen the known over-prediction in the number of low mass galaxies by hierarchical models of galaxy formation.
[91]  [pdf] - 640533
VIPERS: An Unprecedented View of Galaxies and Large-Scale Structure Halfway Back in the Life of the Universe
Comments: 6 pages, 8 figures; introductory paper
Submitted: 2013-03-15
The VIMOS Public Extragalactic Redshift Survey (VIPERS) is an ongoing ESO Large Programme to map in detail the large-scale distribution of galaxies at 0.5 < z <1.2. With a combination of volume and sampling density that is unique for these redshifts, it focuses on measuring galaxy clustering and related cosmological quantities as part of the grand challenge of understanding the origin of cosmic acceleration. VIPERS has also been designed to guarantee a broader legacy, allowing detailed investigations of the properties and evolutionary trends of z~1 galaxies. The survey strategy exploits the specific advantages of the VIMOS spectrograph at the VLT, aiming at a final sample of nearly 100,000 galaxy redshifts to iAB = 22.5 mag, which represents the largest redshift survey ever performed with ESO telescopes. In this introductory article we describe the survey construction, together with early results based on a first sample of ~55,000 galaxies.
[92]  [pdf] - 844753
The VIMOS Public Extragalactic Redshift Survey (VIPERS). An unprecedented view of galaxies and large-scale structure at 0.5<z<1.2
Comments: 20 pages, 19 figures, submitted to A&A. See for version with higher resolution figures and movie
Submitted: 2013-03-11, last modified: 2013-03-13
We describe the construction and general features of VIPERS, the VIMOS Public Extragalactic Redshift Survey. This `Large Programme' has been using the ESO VLT with the aim of building a spectroscopic sample of ~100,000 galaxies with i_{AB}<22.5 and 0.5<z<1.5. The survey covers a total area of ~24 deg^2 within the CFHTLS-Wide W1 and W4 fields. VIPERS is designed to address a broad range of problems in large-scale structure and galaxy evolution, thanks to a unique combination of volume (~ 5 x 10^7 h^{-3} Mpc^3) and sampling rate (~ 40%), comparable to state-of-the-art surveys of the local Universe, together with extensive multi-band optical and near-infrared photometry. Here we present the survey design, the selection of the source catalogue and the development of the spectroscopic observations. We discuss in detail the overall selection function that results from the combination of the different constituents of the project. This includes the masks arising from the parent photometric sample and the spectroscopic instrumental footprint, together with the weights needed to account for the sampling and the success rates of the observations. Using the catalogue of 53,608 galaxy redshifts composing the forthcoming VIPERS Public Data Release 1 (PDR-1), we provide a first assessment of the quality of the spectroscopic data. Benefiting from the combination of size and detailed sampling of this dataset, we conclude by presenting a map showing in unprecedented detail the large-scale distribution of galaxies between 5 and 8 billion years ago. [abridged]
[93]  [pdf] - 1091387
The galaxy stellar mass function and its evolution with time show no dependence on global environment
Comments: final version as published by A&A
Submitted: 2011-11-03, last modified: 2013-01-28
We present the analysis of the galaxy stellar mass function in different environments at intermediate redshift (0.3<z<0.8) for two mass-limited galaxy samples. We use the IMACS Cluster Building Survey (ICBS), at masses M_ast >10^(10.5) M_sun, to study cluster, group, and field galaxies at z=0.3-0.45, and the ESO Distant Cluster Survey (EDisCS), at masses M_ast > 10^(10.2) M_sun, to investigate cluster and group galaxies at z=0.4-0.8. Therefore, in our analysis we include galaxies that are slightly less massive than the Milky Way. Having excluded the brightest cluster galaxies, we show thatthe shape of the mass distribution does not seem to depend on global environment. Our two main results are: (1) Galaxies in the virialized regions of clusters, in groups, and in the field follow a similar mass distribution. (2) Comparing both ICBS and EDisCS mass functions to mass functions in the local Universe, we find evolution from z~0.4-0.6 to z~0.07. The population of low-mass galaxies has proportionally grown with time with respect to that of massive galaxies. This evolution is independent of environment -- the same for clusters and the field. Furthermore, considering only clusters, we find that no differences can be detected neither within the virialized regions, nor when we compare galaxies within and outside the virial radius. Subdividing galaxies in terms of color, in clusters, groups, and field red and blue galaxies are regulated by different mass functions, but comparing separately the blue and red mass functions in different environments, no differences are detected in their shape.
[94]  [pdf] - 1123819
The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology
Comments: 20 pages, 16 figures. MNRAS in press. Minor changes in this version to match the version in press
Submitted: 2012-05-31, last modified: 2013-01-19
We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.
[95]  [pdf] - 616347
A Research Note on the Implementation of Star Formation and Stellar Feedback in Semi-Analytic Models
Comments: 19 pages, 7 figures, Research Note, Comments welcome
Submitted: 2013-01-17
We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.
[96]  [pdf] - 1034778
The effect of thermally pulsating asymptotic giant branch stars on the evolution of the rest-frame near-infrared galaxy luminosity function
Comments: 10 pages, 5 figures, accepted by MNRAS
Submitted: 2010-09-07, last modified: 2012-12-13
We address the fundamental question of matching the rest-frame K-band luminosity function (LF) of galaxies over the Hubble time using semi-analytic models, after modification of the stellar population modelling. We include the Maraston evolutionary synthesis models, that feature a higher contribution by the Thermally Pulsating - Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and the Menci model. We leave all other input physics and parameters unchanged. We find that the modification of the stellar population emission can solve the mismatch between models and the observed rest-frame K-band luminosity from the brightest galaxies derived from UKIDSS data at high redshift. For all explored semi-analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy was recently pointed out. The reason for the success is that at these cosmic epochs the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB phase which makes them redder without the need of changing their mass or age. At the same time, the known overestimation of the faint end is enhanced in the K-band when including the TP-AGB contribution. At lower redshifts (z<2) some of the explored models deviate from the data. This is due to too short merging timescales and inefficient 'radio-mode' AGN feedback. Our results show that a strong evolution in mass predicted by hierarchical models is compatible with no evolution on the bright-end of the K-band LF from z=3 to the local universe. This means that, at high redshifts and contrary to what is commonly accepted, K-band emission is not necessarily a good tracer of galaxy mass.
[97]  [pdf] - 1158499
The evolution of galaxy sizes
Comments: to appear in the proceedings of the IAU S295: The intriguing life of massive galaxies, editors D. Thomas, A. Pasquali & I. Ferreras
Submitted: 2012-12-13
We present a study of galaxy sizes in the local Universe as a function of galaxy environment, comparing clusters and the general field. Galaxies with radii and masses comparable to high-z massive and compact galaxies represent 4.4% of all galaxies more massive than 3 X 10^{10} M_sun in the field. Such galaxies are 3 times more frequent in clusters than in the field. Most of them are early-type galaxies with intermediate to old stellar populations. There is a trend of smaller radii for older luminosity-weighted ages at fixed galaxy mass. We show the relation between size and luminosity-weighted age for galaxies of different stellar masses and in different environments. We compare with high-z data to quantify the evolution of galaxy sizes. We find that, once the progenitor bias due to the relation between galaxy size and stellar age is removed, the average amount of size evolution of individual galaxies between high- and low-z is mild, of the order of a factor 1.6.
[98]  [pdf] - 1157518
Superdense galaxies and the mass-size relation at low redshift
Comments: ApJ, in press
Submitted: 2012-11-05
We search for massive and compact galaxies (superdense galaxies, hereafter SDGs) at z=0.03-0.11 in the Padova-Millennium Galaxy and Group Catalogue, a spectroscopically complete sample representative of the local Universe general field population. We find that compact galaxies with radii and mass densities comparable to high-z massive and passive galaxies represent 4.4% of all galaxies with stellar masses above 3 X 10^10 M_sun, yielding a number density of 4.3 X 10^-4 h^3 Mpc^-3. Most of them are S0s (70%) or ellipticals (23%), are red and have intermediate-to-old stellar populations, with a median luminosity-weighted age of 5.4 Gyr and a median mass-weighted age of 9.2 Gyr. Their velocity dispersions and dynamical masses are consistent with the small radii and high stellar mass estimates. Comparing with the WINGS sample of cluster galaxies at similar redshifts, the fraction of superdense galaxies is three times smaller in the field than in clusters, and cluster SDGs are on average 4 Gyr older than field SDGs. We confirm the existence of a universal trend of smaller radii for older luminosity-weighted ages at fixed galaxy mass. On top of the well known dependence of stellar age on galaxy mass, the luminosity-weighted age of galaxies depends on galaxy compactness at fixed mass, and, for a fixed mass and radius, on environment. This effect needs to be taken into account in order not to overestimate the evolution of galaxy sizes from high- to low-z. Our results and hierarchical simulations suggest that a significant fraction of the massive compact galaxies at high-z have evolved into compact galaxies in galaxy clusters today. When stellar age and environmental effects are taken into account, the average amount of size evolution of individual galaxies between high- and low-z is mild, a factor ~1.6. (abridged)
[99]  [pdf] - 1124926
The VIMOS Public Extragalactic Redshift Survey (VIPERS): spectral classification through Principal Component Analysis
Comments: 15 pages, 20 images, accepted for publication in MNRAS: MN-12-1739-MJ.R1
Submitted: 2012-07-18, last modified: 2012-10-04
We develop a Principal Component Analysis aimed at classifying a sub-set of 27,350 spectra of galaxies in the range 0.4 < z < 1.0 collected by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We apply an iterative algorithm to simultaneously repair parts of spectra affected by noise and/or sky residuals, and reconstruct gaps due to rest-frame transformation, and obtain a set of orthogonal spectral templates that span the diversity of galaxy types. By taking the three most significant components, we find that we can describe the whole sample without contamination from noise. We produce a catalogue of eigen-coefficients and template spectra that will be part of future VIPERS data releases. Our templates effectively condense the spectral information into two coefficients that can be related to the age and star formation rate of the galaxies. We examine the spectrophotometric types in this space and identify early, intermediate, late and starburst galaxies.
[100]  [pdf] - 1123211
Comparison of the VIMOS-VLT Deep Survey with the Munich semi-analytical model. II. The colour-density relation up to z=1.5
Comments: 15 pages, 12 figures, accepted for publication in A&A
Submitted: 2012-05-07, last modified: 2012-09-10
[Abridged] We perform on galaxy mock catalogues the same colour-density analysis made by Cucciati et al. (2006) on a 5 Mpc/h scale using the VVDS-Deep survey, and compare the results from mocks with observed data. We use mocks with the same flux limits (I=24) as the VVDS (CMOCKS), built using the semi- analytic model by De Lucia & Blaizot (2007) applied to the Millennium Simulation. From CMOCKS, we extracted samples of galaxies mimicking the VVDS observational strategy (OMOCKS). We computed the B-band Luminosity Function LF and the colour-density relation (CDR) in the mocks. We find that the LF in mocks roughly agrees with the observed LF, but at z<0.8 the faint-end slope of the model LF is steeper than the VVDS one. Computing the LF for early and late type galaxies, we show that mocks have an excess of faint early-type and of bright late-type galaxies with respect to data. We find that the CDR in OMOCKS is in excellent agreement with the one in CMOCKS. At z~0.7, the CDR in mocks agrees with the VVDS one (red galaxies reside mainly in high densities). Yet, the strength of the CDR in mocks does not vary within 0.2<z<1.5, while the observed relation flattens with increasing z and possibly inverts at z=1.3. We argue that the lack of evolution in the CDR in mocks is not due only to inaccurate prescriptions for satellite galaxies, but that also the treatment of central galaxies has to be revised. The reversal of the CDR can be explained by wet mergers between young galaxies, producing a starburst event. This should be seen on group scales. A residual of this is found in observations at z=1.5 on larger scales, but not in the mocks, suggesting that the treatment of physical processes affecting satellites and central galaxies in models should be revised.
[101]  [pdf] - 1077774
The effect of the environment on the gas kinematics and the structure of distant galaxies
Comments: 26 pages (plus small appendix), 21 Figures. Accepted for publication in MNRAS. v2: typos corrected; references added; and minor changes in tables and plots to include a few missing data points; results and conclusions unchanged; v3: corrected typos in the labels of 3 plots; v4: typo in the label of Figure 18 corrected
Submitted: 2011-07-05, last modified: 2012-08-16
With the aim of distinguishing between possible physical mechanisms acting on galaxies when they fall into clusters, we study the properties of the gas and the stars in a sample of 422 emission-line galaxies from EDisCS in different environments up to z~1. We identify galaxies with kinematical disturbances (from emission-lines in their 2D spectra) and find that they are more frequent in clusters than in the field. The fraction of kinematically-disturbed galaxies increases with cluster velocity dispersion and decreases with distance from the cluster centre, but remains constant with projected galaxy density. We also studied morphological disturbances in the stellar light from HST/F814W images, finding that the fraction of morphologically disturbed galaxies is similar in all environments. Moreover, there is little correlation between the presence of kinematically-disturbed gas and morphological distortions. We also study the dependence of the Tully-Fisher relation, star formation, and extent of the emission on environment, and conclude that the gas disks in cluster galaxies have been truncated, and therefore their star formation is more concentrated than in low-density environments. If spirals transform into S0s, our findings imply that the physical mechanism transforming cluster galaxies efficiently disturbs the star-forming gas and reduces their specific star formation rate. Moreover, this star-forming gas is either removed more efficiently from the outskirts of the galaxies or is driven towards the centre (or both), helping to build the bulges of S0s. These results, in addition to the finding that the transformation mechanism does not seem to induce strong morphological disturbances on the galaxies, suggest that the physical processes involved are related to the intracluster medium, with galaxy-galaxy interactions playing only a limited role in clusters.
[102]  [pdf] - 1125191
On the stellar populations of massive galaxies
Comments: 5 pages, 3 figures; accepted for publication in MNRAS Letters
Submitted: 2012-07-30
In this Letter, we analyse the predicted physical properties of massive galaxies, in the framework of recent semi-analytic models of galaxy formation. All models considered account for winds driven by supernovae explosions and suppression of gas condensation at the centre of relatively massive haloes by active galactic nuclei (AGN). We show that, while these models successfully reproduce the old stellar populations observed for massive galaxies, they fail in reproducing their observed chemical abundances. This problem is alleviate but still present if AGN feedback is completely switched off. Moreover, in this case, model predictions fail in accounting for the old stellar ages of massive galaxies. We argue that the difficulty of semi-analytical models in simultaneously reproducing the observed ages and metallicities of massive galaxies, signals a fundamental problem with the schemes that are currently adopted to model star formation, feedback, and related recycling of gas and metals.
[103]  [pdf] - 1124039
Dark satellites and the morphology of dwarf galaxies
Comments: Submitted to ApJ Letters. 5 pages, 4 figures. A movie showing an encounter between a disky dwarf galaxy and a dark satellite can be found at
Submitted: 2012-06-11
One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.
[104]  [pdf] - 1117503
Milky Way simulations: the Galaxy, its stellar halo and its satellites - insights from a hybrid cosmological approach
Comments: 10 pages, 10 figures, invited plenary review at the Astronomische Gesellschaft meeting (Heidelberg 2011), to appear in Reviews in Modern Astronomy Vol.24 - typos corrected
Submitted: 2012-03-23, last modified: 2012-04-26
Our 'home galaxy' - the Milky Way - is a fairly large spiral galaxy, prototype of the most common morphological class in the local Universe. Although being only a galaxy, it is the only one that can be studied in unique detail: for the MilkyWay and for a number of members of the Local Group, a wealth of observational data is available about the ages and chemical abundances of their stars. Much more information is expected to come in the next few years, from ongoing and planned spectroscopic and astrometric surveys, providing a unique benchmark for modern theories of galaxy formation. In this review, I will summarize recent results on the formation of our Milky Way, its stellar halo, and its satellite galaxies. I will focus, in particular, on results obtained in the framework of hybrid models of galaxy formation, and refer to other reviews in this issue for studies based on hydrodynamical simulations.
[105]  [pdf] - 1091951
The environmental history of group and cluster galaxies in a LambdaCDM Universe
Comments: 18 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2011-11-28, last modified: 2012-03-23
We use publicly available galaxy merger trees, obtained applying semi-analytic techniques to a large high resolution cosmological simulation, to study the environmental history of group and cluster galaxies. Our results highlight the existence of an intrinsic history bias which makes the nature versus nurture (as well as the mass versus environment) debate inherently ill posed. In particular we show that: (i) surviving massive satellites were accreted later than their less massive counterparts, from more massive haloes; (ii) the mixing of galaxy populations is incomplete during halo assembly, which creates a correlation between the time a galaxy becomes satellite and its present distance from the parent halo centre. The weakest trends are found for the most massive satellites, as a result of efficient dynamical friction and late formation times of massive haloes. A large fraction of the most massive group/cluster members are accreted onto the main progenitor of the final halo as central galaxies, while about half of the galaxies with low and intermediate stellar mass are accreted as satellites. Large fractions of group and cluster galaxies (in particular those of low stellar mass) have therefore been `pre-processed' as satellites of groups with mass ~10^13 Msun. To quantify the relevance of hierarchical structure growth on the observed environmental trends, we have considered observational estimates of the passive galaxy fractions, and their variation as a function of halo mass and cluster-centric distance. Comparisons with our theoretical predictions require relatively long times (~5-7 Gyr) for the suppression of star formation in group and cluster satellites. [abridged]
[106]  [pdf] - 1116851
A warm mode of gas accretion on forming galaxies
Comments: 6 pages, 3 figures, accepted for publication in ApJL
Submitted: 2012-02-23
We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}<2.5 \times 10^5$ K, a hot one for $T>10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies.
[107]  [pdf] - 1116365
Simulating the evolution of disc galaxies in a group environment. I. The influence of the global tidal field
Comments: 31 pages, 20 figures, and 3 tables. Accepted for publication in MNRAS
Submitted: 2012-02-02
We present the results of a series of numerical simulations aimed to study the evolution of a disc galaxy within the global tidal field of a group environment. Both the disc galaxy and the group are modelled as multi-component, collision-less, N-body systems, composed by both dark matter and stars. In our simulations, the evolution of disc galaxies is followed as their orbits sink towards the group centre, under the effect of dynamical friction. We explore a broad parameter space, covering several aspects of the galaxy-group interaction that are potentially relevant to galaxy evolution. Namely, prograde and retrograde orbits, orbital eccentricities, disc inclination, role of a central bulge in discs, internal disc kinematics, and galaxy-to-group mass ratios. We find that significant disc transformations occur only after the mean density of the group, measured within the orbit of the galaxy, exceeds ~0.3-1 times the central mean density of the galaxy. The morphological evolution of discs is found to be strongly dependent on the initial inclination of the disc with respect to its orbital plane. That is, discs on face-on and retrograde orbits are shown to retain longer their disc structures and kinematics, in comparison to prograde discs. This suggests that after interacting with the global tidal field alone, a significant fraction of disc galaxies should be found in the central regions of groups. Prominent central bulges are not produced, and pre-existing bulges are not enhanced in discs after the interaction with the group. Assuming that most S0 are formed in group environments, this implies that prominent bulges should be formed mostly by young stars, created only after a galaxy has been accreted by a group. Finally, contrary to some current implementations of tidal stripping in semi-analytical models of galaxy evolution, we find that more massive galaxies suffer more tidal stripping.
[108]  [pdf] - 1091485
Statistics of Substructures in Dark Matter Haloes
Comments: 13 pages, 11 figures and 1 table. MNRAS 2011 in press
Submitted: 2011-11-08
We study the amount and distribution of dark matter substructures within dark matter haloes, using a large set of high-resolution simulations ranging from group size to cluster size haloes, and carried our within a cosmological model consistent with WMAP 7-year data. In particular, we study how the measured properties of subhaloes vary as a function of the parent halo mass, the physical properties of the parent halo, and redshift. The fraction of halo mass in substructures increases with increasing mass. There is, however, a very large halo-to-halo scatter that can be explained only in part by a range of halo physical properties, e.g. concentration. At given halo mass, less concentrated haloes contain significantly larger fractions of mass in substructures because of the reduced strength of tidal disruption. Most of the substructure mass is located at the outskirts of the parent haloes, in relatively few massive subhaloes. This mass segregation appears to become stronger at increasing redshift, and should reflect into a more significant mass segregation of the galaxy population at different cosmic epochs. When haloes are accreted onto larger structures, their mass is significantly reduced by tidal stripping. Haloes that are more massive at the time of accretion (these should host more luminous galaxies) are brought closer to the centre on shorter time-scales by dynamical friction, and therefore suffer of a more significant stripping. The halo merger rate depends strongly on the environment with substructure in more massive haloes suffering more important mergers than their counterparts residing in less massive systems. This should translate into a different morphological mix for haloes of different mass.
[109]  [pdf] - 1076742
Autocorrelations of stellar light and mass at z~0 and ~1: From SDSS to DEEP2
Comments: 10 pages, 4 figures, accepted for publication in Monthly Notices, minor revisions in text
Submitted: 2011-05-19, last modified: 2011-09-13
We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.
[110]  [pdf] - 1084015
What determines the fraction of elliptical galaxies in clusters?
Comments: 8 pages, 3 figures, accepted for publication in MNRAS
Submitted: 2011-09-12
We study the correlation between the morphological mix of cluster galaxies and the assembly history of the parent cluster by taking advantage of two independently developed semi-analytic models for galaxy formation and evolution. In our models, both the number of cluster members and that of elliptical members increase as a function of cluster mass, in such a way that the resulting elliptical fractions are approximately independent of cluster mass. The population of cluster ellipticals exhibit a marked bimodal distribution as a function of galaxy stellar mass, with a dip at masses $\sim 10^{10}\,{\rm M}_{\odot}$. In the framework of our models, this bimodality originates from the combination of a strongly decreasing number of galaxies with increasing stellar mass, and a correspondingly increasing probability of experiencing major mergers. We show that the correlation between the measured elliptical fraction and the assembly history of the parent cluster is weak, and that it becomes stronger in models that adopt longer galaxy merger times. We argue that this results from the combined effect of a decreasing bulge production due to a reduced number of mergers, and an increasing survival probability of pre-existing ellipticals, with the latter process being more important than the former.
[111]  [pdf] - 1083198
On the Radial Stellar Content of Early-Type Galaxies as a Function of Mass and Environment
Comments: 2 figures, accepted for publication on ApJL
Submitted: 2011-08-03
Using optical-optical and optical-NIR colors, we analyze the radial dependence of age and metallicity inside massive (M* > 10^10.5 MSun), low-redshift (z<0.1), early-type galaxies (ETGs), residing in both high-density group regions and the field. On average, internal color gradients of ETGs are mainly driven by metallicity, consistent with previous studies. However, we find that group galaxies feature positive age gradients, Nabla_t, i.e. a younger stellar population in the galaxy center, and steeper metallicity gradients, compared to the field sample, whose Nabla_t ranges from negative in lower mass galaxies, to positive gradients at higher mass. These dependencies yield new constraints to models of galaxy formation and evolution. We speculate that age and metallicity gradients of group ETGs result from (either gas-rich or minor-dry) mergers and/or cold-gas accretion, while field ETGs exhibit the characteristic flatter gradients expected from younger, more metal-rich, stars formed inside--out by later gas-cooling.
[112]  [pdf] - 1078117
Convergence of Galaxy Properties with Merger Tree Temporal Resolution
Comments: 15 pages, 12 figures, submitted to MNRAS. Parameter files and plotting scripts available at and in the anciliary data stored with this arXiv submission
Submitted: 2011-07-20
Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.
[113]  [pdf] - 1052405
The population of Milky Way satellites in the LambdaCDM cosmology
Comments: The paper contains 23 pages, 13 figures. Accepted to MNRAS (2011, June 28)
Submitted: 2011-02-28, last modified: 2011-07-18
We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.
[114]  [pdf] - 1052235
Properties of fossil groups in cosmological simulations and galaxy formation models
Comments: 14 pages, 10 figures; accepted by MNRAS, minor changes to match the accepted version
Submitted: 2011-02-21, last modified: 2011-06-15
It has been a long-standing question whether fossil groups are just sampling the tail of the distribution of ordinary groups, or whether they are a physically distinct class of objects, characterized by an unusual and special formation history. To study this question, we here investigate fossil groups identified in the hydrodynamical simulations of the GIMIC project, which consists of resimulations of five regions in the Millennium Simulation (MS) that are characterized by different large-scale densities, ranging from a deep void to a proto-cluster region. For comparison, we also consider semi-analytic models built on top of the MS, as well as a conditional luminosity function approach. We identify galaxies in the GIMIC simulations as groups of stars and use a spectral synthesis code to derive their optical properties. The X-ray luminosity of the groups is estimated in terms of the thermal bremsstrahlung emission of the gas in the host halos, neglecting metallicity effects. We focus on comparing the properties of fossil groups in the theoretical models and observational results, highlighting the differences between them, and trying to identify possible dependencies on environment for which our approach is particularly well set-up. We find that the optical fossil fraction in all of our theoretical models declines with increasing halo mass, and there is no clear environmental dependence. Combining the optical and X-ray selection criteria for fossil groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over the GIMIC halo mass range we resolve best, 9.0\times1012 \sim 4.0\times1013 h-1 M, the central galaxies in the fossil groups show similar properties as those in ordinary groups, in terms of age, metallicity, color, concentration, and mass-to-light ratio. [abridged]
[115]  [pdf] - 1052129
The other side of Bulge Formation in a Lambda-CDM cosmology: Bulgeless Galaxies in the Local Universe
Comments: Accepted for publication in MNRAS; updated version with additional comparisons with observations
Submitted: 2011-02-15, last modified: 2011-05-11
We study the physical properties, formation histories, and environment of galaxies without a significant "classical" spheroidal component, as predicted by semi-analytical models of galaxy formation and evolution. This work is complementary to the analysis presented in De Lucia et al., (2011), where we focus on the relative contribution of various physical mechanisms responsible for bulge assembly in a Lambda-CDM cosmology. We find that the fraction of bulgeless galaxies is a strong decreasing function of stellar mass: they represent a negligible fraction of the galaxy population with M* > 10^12 Msun, but dominate at M* < 10^10 Msun. We find a clear dichotomy in this galaxy population, between central galaxies of low-mass dark matter haloes, and satellite galaxies in massive groups/clusters. We show that bulgeless galaxies are relatively young systems, that assemble most of their mass at low-redshift, but they can also host very old stellar populations. Since galaxy-galaxy mergers are assumed to lead to the formation of a spheroidal component, in our models these galaxies form preferentially in low-mass haloes that host a small number of satellites galaxies. We show that the adopted modelling for galaxy mergers represents a key ingredient in determining the actual number of bulgeless galaxies. Our results show that these galaxies are not a rare population in theoretical models: at z~0, galaxies with no classical bulge (but often including galaxies with the equivalent of pseudo-bulges) account for up to 14% of the galaxies with 10^11 < M*/Msun < 10^12.
[116]  [pdf] - 1052127
Times, environments and channels of bulge formation in a LambdaCDM cosmology
Comments: 18 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2011-02-15
We analyze predictions from two independently developed galaxy formation models to study the mechanisms, environments, and characteristic times of bulge formation in a LambdaCDM cosmogony. For each model, we test different prescriptions for bulge formation in order to quantify the relative importance of different channels. Our results show that the strong correlation between galaxy and halo mass for central galaxies, and the richer merger history of more massive systems naturally give rise to a strong correlation between galaxy mass and morphology, and between halo mass and morphological type of central galaxies. Large fractions of the bulge mass are acquired through major and minor mergers, but disc instability plays an important role, particularly for intermediate mass galaxies. We find that the modelling of disc instability events, as well as of the galaxy merger times, can affect significantly the timing of bulge formation, and the relative importance of different channels. Bulge dominated galaxies acquire their morphology through major mergers, but this can be modified by cooling of gas from the surrounding hot halo. We find that disc regrowth is a non negligible component of the evolution of bulge dominated galaxies, particularly for low to intermediate masses, and at high redshifts.
[117]  [pdf] - 955389
From dwarf spheroidals to cDs: Simulating the galaxy population in a LCDM cosmology
Comments: Accepted for publication in MNRAS. SQL databases containing the full galaxy data at all redshifts and for both the Millennium and Millennium-II simulations are publicly released at
Submitted: 2010-06-01, last modified: 2011-01-06
We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.
[118]  [pdf] - 1033401
The dependence of AGN activity on stellar and halo mass in Semi-Analytic Models
Comments: 15 pages; 8 figures; 1 table; updated to match MNRAS accepted version
Submitted: 2010-06-29, last modified: 2010-12-17
AGN feedback is believed to play an important role in shaping a variety of observed galaxy properties, as well as the evolution of their stellar masses and star formation rates. In particular, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play a crucial role in regulating the levels of activity in galaxies, in relatively massive halos at low redshift. Only in recent years, however, has detailed statistical information on the dependence of galaxy activity on stellar mass, parent halo mass and hierarchy has become available. In this paper, we compare the fractions of galaxies belonging to different activity classes (star-forming, AGN and radio active) with predictions from four different and independently developed semi-analytical models. We adopt empirical relations to convert physical properties into observables (H_alpha emission lines, OIII line strength and radio power). We demonstrate that all models used in this study reproduce the overall distributions of galaxies belonging to different activity classes as a function of stellar mass and halo mass: star forming galaxies and the strongest radio sources are preferentially associated with low-mass and high-mass galaxies/halos respectively. However, model predictions differ from observational measurements in a number of ways. All models used in our study predict that almost every >1.e12 Msun dark matter halo and/or >1.e11 Msun galaxy should host a bright radio source, while only a small fraction of galaxies belong to this class in the data. In addition, radio brightness is expected to depend strongly on the mass of the parent halo mass in the models, while strong and weak radio galaxies are found in similar environments in data. Our results highlight that the distribution of AGN as a function of stellar mass provides one of the most promising discriminants between different gas accretion schemes.
[119]  [pdf] - 1347409
Modelling the evolution of galaxies as a function of environment
Comments: Invited Review given at the Workshop "Environment and the Formation of Galaxies: 30 years later" held in Lisbon, 6-7 September 2010
Submitted: 2010-12-15
In this review, I provide an overview of theoretical aspects related to the evolution of galaxies as a function of environment. I discuss the main physical processes at play, their characteristic time-scales and environmental dependency, and comment on their treatment in the framework of hierarchical galaxy formation models. I briefly summarize recent results and the main open issues.
[120]  [pdf] - 1041018
Comparison of the VIMOS-VLT Deep Survey with the Munich semi-analytical model - I. Magnitude counts, redshift distribution, colour bimodality, and galaxy clustering
Comments: 15 pages, 14 figures, accepted for publication in A&A
Submitted: 2010-10-02, last modified: 2010-10-06
This paper presents a detailed comparison between high-redshift observations from the VIMOS-VLT Deep Survey (VVDS) and predictions from the Munich semi-analytical model of galaxy formation. In particular, we focus this analysis on the magnitude, redshift, and colour distributions of galaxies, as well as their clustering properties. We constructed 100 quasi-independent mock catalogues, using the output of the semi-analytical model presented in De Lucia & Blaizot (2007).We then applied the same observational selection function of the VVDS-Deep survey, so as to carry out a fair comparison between models and observations. We find that the semi-analytical model reproduces well the magnitude counts in the optical bands. It tends, however, to overpredict the abundance of faint red galaxies, in particular in the i' and z' bands. Model galaxies exhibit a colour bimodality that is only in qualitative agreement with the data. In particular, we find that the model tends to overpredict the number of red galaxies at low redshift and of blue galaxies at all redshifts probed by VVDS-Deep observations, although a large fraction of the bluest observed galaxies is absent from the model. In addition, the model overpredicts by about 14 per cent the number of galaxies observed at 0.2<z<1 with I_AB<24. When comparing the galaxy clustering properties, we find that model galaxies are more strongly clustered than observed ones at all redshift from z=0.2 to z=2, with the difference being less significant above z~1. When splitting the samples into red and blue galaxies, we find that the observed clustering of blue galaxies is well reproduced by the model, while red model galaxies are much more clustered than observed ones, being principally responsible for the strong global clustering found in the model. [abridged]
[121]  [pdf] - 200714
A decline and fall in the future of Italian Astronomy?
Antonelli, Angelo; Antonuccio-Delogu, Vincenzo; Baruffolo, Andrea; Benetti, Stefano; Bianchi, Simone; Biviano, Andrea; Bonafede, Annalisa; Bondi, Marco; Borgani, Stefano; Bragaglia, Angela; Brescia, Massimo; Brucato, John Robert; Brunetti, Gianfranco; Brunino, Riccardo; Cantiello, Michele; Casasola, Viviana; Cassano, Rossella; Cellino, Alberto; Cescutti, Gabriele; Cimatti, Andrea; Comastri, Andrea; Corbelli, Edvige; Cresci, Giovanni; Criscuoli, Serena; Cristiani, Stefano; Cupani, Guido; De Grandi, Sabrina; D'Elia, Valerio; Del Santo, Melania; De Lucia, Gabriella; Desidera, Silvano; Di Criscienzo, Marcella; D'Odorico, Valentina; Dotto, Elisabetta; Fontanot, Fabio; Gai, Mario; Gallerani, Simona; Gallozzi, Stefano; Garilli, Bianca; Gioia, Isabella; Girardi, Marisa; Gitti, Myriam; Granato, Gianluigi; Gratton, Raffaele; Grazian, Andrea; Gruppioni, Carlotta; Hunt, Leslie; Leto, Giuseppe; Israel, Gianluca; Magliocchetti, Manuela; Magrini, Laura; Mainetti, Gabriele; Mannucci, Filippo; Marconi, Alessandro; Marelli, Martino; Maris, Michele; Matteucci, Francesca; Meneghetti, Massimo; Mennella, Aniello; Mercurio, Amata; Molendi, Silvano; Monaco, Pierluigi; Moretti, Alessia; Murante, Giuseppe; Nicastro, Fabrizio; Orio, Marina; Paizis, Adamantia; Panessa, Francesca; Pasian, Fabio; Pentericci, Laura; Pozzetti, Lucia; Rossetti, Mariachiara; Santos, Joana S.; Saro, Alexandro; Schneider, Raffaella; Silva, Laura; Silvotti, Roberto; Smart, Richard; Tiengo, Andrea; Tornatore, Luca; Tozzi, Paolo; Trussoni, Edoardo; Valentinuzzi, Tiziano; Vanzella, Eros; Vazza, Franco; Vecchiato, Alberto; Venturi, Tiziana; Vianello, Giacomo; Viel, Matteo; Villalobos, Alvaro; Viotto, Valentina; Vulcani, Benedetta
Comments: Also available at
Submitted: 2010-07-08
On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.
[122]  [pdf] - 1033580
The colour-magnitude relation of Elliptical and Lenticular galaxies in the ESO Distant Cluster Survey
Comments: 14 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2010-07-08
We study the colour-magnitude relation (CMR) for a sample of 172 morphologically-classified E/S0 cluster galaxies from the ESO Distant Cluster Survey (EDisCS) at 0.4<z<0.8. The intrinsic colour scatter about the CMR is very small (0.076) in rest-frame U-V. Only 7% of the galaxies are significantly bluer than the CMR. The scarcity of blue S0s indicates that, if they are the descendants of spirals, these were already red when they became S0s. We observe no dependence of the CMR scatter with redshift or cluster velocity dispersion. This implies that by the time cluster E/S0s achieve their morphology, the vast majority have already joined the red sequence. We estimate the galaxy formation redshift z_F for each cluster and find that it does not depend on the cluster velocity dispersion. However, z_F increases weakly with cluster redshift. This trend becomes clearer when including higher-z clusters from the literature, suggesting that, at any given z, in order to have a population of fully-formed E and S0s they needed to have formed most of their stars 2-4 Gyr prior to observation. In other words, the galaxies that already have early-type (ET) morphologies also have reasonably-old stellar populations. This is partly a manifestation of the "progenitor bias", but also a consequence of the fact that the vast majority of the ETs in clusters (in particular the massive ones) were already red by the time they achieved their morphology. E and S0 galaxies exhibit very similar colour scatter, implying similar stellar population ages. We also find that fainter ETs finished forming their stars later, consistent with the cluster red sequence being built over time and the brightest galaxies reaching the red sequence earlier than fainter ones. Finally, we find that the ET cluster galaxies must have had their star formation truncated over an extended period of at least 1 Gyr. [abridged]
[123]  [pdf] - 1018432
The VIMOS VLT Deep Survey: the group catalogue
Comments: Submitted to A&A, revised version after referee comments, Table 5 added
Submitted: 2009-11-19, last modified: 2010-05-03
[Abridged] We present a homogeneous and complete catalogue of optical groups identified in the purely flux limited (17.5<=I<=24.0) VIMOS-VLT Deep Survey (VVDS). We use mock catalogues extracted from the MILLENNIUM simulation, to correct for potential systematics that might affect the overall distribution as well as the individual properties of the identified systems. Simulated samples allow us to forecast the number and properties of groups that can be potentially found in a survey with VVDS-like selection functions. We use them to correct for the expected incompleteness and also to asses how well galaxy redshifts trace the line-of-sight velocity dispersion of the underlying mass overdensity. In particular, we train on these mock catalogues the adopted group-finding technique (the Voronoi-Delaunay Method, VDM). The goal is to fine-tune its free parameters, recover in a robust and unbiased way the redshift and velocity dispersion distributions of groups and maximize the level of completeness (C) and purity (P) of the group catalogue. We identify 318 VVDS groups with at least 2 members within 0.2<=z<=1.0, among which 144 (/30) with at least 3 (/5) members. The sample has globally C=60% and P=50%. Nearly 45% of the groups with at least 3 members are still recovered if we run the algorithm with a parameter set which maximizes P (75%). We exploit the group sample to study the redshift evolution of the fraction f_b of blue galaxies (U-B<=1) within 0.2<=z<=1. We find that f_b is significantly lower in groups than in the whole ensemble of galaxies irrespectively of their environment. These quantities increase with redshift, with f_b in groups showing a marginally significant steeper increase. We also confirm that, at any explored redshift, f_b decreases for increasing group richness, and we extend towards fainter luminosities the magnitude range over which this result holds.
[124]  [pdf] - 1018772
Cluster Galaxies Die Hard
Comments: Models of ram-pressure stripping and some extra discussion added, references added. Conclusions unchanged. 20 pages, 15 figures. Accepted for publication in MNRAS.
Submitted: 2009-12-14, last modified: 2010-04-16
We investigate how the specific star formation rates of galaxies of different masses depend on cluster-centric radius and on the central/satellite dichotomy in both field and cluster environments. Recent data from a variety of sources, including the cluster catalogue of von der Linden et al. are compared to the semi-analytic models of De Lucia & Blaizot. We find that these models predict too many passive satellite galaxies in clusters, too few passive central galaxies with low stellar masses, and too many passive central galaxies with high masses. We then outline a series of modifications to the model necessary to solve these problems: a) Instead of instantaneous stripping of the external gas reservoir after a galaxy becomes a satellite, the gas supply is assumed to decrease at the same rate that the surrounding halo loses mass due to tidal stripping, b) The AGN feedback efficiency is lowered to bring the fraction of massive passive centrals in better agreement with the data. We also allow for radio mode AGN feedback in satellite galaxies. c) We assume that satellite galaxies residing in host haloes with masses below 10^12 M_sun do not undergo any stripping. We highlight the fact that in low mass galaxies, the external reservoir is composed primarily of gas that has been expelled from the galactic disk by supernovae driven winds. This gas must remain available as a future reservoir for star formation, even in satellite galaxies. Finally, we present a simple recipe for the stripping of gas and dark matter in satellites that can be used in models where subhalo evolution is not followed in detail.
[125]  [pdf] - 268061
Galactic stellar haloes in the CDM model
Comments: 24 pages, 19 figures; accepted for publication in MNRAS, minor changes; version with high resolution images at
Submitted: 2009-10-19, last modified: 2010-04-14
We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]
[126]  [pdf] - 406802
A semi-analytic model comparison - gas cooling and galaxy mergers
Comments: 21 pages, 14 figures. Accepted for publication in MNRAS.
Submitted: 2010-03-15, last modified: 2010-04-12
We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.
[127]  [pdf] - 1025260
The evolution of the density of galaxy clusters and groups: denser environments at higher redshifts
Comments: accepted for publication in MNRAS
Submitted: 2010-02-18
We show that, observationally, the projected local density distribution in high-z clusters is shifted towards higher values compared to clusters at lower redshift. To search for the origin of this evolution, we analyze a sample of haloes selected from the Millennium Simulation and populated using semi-analytic models, investigating the relation between observed projected density and physical 3D density, using densities computed from the 10 and 3 closest neighbours. Both observationally and in the simulations, we study the relation between number of cluster members and cluster mass, and number of members per unit of cluster mass. We find that the observed evolution of projected densities reflects a shift to higher values of the physical 3D density distribution. In turn, this must be related with the globally higher number of galaxies per unit of cluster volume N/V in the past. We show that the evolution of N/V is due to a combination of two effects: a) distant clusters were denser in dark matter (DM) simply because the DM density within R_{200} (~the cluster virial radius) is defined to be a fixed multiple of the critical density of the Universe, and b) the number of galaxies per unit of cluster DM mass is remarkably constant both with redshift and cluster mass if counting galaxies brighter than a passively evolving magnitude limit. Our results highlight that distant clusters were much denser environments than today's clusters, both in galaxy number and mass, and that the density conditions felt by galaxies in virialized systems do not depend on the system mass.
[128]  [pdf] - 28084
On the nature of the Milky Way satellites
Comments: 18 pages, 16 figures. Minor changes of the presentation to match the published version
Submitted: 2009-09-07, last modified: 2010-01-18
We combine a series of high-resolution simulations with semi-analytic galaxy formation models to follow the evolution of a system resembling the Milky Way and its satellites. The semi-analytic model is based on that developed for the Millennium Simulation, and successfully reproduces the properties of galaxies on large scales, as well as those of the Milky Way. In this model, we are able to reproduce the luminosity function of the satellites around the Milky Way by preventing cooling in haloes with Vvir < 16.7 km/s (i.e. the atomic hydrogen cooling limit) and including the impact of the reionization of the Universe. The physical properties of our model satellites (e.g. mean metallicities, ages, half-light radii and mass-to-light ratios) are in good agreement with the latest observational measurements. We do not find a strong dependence upon the particular implementation of supernova feedback, but a scheme which is more efficient in galaxies embedded in smaller haloes, i.e. shallower potential wells, gives better agreement with the properties of the ultra-faint satellites. Our model predicts that the brightest satellites are associated with the most massive subhaloes, are accreted later (z $\lta$ 1), and have extended star formation histories, with only 1 per cent of their stars made by the end of the reionization. On the other hand, the faintest satellites were accreted early, are dominated by stars with age > 10 Gyr, and a few of them formed most of their stars before the reionization was complete. Objects with luminosities comparable to those of the classical MW satellites are associated with dark matter subhaloes with a peak circular velocity $\gta$ 10 km/s, in agreement with the latest constraints.
[129]  [pdf] - 1024682
Gas cooling in semi-analytic models and SPH simulations: are results consistent?
Comments: Revised version submitted to MNRAS, 15 pages, 9 figures. A High-resolution version of the paper and figures can be found at this
Submitted: 2010-01-18
We present a detailed comparison between the galaxy populations within a massive cluster, as predicted by hydrodynamical SPH simulations and by a semi-analytic model (SAM) of galaxy formation. Both models include gas cooling and a simple prescription of star formation, which consists in transforming instantaneously any cold gas available into stars, while neglecting any source of energy feedback. We find that, in general, galaxy populations from SAMs and SPH have similar statistical properties, in agreement with previous studies. However, when comparing galaxies on an object-by-object basis, we find a number of interesting differences: a) the star formation histories of the brightest cluster galaxies (BCGs) from SAM and SPH models differ significantly, with the SPH BCG exhibiting a lower level of star formation activity at low redshift, and a more intense and shorter initial burst of star formation with respect to its SAM counterpart; b) while all stars associated with the BCG were formed in its progenitors in the semi-analytic model used here, this holds true only for half of the final BCG stellar mass in the SPH simulation, the remaining half being contributed by tidal stripping of stars from the diffuse stellar component associated with galaxies accreted on the cluster halo; c) SPH satellites can loose up to 90 per cent of their stellar mass at the time of accretion, due to tidal stripping, a process not included in the semi-analytic model used in this study; d) in the SPH simulation, significant cooling occurs on the most massive satellite galaxies and this lasts for up to 1 Gyr after accretion. This physical process is not included in the semi-analytic model used in our study, as well as in most of the models discussed in the recent literature.
[130]  [pdf] - 1018702
Metallicity distributions in and around galaxies
Comments: Invited Review given at the Southern Cross Astrophysics Conference on "Galaxy Metabolism" held in Sydney, 22-26 June 2009. Accepted by the Publications of the Astronomical Society of Australia (PASA)
Submitted: 2009-12-09
Metals are found in all baryonic phases and environments, and our knowledge of their distribution `in and around galaxies' has significantly improved over the past few years. Theoretical work has shown that the fraction of metals in different baryonic components can vary significantly when different feedback schemes are adopted. Therefore, studies of element abundances provide important information about all gas-dynamical processes which determine the cosmic evolution of baryons. I give here a brief review of recent observational progress, describe the implications of recent theoretical studies, and discuss briefly future prospects.
[131]  [pdf] - 1018700
Stellar Ages and Metallicities of Central and Satellite Galaxies: Implications for Galaxy Formation and Evolution
Comments: 20 pages, 12 figures, submitted for publication in MNRAS
Submitted: 2009-12-09
Using a large SDSS galaxy group catalogue, we study how the stellar ages and metallicities of central and satellite galaxies depend on stellar mass and halo mass. We find that satellites are older and metal-richer than centrals of the same stellar mass. In addition, the slopes of the age-stellar mass and metallicity-stellar mass relations are found to become shallower in denser environments. This is due to the fact that the average age and metallicity of low mass satellite galaxies increase with the mass of the halo in which they reside. A comparison with the semi-analytical model of Wang et al. (2008) shows that it succesfully reproduces the fact that satellites are older than centrals of the same stellar mass and that the age difference increases with the halo mass of the satellite. This is a consequence of strangulation, which leaves the stellar populations of satellites to evolve passively, while the prolonged star formation activity of centrals keeps their average ages younger. The resulting age offset is larger in more massive environments because their satellites were accreted earlier. The model fails, however, in reproducing the halo mass dependence of the metallicities of low mass satellites, yields metallicity-stellar mass and age-stellar mass relations that are too shallow, and predicts that satellite galaxies have the same metallicities as centrals of the same stellar mass, in disagreement with the data. We argue that these discrepancies are likely to indicate the need to (i) modify the recipes of both supernova feedback and AGN feedback, (ii) use a more realistic description of strangulation, and (iii) include a proper treatment of the tidal stripping, heating and destruction of satellite galaxies. [Abridged]
[132]  [pdf] - 1018127
Photometric redshifts and cluster tomography in the ESO Distant Cluster Survey
Comments: 26 pages, 22 figures, accepted for publication in A&A High-resolution version of the paper and figures can be found at
Submitted: 2009-10-28
This paper reports the results obtained on the photometric redshifts measurement and accuracy, and cluster tomography in the ESO Distant Cluster Survey (EDisCS) fields. Photometric redshifts were computed using two independent codes (Hyperz and G. Rudnick's code). The accuracy of photometric redshifts was assessed by comparing our estimates with the spectroscopic redshifts of ~1400 galaxies in the 0.3<z<1.0 domain. The accuracy for galaxies fainter than the spectroscopic control sample was estimated using a degraded version of the photometric catalog for the spectroscopic sample. The accuracy of photometric redshifts is typically sigma(Delta z/(1+z)) ~ 0.05+/-0.01, depending on the field, the filter set, and the spectral type of the galaxies. The quality of the photometric redshifts degrades by a factor of two in sigma(Delta z/(1+z)) between the brightest (I~22) and the faintest (I~24-24.5) galaxies in the EDisCS sample. The photometric determination of cluster redshifts in the EDisCS fields using a simple algorithm based on photoz is in excellent agreement with the spectroscopic values (delta(z) ~0.03-0.04 in the high-z sample and ~0.05 in the low-z sample). We also developed a method that uses both photz codes jointly to reject interlopers at magnitudes fainter than the spectroscopic limit. When applied to the spectroscopic sample, this method rejects ~50-90% of all spectroscopically confirmed non-members, while retaining at least 90% of all confirmed members. Photometric redshifts are found to be particularly useful for the identification and study of clusters of galaxies in large surveys. They enable efficient and complete pre-selection of cluster members for spectroscopy, allow accurate determinations of the cluster redshifts, and provide a means of determining cluster membership. (Abridged)
[133]  [pdf] - 1017894
Evolution of the Early-Type Galaxy Fraction in Clusters since z = 0.8
Comments: 22 pages, 10 figures, accepted for publication in A&A
Submitted: 2009-10-08
We study the morphological content of a large sample of high-redshift clusters to determine its dependence on cluster mass and redshift. Quantitative morphologies are based on bulge+disk decompositions of cluster and field galaxies on deep VLT/FORS2 images of 18 optically-selected clusters at 0.45 < z < 0.80 from the ESO Distant Cluster Survey (EDisCS). Morphological content is given by the early-type galaxy fraction f_et, and early-type galaxies are selected based on their bulge fraction and image smoothness. A set of 158 SDSS clusters is analyzed exactly as the EDisCS sample to provide a robust local comparison. Our main results are: (1) f_et values for the SDSS and EDisCS clusters exhibit no clear trend as a function of sigma. (2) Mid-z EDisCS clusters around sigma = 500 km/s have f_et ~= 0.5 whereas high-z EDisCS clusters have f_et ~= 0.4 (~25% increase over 2 Gyrs). (3) There is a marked difference in the morphological content of EDisCS and SDSS clusters. None of the EDisCS clusters have f_et greater than 0.6 whereas half of the SDSS clusters lie above this value. This difference is seen in clusters of all velocity dispersions. (4) There is a strong correlation between morphology and star formation in SDSS and EDisCS clusters. This correlation holds independent of sigma and z even though the fraction of [OII] emitters decreases from z~0.8 to z~0.06 in all environments. Our results pose an interesting challenge to structural transformation and star formation quenching processes that strongly depend on the global cluster environment and suggest that cluster membership may be of lesser importance than other variables in determining galaxy properties. (ABRIDGED)
[134]  [pdf] - 1003497
Properties of Long Gamma-Ray Burst Host Galaxies in Cosmological Simulation
Comments: The paper contains 14 figures and 2 tables. Accepted MNRAS
Submitted: 2009-08-17, last modified: 2009-08-19
We use galaxy catalogues constructed by combining high-resolution N-body simulations with semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst (LGRB) host galaxies. We assume that LGRBs originate from the death of massive young stars and analyse how results are affected by different metallicity constraints on the progenitor stars. As expected, the host sample with no metallicity restriction on the progenitor stars provides a perfect tracer of the cosmic star formation history. When LGRBs are required to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star formation rate at lower redshift, as a consequence of the global increase in metallicity. We study the properties of host galaxies up to high redshift (~9), finding that they typically have low-metallicity (Z<0.5 Z_sun) and that they are small (M<10^9 M_sun), bluer and younger than the average galaxy population, in agreement with observational data. They are also less clustered than typical L_* galaxies in the Universe, and their descendents are massive, red and reside in groups of galaxies with halo mass between 10^{13} M_sun to 10^{14} M_sun.
[135]  [pdf] - 1003047
On the role of the post-starburst phase in the build-up of the red-sequence of intermediate redshift clusters
Comments: 11 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2009-07-22
We present new deep spectroscopic observations of 0.05--0.5 L_* galaxies in one cluster (cl1232.5-1250) drawn from the ESO Distant Cluster Survey (EDisCS) sample, at z = 0.54. The new data extend the spectroscopy already available for this cluster by about 1 magnitude. The cluster has a large fraction of passive galaxies and exhibits a well defined and relatively tight colour-magnitude relation. Among spectroscopic members, only six galaxies are classified as `post-starburst' (k+a). For another EDisCS cluster at similar redshift and with as deep spectroscopy, no member is found to have a k+a spectrum. The low measured numbers of post-starburst systems appear to be inadequate to explain the observed increase of faint red galaxies at lower redshift, even when accounting for the infall of new galaxies onto the cluster. Post-starburst galaxies represent a possible channel to move galaxies from the blue star-forming cloud to the red-sequence, but the available data suggest this is not the dominant channel in galaxy clusters. If the galaxies at the faint end of the red-sequence in nearby clusters originate from the blue galaxies observed in distant galaxy clusters, the transformation must have occurred primarily through physical processes that do not lead to a post-starburst phase. In addition, our data exclude a large population of k+a galaxies at faint magnitudes, as found in the Coma cluster.
[136]  [pdf] - 315828
Simulation of the Cosmic Evolution of Atomic and Molecular Hydrogen in Galaxies
Comments: 19 pages, 15 figures
Submitted: 2009-04-14, last modified: 2009-07-15
We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3e7 galaxies, obtained by post-processing the virtual galaxy catalog produced by (De Lucia et al. 2007) on the Millennium Simulation of cosmic structure (Springel et al. 2005). Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local HI- and H2-mass functions, the CO-luminosity function, the cold gas mass--diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local Universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift.
[137]  [pdf] - 17196
On the common mass scale of the Milky Way satellites
Comments: 5 pages, 3 figures, to match the version accepted in MNRAS Letters. One paragraph added to explain the spread in luminosities more explicitly. References updated and Conclusions remain unchanged
Submitted: 2008-10-07, last modified: 2009-06-16
We use a hybrid approach that combines high-resolution simulations of the formation of a Milky Way-like halo with a semi-analytic model of galaxy formation to study the mass content of dwarf galaxies in the concordance $\Lambda$CDM cosmology. We find that the mass within 600 pc of dark matter haloes hosting luminous satellites has a median value of $\sim 3.2\times 10^7$ Msun with very little object-to-object scatter. In contrast, the present day total luminosities of the model satellites span nearly five orders of magnitude. These findings are in very good agreement with the results recently reported in the literature for the dwarf spheroidal galaxies of the Milky Way. In our model, dwarf irregular galaxies like the Small Magellanic Cloud, are predicted to have similar or slightly larger dark matter mass within 600 pc.
[138]  [pdf] - 24697
The Rest-Frame Optical Luminosity Function of Cluster Galaxies at z<0.8 and the Assembly of the Cluster Red Sequence
Comments: Accepted for publication in the Astrophysical Journal. 36 pages, 16 figures, 10 tables
Submitted: 2009-05-27
We present the rest-frame optical luminosity function (LF) of red sequence galaxies in 16 clusters at 0.4<z<0.8 drawn from the ESO Distant Cluster Survey (EDisCS). We compare our clusters to an analogous sample from the Sloan Digital Sky Survey (SDSS) and match the EDisCS clusters to their most likely descendants. We measure all LFs down to M M* + (2.5 - 3.5). At z<0.8, the bright end of the LF is consistent with passive evolution but there is a significant build-up of the faint end of the red sequence towards lower redshift. There is a weak dependence of the LF on cluster velocity dispersion for EDisCS but no such dependence for the SDSS clusters. We find tentative evidence that red sequence galaxies brighter than a threshold magnitude are already in place, and that this threshold evolves to fainter magnitudes toward lower redshifts. We compare the EDisCS LFs with the LF of co-eval red sequence galaxies in the field and find that the bright end of the LFs agree. However, relative to the number of bright red galaxies, the field has more faint red galaxies than clusters at 0.6<z<0.8 but fewer at 0.4<z<0.6, implying differential evolution. We compare the total light in the EDisCS cluster red sequences to the total red sequence light in our SDSS cluster sample. Clusters at 0.4<z<0.8 must increase their luminosity on the red sequence (and therefore stellar mass in red galaxies) by a factor of 1-3 by z=0. The necessary processes that add mass to the red sequence in clusters predict local clusters that are over-luminous as compared to those observed in the SDSS. The predicted cluster luminosities can be reconciled with observed local cluster luminosities by combining multiple previously known effects.
[139]  [pdf] - 1001283
The Many Manifestations of Downsizing: Hierarchical Galaxy Formation Models confront Observations
Comments: MNRAS accepted, 16 pages, 10 figures
Submitted: 2009-01-08, last modified: 2009-05-13
[abridged] It has been widely claimed that several lines of observational evidence point towards a "downsizing" (DS) of the process of galaxy formation over cosmic time. This behavior is sometimes termed "anti-hierarchical", and contrasted with the "bottom-up" assembly of the dark matter structures in Cold Dark Matter models. In this paper we address three different kinds of observational evidence that have been described as DS: the stellar mass assembly, star formation rate and the ages of the stellar populations in local galaxies. We compare a broad compilation of available data-sets with the predictions of three different semi-analytic models of galaxy formation within the Lambda-CDM framework. In the data, we see only weak evidence at best of DS in stellar mass and in star formation rate. We find that, when observational errors on stellar mass and SFR are taken into account, the models acceptably reproduce the evolution of massive galaxies, over the entire redshift range that we consider. However, lower mass galaxies are formed too early in the models and are too passive at late times. Thus, the models do not correctly reproduce the DS trend in stellar mass or the archaeological DS, while they qualitatively reproduce the mass-dependent evolution of the SFR. We demonstrate that these discrepancies are not solely due to a poor treatment of satellite galaxies but are mainly connected to the excessively efficient formation of central galaxies in high-redshift haloes with circular velocities ~100-200 km/s. [abridged]
[140]  [pdf] - 23771
The ESO Distant Cluster Sample: galaxy evolution and environment out to z=1
Comments: to appear in The Messenger, issue June 2009
Submitted: 2009-04-29
The ESO Distant Cluster Survey (EDisCS, P.I. Simon D.M. White, LP 166.A-0162) is an ESO large programme aimed at studying clusters and cluster galaxies at z=0.4-1. How different is the evolution of the star formation activity in clusters, in groups and in the field? Does it depend on cluster mass and/or the local galaxy density? How relevant are starburst and post-starburst galaxies in the different environments? Is there an evolution in the galaxies' structures, and if so, is this related to the changes in their star formation activity? These are some of the main questions that have been investigated using the EDisCS dataset.
[141]  [pdf] - 1001990
"Ab initio" models of galaxy formation: successes and open problems
Comments: 8 pages, to appear in "Probing Stellar Populations out to the Distant Universe", AIP Conf. Proc. Series (Cefalu, Italy, Sep 7-19, 2008)
Submitted: 2009-04-21
In the past decades, different approaches have been developed in order to link the physical properties of galaxies to the dark matter haloes in which they reside. In this review, I give a brief overview of methods, aims, and limits of these techniques, with particular emphasis on semi-analytic models of galaxy formation. For these models, I also provide a brief summary of recent successes and open problems.
[142]  [pdf] - 21706
Frequency and properties of bars in cluster and field galaxies at intermediate redshifts
Comments: 17 pages, accepted for publication in A&A
Submitted: 2009-02-24
We present a study of large-scale bars in field and cluster environments out to redshifts of ~0.8 using a final sample of 945 moderately inclined disk galaxies drawn from the EDisCS project. We characterize bars and their host galaxies and look for relations between the presence of a bar and the properties of the underlying disk. We investigate whether the fraction and properties of bars in clusters are different from their counterparts in the field. The total optical bar fraction in the redshift range z=0.4-0.8 (median z=0.60), averaged over the entire sample, is 25% (20% for strong bars). For the cluster and field subsamples, we measure bar fractions of 24% and 29%, respectively. We find that bars in clusters are on average longer than in the field and preferentially found close to the cluster center, where the bar fraction is somewhat higher (~31%) than at larger distances (~18%). These findings however rely on a relatively small subsample and might be affected by small number statistics. In agreement with local studies, we find that disk-dominated galaxies have a higher optical bar fraction (~45%) than bulge-dominated galaxies (~15%). This result is based on Hubble types and effective radii and does not change with redshift. The latter finding implies that bar formation or dissolution is strongly connected to the emergence of the morphological structure of a disk and is typically accompanied by a transition in the Hubble type. (abridged)
[143]  [pdf] - 1001487
Evolution of Cluster Red-Sequence Galaxies from redshift 0.8 to 0.4: ages, metallicities and morphologies
Comments: Accepted for publication in A&A
Submitted: 2009-02-19
We present a comprehensive analysis of the stellar population properties (age, metallicity and the alpha-element enhancement [E/Fe]) and morphologies of red-sequence galaxies in 24 clusters and groups from z~0.75 to z~0.45. The dataset, consisting of 215 spectra drawn from the ESO Distant Cluster Survey, constitutes the largest spectroscopic sample at these redshifts for which such an analysis has been conducted. Analysis reveals that the evolution of the stellar population properties of red-sequence galaxies depend on their mass: while the properties of most massive are well described by passive evolution and high-redshift formation, the less massive galaxies require a more extended star formation history. We show that these scenarios reproduce the index-sigma relations as well as the galaxy colours. The two main results of this work are (1) the evolution of the line-strength indices for the red-sequence galaxies can be reproduced if 40% of the galaxies with sigma < 175 km/s entered the red-sequence between z=0.75 to z=0.45, in agreement with the fraction derived in studies of the luminosity functions, and (2) the percentage the red-sequence galaxies exhibiting early-type morphologies (E and S0) decreases by 20% from z=0.75 to z=0.45. This can be understood if the red-sequence gets populated at later times with disc galaxies whose star formation has been quenched. We conclude that the processes quenching star formation do not necessarily produce a simultaneous morphological transformation of the galaxies entering the red-sequence.
[144]  [pdf] - 13651
Milky Way type galaxies in a LCDM cosmology
Comments: 10 pages, 7 figures, mne2.cls, MNRAS, replaced with accepted version
Submitted: 2008-06-17, last modified: 2009-01-29
We analyse a sample of 52,000 Milky Way (MW) type galaxies drawn from the publicly available galaxy catalogue of the Millennium Simulation with the aim of studying statistically the differences and similarities of their properties in comparison to our Galaxy. Model galaxies are chosen to lie in haloes with maximum circular velocities in the range 200-250 km/seg and to have bulge-to-disk ratios similar to that of the Milky Way. We find that model MW galaxies formed quietly through the accretion of cold gas and small satellite systems. Only 12 per cent of our model galaxies experienced a major merger during their lifetime. Most of the stars formed in situ, with only about 15 per cent of the final mass gathered through accretion. Supernovae and AGN feedback play an important role in the evolution of these systems. At high redshifts, when the potential wells of the MW progenitors are shallower, winds driven by supernovae explosions blow out a large fraction of the gas and metals. As the systems grow in mass, SN feedback effects decrease and AGN feedback takes over, playing a more important role in the regulation of the star formation activity at lower redshifts. Although model Milky Way galaxies have been selected to lie in a narrow range of maximum circular velocities, they nevertheless exhibit a significant dispersion in the final stellar masses and metallicities. Our analysis suggests that this dispersion results from the different accretion histories of the parent dark matter haloes. Statically, we also find evidences to support the Milky Way as a typical Sb/Sc galaxy in the same mass range, providing a suitable benchmark to constrain numerical models of galaxy formation
[145]  [pdf] - 17446
LCDM predictions for galaxy protoclusters I: the relation between galaxies, protoclusters and quasars at z~6
Comments: Accepted by MNRAS (19 pages, 20 figures, 2 tables). High resolution PDF available at
Submitted: 2008-10-14, last modified: 2008-12-05
Motivated by recent observational studies of the environment of z~6 QSOs, we have used the Millennium Run (MR) simulations to construct a very large (~20 deg^2) mock redshift survey of star-forming galaxies at z~6. We use this simulated survey to study the relation between density enhancements in the distribution of i-dropouts and Lya-emitters, and their relation to the most massive halos and protocluster regions at z~6. Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1 degree, much larger than probed by current surveys. Approximately one third of all z~6 halos hosting i-dropouts brighter than z=26.5 mag (~M* at z=6) become part of z=0 galaxy clusters. i-dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyze two structures of i-dropouts and Lya-emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters-in-formation. In striking contrast, six z~6 QSO fields observed with HST show no significant enhancements in their i-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive halos. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space based instruments that will shed new light on questions related to the large-scale structure at z~6.
[146]  [pdf] - 436203
Quenching of Star Formation
Comments: 7 pages, 3 figures. Conference proceedings in "Classification and Discovery in Large Astronomical Surveys", 2008, C.A.L. Bailer-Jones (ed.)
Submitted: 2008-11-07
In the last decade we have seen an enormous increase in the size and quality of spectroscopic galaxy surveys, both at low and high redshift. New statistical techniques to analyse large portions of galaxy spectra are now finding favour over traditional index based methods. Here we will review a new robust and iterative Principal Component Analysis (PCA) algorithm, which solves several common issues with classic PCA. Application to the 4000AA break region of galaxies in the VIMOS VLT Deep Survey (VVDS) and Sloan Digital Sky Survey (SDSS) gives new high signal-to-noise ratio spectral indices easily interpretable in terms of recent star formation history. In particular, we identify a sample of post-starburst galaxies at z~0.7 and z~0.07. We quantify for the first time the importance of post-starburst galaxies, consistent with being descendants of gas-rich major mergers, for building the red sequence. Finally, we present a comparison with new low and high redshift "mock spectroscopic surveys" derived from a Millennium Run semi-analytic model.
[147]  [pdf] - 18121
The environments of starburst and post-starburst galaxies at z=0.4-0.8
Comments: accepted for publication in ApJ
Submitted: 2008-11-03, last modified: 2008-11-07
Post-starburst (E+A or k+a) spectra, characterized by their exceptionally strong Balmer lines in absorption and the lack of emission lines, belong to galaxies in which the star formation activity ended abruptly sometime during the past Gyr. We perform a spectral analysis of galaxies in clusters, groups, poor groups and the field at z=0.4-0.8 based on the ESO Distant Cluster Survey. The incidence of k+a's at these redshifts depends strongly on environment. K+a's reside preferentially in clusters and, unexpectedly, in a subset of the sigma = 200-400 km/s groups, those that have a low fraction of [OII] emitters. In these environments, 20-30% of the recently star-forming galaxies have had their star formation activity recently truncated. In contrast, there are proportionally fewer k+a's in the field, the poor groups and groups with a high [OII] fraction. The incidence of k+a galaxies correlates with the cluster velocity dispersion: more massive clusters have higher proportions of k+a's. Spectra of dusty starburst candidates, with strong Balmer absorption and emission lines, present a very different environmental dependence from k+a's. They are numerous in all environments at z=0.4-0.8, but they are especially numerous in all types of groups, favoring the hypothesis of triggering by a merger. Our observations are consistent with previous suggestions that cluster k+a galaxies are observed in a transition phase as massive S0 and Sa galaxies, evolving from star-forming later types to passively evolving early-type galaxies. The correlation between k+a fraction and cluster sigma supports the hypothesis that k+a galaxies in clusters originate from processes related to the intracluster medium, while several possibilities are discussed for the origin of the k+a frequency in low-[OII] groups.(abr.)
[148]  [pdf] - 17489
Simulating the formation of a proto-cluster at z~2
Comments: 6 pages, 4 figures, accepted for publication in MNRAS (Letters)
Submitted: 2008-10-15
We present results from two high-resolution hydrodynamical simulations of proto-cluster regions at z~2.1. The simulations have been compared to observational results for the socalled Spiderweb galaxy system, the core of a putative proto-cluster region at z = 2.16, found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with M200~10^14 h-1 Msun (C1) and a rich cluster with M200~2x10^15 h-1 Msun (C2) at z = 0. The simulated proto-clusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy (BCG) of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared to the observed velocities. We argue that the Spiderweb complex resemble the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing AGN feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.
[149]  [pdf] - 11795
The Galaxy and its stellar halo: insights on their formation from a hybrid cosmological approach
Comments: 20 pages, 14 figures. Minor revisions to match version accepted to MNRAS. Version with high-resolution figures available at:
Submitted: 2008-04-15, last modified: 2008-09-02
We use a series of high-resolution simulations of a `Milky-Way' halo coupled to semi-analytic methods to study the formation of our own Galaxy and of its stellar halo. The physical properties of our model Milky Way, as well as the age and metallicity distribution of stars in the different components, are in relatively good agreement with observational measurements. Assuming that the stellar halo builds up from the cores of the satellite galaxies that merged with the Milky Way over its life-time, we are able to study the physical and structural properties of this component. In agreement with previous work, we find that the largest contribution to the stellar halo should come from a few relatively massive (10^8 - 10^10 Msun) satellites, accreted at early times. Our "stellar halo" does not exhibit any clear metallicity gradient, but higher metallicity stars are more centrally concentrated than stars of lower abundance. This indicates that the probability of observing low-metallicity halo stars increases with distance from the Galactic centre. We find that the proposed "dual" nature of the Galactic stellar halo can be explained in our model as a result of a mass-metallicity relation imprinted in the building blocks of this component.
[150]  [pdf] - 15851
The effect of gas-dynamics on semi-analytic modeling of cluster galaxies
Comments: 13 pages, 12 figures, to appear in MNRAS
Submitted: 2008-09-01
We study the degree to which non-radiative gas dynamics affects the merger histories of haloes along with subsequent predictions from a semi-analytic model (SAM) of galaxy formation. To this aim, we use a sample of dark matter only and non-radiative SPH simulations of four massive clusters. The presence of gas-dynamical processes (e.g. ram-pressure from the hot ICM) makes haloes more fragile in the runs which include gas. This results in a 25 per cent decrease in the total number of subhaloes at z = 0. The impact on the galaxy population predicted by SAMs is complicated by the presence of "orphan" galaxies, i.e. galaxies whose parent substructures are reduced below the resolution limit of the simulation. In the model employed in our study, these galaxies survive for a residual merging time that is computed using a variation of the Chandrasekhar formula. Due to ram-pressure stripping, haloes in gas simulations tend to be less massive than their counterparts in the dark matter simulations. The resulting merging times for satellite galaxies are then longer in these simulations. On the other hand, the presence of gas influences the orbits of haloes making them on average more circular and therefore reducing the estimated merging times with respect to the dark matter only simulation. This effect is particularly significant for the most massive satellites and is (at least in part) responsible for the fact that brightest cluster galaxies in runs with gas have stellar masses which are about 25 per cent larger than those obtained from dark matter only simulations. Our results show that gas-dynamics has only a marginal impact on the statistical properties of the galaxy population, but that its impact on the orbits and merging times of haloes strongly influences the assembly of the most massive galaxies.
[151]  [pdf] - 1000813
The Galaxy and its stellar halo - insights from a hybrid cosmological approach
Comments: Proceedings of IAU Symposium 254 "The Galaxy disk in a cosmological context", Copenhagen, June 2008, 6 pages, 4 pages
Submitted: 2008-07-07
We use a series of high-resolution N-body simulations of a `Milky-Way' halo, coupled to semi-analytic techniques, to study the formation of our own Galaxy and of its stellar halo. Our model Milky Way galaxy is a relatively young system whose physical properties are in quite good agreement with observational determinations. In our model, the stellar halo is mainly formed from a few massive satellites accreted early on during the galaxy's lifetime. The stars in the halo do not exhibit any metallicity gradient, but higher metallicity stars are more centrally concentrated than stars with lower abundances. This is due to the fact that the most massive satellites contributing to the stellar halo are also more metal rich, and dynamical friction drags them closer to the inner regions of the host halo.
[152]  [pdf] - 12672
The VIMOS-VLT Deep Survey: Dependence of galaxy clustering on stellar mass
Comments: 4 pages, 2 figures, 43rd Rencontres de Moriond - March 15-22, 2008 - La Thuile (Val d'Aosta, Italy)
Submitted: 2008-05-15
We have investigated the dependence of galaxy clustering on their stellar mass at z~1, using the data from the VIMOS-VLT Deep Survey (VVDS). We have measured the projected two-point correlation function of galaxies, wp(rp) for a set of stellar mass selected samples at an effective redshift <z>=0.85. We have control and quantify all effects on galaxy clustering due to the incompleteness of our low mass samples. We find that more massive galaxies are more clustered. When compared to similar results at z~0.1 in the SDSS, we observed no evolution of the projected correlation function for massive galaxies. These objects present a stronger linear bias at z~1 with respect to low mass galaxies. As expected, massive objects at high redshift are found in the highest pics of the dark matter density field.
[153]  [pdf] - 12436
The relation between star formation, morphology and local density in high redshift clusters and groups
Comments: 21 pages, 14 figures, accepted for publication in ApJ
Submitted: 2008-05-08
We investigate how the [OII] properties and the morphologies of galaxies in clusters and groups at z=0.4-0.8 depend on projected local galaxy density, and compare with the field at similar redshifts and clusters at low-z. In both nearby and distant clusters, higher-density regions contain proportionally fewer star-forming galaxies, and the average [OII] equivalent width of star-forming galaxies is independent of local density. However, in distant clusters the average current star formation rate (SFR) in star-forming galaxies seems to peak at densities ~15-40 galaxies Mpc^{-2}. At odds with low-z results, at high-z the relation between star-forming fraction and local density varies from high- to low-mass clusters. Overall, our results suggest that at high-z the current star formation (SF) activity in star-forming galaxies does not depend strongly on global or local environment, though the possible SFR peak seems at odds with this conclusion. We find that the cluster SFR normalized by cluster mass anticorrelates with mass and correlates with the star-forming fraction. These trends can be understood given a) that the average star-forming galaxy forms about 1 Msun/yr in all clusters; b) that the total number of galaxies scales with cluster mass and c) the dependence of star-forming fraction on cluster mass. We present the morphology-density (MD) relation for our z=0.4-0.8 clusters, and uncover that the decline of the spiral fraction with density is entirely driven by galaxies of types Sc or later. For galaxies of a given Hubble type, we see no evidence that SF properties depend on local environment. In contrast with recent findings at low-z, in our distant clusters the SF-density relation and the MD-relation are equivalent, suggesting that neither of the two is more fundamental than the other.(abr.)
[154]  [pdf] - 12140
SPACE: the SPectroscopic All-sky Cosmic Explorer
Comments: 27 pages, Experimental Astronomy, in press. The SPACE team complete list is available at . The article with full resolution figures is available at
Submitted: 2008-04-28
We describe the scientific motivations, the mission concept and the instrumentation of SPACE, a class-M mission proposed for concept study at the first call of the ESA Cosmic-Vision 2015-2025 planning cycle. SPACE aims to produce the largest three-dimensional evolutionary map of the Universe over the past 10 billion years by taking near-IR spectra and measuring redshifts for more than half a billion galaxies at 0<z<2 down to AB~23 over 3\pi sr of the sky. In addition, SPACE will also target a smaller sky field, performing a deep spectroscopic survey of millions of galaxies to AB~26 and at 2<z<10+. These goals are unreachable with ground-based observations due to the ~500 times higher sky background. To achieve the main science objectives, SPACE will use a 1.5m diameter Ritchey-Chretien telescope equipped with a set of arrays of Digital Micro-mirror Devices (DMDs) covering a total field of view of 0.4 deg2, and will perform large-multiplexing multi-object spectroscopy (e.g. ~6000 targets per pointing) at a spectral resolution of R~400 as well as diffraction-limited imaging with continuous coverage from 0.8mum to 1.8mum.
[155]  [pdf] - 11740
The evolution of the brightest cluster galaxies since z~1 from the ESO Distant Cluster Survey (EDisCS)
Comments: Accepted for publication in MNRAS
Submitted: 2008-04-14
[Abridged] We present K-band data for the brightest cluster galaxies (BCGs) from the ESO Distant Cluster Survey. These data are combined with photometry from Aragon-Salamanca et al. (1998) and a low-redshift comparison sample from von der Linden et al. (2007). The K-band Hubble diagram for BCGs exhibits very low scatter (~0.35mag) since z=1. The colour and $K$-band luminosity evolution of the BCGs are in good agreement with passively-evolving stellar populations formed at z>2. We do not detect any significant change in the stellar mass of the BCG since z~1. These results do not seem to depend on the velocity dispersion of the parent cluster. There is a correlation between the 1D velocity dispersion of the clusters and the K-band luminosity of the BCGs (after correcting for passive evolution). The clusters with large velocity dispersions tend to have brighter BCGs, i.e., BCGs with larger stellar masses. This dependency, although significant, is relatively weak: the stellar mass of the BCGs changes only by ~70% over a two-order-of-magnitude range in cluster mass. This dependency doesn't change significantly with redshift. The models of De Lucia & Blaizot (2007) predict colours which are in reasonable agreement with the observations because the growth in stellar mass is dominated by the accretion of old stars. However, the stellar mass in the model BCGs grows by a factor of 3-4 since z=1, a growth rate which seems to be ruled out by the observations. The models predict a dependency between the BCG's stellar mass and the velocity dispersion of the parent cluster in the same sense as the data, but the dependency is significantly stronger than observed. However, one major difficulty in this comparison is that we have measured fixed metric aperture magnitudes while the models compute total luminosities.
[156]  [pdf] - 11663
The evolution of the colour-magnitude relation and of the star formation activity in galaxy clusters since z~0.8
Comments: 8 pages, 5 figure, to appear in ASP Conference Series (Proceedings of the 1st Subaru International Conference "Panoramic Views of Galaxy Formation and Evolution", held in Japan, 10-15 December 2007)
Submitted: 2008-04-10
We present recent results on the evolution of the colour-magnitude relation and of the star formation activity in galaxy clusters since z~0.8. Results are based on the ESO Distant Cluster Survey (EDisCS) - an ESO large programme aimed at the study of cluster structure and cluster galaxy evolution over a significant fraction of cosmic time - and are discussed in the framework of the current standard paradigm of structure formation.
[157]  [pdf] - 10065
A test of the nature of cosmic acceleration using galaxy redshift distortions
Comments: One PDF file including both main paper and Supplementary Information (28 pages, 3+2 figures). Published version available at
Submitted: 2008-02-13
Observations of distant supernovae indicate that the Universe is now in a phase of accelerated expansion the physical cause of which is a mystery. Formally, this requires the inclusion of a term acting as a negative pressure in the equations of cosmic expansion, accounting for about 75 per cent of the total energy density in the Universe. The simplest option for this "dark energy" corresponds to a cosmological constant, perhaps related to the quantum vacuum energy. Physically viable alternatives invoke either the presence of a scalar field with an evolving equation of state, or extensions of general relativity involving higher-order curvature terms or extra dimensions. Although they produce similar expansion rates, different models predict measurable differences in the growth rate of large-scale structure with cosmic time. A fingerprint of this growth is provided by coherent galaxy motions, which introduce a radial anisotropy in the clustering pattern reconstructed by galaxy redshift surveys. Here we report a measurement of this effect at a redshift of 0.8. Using a new survey of more than 10,000 faint galaxies, we measure the anisotropy parameter b = 0.70 +/- 0.26, which corresponds to a growth rate of structure at that time of f = 0.91 +/- 0.36. This is consistent with the standard cosmological-constant model with low matter density and flat geometry, although the error bars are still too large to distinguish among alternative origins for the accelerated expansion. This could be achieved with a further factor-of-ten increase in the sampled volume at similar redshift.
[158]  [pdf] - 9673
Spectroscopy of clusters in the ESO distant cluster survey (EDisCS).II. Redshifts, velocity dispersions, and substructure for clusters in the last 15 fields
Comments: 31 pages, 23 figures, accepted for publication in A&A
Submitted: 2008-02-01
AIMS. We present spectroscopic observations of galaxies in 15 survey fields as part of the ESO Distant Cluster Survey (EDisCS). We determine the redshifts and velocity dispersions of the galaxy clusters located in these fields, and we test for possible substructure in the clusters. METHODS. We obtained multi-object mask spectroscopy using the FORS2 instrument at the VLT. We reduced the data with particular attention to the sky subtraction. We implemented the method of Kelson for performing sky subtraction prior to any rebinning/interpolation of the data. From the measured galaxy redshifts, we determine cluster velocity dispersions using the biweight estimator and test for possible substructure in the clusters using the Dressler-Shectman test. RESULTS. The method of subtracting the sky prior to any rebinning/interpolation of the data delivers photon-noise-limited results, whereas the traditional method of subtracting the sky after the data have been rebinned/interpolated results in substantially larger noise for spectra from tilted slits. Redshifts for individual galaxies are presented and redshifts and velocity dispersions are presented for 21 galaxy clusters. For the 9 clusters with at least 20 spectroscopically confirmed members, we present the statistical significance of the presence of substructure obtained from the Dressler-Shectman test, and substructure is detected in two of the clusters. CONCLUSIONS. Together with data from our previous paper, spectroscopy and spectroscopic velocity dispersions are now available for 26 EDisCS clusters with redshifts in the range 0.40-0.96 and velocity dispersions in the range 166-1080 km/s.
[159]  [pdf] - 2302
The Dependence of Galaxy Formation on Cosmological Parameters: Can we distinguish the WMAP1 and WMAP3 Parameter Sets?
Comments: 18 pages,17 figures, version accepted by MNRAS, references and subsection on PVD added
Submitted: 2007-06-18, last modified: 2007-12-03
We combine N-body simulations of structure growth with physical modelling of galaxy evolution to investigate whether the shift in cosmological parameters between the 1-year and 3-year results from the Wilkinson Microwave Anisotropy Probe affects predictions for the galaxy population. Structure formation is significantly delayed in the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on the relevant scales. The decrease in dark matter clustering strength is, however, almost entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely altered. In both cosmologies several combinations of physical parameters can reproduce observed, low-redshift galaxy properties; the star formation, supernova feedback and AGN feedback efficiencies can be played off against each other to give similar results. Models which fit observed luminosity functions predict projected 2-point correlation functions which scatter by about 10-20 per cent on large scale and by larger factors on small scale, depending both on cosmology and on details of galaxy formation. Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but careful treatment of the systematics is needed. Given current modelling uncertainties, it is not easy to distinguish the WMAP1 and WMAP3 cosmologies on the basis of low-redshift galaxy properties. Model predictions diverge more dramatically at high redshift. Better observational data at z>2 will better constrain galaxy formation and perhaps also cosmological parameters.
[160]  [pdf] - 7273
Realistic simulations of gravitational lensing by galaxy clusters: extracting arc parameters from mock DUNE images
Comments: 17 pages, 15 figures submitted to A&A
Submitted: 2007-11-21
We present a newly developed code that allows simulations of optical observations of galaxy fields with a variety of instruments. The code incorporates gravitational lensing effects and is targetted at simulating lensing by galaxy clusters. Our goal is to create the tools required for comparing theoretical expectations with observations to obtain a better understanding of how observational noise affects lensing applications such as mass estimates, studies on the internal properties of galaxy clusters and arc statistics. Starting from a set of input parameters, characterizing both the instruments and the observational conditions, the simulator provides a virtual observation of a patch of the sky. It includes several sources of noise such as photon-noise, sky background, seeing, and instrumental noise. Ray-tracing through simulated mass distributions accounts for gravitational lensing. Source morphologies are realistically simulated based on shapelet decompositions of galaxy images retrieved from the GOODS-ACS archive. According to their morphological class, spectral-energy-distributions are assigned to the source galaxies in order to reproduce observations of each galaxy in arbitrary photometric bands. We illustrate our techniques showing virtual observations of a galaxy-cluster core as it would be observed with the space telescope DUNE, which was recently proposed to ESA within its "Cosmic vision" programme. (Abridged)
[161]  [pdf] - 2680
The VIMOS-VLT Deep Survey. The dependence of clustering on galaxy stellar mass at z~1
Comments: 13 pages, 10 figures, accepted in A&A
Submitted: 2007-06-29, last modified: 2007-11-12
Aims: We use the VVDS-Deep first-epoch data to measure the dependence of galaxy clustering on galaxy stellar mass, at z~0.85. Methods: We measure the projected correlation function wp(rp) for sub-samples with 0.5<z<1.2 covering different mass ranges between 10^9 and 10^11 Msun. We quantify in detail the observational selection biases using 40 mock catalogues built from the Millennium run and semi-analytic models. Results: Our simulations indicate that serious incompleteness in mass is present only for log(M/Msun)<9.5. In the mass range log(M/Msun)=[9.0-9.5], the photometric selection function of the VVDS misses 2/3rd of the galaxies. The sample is virtually 100% complete above 10^10 Msun. We present the first direct evidence for a clear dependence of clustering on the galaxy stellar mass at z~0.85. The clustering length increases from r0 ~ 2.76 h^-1 Mpc for galaxies with mass M>10^9 Msun to r0 ~ 4.28 h^-1 Mpc for galaxies more massive than 10^10.5 Msun. At the same time, the slope increases from ~ 1.67 to ~ 2.28. A comparison of the observed wp(rp) to local measurements by the SDSS shows that the evolution is faster for objects less massive than ~10^10.5 Msun. This is interpreted as a higher dependence on redshift of the linear bias b_L for the more massive objects. While for the most massive galaxies b_L decreases from 1.5+/-0.2 at z~0.85 to 1.33+/-0.03 at z~0.15, the less massive population maintains a virtually constant value b_L~1.3. This result is in agreement with a scenario in which more massive galaxies formed at high redshift in the highest peaks of the density field, while less massive objects form at later epochs from the more general population of dark-matter halos.
[162]  [pdf] - 6535
Ram pressure histories of cluster galaxies
Comments: MNRAS accepted
Submitted: 2007-10-30
Ram pressure stripping can remove significant amounts of gas from galaxies that orbit in clusters and massive groups, and thus has a large impact on the evolution of cluster galaxies. In this paper, we reconstruct the present-day distribution of ram-pressure, and the ram pressure histories of cluster galaxies. To this aim, we combine the Millennium Simulation and an associated semi-analytic model of galaxy evolution with analytic models for the gas distribution in clusters. We find that about one quarter of galaxies in massive clusters are subject to strong ram-pressures that are likely to cause an expedient loss of all gas. Strong ram-pressures occur predominantly in the inner core of the cluster, where both the gas density and the galaxy velocity are higher. Since their accretion onto a massive system, more than 64 per cent of galaxies that reside in a cluster today have experienced strong ram-pressures of $>10^{-11}$ dyn cm$^{-2}$ which most likely led to a substantial loss of the gas.
[163]  [pdf] - 88536
The recycling of gas and metals in galaxy formation: predictions of a dynamical feedback model
Comments: 13 pages, 10 figures. Accepted by MNRAS. The galaxy catalogues at z=0 can be downloaded from
Submitted: 2007-01-14, last modified: 2007-05-29
We present results of a new feedback scheme implemented in the Munich galaxy formation model. The new scheme includes a dynamical treatment of galactic winds powered by supernovae explosions and stellar winds in a cosmological context. We find that such a scheme is a good alternative to empirically-motivated recipes for feedback in galaxy formation. Model results are in good agreement with the observed luminosity functions and stellar mass function for galaxies in the local Universe. In particular, the new scheme predicts a number density of dwarfs that is lower than in previous models. This is a consequence of a new feature of the model, that allows an estimate of the amount of mass and metals that haloes can permanently deposit into the IGM. This loss of material leads to the suppression of star formation in small haloes and therefore to the decrease in the number density of dwarf galaxies. The model is able to reproduce the observed mass-stellar metallicity and luminosity-gas metallicity relationships. This demonstrates that our scheme provides a significant improvement in the treatment of the feedback in dwarf galaxies. Despite these successes, our model does not reproduce the observed bimodality in galaxy colours and predicts a larger number of bright galaxies than observed. Finally, we investigate the efficiency of metal injection in winds and in the intergalactic medium. We find that galaxies that reside in haloes with M_vir < 10^12 M_sun/h may deposit most of their metal mass into the intergalactic medium, while groups and clusters at z=0 have lost at most a few percent of their metals before the bulk of the halo mass was accreted.
[164]  [pdf] - 80755
Modelling galaxy clustering in a high resolution simulation of structure formation
Comments: 12 pages, 11 figures, published version
Submitted: 2006-03-20, last modified: 2007-01-24
We use the Millennium Simulation, a 10 billion particle simulation of the growth of cosmic structure, to construct a new model of galaxy clustering. We adopt a methodology that falls midway between the traditional semi-analytic approach and the halo occupation distribution(HOD) approach. In our model, we adopt the positions and velocities of the galaxies that are predicted by following the orbits and merging histories of the substructures in the simulation. Rather than using star formation and feedback `recipes' to specify the physical properties of the galaxies, we adopt parametrized functions to relate these properties to the quantity M_{infall}, defined as the mass of the halo at the epoch when the galaxy was last the central dominant object in its own halo. We test whether these parametrized relations allow us to recover the basic statistical properties of galaxies in the semi-analytic catalogues, including the luminosity function, the stellar mass function and the shape and amplitude of the two-point correlation function evaluated in different stellar mass and luminosity ranges. We then use our model to interpret recent measurements of these quantities from Sloan Digital Sky Survey data. We derive relations between the luminosities and the stellar masses of galaxies in the local Universe and their host halo masses. Our results are in excellent agreement with recent determinations of these relations by Mandelbaum et al using galaxy-galaxy weak lensing measurements from the SDSS.
[165]  [pdf] - 88811
Modelling and interpreting the dependence of clustering on the spectral energy distributions of galaxies
Comments: 14 pages, 11 figures, submitted to MNRAS
Submitted: 2007-01-24
We extend our previous physically-based halo occupation distribution models to include the dependence of clustering on the spectral energy distributions of galaxies. The high resolution Millennium Simulation is used to specify the positions and the velocities of the model galaxies. The stellar mass of a galaxy is assumed to depend only on M_{infall}, the halo mass when the galaxy was last the central dominant object of its halo. Star formation histories are parametrized using two additional quantities that are measured from the simulation for each galaxy: its formation time (t_{form}), and the time when it first becomes a satellite (t_{infall}). Central galaxies begin forming stars at time t_{form} with an exponential time scale tau_c. If the galaxy becomes a satellite, its star formation declines thereafter with a new time scale tau_s. We compute 4000 \AA break strengths for our model galaxies using stellar population synthesis models. By fitting these models to the observed abundances and projected correlations of galaxies as a function of break strength in the Sloan Digital Sky Survey, we constrain tau_c and tau_s as functions of galaxy stellar mass. We find that central galaxies with large stellar masses have ceased forming stars. At low stellar masses, central galaxies display a wide range of different star formation histories, with a significant fraction experiencing recent starbursts. Satellite galaxies of all masses have declining star formation rates, with similar e-folding times, tau_s ~ 2.5Gyr. One consequence of this long e-folding time is that the colour-density relation is predicted to flatten at redshifts > 1.5, because star formation in the majority of satellites has not yet declined by a significant factor. This is consistent with recent observational results from the DEEP and VVDS surveys.
[166]  [pdf] - 81483
Substructure in lensing clusters and simulations
Comments: 15 pages, 5 figures, 2 tables, in press MNRAS
Submitted: 2006-04-19, last modified: 2006-12-19
We present high-resolution mass reconstructions for five massive cluster-lenses spanning a redshift range from $z = 0.18$--0.57 utilising archival {\it Hubble Space Telescope} ({\it HST}) data and applying galaxy-galaxy lensing techniques. These detailed mass models were obtained from the observations by combining constraints from the strong and weak lensing regimes. We ascribe local weak distortions in the shear maps to perturbations induced by the presence of galaxy haloes around individual bright early-type cluster member galaxies. This technique constrains the mass enclosed within an aperture for these subhaloes. We are sensitive to a specific mass range for these subhaloes, $10^{11}$ -- $10^{12.5} \msun$, which we associate with galaxy-scale subhaloes. Adopting a parametric model for the subhaloes, we also derive their velocity dispersion function and the aperture radius function. The mass spectrum of substructure in the inner regions of the observed clusters is directly compared with that in simulated clusters extracted from the {\it Millennium Simulation}. The massfunction, aperture radii and velocity dispersion function are compared in detail. Overall, we find good agreement between the distribution of substructure properties retrieved using the lensing analysis and those obtained from the simulation (truncated abstract).
[167]  [pdf] - 82959
The hierarchical formation of the brightest cluster galaxies
Comments: 13 pages, 17 figures, MNRAS accepted version
Submitted: 2006-06-21, last modified: 2006-11-14
We use semi-analytic techniques to study the formation and evolution of brightest cluster galaxies (BCGs). We show the extreme hierarchical nature of these objects and discuss the limits of simple ways to capture their evolution. In a model where cooling flows are suppressed at late times by AGN activity, the stars of BCGs are formed very early (50 per cent at z~5, 80 per cent at z~3) and in many small galaxies. The high star formation rates in these high-z progenitors are fuelled by rapid cooling, not by merger-triggered starbursts. We find that model BCGs assemble surprisingly late: half their final mass is typically locked-up in a single galaxy after z~0.5. Because most of the galaxies accreted onto BCGs have little gas content and red colours, late mergers do not change the apparent age of BCGs. It is this accumulation of a large number of old stellar populations -- driven mainly by the merging history of the dark matter halo itself -- that yields the observed homogeneity of BCG properties. In the second part of the paper, we discuss the evolution of BCGs to high redshifts, from both observational and theoretical viewpoints. We show that our model BCGs are in qualitative agreement with high-z observations. We discuss the hierarchical link between high-z BCGs and their local counter-parts. We show that high-z BCGs belong to the same population as the massive end of local BCG progenitors, although they are not in general the same galaxies. Similarly, high-z BCGs end-up as massive galaxies in the local Universe, although only a fraction of them are actually BCGs of massive clusters.
[168]  [pdf] - 86736
How `heredity' and `environment' shape galaxy properties
Comments: 8 pages, 4 figure, to appear in ASP Conference Series (Proceedings of the "Cosmic Frontiers" conference held in Durham, 31st July - 4th August 2006)
Submitted: 2006-11-10
In this review, I give a brief summary of galaxy evolution processes in hierarchical cosmologies and of their relative importance at different masses, times, and environments. I remind the reader of the processes that are commonly included in modern semi-analytic models of galaxy formation, and I comment on recent results and open issues.
[169]  [pdf] - 85797
The build-up of the colour-magnitude relation in galaxy clusters since z~0.8
Comments: 15 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2006-10-12, last modified: 2006-10-18
Using galaxy clusters from the ESO Distant Cluster Survey, we study how the distribution of galaxies along the colour-magnitude relation has evolved since z~0.8. While red-sequence galaxies in all these clusters are well described by an old, passively evolving population, we confirm our previous finding of a significant evolution in their luminosity distribution as a function of redshift. When compared to galaxy clusters in the local Universe, the high redshift EDisCS clusters exhibit a significant "deficit" of faint red galaxies. Combining clusters in three different redshift bins, and defining as `faint' all galaxies in the range 0.4 > L/L* > 0.1, we find a clear decrease in the luminous-to-faint ratio of red galaxies from z~0.8 to z~0.4. The amount of such a decrease appears to be in qualitative agreement with predictions of a model where the blue bright galaxies that populate the colour-magnitude diagram of high redshift clusters, have their star formation suppressed by the hostile cluster environment. Although model results need to be interpreted with caution, our findings clearly indicate that the red-sequence population of high-redshift clusters does not contain all progenitors of nearby red-sequence cluster galaxies. A significant fraction of these must have moved onto the red-sequence below z~0.8.
[170]  [pdf] - 74935
The many lives of active galactic nuclei: cooling flows, black holes and the luminosities and colours of galaxies
Comments: 21 pages, 12 figures, accepted MNRAS, Millennium Run semi-analytic galaxy catalogues are available for download at (~9 million galaxies brighter than -17.4 in r), mini Millennium Run catalogues also available for testing purposes
Submitted: 2005-08-02, last modified: 2006-06-14
We simulate the growth of galaxies and their central supermassive black holes by implementing a suite of semi-analytic models on the output of the Millennium Run, a very large simulation of the concordance LCDM cosmogony. Our procedures follow the detailed assembly history of each object and are able to track the evolution of all galaxies more massive than the Small Magellanic Cloud throughout a volume comparable to that of large modern redshift surveys. In this first paper we supplement previous treatments of the growth and activity of central black holes with a new model for `radio' feedback from those AGN that lie at the centre of a quasistatic X-ray emitting atmosphere in a galaxy group or cluster. We show that for energetically and observationally plausible parameters such a model can simultaneously explain: (i) the low observed mass drop-out rate in cooling flows; (ii) the exponential cut-off at the bright end of the galaxy luminosity function; and (iii) the fact that the most massive galaxies tend to be bulge-dominated systems in clusters and to contain systematically older stars than lower mass galaxies. This success occurs because static hot atmospheres form only in the most massive structures, and radio feedback (in contrast, for example, to supernova or starburst feedback) can suppress further cooling and star formation without itself requiring star formation. We discuss possible physical models which might explain the accretion rate scalings required for our phenomenological `radio mode' model to be successful.
[171]  [pdf] - 79641
Erratum - The many lives of AGN: cooling flows, black holes and the luminosities and colours of galaxies
Comments: 1 page, accepted to MNRAS, erratum for MNRAS, 365, 11
Submitted: 2006-02-03
In Figure 6 we inadvertently labeled the proxy circular velocity as the virial velocity of the dark matter halo instead of what is actually plotted, the maximum circular velocity of the dark matter halo. The maximum halo circular velocity is a much better estimate of the disk V_c than is V_vir. This confusion influenced the discussion of the Tully-Fisher relation in our paper. In fact, Figure 6 demonstrates that it is possible to simultaneously reproduce both the local Tully-Fisher relation and luminosity function using semi-analytic techniques applied to the standard LCDM cosmology, thus contradicting previous studies of this issue and our own discussion in Section 3.6.
[172]  [pdf] - 78591
The evolution of the star formation activity in galaxies and its dependence on environment
Comments: 28 pages, 19 figures, ApJ in press
Submitted: 2005-12-15
We study how the proportion of star-forming galaxies evolves between z=0.8 and z=0 as a function of galaxy environment, using the [OII] line in emission as a signature of ongoing star formation. Our high-z dataset comprises 16 clusters, 10 groups and another 250 galaxies in poorer groups and the field at z=0.4-0.8 from the ESO Distant Cluster Survey, plus another 9 massive clusters at similar redshifts. As a local comparison, we use samples of galaxy systems selected from the Sloan Digital Sky Survey at 0.04< z < 0.08. At high-z most systems follow a broad anticorrelation between the fraction of star-forming galaxies and the system velocity dispersion. At face value, this suggests that at z=0.4-0.8 the mass of the system largely determines the proportion of galaxies with ongoing star formation. At these redshifts the strength of star formation (as measured by the [OII] equivalent width) in star-forming galaxies is also found to vary systematically with environment. Sloan clusters have much lower fractions of star-forming galaxies than clusters at z=0.4-0.8 and, in contrast with the distant clusters, show a plateau for velocity dispersions $ \ge 550 km s^-1$, where the fraction of galaxies with [OII] emission does not vary systematically with velocity dispersion. We quantify the evolution of the proportion of star-forming galaxies as a function of the system velocity dispersion and find it is strongest in intermediate-mass systems (sigma ~ 500-600 km s^-1 at z=0). To understand the origin of the observed trends, we use the Press-Schechter formalism and the Millennium Simulation and show that galaxy star formation histories may be closely related to the growth history of clusters and groups. We propose a scheme that is able to account for the observed relations between the star-forming fraction and \sigma [abridged].
[173]  [pdf] - 78104
Weak lensing mass reconstructions of the ESO Distant Cluster Survey
Comments: 35 pages, 24 figures, accepted to A&A, a version with better figure resolution can be found at
Submitted: 2005-11-28
We present weak lensing mass reconstructions for the 20 high-redshift clusters i n the ESO Distant Cluster Survey. The weak lensing analysis was performed on deep, 3-color optical images taken with VLT/FORS2, using a composite galaxy catalog with separate shape estimators measured in each passband. We find that the EDisCS sample is composed primarily of clusters that are less massive than t hose in current X-ray selected samples at similar redshifts, but that all of the fields are likely to contain massive clusters rather than superpositions of low mass groups. We find that 7 of the 20 fields have additional massive structures which are not associated with the clusters and which can affect the weak lensing mass determination. We compare the mass measurements of the remaining 13 clusters with luminosity measurements from cluster galaxies selected using photometric redshifts and find evidence of a dependence of the cluster mass-to-light ratio with redshift. Finally we determine the noise level in the shear measurements for the fields as a function of exposure time and seeing and demonstrate that future ground-based surveys which plan to perform deep optical imaging for use in weak lensing measurements must achieve point-spread functions smaller than a median of 0.6" FWHM.
[174]  [pdf] - 76306
The formation history of elliptical galaxies
Comments: 12 pages, 11 figures, minor revisions, MNRAS accepted
Submitted: 2005-09-23, last modified: 2005-11-21
We take advantage of the largest high-resolution simulation of cosmic structure growth ever carried out -- the Millennium Simulation of the concordance LambdaCDM cosmogony -- to study how the star formation histories, ages and metallicities of elliptical galaxies depend on environment and on stellar mass. We concentrate on a galaxy formation model which is tuned to fit the joint luminosity/colour/morphology distribution of low redshift galaxies. Massive ellipticals in this model have higher metal abundances, older luminosity-weighted ages, shorter star formation timescales, but lower assembly redshifts than less massive systems. Within clusters the typical masses, ages and metal abundances of ellipticals are predicted to decrease, on average, with increasing distance from the cluster centre. We also quantify the effective number of progenitors of ellipticals as a function of present stellar mass, finding typical numbers below 2 for M* < 10^{11} Msun, rising to about 5 for the most massive systems. These findings are consistent with recent observational results that suggest ``down-sizing'' or ``anti-hierarchical'' behaviour for the star formation history of the elliptical galaxy population, despite the fact that our model includes all the standard elements of hierarchical galaxy formation and is implemented on the standard, LambdaCDM cosmogony.
[175]  [pdf] - 76901
Gas Infall and Stochastic Star Formation in Galaxies in the Local Universe
Comments: Submitted to MNRAS. Note that galaxy physical parameters for DR4 (masses, metallicities, SFR, AGN catalog) are now available at
Submitted: 2005-10-13, last modified: 2005-10-17
We study the recent star formation histories of local galaxies by analyzing the scatter in their colours and spectral properties. We present evidence that the distribution of star formation histories changes qualitatively above a characteristic stellar surface mass density of 3x10^8 M_sol kpc^-2, corresponding to the transition between disk-dominated and bulge-dominated galaxies. When we average over subpopulations of galaxies with densities below this value, we find that subpopulations of all masses form their stars at the same average rate per unit stellar mass. However, the scatter in galaxy colours and emission line strengths is larger for more compact galaxies of a given mass. This suggests that star formation occurs in shorter, higher amplitude events in galaxies with smaller sizes. Above the characteristic density, galaxy growth through star formation shuts down. We propose that star formation events are triggered when cold gas is accreted onto a galaxy. We have used the new high resolution Millennium Simulation of structure formation in a concordance LCDM Universe to quantify the incidence of these accretion events and we show that the observational data are well fit by a model in which the consumption time of accreted gas decreases with the surface density of the galaxy as t_cons \propto mu*^-1. The dark matter halos hosting massive galaxies with high densities are also expected to grow through accretion, but the observations indicate that star formation in bulge- dominated galaxies is no longer coupled to the hierarchical build-up of their surrounding dark matter halos.
[176]  [pdf] - 75240
EDisCS -- the ESO Distant Cluster Survey -- Sample Definition and Optical Photometry
Comments: 21 pages, 6 figures, A&A in press. A version with full resolution figures is available at
Submitted: 2005-08-16
We present the ESO Distant Cluster Survey (EDisCS) a survey of 20 fields containing distant galaxy clusters with redshifts ranging from 0.4 to almost 1.0. Candidate clusters were chosen from among the brightest objects identified in the Las Campanas Distant Cluster Survey, half with estimated redshift z_est~0.5 and half with z_est~0.8. They were confirmed by identifying red sequences in moderately deep two colour data from VLT/FORS2. For confirmed candidates we have assembled deep three-band optical photometry using VLT/FORS2, deep near-infrared photometry in one or two bands using NTT/SOFI, deep optical spectroscopy using VLT/FORS2, wide field imaging in two or three bands using the ESO Wide Field Imager, and HST/ACS mosaic images for 10 of the most distant clusters. This first paper presents our sample and the VLT photometry we have obtained. We present images, colour-magnitude diagrams and richness estimates for our clusters, as well as giving redshifts and positions for the brightest cluster members. Subsequent papers will present our infrared photometry, spectroscopy, HST and wide-field imaging, as well as a wealth of further analysis and science results. Our reduced data become publicly available as these papers are accepted.
[177]  [pdf] - 1456389
The Capodimonte Deep Field I - Presentation of the survey and first follow-up studies
Comments: 24 pages, accepted for pubblication on A&A
Submitted: 2004-08-12
We present the Capodimonte Deep Field (OACDF), a deep field covering an area of 0.5 square degrees in the B, V, R optical bands plus six medium-band filters in the wavelength range 773--913 nm. The field reaches the following limiting AB magnitudes: B~25.3, V~24.8 and R~25.1 and contains about 50000 extended sources in the AB magnitude range 18 < R < 25 . Hence, it is intermediate between deep pencil beam surveys and very wide but shallow surveys. The main scientific goal of the OACDF is the identification and characterization of early-type field galaxies at different look-back times in order to study different scenarios of galaxy formation. Parallel goals include the search for groups and clusters of galaxies and the search for rare and peculiar objects (gravitational lenses, QSOs, halo White Dwarfs). In this paper we describe the OACDF data reduction, the methods adopted for the extraction of the photometric catalogs, the photometric calibration and the quality assessment of the catalogs by means of galaxy number counts, spectroscopic and photometric redshifts and star colors. We also present the first results of the search for galaxy overdensities. The depth of the OACDF and its relatively large spatial coverage with respect to pencil beam surveys make it a good tool for further studies of galaxy formation and evolution in the redshift range 0--1, as well as for stellar studies.
[178]  [pdf] - 66562
Spectroscopy of clusters in the ESO Distant Cluster Survey (EDisCS)
Comments: 18 pages, 27 figures, accepted for publication in A&A, Table 4 is available ahead of journal publication by downloading the source files for this astro-ph submission or from first author on request (
Submitted: 2004-08-04
We present spectroscopic observations of galaxies in 4 clusters at z = 0.7-0.8 and in one cluster at z~0.5 obtained with the FORS2 spectrograph on the VLT as part of the ESO Distant Cluster Survey (EDisCS), a photometric and spectroscopic survey of 20 intermediate to high redshift clusters. We describe our target selection, mask design, observation and data reduction procedures, using these first 5 clusters to demonstrate how our strategies maximise the number of cluster members for which we obtain spectroscopy. We present catalogues containing positions, I-band magnitudes and spectroscopic redshifts for galaxies in the fields of our 5 clusters. These contain 236 cluster members, with the number of members per cluster ranging from 30 to 67. Our spectroscopic success rate, i.e. the fraction of spectroscopic targets which are cluster members, averages 50% and ranges from 30% to 75%. We use a robust biweight estimator to measure cluster velocity dispersions from our spectroscopic redshift samples. We also make a first assessment of substructure within our clusters. The velocity dispersions range from 400 to 1100 km s-1. Some of the redshift distributions are significantly non-Gaussian and we find evidence for significant substructure in two clusters, one at z~0.79 and the other at z~0.54. Both have velocity dispersions exceeding 1000 km s-1 but are clearly not fully virialised; their velocity dispersions may thus be a poor indicator of their masses. The properties of these first 5 EDisCS clusters span a wide range in redshift, velocity dispersion, richness and substructure, but are representative of the sample as a whole. Spectroscopy for the full dataset will allow a comprehensive study of galaxy evolution as a function of cluster environment and redshift.
[179]  [pdf] - 65611
The build-up of the red-sequence in galaxy clusters since z~0.8
Comments: 4 pages; 3 figures; ApJL accepted
Submitted: 2004-06-21
We study the rest-frame (U-V) color-magnitude relation in 4 clusters at redshifts 0.7-0.8, drawn from the ESO Distant Cluster Survey. We confirm that red-sequence galaxies in these clusters can be described as an old, passively-evolving population and we demonstrate, by comparison with the Coma cluster, that there has been significant evolution in the stellar mass distribution of red-sequence galaxies since z~0.75. The EDisCS clusters exhibit a deficiency of low luminosity passive red galaxies. Defining as `faint' all galaxies in the passive evolution corrected range 0.4>~ L/L*>~0.1, the luminous-to-faint ratio of red-sequence galaxies varies from 0.34+/-0.06 for the Coma cluster to 0.81+/-0.18 for the high redshift clusters. These results exclude a synchronous formation of all red-sequence galaxies and suggest that a large fraction of the faint red galaxies in current clusters moved on to the red-sequence relatively recently. Their star formation activity presumably came to an end at z<~0.8.
[180]  [pdf] - 64530
Galaxies and subhaloes in LCDM galaxy clusters
Comments: 6 pages, 2 figures, accepted for publication in Monthly Notices Letters. Final version, figures updated, minor changes in text
Submitted: 2004-05-03, last modified: 2004-06-04
We combine 10 high resolution resimulations of cluster--sized dark haloes with semi--analytic galaxy formation modelling in order to compare the number density and velocity dispersion profiles of cluster galaxies with those of dark matter substructures (subhaloes). While the radial distribution of galaxies follows closely that of the dark matter, the distribution of dark matter subhaloes is much less centrally concentrated. The velocity dispersion profiles of galaxies are also very similar to those of the dark matter, while those for subhaloes are biased high, particularly in the inner regions of clusters. We explain how these differences, already clearly visible in earlier work, are a consequence of the formation of galaxies at the centres of dark matter haloes. Galaxies and subhaloes represent different populations and are not directly comparable. Evolution produces a complex and strongly position--dependent relation between galaxies and the subhaloes in which they reside. This relation can be properly modelled only by appropriate physical representation of the galaxy formation process.
[181]  [pdf] - 63992
Observing the build-up of the colour-magnitude relation at redshift ~0.8
Comments: 5 pages, 2 figures, to appear in Proc. of IAU Colloq. 195: "Outskirts of Galaxy Clusters: Intense Life in the Suburbs" -- minor typos corrected
Submitted: 2004-04-05, last modified: 2004-04-06
We analyse the rest-frame (U-V) colour-magnitude relation for 2 clusters at redshift 0.7 and 0.8, drawn from the ESO Distant Cluster Survey. By comparing with the population of red galaxies in the Coma cluster, we show that the high redshift clusters exhibit a deficit of passive faint red galaxies. Our results show that the red-sequence population cannot be explained in terms of a monolithic and synchronous formation scenario. A large fraction of faint passive galaxies in clusters today has moved onto the red sequence relatively recently as a consequence of the fact that their star formation activity has come to an end at z<0.8.
[182]  [pdf] - 62816
Effects of asphericity and substructure on the determination of cluster mass with weak gravitational lensing
Comments: 10 pages, 10 figures, accepted by MNRAS
Submitted: 2004-02-12
Weak gravitational lensing can be used to directly measure the mass along a line-of-sight without any dependence on the dynamical state of the mass, and thus can be used to measure the masses of clusters even if they are not relaxed. One common technique used to measure cluster masses is fitting azimuthally-averaged gravitational shear profiles with a spherical mass model. In this paper we quantify how asphericity and projected substructure in clusters can affect the virial mass and concentration measured with this technique by simulating weak lensing observations on 30 independent lines-of-sights through each of four high-resolution N-body cluster simulations. We find that the variations in the measured virial mass and concentration are of a size similar to the error expected in ideal weak lensing observations and are correlated, but that the virial mass and concentration of the mean shear profile agree well with that measured in three dimensional models of the clusters. The dominant effect causing the variations is the proximity of the line-of-sight to the major axis of the 3-D cluster mass distribution, with projected substructure only causing minor perturbations in the measured concentration. Finally we find that the best-fit "universal" CDM models used to fit the shear profiles over-predict the surface density of the clusters due to the cluster mass density falling off faster than the r^{-3} model assumption.
[183]  [pdf] - 59931
Chemical enrichment of the intra-cluster and intergalactic medium in a hierarchical galaxy formation model
Comments: Accepted for publication in MNRAS. 17 pages, 13 figures. Minor changes to submitted version
Submitted: 2003-10-09, last modified: 2004-01-08
(Abridged) We use a combination of high resolution N-body simulations and semi-analytic techniques to follow the formation, the evolution and the chemical enrichment of cluster galaxies in a Lambda-CDM Universe. We model the transport of metals between the stars, the cold gas in galaxies, the hot gas in dark matter haloes, and the intergalactic gas outside virialized haloes. We have compared three different feedback schemes. The `retention' model assumes that material reheated by supernova explosions is able to leave the galaxy, but not the dark matter halo. The `ejection' model assumes that this material leaves the halo and is then re-incorporated when structure collapses on larger scales. The `wind' model uses prescriptions that are motivated by observations of local starburst galaxies. We require that our models reproduce the cluster LF from the 2dF survey, the relations between stellar mass, gas mass and metallicity inferred from new SDSS data, and the observed amount of metals in the ICM. With suitable adjustment of the free parameters in the model, a reasonable fit to the observational results at redshift zero can be obtained for all three feedback schemes. All three predict that the chemical enrichment of the ICM occurs at high redshift with 60-80 per cent of the metals currently in the ICM ejected at redshifts larger than 1. Massive galaxies are important contributors to the chemical pollution. The observed decline in baryon fraction from rich clusters to galaxy groups is reproduced only in an `extreme' ejection scheme, where material ejected from dark matter haloes is re-incorporated on a timescale comparable to the age of the Universe. We explore how the metal abundance in the intergalactic medium as a function of redshift can constraint how and when galaxies ejected their metals.
[184]  [pdf] - 57286
Substructures in Cold Dark Matter Haloes
Comments: 13 pages, 15 figure. Accepted to MNRAS
Submitted: 2003-06-10, last modified: 2003-11-05
We analyse the properties of substructures within dark matter halos (subhalos) using a set of high-resolution numerical simulations of the formation of structure in a Lambda-CDM Universe. Our simulation set includes 11 high-resolution simulations of massive clusters as well as a region of mean density, allowing us to study the spatial and mass distribution of substructures down to a mass resolution limit of 10^9 h^(-1)Mo. We also investigate how the properties of substructures vary as a function of the mass of the `parent' halo in which they are located. We find that the substructure mass function depends at most weakly on the mass of the parent halo and is well described by a power-law. The radial number density profiles of substructures are steeper in low mass halos than in high mass halos. More massive substructures tend to avoid the centres of halos and are preferentially located in the external regions of their parent halos. We also study the mass accretion and merging histories of substructures, which we find to be largely independent of environment. We find that a significant fraction of the substructures residing in clusters at the present day were accreted at redshifts z < 1. This implies that a significant fraction of present-day `passive' cluster galaxies should have been still outside the cluster progenitor and more active at z~1.
[185]  [pdf] - 59140
Evolution of galaxies in clusters: semi-analytic models and observations
Comments: 2 pages, 2 figures, to appear in "Multiwavelength Cosmology", Mykonos Conference Proceedings
Submitted: 2003-09-11
We present two different projects. The first one is related to the development of a semi-analytic model to follow the formation, the evolution and chemical enrichment of cluster galaxies in a hierarchical dark matter model. The second project concerns the comparison between high resolution N-body simulations of clusters of galaxies with the EDisCS observational dataset.
[186]  [pdf] - 45306
The puzzle of the polar structure in NGC 4650A
Comments: 26 pages, 10 figures, accepted for publication to the Astronomical Journal
Submitted: 2001-10-10
This work presents new surface photometry and two-dimensional modeling of the light distribution of the Polar Ring Galaxy NGC 4650A, based on near-infrared (NIR) observations and high resolution optical imaging acquired during the Hubble Heritage program. The NIR and optical integrated colors of the S0 and the polar ring, and their scale parameters, are compared with those for standard galaxy morphological types. The polar structure appears to be a disk of a very young age, while the colors and light distribution of the host galaxy do not resemble that of a typical early-type system. We compare these observational results with the predictions from different formation scenarios for polar ring galaxies. The peculiarities of the central S0 galaxy, the polar disk structure and stellar population ages suggest that the polar ring galaxy NGC 4650A may be the result of a dissipative merger event, rather than of an accretion process.