Search

Context: New spectroscopic surveys will increase the number of astronomical objects requiring characterization by over tenfold.. Machine learning tools are required to address this data deluge in a fast and accurate fashion. Most machine learning algorithms can not estimate error directly, making them unsuitable for reliable science. Aims: We aim to train a supervised deep-learning algorithm tailored for high-resolution observational stellar spectra. This algorithm accurately infer precise estimates while providing coherent estimates of uncertainties by leveraging information from both the neural network and the spectra. Methods: We train a conditional Invertible Neural Network (cINN) on observational spectroscopic data obtained from the GIRAFFE spectrograph (HR10 and HR21 setups) within the Gaia-ESO survey. A key features of cINN is its ability to produce the Bayesian posterior distribution of parameters for each spectrum. By analyzing this distribution, we inferred parameters and their uncertainties. Several tests have been applied to study how parameters and errors are estimated. Results: We achieved an accuracy of 28K in $T_{\text{eff}}$, 0.06 dex in $\log g$, 0.03 dex in $[\text{Fe/H}]$, and between 0.05 dex and 0.17 dex for the other abundances for high quality spectra. Accuracy remains stable with low signal-to-noise ratio spectra. The uncertainties obtained are well within the same order of magnitude. The network accurately reproduces astrophysical relationships both on the scale of the Milky Way and within smaller star clusters. We created a table containing the new parameters generated by our cINN. Conclusion: This neural network represents a compelling proposition for future astronomical surveys. These coherent derived uncertainties make it possible to reuse these estimates in other works as Bayesian priors and thus present a solid basis for future work.

We present spatially resolved, simple stellar population equivalent ages, stellar metallicities, and abundance ratios for 456 massive ($10.3\lesssim\log(\mathrm{M}_*/\mathrm{M}_\odot)\lesssim11.8$) quiescent galaxies at $0.6\lesssim z\lesssim1.0$ from the Large Early Galaxy Astrophysics Census, derived using full-spectrum models. Typically, we find flat age and [Mg/Fe] gradients, and negative [Fe/H] gradients, implying iron-rich cores. We also estimate intrinsic [Fe/H] gradients via forward modelling. We examine the observed gradients in three age bins. Younger quiescent galaxies typically have negative [Fe/H] gradients and positive age gradients, possibly indicating a recent central starburst. Additionally, this finding suggests that photometrically measured flat colour gradients in young quiescent galaxies are the result of the positive age and negative metallicity gradients cancelling each other. For older quiescent galaxies, the age gradients become flat and [Fe/H] gradients weaken, though remain negative. Thus, negative colour gradients at older ages are likely driven by metallicity gradients. The diminishing age gradient may result from the starburst fading. Furthermore, the persistence of the [Fe/H] gradients may suggest that the outskirts are simultaneously built up by mergers with lower metallicity satellites. On the other hand, the gradients could be inherited from the star-forming phase, in which case mergers may not be needed to explain our findings. This work illustrates the need for resolved spectroscopy, instead of just photometry, to measure stellar population gradients. Extending these measurements to higher redshift is imperative for understanding how stellar populations in quiescent galaxies are assembled over cosmic time., Comment: Accepted for publication in MNRAS; minor typesetting corrections after copyediting

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6

