Search

We present the integral-field kinematics and stellar population properties of 64 galaxies (61 are Early-Type galaxies, ETGs) with Counter-Rotating stellar Disks (CRD) selected from about 4000 galaxies in the MaNGA survey, based on evidence of counter-rotation or two velocity dispersion peaks in the kinematic maps. For 17 CRDs, the counter-rotating components can also be separated spectroscopically. The frequency of CRDs in MaNGA is, 26 pages, 22 figures, 1 table. Accepted for publication on December 16, 2021; published as MNRAS 511, 139 (2022)

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5�� allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar), 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

This work is based on deep multi-band (g, r, i) data from the Fornax Deep Survey with VST. We analyse the surface brightness profiles of the 19 bright ETGs inside the virial radius of the Fornax cluster. The main aim of this work is to identify signatures of accretion onto galaxies by studying the presence of outer stellar halos, and understand their nature and occurrence. Our analysis also provides a new and accurate estimate of the intra-cluster light inside the virial radius of Fornax. We performed multi-component fits to the azimuthally averaged surface brightness profiles available for all sample galaxies. This allows to quantify the relative weight of all components in the galaxy structure that contribute to the total light. In addition, we derived the average g-i colours in each component identified by the fit, as well as the azimuthally averaged g-i colour profiles, to correlate them with the stellar mass of each galaxy and the location inside the cluster. We find that in the most massive and reddest ETGs the fraction of light in, probably accreted, halos is much larger than in the other galaxies. Less-massive galaxies have an accreted mass fraction lower than 30%, bluer colours and reside in the low-density regions of the cluster. Inside the virial radius of the cluster, the total luminosity of the intra-cluster light, compared with the total luminosity of all cluster members, is about 34%. Inside the Fornax cluster there is a clear correlation between the amount of accreted material in the stellar halos of galaxies and the density of the environment in which those galaxies reside. By comparing this quantity with theoretical predictions and previous observational estimates, there is a clear indication that the driving factor for the accretion process is the total stellar mass of the galaxy, in agreement with the hierarchical accretion scenario., 18 pages, 10 figures. Accepted for publication in A&A

In this work we analyse the structural and photometric properties of 21 barred simulated galaxies from the Auriga Project. These consist of Milky Way-mass magnetohydrodynamical simulations in a Λ cold dark matter (ΛCDM) cosmological context. In order to compare with observations, we generate synthetic SDSS-like broad-band images from the numerical data at z = 0 with different inclinations (from face-on to edge-on). Ellipse fits are used to determine the bar lengths, and 2D bulge/disc/bar decompositions with galfit are also performed, modelling the bar component with the modified Ferrer profile. We find a wide range of bar sizes and luminosities in the sample, and their structural parameters are in good agreement with the observations. All bulges present low Sérsic indexes, and are classified as pseudobulges. In regard to the discs, the same breaks in the surface brightness profiles observed in real galaxies are found, and the radii at which these take place are in agreement with the observations. Also, from edge-on unsharp-masked images at z = 0, boxy or peanut-shaped (B/P) structures are clearly identified in the inner part of four bars, and also two more bars are found in buckling phase. The sizes of the B/P match fairly well with those obtained from observations. We thus conclude that the observed photometric and structural properties of galaxies with bars, which are the main drivers of secular evolution, can be developed in present state-of-the-art ΛCDM cosmological simulations.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh–Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work, we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh–Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh–Mstar relation requires a mean radiative efficiency ε ≳ 0.15, in line with theoretical expectations for accretion on to spinning black holes. However, matching the ‘raw’ observed relation for inactive black holes requires ε ∼ 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve ε ∼ 0.12–0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized., FS acknowledges partial support from a Leverhulme Trust Research Fellowship. RC acknowledges financial support from CONICYT Doctorado Nacional N° 21161487 and CONICYT PIA ACT172033. DMA thanks the Science and Technology Facilities Council (STFC) for support from grant no. ST/L00075X/1. ID is supported by the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 788679. MM acknowledges support from the Beatriu de Pinos fellowship (2017-BP-00114).

