Source author record

Stefano Zibetti

Stefano Zibetti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

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.

preprint2022arXiv

Stellar mass as the "glocal" driver of galaxies' stellar population properties

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 <mu>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, <mu>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.

preprint2021arXiv

Galaxy evolution across environments as probed by the ages, stellar metallicities and [alpha/Fe] of central and satellite galaxies

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.

preprint2020arXiv

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

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δ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $Hδ_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ Å/$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_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.

preprint2020arXiv

The Fundamental Plane of Massive Quiescent Galaxies at z~2

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.5<z<2.5$ to date. The FP ($r_{e}, σ_{e}, I_{e}$) is established using $19$ $UVJ$ quiescent galaxies from COSMOS with $Hubble$ $Space$ $Telescope$ $(HST)$ $H_{F160W}$ rest-frame optical sizes and X-shooter absorption line measured stellar velocity dispersions. For a very massive, ${\rm{log}}(M_{\ast}/M_{\odot})>11.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.

preprint2019arXiv

X-shooter Spectroscopy and HST Imaging of 15 Ultra Massive Quiescent Galaxies at $z\gtrsim2$

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.

preprint2016arXiv

Asymmetric mass models of disk galaxies - I. Messier 99

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.

preprint2016arXiv

Variations of the initial mass function in semi-analytical models: implications for the mass assembly and the chemical enrichment of galaxies in the GAEA model

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 [$α$/Fe]-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of $α$-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 $α$-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 [$α$/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.

preprint2015arXiv

The Science Case for Multi-Object Spectroscopy on the European ELT

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.

preprint2014arXiv

Dynamical signatures of a LCDM-halo and the distribution of the baryons in M33

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.

preprint2014arXiv

Kinematics and simulations of the stellar stream in the halo of the Umbrella Galaxy

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.

preprint2014arXiv

The Herschel Virgo Cluster Survey. XVII. SPIRE point-source catalogs and number counts

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.

preprint2014arXiv

The influence of the environmental history on quenching star formation in a $Λ$CDM universe

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.

preprint2013arXiv

Cannibalization and Rebirth in the NGC5387 System. I. The Stellar Stream & Star Forming Region

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.

preprint2013arXiv

More than just halo mass: Modelling how the red galaxy fraction depends on multiscale density in a HOD framework

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.

preprint2013arXiv

The Herschel Virgo Cluster Survey - XIII. Dust in early-type galaxies

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.

preprint2013arXiv

The Herschel Virgo Cluster Survey XIV: transition-type dwarf galaxies in the Virgo cluster

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)

preprint2012arXiv

Direct Constraints on the Impact of TP-AGB Stars on the SED of Galaxies from Near-Infrared Spectroscopy

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.

preprint2012arXiv

Dwarfs Gobbling Dwarfs: A Stellar Tidal Stream Around NGC 4449 and Hierarchical Galaxy Formation on Small Scales

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.

preprint2012arXiv

Measures of Galaxy Environment - I. What is "Environment"?

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.

preprint2012arXiv

Near Infrared spectroscopy of post-starburst galaxies: a limited impact of TP-AGB stars on galaxy SEDs

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.

preprint2012arXiv

Shape Evolution of Massive Early-Type Galaxies: Confirmation of Increased Disk Prevalence at z>1

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.6<z<1.8 are accurate to 0.1 or better. The projected axis-ratio distributions are then compared with lower redshift samples. We find strong evidence for evolution of the population properties: early-type galaxies at z>1 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.

preprint2012arXiv

The Herschel Virgo Cluster Survey. XI. Environmental effects on molecular gas and dust in spiral disks

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 μ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.

preprint2011arXiv

Resolved optical-infrared SEDs of galaxies: universal relations and their break-down on local scales

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.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

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.

preprint2011arXiv

The formation of shell galaxies similar to NGC 7600 in the cold dark matter cosmogony

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.