We employ a recently developed population-orbit superposition technique to simultaneously fit the stellar kinematic and age maps of 82 CALIFA spiral galaxies and obtain the ages of stars in different dynamical structures. We first evaluated the capabilities of this method on CALIFA-like mock data created from the Auriga simulations. The recovered mean ages of dynamically cold, warm, and hot components match the true values well, with an observational error of up to $20\%$ in the mock age maps. For CALIFA spiral galaxies, we find that the stellar ages of the cold, warm, and hot components all increase with the stellar mass of the galaxies, from $\overline{t_{\rm cold}}\sim2.2$ Gyr, $\overline{t_{\rm warm}}\sim2.3$ Gyr, and $\overline{t_{\rm hot}}\sim2.6$ Gyr for galaxies with stellar mass $M_*<10^{10}\,\rm M_{\odot}$, to $\overline{t_{\rm cold}}\sim4.0$ Gyr, $\overline{t_{\rm warm}}\sim5.1$ Gyr, and $\overline{t_{\rm hot}}\sim5.9$ Gyr for galaxies with $M_*>10^{11}\,\rm M_{\odot}$. About $80\%$ of the galaxies in our sample have $t_{\rm hot}>t_{\rm cold}$, and the mean values of $t_{\rm hot}-t_{\rm cold}$ also increase with stellar mass, from $0.7_{-0.2}^{+0.6}$ Gyr in low-mass galaxies ($10^{8.9}\,\rm M_{\odot}

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe., Comment: 28 pages, 12 figures, 4 tables, Accepted to ApJ

We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $\sigma$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $\sigma$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent., Comment: 19 pages, 9 figures Accepted for publication in MNRAS

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity., Comment: 26 pages, 15 figures Accepted for publication in MNRAS

We present a nonparametric morphology analysis of the stellar continuum and nebular emission lines for a sample of local galaxies. We explore the dependence of the various morphological parameters on wavelength and morphological type. Our goal is to quantify the difference in morphology between the stellar and nebular components. We derive the nonparametric morphological indicators of 364 galaxies from the CALIFA Survey. To calculate those indicators, we apply the StatMorph package on the high-quality integral field spectroscopic data cubes, as well as to the most prominent nebular emission-line maps, namely [OIII]$\lambda$5007, H$\alpha$, and [NII]$\lambda$6583. We show that the physical size of galaxies, M$_{20}$ index, and concentration have a strong gradient from blue to red optical wavelengths. We find that the light distribution of the nebular emission is less concentrated than the stellar continuum. A comparison between the nonparametric indicators and the galaxy physical properties revealed a very strong correlation of the concentration with the specific star-formation rate and morphological type. Furthermore, we explore how the galaxy inclination affects our results. We find that edge-on galaxies show a more rapid change in physical size and concentration with increasing wavelength due to the increase in optical free path. We conclude that the apparent morphology of galaxies originates from the pure stellar distribution, but the morphology of the ISM presents differences with respect to the morphology of the stellar component. Our analysis also highlights the importance of dust attenuation and galaxy inclination in the measurement of nonparametric morphological indicators, especially in the the wavelength range 4000-5000 \r{A}.

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator., Comment: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

The properties of the stellar populations in a galaxy are known to correlate with the amount and the distribution of stellar mass. We take advantage of the maps of light-weighted mean stellar age Agewr and metallicity Z*wr for a sample of 362 galaxies from the integral-field spectroscopic survey CALIFA (summing up to >600,000 individual regions of approximately 1 kpc linear size), produced in our previous works, to investigate how these local properties react to the local stellar-mass surface density mu* and to the global total stellar mass M* and mean stellar-mass surface density e. We establish the existence of i) a dual mu*-Agewr relation, resulting in a young sequence and an old ridge, and ii) a mu*-Z*wr relation, overall independent of the age of the regions. The global mass parameters (M* and, possibly secondarily, e) determine the distribution of mu* in a galaxy and set the maximum attainable mu*, which increases with M*. M* affects the shape and normalization of the local relations up to a threshold mass of $\sim 10^{10.3}$ MSun, above which they remain unchanged. We conclude that stellar mass is a "glocal" (i.e. simultaneously global and local) driver of the stellar population properties. We consider how the local and global mass-age and mass-metallicity relations are connected, and in particular discuss how it is possible, from a single local relation, to produce two different global mass-metallicity relations for quiescent and star-forming galaxies respectively, as reported in the literature. Structural differences in these two classes of galaxies are key to explain the duality in global scaling relations and appear as essential in modelling the baryonic cycle of galaxies., Comment: Accepted for publication on MNRAS

We present the first comparison of observed stellar continuum spectra of high-redshift galaxies and mock galaxy spectra generated from hydrodynamical simulations. The mock spectra are produced from the IllustrisTNG TNG100 simulation combined with stellar population models and take into account dust attenuation and realistic observational effects (aperture effects and noise). We compare the simulated $D_n4000$ and EW(H$\delta$) of galaxies with $10.5 \leq \log(M_\ast/M_\odot) \leq 11.5$ at $0.6 \leq z \leq 1.0$ to the observed distributions from the LEGA-C survey. TNG100 globally reproduces the observed distributions of spectral indices, implying that the age distribution of galaxies in TNG100 is generally realistic. Yet there are small but significant differences. For old galaxies, TNG100 shows small $D_n4000$ when compared to LEGA-C, while LEGA-C galaxies have larger EW(H$\delta$) at fixed $D_n4000$. There are several possible explanations: 1) LEGA-C galaxies have overall older ages combined with small contributions (a few percent in mass) from younger ($<1$~Gyr) stars, while TNG100 galaxies may not have such young sub-populations; 2) the spectral mismatch could be due to systematic uncertainties in the stellar population models used to convert stellar ages and metallicities to observables. In conclusion, the latest cosmological galaxy formation simulations broadly reproduce the global age distribution of galaxies at $z\sim1$ and, at the same time, the high quality of the latest observed and simulated datasets help constrain stellar population synthesis models as well as the physical models underlying the simulations., Comment: 17 pages, 12 figures, accepted by AJ

We present the third and final data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO/VLT public spectroscopic survey targeting $0.6 < z < 1.0$, Ks-selected galaxies. The data release contains 3528 spectra with measured stellar velocity dispersions and stellar population properties, a 25-fold increase in sample size compared to previous work. This $K_s$-selected sample probes the galaxy population down to $\sim0.3 L^*$, for all colors and morphological types. Along with the spectra we publish a value-added catalog with stellar and ionized gas velocity dispersions, stellar absorption line indices, emission line fluxes and equivalent widths, complemented with structural parameters measured from HST/ACS imaging. With its combination of high precision and large sample size, LEGA-C provides a new benchmark for galaxy evolution studies., Comment: Accepted for publication in ApJS. Data available at http://archive.eso.org/cms/eso-archive-news/Third-and-final-release-of-the-Large-Early-Galaxy-Census-LEGA-C-Spectroscopic-Public-Survey-published.html or https://users.ugent.be/\string~avdrwel/research.html\#legac