We analyse the sizes, colour gradients, and resolved stellar mass distributions for 36 massive and passive galaxies in the cluster XMMUJ2235-2557 at z=1.39 using optical and near-infrared Hubble Space Telescope imaging. We derive light-weighted S\'ersic fits in five HST bands ($i_{775},z_{850},Y_{105},J_{125},H_{160}$), and find that the size decreases by ~20% going from $i_{775}$ to $H_{160}$ band, consistent with recent studies. We then generate spatially resolved stellar mass maps using an empirical relationship between $M_{*}/L_{H_{160}}$ and $(z_{850}-H_{160})$ and use these to derive mass-weighted S\'ersic fits: the mass-weighted sizes are ~41% smaller than their rest-frame $r$-band counterparts compared with an average of ~12% at z~0. We attribute this evolution to the evolution in the $M_{*}/L_{H_{160}}$ and colour gradient. Indeed, as expected, the ratio of mass-weighted to light-weighted size is correlated with the $M_{*}/L$ gradient, but is also mildly correlated with the mass surface density and mass-weighted size. The colour gradients $(\nabla_{z-H})$ are mostly negative, with a median value of $\sim0.45$ mag dex$^{-1}$, twice the local value. The evolution is caused by an evolution in age gradients along the semi-major axis ($a$), with $\nabla_{age} = d \log(age) / d \log(a)$ $\sim-0.33$, while the survival of weaker colour gradients in old, local galaxies implies that metallicity gradients are also required, with $\nabla_{Z} = d \log(Z) / d \log(a)$ $\sim-0.2$. This is consistent with recent observational evidence for the inside-out growth of passive galaxies at high redshift, and favours a gradual mass growth mechanism, such as minor mergers., Comment: 30 pages, 25 figures, 4 tables, accepted for publication in MNRAS

Context. The upcoming new generation of optical spectrographs on four-meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. Aims. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. Methods. We produced spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assessed our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise ratios. We then analyzed such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. Results. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ∼0.1 Gyr for stellar populations older than ∼1.5 Gyr, pushing up to an age of ∼1 Gyr for stellar populations older than ∼5 Gyr. Conclusions. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.© ESO 2019., L. C. wish to acknowledge financial support from Premiale 2015 MITiC, program 1.05.06.10 and Comunidad de Madrid under Atraccion de Talento grant 2018-T2/TIC-11612. A. I., S. Z., M. L., A. G., A. M., and C. T. acknowledge the financial support from the INAF PRIN-SKA 2017 program 1.05.01.88.04.ESKAPE-HI. I. L. acknowledges financial support from INAF-WEAVE funds, program 1.05.03.04.05. R. G. B. acknowledges support from the Spanish Ministry of Economy and Competitiveness through grant 205 AYA2016-77846-P and the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award to 2017 the Instituto de Astrofisica de Andalucia (SEV-2017-0709). A. F. M. has received financial support through the Postdoctoral Junior Leader Fellowship Programme from >La Caixa> Banking Foundation (LCF/BQ/LI18/11630007). B. V. acknowledges the financial support from INAF Main Stream 2018 (P. I.: B. Vulcani).

Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range $10^9-10^{11.5}~\mathrm{M}_{\odot}$. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile ($C_{82}$). S/T is correlated with stellar mass and star-formation activity, while $C_{82}$ correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and $C_{82}$. Low and high mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies ($M_{\star}\approx10^{10}-10^{10.5}~\mathrm{M}_{\odot}$) are dominated by ordered rotation. Whereas higher-mass galaxies are typical spheroids with high concentrations, lower-mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star-formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and $C_{82}$) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier-forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in-situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with $M_{\star}>10^{11}~\mathrm{M}_{\odot}$. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent., 20 pages with 14 figures (+ appendix). Accepted by MNRAS

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4, Comment: Submitted to A&A

The surface-brightness profiles of galaxies I(R) and the density profiles of dark-matter halos rho(r) are well represented by the same analytic function, named after either S\'ersic, I~exp[-(R/R*)^(1/m)], or Einasto, rho~[exp[-(r/r*)^alpha], where R* and r* are characteristic radii. Systems with high S\'ersic index m (or low Einasto index alpha) have steep central profiles and shallow outer profiles, while systems with low m (or high alpha) have shallow central profiles and steep profiles in the outskirts. We present the results of idealized numerical experiments which suggest that the origin of these profiles can be traced back to the initial density fluctuation field: high-alpha (low-m) systems form in smooth regions via few mergers, while low-alpha (high-m) systems form in clumpy regions via several mergers., Comment: 3 pages, 2 figures, Proceedings of IAU Symposium 321, "Formation and evolution of galaxy outskirts", Eds. A. Gil de Paz, J. C. Lee & J. H. Knapen, Cambridge University Press, Cambridge

We demonstrate that the existence of a Noether symmetry in f ( R ) theories of gravity gives rise to an additional gravitational radius, besides the standard Schwarzschild one, determining the dynamics at galactic scales. By this feature, it is possible to explain the baryonic Tully-Fisher relation and the rotation curve of gas-rich galaxies without the dark matter hypothesis. Furthermore, under the same standard, the Fundamental Plane of elliptical galaxies can be addressed.