preprint2010arXiv

A multiscale approach to environment

Physical processes influencing the properties of galaxies can be traced by the dependence and evolution of galaxy properties on their environment. A detailed understanding of this dependence can only be gained through comparison of observations with models, with an appropriate quantification of the rich parameter space describing the environment of the galaxy. We present a new, multiscale parameterization of galaxy environment which retains an observationally motivated simplicity whilst utilizing the information present on different scales. We examine how the distribution of galaxy (u-r) colours in the Sloan Digital Sky Survey (SDSS), parameterized using a double gaussian (red plus blue peak) fit, depends upon multiscale density. This allows us to probe the detailed dependence of galaxy properties on environment in a way which is independent of the halo model. Nonetheless, cross-correlation with the group catalogue constructed by Yang et al, 2007 shows that galaxy properties trace environment on different scales in a way which mimics that expected within the halo model. This provides independent support for the existence of virialized haloes, and important additional clues to the role played by environment in the evolution of the galaxy population. This work is described in full by Wilman et al., 2010, MNRAS, accepted

preprint2010arXiv

A multiscale approach to environment and its influence on the colour distribution of galaxies

We present a multiscale approach to measurements of galaxy density, applied to a volume-limited sample constructed from SDSS DR5. We populate a rich parameter space by obtaining independent measurements of density on different scales for each galaxy, avoiding the implicit assumptions involved, e.g., in the construction of group catalogues. As the first application of this method, we study how the bimodality in galaxy colour distribution (u-r) depends on multiscale density. The u-r galaxy colour distribution is described as the sum of two gaussians (red and blue) with five parameters: the fraction of red galaxies (f_r) and the position and width of the red and blue peaks (mu_r, mu_b, sigma_r and sigma_b). Galaxies mostly react to their smallest scale (< 0.5 Mpc) environments: in denser environments red galaxies are more common (larger f_r), redder (larger mu_r) and with a narrower distribution (smaller sigma_r), while blue galaxies are redder (larger mu_b) but with a broader distribution (larger sigma_b). There are residual correlations of f_r and mu_b with 0.5 - 1 Mpc scale density, which imply that total or partial truncation of star formation can relate to a galaxy's environment on these scales. Beyond 1 Mpc (0.5 Mpc for mu_r) there are no positive correlations with density. However f_r (mu_r) anti-correlates with density on >2 (1) Mpc scales at fixed density on smaller scales. We examine these trends qualitatively in the context of the halo model, utilizing the properties of haloes within which the galaxies are embedded, derived by Yang et al, 2007 and applied to a group catalogue. This yields an excellent description of the trends with multiscale density, including the anti-correlations on large scales, which map the region of accretion onto massive haloes. Thus we conclude that galaxies become red only once they have been accreted onto haloes of a certain mass.

preprint2010arXiv

A simple connection between the near- and mid-infrared emission of galaxies and their star-formation rates

We have measured the near-infrared colors and the fluxes of individual pixels in 68 galaxies common to the Spitzer Infrared Nearby Galaxies Survey and the Large Galaxy Atlas Survey. Each galaxy was separated into regions of increasingly red near-infrared colors. In the absence of dust extinction and other non-stellar emission, stellar populations are shown to have relatively constant NIR colors, independent of age. In regions of high star formation, the average intensity of pixels in red-excess regions (at 1.25, 3.6, 4.5, 5.6, 8.0 and 24 micron) scales linearly with the intrinsic intensity of Halpha emission, and thus with the star-formation rate within the pixel. This suggests that most NIR-excess regions are not red because their light is being depleted by absorption. Instead, they are red because additional infrared light is being contributed by a process linked to star-formation. This is surprising because the shorter wavelength bands in our study (1.25 micron-5.6 micron) do not probe emission from cold (10-20 K) and warm (50-100 K) dust associated with star-formation in molecular clouds. However, emission from hot dust (700-1000 K) and/or Polycyclic Aromatic Hydrocarbon molecules can explain the additional emission seen at the shorter wavelengths in our study. The contribution from hot dust and/or PAH emission at 2-5micron and PAH emission at 5.6 and 8.0 micron scales linearly with warm dust emission at 24 micron and the intrinsic Halpha emission. Since both are tied to the star-formation rate, our analysis shows that the NIR excess continuum emission and PAH emission at ~1-8 micron can be added to spectral energy distribution models in a very straight-forward way, by simply adding an additional component to the models that scales linearly with star-formation rate.