This work presents an in-depth analysis of four gravitationally lensed red galaxies at z = 1.6-3.2. The sources are magnified by factors of 2.7-30 by foreground clusters, enabling spectral and morphological measurements that are otherwise challenging. Our sample extends below the characteristic mass of the stellar mass function and is thus more representative of the quiescent galaxy population at z > 1 than previous spectroscopic studies. We analyze deep VLT/X-SHOOTER spectra and multi-band Hubble Space Telescope photometry that cover the rest-frame UV-to-optical regime. The entire sample resembles stellar disks as inferred from lensing-reconstructed images. Through stellar population synthesis analysis we infer that the targets are young (median age = 0.1-1.2 Gyr) and formed 80% of their stellar masses within 0.07-0.47 Gyr. Mg II $\lambda\lambda 2796,2803$ absorption is detected across the sample. Blue-shifted absorption and/or redshifted emission of Mg II is found in the two youngest sources, indicative of a galactic-scale outflow of warm ($T\sim10^{4}$ K) gas. The [O III] $\lambda5007$ luminosity is higher for the two young sources (median age less than 0.4 Gyr) than the two older ones, perhaps suggesting a decline in nuclear activity as quenching proceeds. Despite high-velocity ($v\approx1500$ km s$^{-1}$) galactic-scale outflows seen in the most recently quenched galaxies, warm gas is still present to some extent long after quenching. Altogether our results indicate that star formation quenching at high redshift must have been a rapid process (< 1 Gyr) that does not synchronize with bulge formation or complete gas removal. Substantial bulge growth is required if they are to evolve into the metal-rich cores of present-day slow-rotators., Comment: Accepted for publication in the Astrophysical Journal. 37 pages, 20 figures, 10 tables

We examine the Fundamental Plane (FP) and mass-to-light ratio ($M/L$) scaling relations using the largest sample of massive quiescent galaxies at $1.511.26$, subset of 8 quiescent galaxies at $z>2$, from Stockmann et al. (2020), we show that they cannot passively evolve to the local Coma cluster relation alone and must undergo significant structural evolution to mimic the sizes of local massive galaxies. The evolution of the FP and $M/L$ scaling relations, from $z=2$ to present-day, for this subset are consistent with passive aging of the stellar population and minor merger structural evolution into the most massive galaxies in the Coma cluster and other massive elliptical galaxies from the MASSIVE Survey. Modeling the luminosity evolution from minor merger added stellar populations favors a history of merging with "dry" quiescent galaxies., Comment: 17 pages, 6 figures