We study the vertical structure of polytropic, $P\propto \rho^\Gamma$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by the sound speed and surface density, $c_s^2/G\Sigma$, in SG discs, is smaller than either of these scales. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM halos with roughly flat circular velocity curves generally "flare" outwards. For the polytropic equation of state of the EAGLE simulations, discs whose mass and size match observational constraints are stable (NSG) for $M_d< 3\times 10^9\, M_\odot$ and unstable (SG) at higher masses, if fully gaseous. We test these analytic results using a set of idealized SPH simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are non-self-gravitating so the softening plays no role in their vertical structure. High-mass discs, on the other hand, are expected to be self-gravitating and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true: resolving instabilities does not guarantee that the vertical structure is resolved., Comment: Submitted to MNRAS

Observations of diffuse starlight in the outskirts of galaxies are thought to be a fundamental source of constraints on the cosmological context of galaxy assembly in the Λ CDM model. Such observations are not trivial because of the extreme faintness of such regions. In this work, we investigated the photometric properties of six massive early-type galaxies (ETGs) in the VST Elliptical GAlaxies Survey (VEGAS) sample (NGC 1399, NGC 3923, NGC 4365, NGC 4472, NGC 5044, and NGC 5846) out to extremely low surface brightness levels with the goal of characterizing the global structure of their light profiles for comparison to state-of-the-art galaxy formation models. We carried out deep and detailed photometric mapping of our ETG sample taking advantage of deep imaging with VST/OmegaCAM in the g and i bands. By fitting the light profiles, and comparing the results to simulations of elliptical galaxy assembly, we have identified signatures of a transition between relaxed and unrelaxed accreted components and can constrain the balance between in situ and accreted stars. The very good agreement of our results with predictions from theoretical simulations demonstrates that the full VEGAS sample of ∼ 100 ETGs will allow us to use the distribution of diffuse light as a robust statistical probe of the hierarchical assembly of massive galaxies.

Studying the coupling between the energy output produced by the central quasar and the host galaxy is fundamental to fully understand galaxy evolution. Quasar feedback is indeed supposed to dramatically affect the galaxy properties by depositing large amounts of energy and momentum into the ISM. In order to gain further insights on this process, we study the SEDs of sources at the brightest end of the quasar luminosity function, for which the feedback mechanism is supposed to be at its maximum. We model the rest-frame UV-to-FIR SEDs of 16 WISE-SDSS Selected Hyper-luminous (WISSH) quasars at 1.8 < z < 4.6 disentangling the different emission components and deriving physical parameters of both the nuclear component and the host galaxy. We also use a radiative transfer code to account for the contribution of the quasar-related emission to the FIR fluxes. Most SEDs are well described by a standard combination of accretion disk+torus and cold dust emission. However, about 30% of them require an additional emission component in the NIR, with temperatures peaking at 750K, which indicates the presence of a hotter dust component in these powerful quasars. We measure extreme values of both AGN bolometric luminosity (LBOL > 10^47 erg/s) and SFR (up to 2000 Msun/yr). A new relation between quasar and star-formation luminosity is derived (LSF propto LQSO^(0.73)) by combining several Herschel-detected quasar samples from z=0 to 4. Future observations will be crucial to measure the molecular gas content in these systems, probe the impact between quasar-driven outflows and on-going star-formation, and reveal the presence of merger signatures in their host galaxies., 19 pages, 12 figures; Accepted for publication in Astronomy & Astrophysics on June 13, 2017

Observations of diffuse starlight in the outskirts of galaxies are thought to be a fundamental source of constraints on the cosmological context of galaxy assembly in the $\Lambda$CDM model. Such observations are not trivial because of the extreme faintness of such regions. In this work, we investigate the photometric properties of six massive early type galaxies (ETGs) in the VEGAS sample (NGC 1399, NGC 3923, NGC 4365, NGC 4472, NGC 5044, and NGC 5846) out to extremely low surface brightness levels, with the goal of characterizing the global structure of their light profiles for comparison to state-of-the-art galaxy formation models. We carry out deep and detailed photometric mapping of our ETG sample taking advantage of deep imaging with VST/OmegaCAM in the g and i bands. By fitting the light profiles, and comparing the results to simulations of elliptical galaxy assembly, we identify signatures of a transition between "relaxed" and "unrelaxed" accreted components and can constrain the balance between in situ and accreted stars. The very good agreement of our results with predictions from theoretical simulations demonstrates that the full VEGAS sample of $\sim 100$ ETGs will allow us to use the distribution of diffuse light as a robust statistical probe of the hierarchical assembly of massive galaxies., Comment: Accepted for publication in Astronomy & Astrophysics

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process., Comment: 8 pages, 3 figures, MNRAS submitted