preprint2010arXiv

Stellar Tidal Streams in Spiral Galaxies of the Local Volume: A Pilot Survey with Modest Aperture Telescopes

[Abridged] Within the hierarchical framework for galaxy formation, minor merging and tidal interactions are expected to shape all large galaxies to the present day. As a consequence, most seemingly normal disk galaxies should be surrounded by spatially extended stellar 'tidal features' of low surface brightness. As part of a pilot survey for such interaction signatures, we have carried out ultra deep, wide field imaging of 8 isolated spiral galaxies in the Local Volume, with data taken at small (D=0.1-0.5m) robotic telescopes that provide exquisite surface brightness sensitivity (mu_V)~28.5$ mag/arcsec^2). This initial observational effort has led to the discovery of six previously undetected extensive (to ~30 kpc) stellar structures in the halos surrounding these galaxies, likely debris from tidally disrupted satellites. In addition, we confirm and clarify several enormous stellar over-densities previously reported in the literature, but never before interpreted as tidal streams. Even this pilot sample of galaxies exhibits strikingly diverse morphological characteristics of these extended stellar features: great circle-like features that resemble the Sagittarius stream surrounding the Milky Way, remote shells and giant clouds of presumed tidal debris far beyond the main stelar body, as well as jet-like features emerging from galactic disks. A qualitative comparison with available simulations set in a Lambda-Cold Dark Matter cosmology shows that the extraordinary variety of stellar morphologies detected in this pilot survey matches that seen in those simulations. The common existence of these tidal features around 'normal' disk galaxies and the morphological match to the simulations constitutes new evidence that these theoretical models also apply to a large number of other Milky Way-mass disk galaxies in the Local Volume.

preprint2009arXiv

Probing star formation across cosmic time with absorption line systems

We present an empirical connection between cold gas in galactic halos and star formation. Using a sample of more than 8,500 MgII absorbers from SDSS quasar spectra, we report the detection of a 15 sigma correlation between the rest equivalent width W0 of MgII absorbers and the associated OII luminosity, an estimator of star formation rate. This correlation has interesting implications: using only observable quantities we show that MgII absorbers trace a substantial fraction of the global OII luminosity density and recover the overall star formation history of the Universe derived from classical emission estimators up to z~2. We then show that the distribution function of MgII rest equivalent widths, dN/dW0 inherits both its shape and amplitude from the OII luminosity function Phi(L). These distributions can be naturally connected, without any free parameter. Our results imply a high covering factor of cold gas around star forming galaxies: C>0.5, favoring outflows as the mechanism responsible for MgII absorption. We then argue that intervening MgII absorbers and blue-shifted MgII absorption seen in the spectra of star forming galaxies are essentially the same systems. These results not only shed light on the nature of MgII absorbers but also provide us with a new probe of star formation, in absorption, i.e. in a way which does not suffer from dust extinction and with a redshift-independent sensitivity. As shown in this analysis, such a tool can be applied in a noise-dominated regime, i.e. using a dataset for which emission lines are not detected in individual objects. This is of particular interest for high redshift studies.

preprint2009arXiv

Resolved maps of stellar mass and SED of galaxies from optical/NIR imaging and SPS models