We explore how the star formation and metal enrichment histories of present-day galaxies have been affected by environment combining stellar population parameter estimates and group environment characterization for SDSS DR7. We compare stellar ages, stellar metallicities and element abundance ratios [alpha/Fe] of satellite and central galaxies, as a function of their stellar and host group halo mass, controlling for the current star formation rate and for the infall epoch. We confirm that below log(Mstar/Msun)=10.5 satellites are older and metal-richer than equally-massive central galaxies. On the contrary, we do not detect any difference in their [alpha/Fe]: this depends primarily on stellar mass and not on group hierarchy nor host halo mass. We also find that the differences in the median age and metallicity of satellites and centrals at stellar mass below 10^{10.5}Msun are largely due to the higher fraction of passive galaxies among satellites and as a function of halo mass. We argue that the observed trends at low masses reveal the action of satellite-specific environmental effects in a `delayed-then-rapid' fashion. When accounting for the varying quiescent fraction, small residual excess in age, metallicity and [alpha/Fe] emerge for satellites dominated by old stellar populations and residing in halos more massive than 10^{14}Msun, compared to equally-massive central galaxies. This excess in age, metallicity and [alpha/Fe] pertain to ancient infallers, i.e. satellites that have accreted onto the current halo more than 5 Gyr ago. This result points to the action of environment in the early phases of star formation in galaxies located close to cosmic density peaks., Comment: accepted for publication on MNRAS, minor revisions after referee report, 24 pages, 12 figures, 3 Appendices

We present the Phase A Science Case for the Multi-conjugate Adaptive-optics Visible Imager-Spectrograph (MAVIS), planned for the Adaptive Optics Facility (AOF) of the Very Large Telescope (VLT). MAVIS is a general-purpose instrument for exploiting the highest possible angular resolution of any single optical telescope available in the next decade, either on Earth or in space, and with sensitivity comparable to (or better than) larger aperture facilities. MAVIS uses two deformable mirrors in addition to the deformable secondary mirror of the AOF, providing a mean V-band Strehl ratio of >10% (goal >15%) across a relatively large (30 arc second) science field. This equates to a resolution of <20mas at 550nm - comparable to the K-band diffraction limit of the next generation of extremely large telescopes, making MAVIS a genuine optical counterpart to future IR-optimised facilities like JWST and the ELT. Moreover, MAVIS will have unprecedented sky coverage for a high-order AO system, accessing at least 50% of the sky at the Galactic Pole, making MAVIS a truly general purpose facility instrument. As such, MAVIS will have both a Nyquist-sampled imager (30x30 arcsec field), and a powerful integral field spectrograph with multiple spatial and spectral modes spanning 370-1000nm. This science case presents a distilled set of thematically linked science cases drawn from the MAVIS White Papers (www.mavis-ao.org/whitepapers), selected to illustrate the driving requirements of the instrument resulting from the recent MAVIS Phase A study., Comment: 141 pages, 79 figures, 10Mb. Modified version of the Phase A Science Case submitted to ESO in April 2020 as part of the MAVIS Phase A design study. Full resolution version available upon request

We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $H{\delta}_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $H{\delta}_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $H{\delta}_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ \r{A}/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ \r{A}$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius., Comment: 18 pages, 9 figures

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years., Comment: 30 pages, 10 figures, accepted in ApJ

We perform spatially resolved stellar population analysis for a sample of 69 early-type galaxies (ETGs) from the CALIFA integral field spectroscopic survey, including 48 ellipticals and 21 S0's. We generate and quantitatively characterize profiles of light-weighted mean stellar age and metallicity within $\lesssim 2R_e$, as a function of radius and stellar-mass surface density $\mu_*$. We study in detail the dependence of profiles on galaxies' global properties, including velocity dispersion $\sigma_e$, stellar mass, morphology. ETGs are universally characterized by strong, negative metallicity gradients ($\sim -0.3\,\text{dex}$ per $R_e$) within $1\,R_e$, which flatten out moving towards larger radii. A quasi-universal local $\mu_*$-metallicity relation emerges, which displays a residual systematic dependence on $\sigma_e$, whereby higher $\sigma_e$ implies higher metallicity at fixed $\mu_*$. Age profiles are typically U-shaped, with minimum around $0.4\,R_e$, asymptotic increase to maximum ages beyond $\sim 1.5\,R_e$, and an increase towards the centre. The depth of the minimum and the central increase anti-correlate with $\sigma_e$. A possible qualitative interpretation of these observations is a two-phase scenario. In the first phase, dissipative collapse occurs in the inner $1\,R_e$, establishing a negative metallicity gradient. The competition between the outside-in quenching due to feedback-driven winds and some form of inside-out quenching, possibly caused by central AGN feedback or dynamical heating, determines the U-shaped age profiles. In the second phase, the accretion of ex-situ stars from quenched and low-metallicity satellites shapes the flatter stellar population profiles in the outer regions., Comment: Accepted version to appear in MNRAS. Improved discussion with respect to original submission and additional tests included

We analyze radial stellar metallicity and kinematic profiles out to 1Re in 244 CALIFA galaxies ranging from morphological type E to Sd, to study the evolutionary mechanisms of stellar population gradients. We find that linear metallicity gradients exhibit a clear correlation with galaxy morphological type - with early type galaxies showing the steepest gradients. We show that the metallicity gradients simply reflect the local mass-metallicity relation within a galaxy. This suggests that the radial stellar population distribution within a galaxys effective radius is primarily a result of the \emph{in-situ} local star formation history. In this simple picture, the dynamically derived stellar surface mass density gradient directly predicts the metallicity gradient of a galaxy. We show that this correlation and its scatter can be reproduced entirely by using independent empirical galaxy structural and chemical scaling relations. Using Schwarzschild dynamical models, we also explore the link between a galaxys local stellar populations and their orbital structures. We find that galaxies angular momentum and metallicity gradients show no obvious causal link. This suggests that metallicity gradients in the inner disk are not strongly shaped by radial migration, which is confirmed by the lack of correlation between the metallicity gradients and observable probes of radial migration in the galaxies, such as bars and spiral arms. Finally, we find that galaxies with positive metallicity gradients become increasingly common towards low mass and late morphological types - consistent with stellar feedback more efficiently modifying the baryon cycle in the central regions of these galaxies., Comment: 20 pages, 13 figures

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context., Comment: Published in Nature Astronomy, 1 January 2018; doi:10.1038/s41550-017-0348-1; 22pages, 8 figures, 2 tables

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies., Comment: Nature (June 22, 2017)

Galaxies in the local Universe are known to follow bimodal distributions in the global stellar populations properties. We analyze the distribution of the local average stellar-population ages of 654,053 sub-galactic regions resolved on ~1-kpc scales in a volume-corrected sample of 394 galaxies, drawn from the CALIFA-DR3 integral-field-spectroscopy survey and complemented by SDSS imaging. We find a bimodal local-age distribution, with an old and a young peak primarily due to regions in early-type galaxies and star-forming regions of spirals, respectively. Within spiral galaxies, the older ages of bulges and inter-arm regions relative to spiral arms support an internal age bimodality. Although regions of higher stellar-mass surface-density, mu*, are typically older, mu* alone does not determine the stellar population age and a bimodal distribution is found at any fixed mu*. We identify an "old ridge" of regions of age ~9 Gyr, independent of mu*, and a "young sequence" of regions with age increasing with mu* from 1-1.5 Gyr to 4-5 Gyr. We interpret the former as regions containing only old stars, and the latter as regions where the relative contamination of old stellar populations by young stars decreases as mu* increases. The reason why this bimodal age distribution is not inconsistent with the unimodal shape of the cosmic-averaged star-formation history is that i) the dominating contribution by young stars biases the age low with respect to the average epoch of star formation, and ii) the use of a single average age per region is unable to represent the full time-extent of the star-formation history of "young-sequence" regions., Comment: 17 pages, 11 figures, MNRAS accepted

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., Comment: 14 pages, 8 figures, 2 tables, submitted to MNRAS

Mass models of galactic disks traditionally rely on axisymmetric density and rotation curves, paradoxically acting as if their most remarkable asymmetric features, such as lopsidedness or spiral arms, were not important. In this article, we relax the axisymmetry approximation and introduce a methodology that derives 3D gravitational potentials of disk-like objects and robustly estimates the impacts of asymmetries on circular velocities in the disk midplane. Mass distribution models can then be directly fitted to asymmetric line-of-sight velocity fields. Applied to the grand-design spiral M99, the new strategy shows that circular velocities are highly nonuniform, particularly in the inner disk of the galaxy, as a natural response to the perturbed gravitational potential of luminous matter. A cuspy inner density profile of dark matter is found in M99, in the usual case where luminous and dark matter share the same center. The impact of the velocity nonuniformity is to make the inner profile less steep, although the density remains cuspy. On another hand, a model where the halo is core dominated and shifted by 2.2-2.5 kpc from the luminous mass center is more appropriate to explain most of the kinematical lopsidedness evidenced in the velocity field of M99. However, the gravitational potential of luminous baryons is not asymmetric enough to explain the kinematical lopsidedness of the innermost regions, irrespective of the density shape of dark matter. This discrepancy points out the necessity of an additional dynamical process in these regions: possibly a lopsided distribution of dark matter., Comment: Accepted for publication in Astronomy & Astrophysics, 20 pages, 17 figures. 2nd version after language editing

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

We present three independent catalogs of point-sources extracted from SPIRE images at 250, 350, and 500 micron as a part of the Herschel Virgo Cluster Survey (HeViCS). The source positions are determined by estimating the likelihood to be a real source for each peak on the maps and the flux densities are estimated using the sourceExtractorTimeline, a timeline-based point source fitter. Afterwards, each source is subtracted from the maps, removing a Gaussian function in every position with the full width half maximum equal to that estimated in sourceExtractorTimeline. This procedure improves the robustness of our algorithm in terms of source identification. The HeViCS catalogs contain about 52000, 42200, and 18700 sources selected at 250, 350, and 500 micron above 3sigma and are ~ 75%, 62%, and 50% complete at flux densities of 20 mJy at 250, 350, 500 micron, respectively. We then measured source number counts at 250, 350, and 500 micron and we also cross-correlated the catalogs with the Sloan Digital Sky Survey to investigate the redshift distribution of the nearby sources. From this cross-correlation, we select ~2000 sources with reliable fluxes and a high signal-to-noise ratio, finding an average redshift z~0.3+/-0.22. The number counts at 250, 350, and 500 micron show an increase in the slope below 200 mJy, indicating a strong evolution in number of density for galaxies at these fluxes. In general, models tend to overpredict the counts at brighter flux densities, underlying the importance of studying the Rayleigh-Jeans part of the spectral energy distribution to refine the theoretical recipes of the models. Our iterative method for source identification allowed the detection of a family of 500 micron sources that are not foreground objects belonging to Virgo and not found in other catalogs., Comment: Accepted for publication in Astronomy and Astrophysics

We map the neutral atomic gas content of M33 using high resolution VLA and GBT observations and fit a tilted ring model to determine the orientation of the extended gaseous disk and its rotation curve. The disk of M33 warps from 8 kpc outwards without substantial change of its inclination with respect to the line of sight. Rotational velocities rise steeply with radius in the inner disk, reaching 100 km/s in 4 kpc, then the rotation curve becomes more perturbed and flatter with velocities as high as 120-130 km/s out to 23 kpc. We derive the stellar mass surface density map of M33's optical disk, via pixel -SED fitting methods based on population synthesis models, which highlights variations in the mass-to-light ratio. The stellar mass surface further out is estimated from deep images of outer disk fields. Stellar and gas maps are then used in the dynamical analysis of the rotation curve to constrain the dark matter distribution which is relevant at all radii. A dark matter halo with a Navarro-Frenk-White density profile in a LCDM cosmology, provides the best fit to the rotation curve for a dark halo concentration C=10 and a total halo mass of 4.3 10^{11}Msun. This imples a baryonic fraction of order 0.02 and the evolutionary history of this galaxy should account for loss of a large fraction of its original baryonic content., Comment: 19 pages, 18 figures, Accepted for publication in A&A (replacement for corrected references to two papers)

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., Comment: 24 pages, 14 figures, accepted for publication in MNRAS

We study the dynamics of faint stellar substructures around the Umbrella Galaxy, NGC 4651, which hosts a dramatic system of streams and shells formed through the tidal disruption of a nucleated dwarf elliptical galaxy. We elucidate the basic characteristics of the system (colours, luminosities, stellar masses) using multi-band Subaru/Suprime-Cam images. The implied stellar mass-ratio of the ongoing merger event is about 1:50. We identify candidate kinematic tracers (globular clusters, planetary nebulae, H ii regions), and follow up a subset with Keck/DEIMOS spectroscopy to obtain velocities. We find that 15 of the tracers are likely associated with halo substructures, including the probable stream progenitor nucleus. These objects delineate a kinematically cold feature in position-velocity phase space. We model the stream using single test-particle orbits, plus a rescaled pre-existing N-body simulation. We infer a very eccentric orbit with a period of roughly 0.35 Gyr and turning points at approximately 2-4 and 40 kpc, implying a recent passage of the satellite through the disc, which may have provoked the visible disturbances in the host galaxy. This work confirms that the kinematics of low surface brightness substructures can be recovered and modeled using discrete tracers - a breakthrough that opens up a fresh avenue for unraveling the detailed physics of minor merging., Comment: 23 pages, 10 figures, 7 tables, accepted for publication in MNRAS 3D animated gif N-body model: http://sages.ucolick.org/downloads/Umbrella_3d.gif

The fraction of galaxies with red colours depends sensitively on environment, and on the way in which environment is measured. To distinguish competing theories for the quenching of star formation, a robust and complete description of environment is required, to be applied to a large sample of galaxies. The environment of galaxies can be described using the density field of neighbours on multiple scales - the multiscale density field. We are using the Millennium simulation and a simple HOD prescription which describes the multiscale density field of Sloan Digital Sky Survey DR7 galaxies to investigate the dependence of the fraction of red galaxies on the environment. Using a volume limited sample where we have sufficient galaxies in narrow density bins, we have more dynamic range in halo mass and density for satellite galaxies than for central galaxies. Therefore we model the red fraction of central galaxies as a constant while we use a functional form to describe the red fraction of satellites as a function of halo mass which allows us to distinguish a sharp from a gradual transition. While it is clear that the data can only be explained by a gradual transition, an analysis of the multiscale density field on different scales suggests that colour segregation within the haloes is needed to explain the results. We also rule out a sharp transition for central galaxies, within the halo mass range sampled., Comment: 24 pages, 21 figures, accepted for publication by MNRAS

We have identified a low surface brightness stellar stream from visual inspection of SDSS imaging for the edge-on, spiral galaxy NGC5387. A blue overdensity was also identified in SDSS coincident with the stream intersection with the NGC5387 disk. The overdensity was also detected in the GALEX Deep Imaging Survey and found to contribute 38% of the total FUV integrated flux from NGC5387, which suggests that the region is actively forming stars. Deeper imaging was acquired with the Vatican Advanced Technology Telescope (VATT) in the B, V, and R filters that confirmed the presence of both the stellar stream and the blue overdensity. Analysis of the VATT photometry indicates the stellar stream is red in color, B-V = 0.7, and has a stellar mass of 6x10^8 M$_{\odot}$, which implies a merger ratio of 1:50. Assessment of the stream morphology suggests that the merger event has a circular radius, R~16 kpc, the stream formed in ~400 Myr, and the progenitor had a total mass of ~2x10 M$_{\odot}$. Spectroscopy from LBT+MODS1 was used to determine that the blue overdensity is at the same redshift as NGC5387, consists of young stellar populations (~10 Myr), is metal-poor (12 + log(O/H) = 8.03), and forming stars at an enhanced rate (~1-3 M$_{\odot}$/year) given its total stellar mass (2x10^7 M$_{\odot}$). Several interpretations are posited to explain these observational data, of which the most likely are (i) that the blue overdensity is a region of enhanced star formation in the outer disk of NGC5387 induced by the minor accretion event, and (ii) that the blue overdensity is the progenitor of the stellar stream undergoing a period of enhanced star formation as a result of its interaction with NGC5387. Confirmation and theoretical exploration of these scenarios are presented in a companion paper., Comment: 16 pages, 6 figures, submitted to ApJ

We use dust scaling relations to investigate the hypothesis that Virgo cluster transition-type dwarfs are infalling star-forming field galaxies, which is argued based on their optical features (e.g. disks, spiral arms, bars) and kinematic properties similar to late-type galaxies. After their infall, environmental effects gradually transform them into early-type galaxies through the removal of their interstellar medium and quenching of all star formation activity. In this paper, we aim to verify whether this hypothesis holds using far-infrared diagnostics based on Herschel observations of the Virgo cluster taken as part of the Herschel Virgo Cluster Survey (HeViCS). We select transition-type objects in the nearest cluster, Virgo, based on spectral diagnostics indicative for their residual or ongoing star formation. We detect dust Md ~ 10^{5-6} Msun in 36% of the transition-type dwarfs located on the high end of the stellar mass distribution. This suggests that the dust reservoirs present in non-detections fall just below the Herschel detection limit (< 1.1x10^5 Msun). Dust scaling relations support the hypothesis of a transformation between infalling late-type galaxies to quiescent low-mass spheroids governed by environmental effects, with dust-to-stellar mass fractions for transition-type dwarfs in between values characteristic for late-type objects and the lower dust fractions observed in early-type galaxies. Several transition-type dwarfs demonstrate blue central cores, hinting at the radially outside-in removal of gas and quenching of star formation activity. The fact that dust is also confined to the inner regions suggests that metals are stripped in the outer regions along with the gas. In the scenario of most dust being stripped from the galaxy along with the gas, we argue that... (abridged), Comment: 18 pages, 7 figues, accepted for publication in MNRAS

Aims. We study the dust content of a large optical input sample of 910 early-type galaxies (ETG) in the Virgo cluster, extending also to the dwarf ETG, and examine the results in relation with those on the other cold ISM components. Methods. We searched for far-infrared emission in all galaxies of the input sample using the 250 micron image of the Herschel Virgo Cluster Survey (HeViCS). This image covers a large fraction of the cluster. For the detected ETG we measured fluxes in 5 bands from 100 to 500 micron, and estimated the dust mass and temperature with modified black-body fits. Results. Dust is detected above the completeness limit of 25.4 mJy at 250 micron in 46 ETG, 43 of which are in the optically complete part of the input sample. In addition dust is present at fainter levels in another 6 ETG. We detect dust in the 4 ETG with synchrotron emission, including M 87. Dust appears to be much more concentrated than stars and more luminous ETG have higher dust temperatures. Dust detection rates down to the 25.4 mJy limit are 17% for ellipticals, about 40% for lenticulars (S0 + S0a) and around 3% for dwarf ETG. Dust mass does not correlate clearly with stellar mass and is often much more than that expected for a passive galaxy in a closed-box model. The dust-to-stars mass ratio anticorrelates with galaxy luminosity, and for some dwarf ETG reaches values as high as for dusty late-type galaxies. In the Virgo cluster slow rotators appear more likely to contain dust than fast ones. Comparing the dust results with those on HI from ALFALFA, there are only 8 ETG detected both in dust and in HI in the HeViCS area; 39 have dust but only an upper limit on HI, and 8 have HI but only an upper limit on dust. The locations of these galaxies in the cluster are different, with the dusty ETG concentrated in the densest regions, while the HI rich ETG are at the periphery., Comment: Accepted by Astronomy and Astrophysics; modified to reflect the on-line forthcoming version on the A&A web site

We use high-resolution K-band VLT/HAWK-I imaging over 0.25 square degrees to study the structural evolution of massive early-type galaxies since z~1. Mass-selected samples, complete down to log(M/M_sun)~10.7 such that `typical' L* galaxies are included at all redshifts, are drawn from pre-existing photometric redshift surveys. We then separated the samples into different redshift slices and classify them as late- or early-type galaxies on the basis of their specific star-formation rate. Axis-ratio measurements for the ~400 early-type galaxies in the redshift range 0.61 are, on average, flatter than at z<1 and the median projected axis ratio at a fixed mass decreases with redshift. However, we also find that at all epochs z<~2 the very most massive early-type galaxies (log(M/M_sun)>11.3) are the roundest, with a pronounced lack among them of galaxies that are flat in projection. Merging is a plausible mechanism that can explain both results: at all epochs merging is required for early-type galaxies to grow beyond log(M/M_sun)~11.3, and all early types over time gradually and partially loose their disk-like characteristics., Comment: accepted for publication in ApJ; 11 pages; 10 figures; 4 tables