Context. The determination of the relative frequency of active galactic nuclei (AGN) versus other spectral classes, for example, HII region-like (HII), transition objects (TRAN), passive (PAS), and retired (RET), in a complete set of galaxies in the local Universe is of primary importance to discriminate the source of ionization in the nuclear region of galaxies. Aims. Here we aim to provide a spectroscopic characterization of the nuclei of galaxies belonging to the Herschel Reference Survey (HRS), a volume and magnitude limited sample representative of the local Universe, which has become a benchmark for local and high-z studies, for semianalytical models and cosmological simulations. The comparison between the nuclear spectral classification and the one determined on the global galactic scale provides information about how galaxy properties change from the nuclear to the outer regions. Moreover, the extrapolation of the global star formation (SF) properties from the SDSS fiber spectroscopy compared to the one computed by Halpha photometry can be useful for testing the method based on aperture correction for determining the global star formation rate (SFR) for local galaxies. Methods. By collecting the existing nuclear spectroscopy available from the literature, complemented with new observations obtained using the Loiano 1.52m telescope, we analyze the 322 nuclear spectra of HRS galaxies. Results. Using two diagnostic diagrams (the BPT and the WHAN) we provide a nuclear and an integrated spectral classification for the HRS galaxies. Conclusions. We find that the fraction of AGNs increases with stellar mass, such that at logM > 10.0 M\odot or 66% of the LTGs are AGNs or TRAN., 21 pages

We present the serendipitous ALMA detection of a faint submillimeter galaxy (SMG) lensed by a foreground z~1 galaxy. By optimizing the source detection to deblend the system, we accurately build the full spectral energy distribution of the distant galaxy from the I814 band to radio wavelengths. It is extremely red, with a I-K colour larger than 2.5. We estimate a photometric redshift of 3.28 and determine the physical parameters. The distant galaxy turns out to be magnified by the foreground lens by a factor of ~1.5, which implies an intrinsic Ks-band magnitude of ~24.5, a submillimeter flux at 870um of ~2.5 mJy and a SFR of ~150-300Msun/yr, depending on the adopted tracer. These values place our source towards the faint end of the distribution of observed SMGs, and in particular among the still few faint SMGs with a fully characterized spectral energy distribution, which allows us not only to accurately estimate its redshift but also to measure its stellar mass and other physical properties. The galaxy studied in this work is a representative of the population of faint SMGs, of which only few objects are known to date, that are undetected in optical and therefore are not typically accounted for when measuring the cosmic star formation history (SFH). This faint galaxy population thus likely represents an important and missing piece in our understanding of the cosmic SFH. Its observation and characterization is of major importance to achieve a solid picture of galaxy evolution., A&A, in press. Very minor changes to match the printed version

We present an automatic procedure to perform reliable photometry of galaxies on SDSS images. We selected a sample of 5853 galaxies in the Coma and Virgo superclusters. For each galaxy, we derive Petrosian g and i magnitudes, surface brightness profiles and color profiles. Unlike the SDSS pipeline, our procedure is not affected by the well known shredding problem and efficiently extracts Petrosian magnitudes for all galaxies. Hence we derived magnitudes even from the population of galaxies missed by the SDSS that represents 25% of all Local supercluster galaxies and ~95% of galaxies with g < 11 mag. After correcting the g and i magnitudes for Galactic and internal extinction, the blue and red sequences in the color magnitude diagram are well separated, with similar slopes. In addition, we study (i) the color-magnitude diagrams in different galaxy regions, the inner (r =0.35R{Pet})) zone, and (ii), we compute template color profiles, discussing the dependences of the templates on the galaxy masses and on their morphological type. The two analyses consistently lead to a picture where elliptical galaxies show no color gradients, irrespective of their masses. Spirals, instead, display a steeper gradient in their color profiles with increasing mass, consistently with the growing relevance of a bulge and/or a bar component above 10^{10} Modot., 18 pages, 12 figures, Accepted for publication by Astronomy & Astrophysics