We report on the method developed by Zibetti, Charlot and Rix (2009) to construct resolved stellar mass maps of galaxies from optical and NIR imaging. Accurate pixel-by-pixel colour information (specifically g-i and i-H) is converted into stellar mass-to-light ratios with typical accuracy of 30%, based on median likelihoods derived from a Monte Carlo library of 50,000 stellar population synthesis models that include dust and updated TP-AGB phase prescriptions. Hence, surface mass densities are computed. In a pilot study, we analyze 9 galaxies spanning a broad range of morphologies. Among the main results, we find that: i) galaxies appear much smoother in stellar mass maps than at any optical or NIR wavelength; ii) total stellar mass estimates based on unresolved photometry are biased low with respect to the integral of resolved stellar mass maps, by up to 40%, due to dust obscured regions being under-represented in global colours; iii) within a galaxy, on local scales colours correlate with surface stellar mass density; iv) the slope and tightness of this correlation reflect/depend on the morphology of the galaxy.

preprint2009arXiv

Resolved stellar mass maps of galaxies. I: method and implications for global mass estimates

(Abridged) We introduce a novel technique to construct spatially resolved maps of stellar mass surface density in galaxies based on optical and near IR imaging. We use optical/NIR colour(s) to infer effective stellar mass-to-light ratios (M/L) at each pixel, which are then multiplied by the surface brightness to obtain the local stellar surface mass density. We build look-up tables to express M/L as a function of colour(s) by marginalizing over a Monte Carlo library of 50,000 stellar population synthesis (SPS) models by Charlot & Bruzual (2007), which include a revised prescription for the TP-AGB stellar evolutionary phase, with a wide range of dust exinctions. In order to extract reliable flux and colour information at any position in the galaxy, we perform a median adaptive smoothing of the images that preserves the highest possible spatial resolution. As the most practical and robust, and hence fiducial method, we express the M/L in the H band as a function of (g-i) and (i-H). Stellar mass maps computed in this way have a typical accuracy of 30 per cent or less at any given pixel, determined from the scatter in the models. We compare maps obtained with our fiducial method with those derived using other combinations of bandpasses and the old BC03 TP-AGB prescription. Finally, we compare total stellar mass estimates obtained by integrating resolved mass maps with those obtained with unresolved photometry. In galaxies with evident dust lanes, unresolved estimates may miss up to 40 per cent of the total stellar mass because dusty regions are strongly under-represented in the luminous fluxes.

preprint2007arXiv

Optical properties and spatial distribution of MgII absorbers from SDSS image stacking

We present a statistical analysis of the photometric properties and spatial distribution of more than 2,800 MgII absorbers with 0.37<z<1 and rest equivalent width W_0(\lambda2796)>0.8Ådetected in SDSS quasar spectra. Using an improved image stacking technique, we measure the cross-correlation between MgII gas and light (in the g, r, i and z-bands) from 10 to 200 kpc and infer the light-weighted impact parameter distribution of MgII absorbers. Such a quantity is well described by a power-law with an index that strongly depends on W_0, ranging from ~-1 for W_0<~1Åto ~-2 for W_0>~ 1.5Å. At redshift 0.37<z<0.55, we find the average luminosity enclosed within 100 kpc around MgII absorbers to be M_g=-20.65+-0.11 mag, which is ~0.5 L_g*. The global luminosity-weighted colors are typical of present-day intermediate type galaxies. However, while the light of weaker absorbers originates mostly from red passive galaxies, stronger systems display the colors of blue star-forming galaxies. Based on these observations, we argue that the origin of strong MgII absorber systems might be better explained by models of metal-enriched gas outflows from star-forming/bursting galaxies. Our analysis does not show any redshift dependence for both impact parameter and rest-frame colors up to z=1. However, we do observe a brightening of the absorbers related light at high redshift (~50% from z~0.4 to 1). We argue that MgII absorbers are a phenomenon typical of a given evolutionary phase that more massive galaxies experience earlier than less massive ones, in a downsizing fashion. (abridged)