We investigate the dust-to-gas mass ratio and the environmental effects on the various components of the interstellar medium for a spatially resolved sample of Virgo spirals. We have used the IRAM-30m telescope to map over their full extent NGC 4189, NGC 4298, NGC 4388, and NGC 4299 in the 12CO(1-0) and the 12CO(2-1) lines. We observed the same lines in selected regions of NGC 4351, NGC 4294, and NGC 4424. The CO observations are combined with Herschel maps in 5 bands between 100-500 {\mu}m from the HeViCS survey, and with HI data from the VIVA survey, to obtain spatially resolved dust and gas distributions. We studied the environmental dependencies by adding to our sample eight galaxies with 12CO(1-0) maps from the literature. We estimate the integrated mass of molecular hydrogen for the galaxies observed in the CO lines. We find molecular-to-total gas mass fractions between 0.04 \leq fmol \leq 0.65, with the lowest values for the dimmest galaxy in the B-band. The integrated dust-to-gas ratio ranges between 0.011 and 0.004. For the 12 mapped galaxies we derive the radial distributions of the atomic gas, molecular gas, and dust. We also study the effect of different CO-to-H2 conversion factors. Both the molecular gas and the dust distributions show steeper radial profiles for HI-deficient galaxies and the average dust-to-gas ratio for these galaxies increases or stays radially constant. On scales of \sim 3 kpc, we find a strong correlation between the molecular gas and the 250 micron surface brightness that is tighter than average for non-deficient galaxies. The correlation becomes linear if we consider the total gas surface mass density. However, the inclusion of atomic hydrogen does not improve the statistical significance of the correlation. The environment can modify the distributions of molecules and dust within a galaxy, although these components are more tightly bound than the atomic gas., Comment: 17pp, 14 fig., accepted for publications in Astronomy and Astrophysics