We present a new and simple method to measure the instantaneous mass and radial growth rates of the stellar discs of spiral galaxies, based on their star formation rate surface density (SFRD) profiles. Under the hypothesis that discs are exponential with time-varying scalelengths, we derive a universal theoretical profile for the SFRD, with a linear dependence on two parameters: the specific mass growth rate $\nu_\textrm{M} \equiv \dot{M_\star}/M_\star$ and the specific radial growth rate $\nu_\textrm{R} \equiv \dot{R}_\star/R_\star$ of the disc. We test our theory on a sample of 35 nearby spiral galaxies, for which we derive a measurement of $\nu_\textrm{M}$ and $\nu_\textrm{R}$. 32/35 galaxies show the signature of ongoing inside-out growth ($\nu_\textrm{R} > 0$). The typical derived e-folding timescales for mass and radial growth in our sample are ~ 10 Gyr and ~ 30 Gyr, respectively, with some systematic uncertainties. More massive discs have a larger scatter in $\nu_\textrm{M}$ and $\nu_\textrm{R}$, biased towards a slower growth, both in mass and size. We find a linear relation between the two growth rates, indicating that our galaxy discs grow in size at ~ 0.35 times the rate at which they grow in mass; this ratio is largely unaffected by systematics. Our results are in very good agreement with theoretical expectations if known scaling relations of disc galaxies are not evolving with time., Comment: MNRAS, accepted. 14 pages, 4 figures, 3 tables. Additional material (Atlas.pdf) available at http://www.filippofraternali.com/downloads/index.html

We present the VST Early-type GAlaxy Survey (VEGAS), which is designed to obtain deep multiband photometry in g, r, i, of about one hundred nearby galaxies down to 27.3, 26.8, and 26 mag/arcsec^2 respectively, using the ESO facility VST/OmegaCAM. The goals of the survey are 1) to map the light distribution up to ten effective radii, r_e, 2) to trace color gradients and surface brightness fluctuation gradients out to a few r_e for stellar population characterization, and 3) to obtain a full census of the satellite systems (globular clusters and dwarf galaxies) out to 20% of the galaxy virial radius. The external regions of galaxies retain signatures of the formation and evolution mechanisms that shaped them, and the study of nearby objects enables a detailed analysis of their morphology and interaction features. To clarify the complex variety of formation mechanisms of early-type galaxies (ETGs), wide and deep photometry is the primary observational step, which at the moment has been pursued with only a few dedicated programs. The VEGAS survey has been designated to provide these data for a volume-limited sample with exceptional image quality. In this commissioning photometric paper we illustrate the capabilities of the survey using g- and i-band VST/OmegaCAM images of the nearby galaxy NGC 4472 and of smaller ETGs in the surrounding field. Our surface brightness profiles reach rather faint levels and agree excellently well with previous literature. Genuine new results concern the detection of an intracluster light tail in NGC 4472 and of various substructures at increasing scales. We have also produced extended (g-i) color profiles. The VST/OmegaCAM data that we acquire in the context of the VEGAS survey provide a detailed view of substructures in the optical emission from extended galaxies, which can be as faint as a hundred times below the sky level., Accepted for publication on A&A. 36 pages, 27 figures

The surface-brightness profiles of galaxies are well described by the S\'ersic law: systems with high S\'ersic index m have steep central profiles and shallow outer profiles, while systems with low m have shallow central profiles and steep outer profiles. R. Cen (2014, ApJL, 790, L24) has conjectured that these profiles arise naturally in the standard cosmological model with initial density fluctuations represented by a Gaussian random field (GRF). We explore and confirm this hypothesis with N-body simulations of dissipationless collapses in which the initial conditions are generated from GRFs with different power spectra. The numerical results show that GRFs with more power on small scales lead to systems with higher m. In our purely dissipationless simulations the S\'ersic index is in the range 2, Comment: 6 pages, 2 figures, accepted for publication in the ApJL

International audience; A growing body of evidence indicates that the star formation rate per unit stellar mass (sSFR) decreases with increasing mass in normal main-sequence star-forming galaxies. Many processes have been advocated as being responsible for this trend (also known as mass quenching), e.g., feedback from active galactic nuclei (AGNs), and the formation of classical bulges. In order to improve our insight into the mechanisms regulating the star formation in normal star-forming galaxies across cosmic epochs, we determine a refined star formation versus stellar mass relation in the local Universe. To this end we use the H alpha narrow-band imaging follow-up survey (H alpha 3) of field galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) in the Coma and Local superclusters. By complementing this local determination with high-redshift measurements from the literature, we reconstruct the star formation history of main-sequence galaxies as a function of stellar mass from the present epoch up to z = 3. In agreement with previous studies, our analysis shows that quenching mechanisms occur above a threshold stellar mass M-knee that evolves with redshift as proportional to(1 + z)(2). Moreover, visual morphological classification of individual objects in our local sample reveals a sharp increase in the fraction of visually classified strong bars with mass, hinting that strong bars may contribute to the observed downturn in the sSFR above Mknee. We test this hypothesis using a simple but physically motivated numerical model for bar formation, finding that strong bars can rapidly quench star formation in the central few kpc of field galaxies. We conclude that strong bars contribute significantly to the red colors observed in the inner parts of massive galaxies, although additional mechanisms are likely required to quench the star formation in the outer regions of massive spiral galaxies. Intriguingly, when we extrapolate our model to higher redshifts, we successfully recover the observed redshift evolution for M-knee. Our study highlights how the formation of strong bars in massive galaxies is an important mechanism in regulating the redshift evolution of the sSFR for field main-sequence galaxies.