We present VLT-ISAAC NIR spectro-photometric observations of 16 post-starburst (PSB) galaxies aimed at constraining the debated influence of TP-AGB stars on the SED of galaxies with stellar ages of 0.5-2 Gyr, critical for high-redshift studies. PSB galaxies have negligible on-going star formation and a SED dominated by the stars formed in a recent (<2 Gyr) burst. By spectroscopically selecting PSB galaxies with mean luminosity-weighted ages between 0.5 and 1.5 Gyr and a broad range of metallicities, we explore the parameter space around the peak of the relative energy output of TP-AGB stars. Crucially, we target galaxies at z~0.2, so that two main spectral features of TP-AGB stars (C-molecule band-head drops at 1.41 and 1.77mum), which are heavily absorbed by our atmosphere for targets at z~0, move inside the H and K atmospheric windows and can be constrained for the first time to high accuracy. Our findings provide key constraints to stellar population synthesis models: i) The NIR regions around 1.41 and 1.77mum (rest-frame) are featureless for all galaxies in our sample over the whole range of relevant ages and metallicities at variance with the Maraston (2005) TP-AGB heavy models, which exhibit marked drops there; ii) No flux boosting is observed in the NIR. The optical-NIR SEDs of most of our PSB galaxies can be consistently reproduced with Bruzual & Charlot (2003) models, using either simple stellar populations of corresponding light-weighted ages and metallicities, or a more realistic burst plus an old population containing <~60% of the total stellar mass. In contrast, all combinations of this kind based on the Maraston (2005) models are unable to simultaneously reproduce the smoothness of the NIR spectra and the relatively blue optical-NIR colours in the observations. The data collected in this study appear to disfavour TP-AGB heavy models with respect to TP-AGB light ones., Comment: MNRAS accepted. 21 pages, 11 figures. New version improves mainly on clarity of figures and discussion of burst fraction/SFH

We present new spectro-photometric NIR observations of 16 post-starburst galaxies especially designed to test for the presence of strong carbon features of thermally pulsing AGB (TP-AGB) stars, as predicted by recent models of stellar population synthesis. Selection based on clear spectroscopic optical features indicating the strong predominance of stellar populations with ages between 0.5 and 1.5 Gyr and redshift around 0.2 allows us to probe the spectral region that is most affected by the carbon features of TP-AGB stars (unaccessible from the ground for z~0 galaxies) in the evolutionary phase when their impact on the IR luminosity is maximum. Nevertheless, none of the observed galaxies display such features. Moreover the NIR fluxes relative to optical are consistent with those predicted by the original Bruzual & Charlot (2003) models, where the impact of TP-AGB stars is much lower than has been recently advocated., Comment: 3 pages, to appear in the Proceedings of the IAU Symposium No. 284, 2011, "The Spectral Energy Distribution of Galaxies", R.J. Tuffs & C.C.Popescu, eds

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf starburst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique combination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between luminosity or stellar-mass between the two galaxies is ~1:50, while the indirectly measured dynamical mass-ratio, when including dark matter, may be ~1:10-1:5. This system may thus represent a "stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and possibly in triggering their starbursts., Comment: ApJ Letters, in press. Minor changes. 6 pages, 3 figures. A high resolution version of the paper and Subaru image are available at http://www.cosmotography.com/images/small_ngc4449.html