We present new Halpha+[NII] imaging data of late-type galaxies in the Herschel Reference Survey aimed at studying the star formation properties of a K-band-selected, volume-limited sample of nearby galaxies. The Halpha+[NII] data are corrected for [NII] contamination and dust attenuation using different recipes based on the Balmer decrement and the 24mic luminosities. We show that the L(Halpha) derived with different corrections give consistent results only whenever the uncertainty on the estimate of the Balmer decrement is, Accepted for publication on A&A

We studied the chronology of galactic bulge and disc formation by analysing the relative contributions of these components to the B-band rest-frame luminosity density at different epochs. We present the first estimate of the evolution of the fraction of rest-frame B-band light in galactic bulges and discs since redshift z~0.8. We performed a bulge-to-disc decomposition of HST/ACS images of 3266 galaxies in the zCOSMOS-bright survey with spectroscopic redshifts in the range 0.7 < z < 0.9. We find that the fraction of B-band light in bulges and discs is $(26 \pm 4)%$ and $(74 \pm 4)%$, respectively. When compared with rest-frame B-band measurements of galaxies in the local Universe in the same mass range ($10^{9} M_{\odot}\lessapprox M \lessapprox 10^{11.5} M_{\odot}$), we find that the B-band light in discs decreases by ~30% from z~0.7-0.9 to z~0, while the light from the bulge increases by ~30% over the same period of time. We interpret this evolution as the consequence of star formation and mass assembly processes, as well as morphological transformation, which gradually shift stars formed at half the age of the Universe from star-forming late-type/irregular galaxies toearlier types and ultimately into spheroids., Letter to the Editor, 4 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (, Letter in A&A 568, L8 (2014). This replacement matches the published version

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z, Accepted for publication in A&A

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., 15 pages, 20 images, accepted for publication in MNRAS: MN-12-1739-MJ.R1

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy., ApJ in press, 16 pages, 7 figures. Updated to address referee comments

The chemical evolution of galaxies on a cosmological timescale is still a matter of debate despite the increasing number of available data provided by spectroscopic surveys of star-forming galaxies at different redshifts. The fundamental relations involving metallicity, such as the mass-metallicity relation (MZR) or the fundamental-metallicity relation, give controversial results about the reality of evolution of the chemical content of galaxies at a given stellar mass. In this work we shed some light on this issue using the completeness reached by the 20k bright sample of the zCOSMOS survey and using for the first time the nitrogen-to-oxygen ratio (N/O) as a tracer of the gas phase chemical evolution of galaxies that is independent of the star formation rate. Emission-line galaxies both in the SDSS and 20k zCOSMOS bright survey were used to study the evolution from the local Universe of the $MZR up to a redshift of 1.32 and the relation between stellar mass and nitrogen-to-oxygen ratio (MNOR) up to a redshift of 0.42 using the N2S2 parameter. All the physical properties derived from stellar continuum and gas emission-lines, including stellar mass, star formation rates, metallicity and N/O, were calculated in a self-consistent way over the full redshift range. We confirm the trend to find lower metallicities in galaxies of a given stellar mass in a younger Universe. This trend is even observed when taking possible selection effects into account that are due to the observed larger median star formation rates for galaxies at higher redshifts. We also find a significant evolution of the MNOR up to z = 0.4. Taking the slope of the O/H vs. N/O relation into account for the secondary-nitrogen production regime, the observed evolution of the MNOR is consistent with the trends found for both the MZR and its equivalent relation using new expressions to reduce its dependence on SFR., Comment: 15 pages, 13 figures. Accepted for publication in Astronomy & Astrophysics. Replaced to match published version and references corrected

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7, Comment: 25 pages, 15 figures, accepted to ApJ; v3 include many minor modifications following the referee report, added a section comparing various z~2 IR template calibrations

X-ray surveys contain sizable numbers of star forming galaxies, beyond the AGN which usually make the majority of detections. Many methods to separate the two populations are used in the literature, based on X-ray and multiwavelength properties. We aim at a detailed test of the classification schemes and to study the X-ray properties of the resulting samples. We build on a sample of galaxies selected at 1.4 GHz in the VLA-COSMOS survey, classified by Smolcic et al. (2008) according to their optical colours and observed with Chandra. A similarly selected control sample of AGN is also used for comparison. We review some X-ray based classification criteria and check how they affect the sample composition. The efficiency of the classification scheme devised by Smolcic et al. (2008) is such that ~30% of composite/misclassified objects are expected because of the higher X-ray brightness of AGN with respect to galaxies. The latter fraction is actually 50% in the X-ray detected sources, while it is expected to be much lower among X-ray undetected sources. Indeed, the analysis of the stacked spectrum of undetected sources shows, consistently, strongly different properties between the AGN and galaxy samples. X-ray based selection criteria are then used to refine both samples. The radio/X-ray luminosity correlation for star forming galaxies is found to hold with the same X-ray/radio ratio valid for nearby galaxies. Some evolution of the ratio may be possible for sources at high redshift or high luminosity, tough it is likely explained by a bias arising from the radio selection. Finally, we discuss the X-ray number counts of star forming galaxies from the VLA- and C-COSMOS surveys according to different selection criteria, and compare them to the similar determination from the Chandra Deep Fields. The classification scheme proposed here may find application in future works and surveys., Comment: 15 pages, 7 figures, 3 tables

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

Aims. The aim of this work is to study the contribution of the Ly-a emitters (LAE) to the star formation rate density (SFRD) of the Universe in the interval 21.5x10^18 erg/s/cm^2 are 33 galaxies/arcmin^2 and 4x10^-2 Mpc^-3, respectively. We find that the the observed luminosity function of LAE does not evolve from z=2 to z=6. This implies that, after correction for the redshift-dependent IGM absorption, the intrinsic LF must have evolved significantly over 3 Gyr. The SFRD from LAE contributes to about 20% of the SFRD at z =2-3, while the LAE appear to be the dominant source of star formation producing ionizing photons in the early universe z>5-6, becoming equivalent to that of Lyman Break galaxies., Comment: 17 pages, 20 figures, accepted for publication to A&A

We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 \times 10^9 M\odot - 2 \times 10^11 M\odot with a mean value of ~ 2 \times 10^10 M\odot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically "hot", with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet "in place" at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role., A&A accepted

With the aim of quantifying the contribution of the environment on the evolution of galaxies at z=0 we have used the DR7 catalogue of the Sloan Digital Sky Survey (SDSS) to reconstruct the 3-D distribution of 4132 galaxies in 420 square degrees of the Coma supercluster, containing two rich clusters (Coma and A1367), several groups, and many filamentary structures belonging to the "Great Wall", at the approximate distance of 100 Mpc. At this distance the galaxy census is complete to Mi=-17.5 mag, i.e. approx 4 mag fainter than M*. The morphological classification of galaxies into early- (ellipticals) and late-types (spirals) was carried out by inspection of individual SDSS images and spectra. The density around each galaxies was determined in cylinders of 1 Mpc radius and 1000 km s^-1 half length. The color-luminosity relation was derived for galaxies in bins morphological type and in four thresholds of galaxy density-contrast, ranging from delta{1,1000} 20 (UH = the cluster's cores). The fraction of early-type galaxies increases with the log of the over-density. A well defined "red sequence" composed of early-type galaxies exists in all environments at high luminosity, but it lacks of low luminosity (dwarf) galaxies in the lowest density environment. Conversely low luminosity isolated galaxies are predominantly of late-type. In other words the low luminosity end of the distribution is dominated by red dE galaxies in clusters and groups and by dwarf blue amorphous systems in the lowest density regions. At z=0 we find evidence for strong evolution induced by the environment (Nurture). Transformations take place mostly at low luminosity when star forming dwarf galaxies inhabiting low density environments migrate into amorphous passive dwarf ellipticals in their infall into denser regions. The mechanism involves suppression of the star formation due to gas stripping, without significant mass growth, as proposed by Boselli et al. (2008a). This process is more efficient and fast in ambients of increasing density. In the highest density environments (around clusters) the truncation of the star formation happens fast enough (few 100 Myr) to produce the signature of post-star-burst in galaxy spectra. PSB galaxies, that are in fact found significantly clustered around the largest dynamical units, represent the remnants of star forming isolated galaxies that had their star formation violently suppressed during their infall in clusters in the last 0.5-1.5 Gyrs, and the progenitors of future dEs., 14 pages, 14 figures, Astronomy and Astrophysics, in press

We report on a complete CCD imaging survey of 226 elliptical galaxies in the North-East quadrant of the Virgo cluster, representative of the properties of giant and dwarf elliptical galaxies in this cluster. We fit their radial light profiles with the Sersic r^(1/n) model of light distribution. We confirm the result of Graham & Guzman (2003) that the apparent dichotomy between E and dE galaxies in the luminosity-_e plane no longer appears when other structural parameters are considered and can be entirely attributed to the onset of "core" galaxies at B_T sim -20.5 mag. When "core" galaxies are not considered, E and dE form a unique family with n linearly increasing with the luminosity. For 90 galaxies we analyze the B-I color indices, both in the nuclear and in the outer regions. Both indices are bluer toward fainter luminosities. We find also that the outer color gradients do not show any significant correlation with the luminosity. The scatter in all color indicators increases significantly toward lower luminosities, e.g. galaxies fainter than B_T sim -15 have a B-I spread 0.5 mag., 21 pages, 14 figures, 2 tables. Accepted for publication on A&A. High resolution figures at http://goldmine.mib.infn.it/papers/WFS_04/WFS_04-frame.html

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2, Comment: 12 pages, 6 figures, A&A in press

Recently, the temperature T and luminosity L_X of the hot gas halos of early type galaxies have been derived with unprecedented accuracy from Chandra data, for 30 galaxies covering a wider range of galactic luminosity (and central velocity dispersion sigma_c) than before. This work investigates the origin of the observed temperatures, by examining the relationship between them and the galaxy structure, the gas heating due to Type Ia supernovae (SNIa's) and the gravitational potential, and the dynamical status of the gas flow. In galaxies with sigma_c, Comment: 16 pages, accepted for publication in ApJ

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z ∼ 2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5 < z < 2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera and Spectrometer (PACS) onboard Herschel, as part of the PACS evolutionary probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 μm and UV SFRs. We find that SFRs from 24 μm alone are higher by a factor of ∼4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24 μm) > 12.2 L ⊙. The SFGs and AGNs tend to exhibit the same 24 μm excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions. © 2010 ESO.

We study the evolution of the observed correlation between central stellar velocity dispersion $\sigma_\mathrm{e}$ and stellar mass $M_*$ of massive ($M_*\gtrsim 3\times 10^{10}\,\mathrm{M_\odot}$) early-type galaxies (ETGs) out to redshift $z\approx 2.5$, exploiting a Bayesian hierarchical inference formalism. Collecting ETGs from state-of-the-art literature samples, we build a $fiducial$ sample ($0\lesssim z\lesssim 1$), which is obtained with homogeneous selection criteria, but also a less homogeneous $extended$ sample ($0\lesssim z\lesssim 2.5$). Based on the fiducial sample, we find that the $M_*$-$\sigma_\mathrm{e}$ relation is well represented by $\sigma_\mathrm{e}\propto M_*^{\beta}(1+z)^{\zeta}$, with $\beta\simeq 0.18$ independent of redshift and $\zeta\simeq 0.4$ (at given $M_*$, $\sigma_\mathrm{e}$ decreases for decreasing $z$, for instance by a factor of $\approx1.3$ from $z=1$ to $z=0$). When the slope $\beta$ is allowed to evolve, we find it increasing with redshift: $\beta(z)\simeq 0.16+0.26\log(1+z)$ describes the data as well as constant $\beta\simeq 0.18$. The intrinsic scatter of the $M_*$-$\sigma_\mathrm{e}$ relation is $\simeq0.08$ dex in $\sigma_\mathrm{e}$ at given $M_*$, independent of redshift. Our results suggest that, on average, the velocity dispersion of $individual$ massive ($M_*\gtrsim 3\times 10^{11}\,M_\odot$) ETGs decreases with time while they evolve from $z\approx 1$ to $z\approx 0$. The analysis of the extended sample leads to results similar to that of the fiducial sample, with slightly stronger redshift dependence of the normalisation ($\zeta\simeq 0.5$) and weaker redshift dependence of the slope (${\rm d} \beta/{\rm d} \log (1+z)\simeq 0.18$) when $\beta$ varies with time. At $z=2$ ETGs with $M_*\approx 10^{11}\,M_\odot$ have, on average, $\approx1.7$ higher $\sigma_\mathrm{e}$ than ETGs of similar stellar mass at $z=0$., Comment: 26 pages, 15 figures, accepted for publication in MNRAS