We present new deep observations of 'shell' structures in the halo of the nearby elliptical galaxy NGC 7600, alongside a movie of galaxy formation in a cold dark matter universe (available at http://www.virgo.dur.ac.uk/shell-galaxies). The movie, based on an ab initio cosmological simulation, shows how continuous accretion of clumps of dark matter and stars creates a swath of diffuse circumgalactic structures. The disruption of a massive clump on a near-radial orbit creates a complex system of transient concentric shells which bare a striking resemblance to those of NGC 7600. With the aid of the simulation we interpret NGC 7600 in the context of the CDM model., Comment: Accepted for publication in ApJ Letters. 6 pages, 3 figures. For the movies discussed in the paper, high resolution figures and more, see http://www.virgo.dur.ac.uk/shell-galaxies

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed., Comment: 14 pages, 9 figures, 1 table, published in MNRAS

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars., Comment: 19 pages, 12 figures; published in ApJ

A large body of evidence has demonstrated that the global rest-frame optical and IR colours of galaxies correlate well with each other, as well as with other galactic properties such as surface brightness and morphology. However the processes that lead to the observed correlations are contrary; the stellar light that contributes to the optical is readily absorbed by dust which emits in the IR. Thus on small scales we expect these correlations to break down. We examine seven nearby galaxies ranging from early- to late-types, on a pixel-by-pixel basis and we demonstrate that there is disconnect between the optical and IR when normalized to the near-IR (H-band). We can decompose this disconnect into two distinct components through a Principal Component Analysis of the H-band normalized SED of the pixels: one mainly correlated with variations in the IR, the other correlated with variations in the optical. By mapping these two components it is clear they arise from distinct spatial regions. The IR dominated component is strongly associated with the specific star-formation rate, while the optical-dominated component is broadly associated with the stellar mass density. However, when the pixels of all galaxies are compared, the well known optical-IR colour correlations return, demonstrating that the variance observed within galaxies is around a mean which follows the well-known trend. We also examine the extremely strong correlations between the IRAC-NIR colours and demonstrate that they are tight enough to use a single IRAC-NIR colour (i.e. 8mum-H) to determine the fluxes in the other IRAC bands. These correlations arise from the differing contribution of stellar light and dust to the IRAC bands, enabling us to determine pure stellar colours for these bands, but still demonstrating the need for dust (or stellar) corrections in these bands when being used as stellar (dust) tracers., Comment: 20 pages plus 25 pages appendix. Resubmitted to MNRAS after referee's comments