Source author record

Stephane Courteau

Stephane Courteau 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

26works
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

26 published item(s)

preprint2020arXiv

Abundance matching tested on small scales with galaxy dynamics

We present a comprehensive test of the relation between stellar and total mass in galaxies as predicted by popular models based on abundance matching (AM) techniques. We use the ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey with photometric and dynamical profiles for 190 Virgo cluster galaxies to establish a relation between the stellar and dynamical masses measured within the isophotal radius $r_{23.5}$. Various dark matter and galaxy scaling relations are combined with results from the NIHAO suite of hydrodynamical simulations to recast AM predictions in terms of these observed quantities. Our results are quite insensitive to the exact choice of dark matter profile and halo response to baryon collapse. We find that theoretical models reproduce the slope and normalization of the observed stellar-to-halo mass relation (SHMR) over more than three orders of magnitude in stellar mass $(10^8 < M_*/M_{\odot} < 2 \times 10^{11})$. However, the scatter of the observed SHMR exceeds that of AM predictions by a factor $\sim$5. For systems with stellar masses exceeding $5 \times 10^{10}~M_{\odot}$, AM overpredicts the observed stellar masses for a given dynamical mass. The latter offset may support previous indications of a different stellar initial mass function in these massive galaxies. Overall our results support the validity of AM predictions on a wide dynamical range.

preprint2020arXiv

The Red Supergiant Content of M31 and M33

We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of log L/Lo=4.0. Our archival survey data cover 5 deg^2 of M31, and 3 deg^2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in the color-magnitude diagram. The photometry is used to derive effective temperatures and bolometric luminosities via MARCS stellar atmosphere models. The resulting H-R diagrams show superb agreement with the evolutionary tracks of the Geneva evolutionary group. Our census includes 6400 RSGs in M31 and 2850 RSGs in M33 within their Holmberg radii; by contrast, only a few hundred RSGs are known so far in the Milky Way. Our catalog serves as the basis for a study of the RSG binary frequency being published separately, as well as future studies relating to the evolution of massive stars. Here we use the matches between the NIR-selected RSGs and their optical counterparts to show that the apparent similarity in the reddening of OB stars in M31 and M33 is the result of Malmquist bias; the average extinction in M31 is likely higher than that of M33. As expected, the distribution of RSGs follows that of the spiral arms, while the much older AGB population is more uniformly spread across each galaxy's disk.

preprint2016arXiv

The Evolution of Star Formation Activity in Cluster Galaxies Over $0.15<z<1.5$

We explore 7.5 billion years of evolution in the star formation activity of massive ($M_{\star}>10^{10.1}\,M_{\odot}$) cluster galaxies using a sample of 25 clusters over $0.15<z<1$ from the Cluster Lensing And Supernova survey with Hubble and 11 clusters over $1<z<1.5$ from the IRAC Shallow Cluster Survey. Galaxy morphologies are determined visually using high-resolution Hubble Space Telescope images. Using the spectral energy distribution fitting code CIGALE, we measure star formation rates, stellar masses, and 4000 Å break strengths. The latter are used to separate quiescent and star-forming galaxies (SFGs). From $z\sim1.3$ to $z\sim0.2$, the specific star formation rate (sSFR) of cluster SFGs and quiescent galaxies decreases by factors of three and four, respectively. Over the same redshift range, the sSFR of the entire cluster population declines by a factor of 11, from $0.48\pm0.06\;\mathrm{Gyr}^{-1}$ to $0.043\pm0.009\;\mathrm{Gyr}^{-1}$. This strong overall sSFR evolution is driven by the growth of the quiescent population over time; the fraction of quiescent cluster galaxies increases from $28^{+8}_{-19}\%$ to $88^{+5}_{-4}\%$ over $z\sim1.3\rightarrow0.2$. The majority of the growth occurs at $z\gtrsim0.9$, where the quiescent fraction increases by 0.41. While the sSFR of the majority of star-forming cluster galaxies is at the level of the field, a small subset of cluster SFGs have low field-relative star formation activity, suggestive of long-timescale quenching. The large increase in the fraction of quiescent galaxies above $z\sim0.9$, coupled with the field-level sSFRs of cluster SFGs, suggests that higher redshift cluster galaxies are likely being quenched quickly. Assessing those timescales will require more accurate stellar population ages and star formation histories.

preprint2016arXiv

The Next Generation Virgo Cluster Survey (NGVS). XIII. The Luminosity and Mass Function of Galaxies in the Core of the Virgo Cluster and the Contribution from Disrupted Satellites

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $α=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $α= -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $α=-1.21\pm0.05$.

preprint2015arXiv

Bayesian Investigation of Isochrone Consistency Using the Old Open Cluster NGC 188

This paper provides a detailed comparison of the differences in parameters derived for a star cluster from its color-magnitude diagrams depending on the filters and models used. We examine the consistency and reliability of fitting three widely-used stellar evolution models to fifteen combinations of optical and near-IR photometry for the old open cluster NGC 188. The optical filter response curves match those of the theoretical systems and are thus not the source of fit inconsistencies. NGC 188 is ideally suited to the present study thanks to a wide variety of high-quality photometry and available proper motions and radial velocities which enable us to remove non-cluster members and many binaries. Our Bayesian fitting technique yields inferred values of age, metallicity, distance modulus, and absorption as a function of the photometric band combinations and stellar models. We show that the historically-favored three band combinations of UBV and VRI can be meaningfully inconsistent with each other and with longer baseline datasets such as UBVRIJHKs. Differences among model sets can also be substantial. For instance, fitting Yi et al. (2001) and Dotter et al. (2008) models to UBVRIJHKs photometry for NGC 188 yields the following cluster parameters: age={5.78+ 0.03, 6.45+-0.04} Gyr, [Fe/H]={+0.125+-0.003, -0.077+-0.003} dex, m-M={11.441+-0.007, 11.525+-0.005} mag, and Av={0.162+-0.003, 0.236+-0.003} mag, respectively. Within the formal fitting errors, these two fits are substantially and statistically different. Such differences amongst fits using different filters and models are a cautionary tale regarding our current ability to fit star cluster color-magnitude diagrams. Additional modeling of this kind, with more models and star clusters, and future GAIA parallaxes are critical for isolating and quantifying the most relevant uncertainties in stellar evolutionary models.

preprint2015arXiv

On the Global Mass Distribution in Disk Galaxies

We present a summary of recent assessments of the mass distribution in disk galaxies. Of issue in order to characterize galaxy formation models is to determine the relative fraction of baryons and dark matter at all radii in galaxies. For disk galaxies, various measurements of the mass distribution in galaxies based on vertical kinematics, strong lensing, residuals of scaling relations, fluid dynamical modeling, bar strength and pattern speed, warps, and others, have called for either a maximal or sub-maximal contribution of the baryons in the inner parts of the disk. We propose a global picture whereby all galaxies are typically baryon-dominated (maximal) at the center and dark-matter dominated (sub-maximal) in their outskirts. Using as a fiducial radius the peak of the rotation curve of a pure baryonic exponential disk at R_2.2=2.2 Rd, where Rd is the disk scale length, the transition from maximal to sub-maximal baryons occurs within or near R_2.2 for low-mass disk galaxies (with V_tot < 200 km/s) and beyond R_2.2 for more massive systems. The mass fraction's inverse dependence on circular velocity in disk galaxies is largely insensitive to the presence of a bar at R_2.2. The mean mass fractions of late- and early-type galaxies are shown to be different at the same fiducial radius and circular velocity, pointing to different galaxy formation mechanisms. Feedback, dynamical friction, size, and IMF variations are key ingredients for understanding these differences.

preprint2015arXiv

On the Uncertainties of Stellar Mass Estimates via Colour Measurements

Mass-to-light versus colour relations (MLCRs), derived from stellar population synthesis models, are widely used to estimate galaxy stellar masses (M$_*$) yet a detailed investigation of their inherent biases and limitations is still lacking. We quantify several potential sources of uncertainty, using optical and near-infrared (NIR) photometry for a representative sample of nearby galaxies from the Virgo cluster. Our method for combining multi-band photometry with MLCRs yields robust stellar masses, while errors in M$_*$ decrease as more bands are simultaneously considered. The prior assumptions in one's stellar population modelling dominate the error budget, creating a colour-dependent bias of up to 0.6 dex if NIR fluxes are used (0.3 dex otherwise). This matches the systematic errors associated with the method of spectral energy distribution (SED) fitting, indicating that MLCRs do not suffer from much additional bias. Moreover, MLCRs and SED fitting yield similar degrees of random error ($\sim$0.1-0.14 dex) when applied to mock galaxies and, on average, equivalent masses for real galaxies with M$_* \sim$ 10$^{8-11}$ M$_{\odot}$. The use of integrated photometry introduces additional uncertainty in M$_*$ measurements, at the level of 0.05-0.07 dex. We argue that using MLCRs, instead of time-consuming SED fits, is justified in cases with complex model parameter spaces (involving, for instance, multi-parameter star formation histories) and/or for large datasets. Spatially-resolved methods for measuring M$_*$ should be applied for small sample sizes and/or when accuracies less than 0.1 dex are required. An Appendix provides our MLCR transformations for ten colour permutations of the $grizH$ filter set.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. IX. Estimating the Efficiency of Galaxy Formation on the Lowest-Mass Scales

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$.

preprint2014arXiv

Revised Mass-to-Light Ratios For Nearby Galaxy Groups and Clusters

We present a detailed investigation of the cluster stellar mass-to-light (M*/L) ratio and cumulative stellar masses, derived on a galaxy-by-galaxy basis, for 12 massive (M500 ~ 10^14 - 10^15 Msun), nearby clusters with available optical imaging data from the Sloan Digital Sky Survey Data Release 10 and X-ray data from the Chandra X-ray Observatory. Our method involves a statistical cluster membership using both photometric and spectroscopic redshifts when available to maximize completeness whilst minimizing contamination effects. We show that different methods of estimating the stellar mass-to-light ratio from observed photometry result in systematic discrepancies in the total stellar masses and average mass-to-light ratios of cluster galaxies. Nonetheless, all conversion methodologies point to a lack of correlation between M*/Li and total cluster mass, even though low-mass groups contain relatively more blue galaxies. We also find no statistically significant correlation between M*/Li and the fraction of blue galaxies. For the mass range covered by our sample, the assumption of a Chabrier IMF yields an integrated M*/Li = 1.7 +/- 0.2 Msun/Lsun, a lower value than used in most similar studies, though consistent with the study of low-mass galaxy groups by Leauthaud et al. (2012). A light (diet) Salpeter IMF would imply a ~60% increase in M*/Li.

preprint2014arXiv

The Stellar Mass of M31 as inferred by the Andromeda Optical & Infrared Disk Survey

Our proximity and external vantage point make M31 an ideal testbed for understanding the structure of spiral galaxies. The Andromeda Optical and Infrared Disk Survey (ANDROIDS) has mapped M31's bulge and disk out to R=40 kpc in $ugriJK_s$ bands with CFHT using a careful sky calibration. We use Bayesian modelling of the optical-infrared spectral energy distribution (SED) to estimate profiles of M31's stellar populations and mass along the major axis. This analysis provides evidence for inside-out disk formation and a declining metallicity gradient. M31's $i$-band mass-to-light ratio ($M/L_i^*$) decreases from 0.5 dex in the bulge to $\sim 0.2$ dex at 40 kpc. The best-constrained stellar population models use the full $ugriJK_s$ SED but are also consistent with optical-only fits. Therefore, while NIR data can be successfully modelled with modern stellar population synthesis, NIR data do not provide additional constraints in this application. Fits to the $gi$-SED alone yield $M/L_i^*$ that are systematically lower than the full SED fit by 0.1 dex. This is still smaller than the 0.3 dex scatter amongst different relations for $M/L_i$ via $g-i$ colour found in the literature. We advocate a stellar mass of $M_*(30\mathrm{kpc})=10.3^{+2.3}_{-1.7}\times 10^{10}\mathrm{M}_\odot$ for the M31 bulge and disk.

preprint2014arXiv

Validation of optimised population synthesis through mock spectra and Galactic globular clusters

Optimised population synthesis provides an empirical method to extract the relative mix of stellar evolutionary stages and the distribution of atmospheric parameters within unresolved stellar systems, yet a robust validation of this method is still lacking. We here provide a calibration of population synthesis via non-linear bound-constrained optimisation of stellar populations based upon optical spectra of mock stellar systems and observed Galactic Globular Clusters (GGCs). The MILES stellar library is used as a basis for mock spectra as well as templates for the synthesis of deep GGC spectra from Schiavon et al. (2005). Optimised population synthesis applied to mock spectra recovers mean light-weighted stellar atmospheric parameters to within a mean uncertainty of 240 K, 0.04 dex, and 0.03 dex for T_eff, log(g), and [Fe/H], respectively. Decompositions of both mock and GGC spectra confirm the method's ability to recover the expected mean light-weighted metallicity in dust-free conditions (E[B-V] < 0.15) with uncertainties comparable to evolutionary population synthesis methods. Dustier conditions require either appropriate dust-modelling when fitting to the full spectrum, or fitting only to select spectral features. We derive light-weighted fractions of stellar evolutionary stages from our population synthesis fits to GGCs, yielding on average a combined 25+/-6 per cent from main sequence and turnoff dwarfs, 64+/-7 per cent from subgiant, red giant and asymptotic giant branch stars, and 15+/-7 per cent from horizontal branch stars and blue stragglers. Excellent agreement is found between these fractions and those estimated from deep HST/ACS CMDs. Overall, optimised population synthesis remains a powerful tool for understanding the stellar populations within the integrated light of galaxies and globular clusters.

preprint2013arXiv

Constraining stellar population models - I. Age, metallicity, and abundance pattern compilation for Galactic globular clusters

We present an extenstive literature compilation of age, metallicity, and chemical abundance pattern information for the 41 Galactic globular clusters (GGCs) studied by Schiavon et al. (2005). Our compilation constitutes a notable improvement over previous similar work, particularly in terms of chemical abundances. Its primary purpose is to enable detailed evaluations of and refinements to stellar population synthesis models designed to recover the above information for unresolved stellar systems based on their integrated spectra. However, since the Schiavon sample spans a wide range of the known GGC parameter space, our compilation may also benefit investigations related to a variety of astrophysical endeavours, such as the early formation of the Milky Way, the chemical evolution of GGCs, and stellar evolution and nucleosynthesis. For instance, we confirm with our compiled data that the GGC system has a bimodal metallicity distribution and is uniformly enhanced in the alpha-elements. When paired with the ages of our clusters, we find evidence that supports a scenario whereby the Milky Way obtained its globular clusters through two channels, in situ formation and accretion of satellite galaxies. The distributions of C, N, O, and Na abundances and the dispersions thereof per cluster corroborate the known fact that all GGCs studied so far with respect to multiple stellar populations have been found to harbour them. Finally, using data on individual stars, we also confirm that the atmospheres of stars become progressively polluted by CN(O)-processed material after they leave the main sequence and uncover evidence which suggests the alpha-elements Mg and Ca may originate from more than one nucleosynthetic production site. [abridged]

preprint2013arXiv

Star Clusters in M31. V. Evidence for Self-Enrichment in Old M31 Clusters from Integrated Spectroscopy

In the past decade, the notion that globular clusters (GCs) are composed of coeval stars with homogeneous initial chemical compositions has been challenged by growing evidence that they host an intricate stellar population mix, likely indicative of a complex history of star formation and chemical enrichment. Several models have been proposed to explain the existence of multiple stellar populations in GCs, but no single model provides a fully satisfactory match to existing data. Correlations between chemistry and global parameters such as cluster mass or luminosity are fundamental clues to the physics of GC formation. In this Letter, we present an analysis of the mean abundances of Fe, Mg, C, N, and Ca for 72 old GCs from the Andromeda galaxy. We show for the first time that there is a correlation between the masses of GCs and the mean stellar abundances of nitrogen, spanning almost two decades in mass. This result sheds new light on the formation of GCs, providing important constraints on their internal chemical evolution and mass loss history.

preprint2012arXiv

Stellar Populations and Radial Migrations in Virgo Disk Galaxies

We present stellar age profiles for 64 Virgo cluster disk galaxies whose analysis poses a challenge for current galaxy formation models. Our results can be summarized as follows: first, and contrary to observations of field galaxies, these cluster galaxies are distributed almost equally amongst the three main types of disk galaxy luminosity profiles (I/II/III), indicating that the formation and/or survival of Type II breaks is suppressed within the cluster environment. Second, we find examples of statistically-significant inversions ("U-shapes") in the age profiles of all three disk galaxy types, reminescent of predictions from high-resolution simulations of classically-truncated Type II disks in the field. These features characterize the age profiles for only about a third (<36%) of each disk galaxy type in our sample. An even smaller fraction of cluster disks (~11% of the total sample), exhibit age profiles which decrease outwards (i.e., negative age gradients). Instead, flat and/or positive age gradients prevail (>50%) within our Type I, II and III sub-samples. These observations thus suggest that while stellar migrations and inside-out growth can play a significant role in the evolution of all disk galaxy types, other factors contributing to the evolution of galaxies can overwhelm the predicted signatures of these processes. We interpret our observations through a scenario whereby Virgo cluster disk galaxies formed initially like their bretheren in the field but which, upon falling into the cluster, were transformed into their present state through external processes linked to the environment. Current disk galaxy formation models fail to reproduce these results, thus calling for adequate hydrodynamical simulations of dense galaxy environments, for which the current paper provides many constraints. [Abridged]

preprint2011arXiv

A Survey of 286 Virgo Cluster Galaxies at Optical griz and Near-IR H-band: Surface Brightness Profiles and Bulge-Disk Decompositions

We present griz and H-band surface brightness profiles and bulge-disk decompositions for a morphologically-broad sample of 286 Virgo cluster galaxies. The H-band data come from a variety of sources including our survey of 171 galaxies at the UH 2.2-m, CFHT and UKIRT telescopes, and another 115 galaxies from the Two-Micron All-Sky Survey and GOLDMine archives. The optical data for all 286 Virgo galaxies were extracted from Sloan Digital Sky Survey (SDSS) images. The H-band and SDSS griz data were analyzed in a homogeneous manner using our own software, yielding a consistent set of deep, multi-band surface brightness profiles for each galaxy. Average surface brightness profiles per morphological bin were created in order to characterize the variety of galaxy light profiles across the Hubble sequence. The 1D bulge-disk decomposition parameters, as well as non-parametric galaxy measures, such as effective radius, effective surface brightness and light concentration, are presented for all 286 galaxies in each of the five optical/near-infrared wavebands. The profile decompositions account for bulge and disk components, spiral arms, nucleus and atmospheric blurring. The Virgo spiral galaxy bulges typically have a Sersic index n ~ 1, while elliptical galaxies prefer n ~ 2. No galaxy spheroid requires n > 3. The light profiles for 70% of the Virgo elliptical galaxies reveal the presence of both a spheroid and disk component. A more in-depth discussion of the structural parameter trends can be found in McDonald et al (2009b). The data provided here should serve as a base for studies of galaxy structure and stellar populations in the cluster environment. The galaxy light profiles and bulge-disk decomposition results are available at the Centre de Donnees astronomiques de Strasbourg (CDS; http://cds.u-strasbg.fr) and the author's own website (http://www.astro.queensu.ca/virgo).

preprint2011arXiv

An Investigation of Sloan Digital Sky Survey Imaging Data and Multi-Band Scaling Relations of Spiral Galaxies (with Dynamical Information)

We have compiled a sample of 3041 spiral galaxies with multi-band gri imaging from the Sloan Digital Sky Survey (SDSS) Data Release 7 and available galaxy rotational velocities derived from HI line widths. We compare the data products provided through the SDSS imaging pipeline with our own photometry of the SDSS images, and use the velocities (V) as an independent metric to determine ideal galaxy sizes (R) and luminosities (L). Our radial and luminosity parameters improve upon the SDSS DR7 Petrosian radii and luminosities through the use of isophotal fits to the galaxy images. This improvement is gauged via VL and RV relations whose respective scatters are reduced by ~8% and ~30% compared to similar relations built with SDSS parameters. The tightest VRL relations are obtained with the i-band radius, R235i, measured at 23.5 mag/arcsec^-2, and the luminosity L235i, measured within R235i. Our VRL scaling relations compare well, both in scatter and slope, with similar studies (such comparisons however depend sensitively on the nature and size of the compared samples). The typical slopes, b, and observed scatters, sigma, of the i-band VL, RL and RV relations are bVL=0.27+/-0.01, bRL=0.41+/-0.01, bRV=1.52+/-0.07, and sigmaVL=0.074, sigmaRL=0.071, sigmaRV=0.154 dex. Similar results for the SDSS g and r bands are also provided. Smaller scatters may be achieved for more pruned samples. We also compute scaling relations in terms of the baryonic mass (stars + gas), Mbar, ranging from 10^8.7 Msol to 10^11.6 Msol. Our baryonic velocity-mass (VM) relation has slope 0.29+/-0.01 and a measured scatter sigma_meas = 0.076 dex. While the observed VL and VM relations have comparable scatter, the stellar and baryonic VM relations may be intrinsically tighter, and thus potentially more fundamental, than other VL relations of spiral galaxies.

preprint2011arXiv

Dark halo response and the stellar initial mass function in early-type and late-type galaxies

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

preprint2011arXiv

Star Clusters in M 31. IV. A Comparative Analysis of Absorption Line Indices in Old M 31 and Milky Way Clusters

We present absorption line indices measured in the integrated spectra of globular clusters both from the Galaxy and from M 31. Our samples include 41 Galactic globular clusters, and more than 300 clusters in M 31. The conversion of instrumental equivalent widths into the Lick system is described, and zero-point uncertainties are provided. Comparison of line indices of old M 31 clusters and Galactic globular clusters suggests an absence of important differences in chemical composition between the two cluster systems. In particular, CN indices in the spectra of M 31 and Galactic clusters are essentially consistent with each other, in disagreement with several previous works. We reanalyze some of the previous data, and conclude that reported CN differences between M 31 and Galactic clusters were mostly due to data calibration uncertainties. Our data support the conclusion that the chemical compositions of Milky Way and M 31 globular clusters are not substantially different, and that there is no need to resort to enhanced nitrogen abundances to account for the optical spectra of M 31 globular clusters.

preprint2011arXiv

The Formation and Evolution of Virgo Cluster Galaxies - I. Broadband Optical & Infrared Colours

We use a combination of deep optical (gri) and near-infrared (H) photometry to study the radially-resolved colours of a broad sample of 300 Virgo cluster galaxies. For most galaxy types, we find that the median g-H colour gradient is either flat (gas-poor giants and gas-rich dwarfs) or negative (i.e., colours become bluer with increasing radius; gas-poor dwarfs, spirals, and gas-poor peculiars). Later-type galaxies typically exhibit more negative gradients than early-types. Given the lack of a correlation between the central colours and axis ratios of Virgo spiral galaxies, we argue that dust likely plays a small role, if at all, in setting those colour gradients. We search for possible correlations between galaxy colour and photometric structure or environment and find that the Virgo galaxy colours become redder with increasing concentration, luminosity and surface brightness, while no dependence with cluster-centric radius or local galaxy density is detected (over a range of ~2 Mpc and ~3-16 Mpc^-2, respectively). However, the colours of gas-rich Virgo galaxies do correlate with neutral gas deficiency, such that these galaxies become redder with higher deficiencies. Comparisons with stellar population models suggest that these colour gradients arise principally from variations in stellar metallicity within these galaxies, while age variations only make a significant contribution to the colour gradients of Virgo irregulars. A detailed stellar population analysis based on this material is presented in Roediger et al (2011b; arXiv:1011.3511).

preprint2011arXiv

The Formation and Evolution of Virgo Cluster Galaxies - II. Stellar Populations

We use a combination of deep optical and near-infrared light profiles for a morphologically diverse sample of Virgo cluster galaxies to study the radially-resolved stellar populations of cluster galaxies over a wide range of galaxy structure. We find that, in the median, the age gradients of Virgo galaxies are either flat (lenticulars and Sa-Sb spirals) or positive (ellipticals, Sbc+Sc spirals, gas-rich dwarfs, and irregulars), while all galaxy types have a negative median metallicity gradient. Comparison of the galaxy stellar population diagnostics (age, metallicity, and gradients thereof) against structural and environmental parameters also reveals that the ages of gas-rich systems depend mainly on their atomic gas deficiencies. Conversely, the metallicities of Virgo gas-poor galaxies depend on their concentrations, luminosities, and surface brightnesses. The stellar population gradients of all Virgo galaxies exhibit no dependence on either their structure or environment. We interpret these stellar population data for Virgo galaxies in the context of popular formation and evolution scenarios, and suggest that gas-poor giants grew hierarchically (through dissipative starbursts), gas-poor dwarfs have descended from at least two different production channels (e.g., environmental transformation and merging), while spirals formed inside-out, but with star formation in the outskirts of a significant fraction of the population having been quenched due to ram pressure stripping. (Abridged)

preprint2011arXiv

The Luminosity Profile and Structural Parameters of the Andromeda Galaxy

We have constructed an extended composite luminosity profile for the Andromeda galaxy, M31, and have decomposed it into three basic luminous structural components: a bulge, a disk and a halo. The dust-free Spitzer/IRAC imaging and extended spatial coverage of ground-based optical imaging and deep star counts allow us to map M31's structure from its center to 22 kpc along the major axis. We apply different decomposition methods for the 1D luminosity profiles and 2D images. These include non-linear least-squares and Bayesian Monte-Carlo Markov-chain analyses. The basic photometric model for M31 has a Sersic bulge with shape index n = 2.2 +/- 0.3 and effective radius R_e = 1.0 +/- 0.2 kpc, a dust-free exponential disk of scale length R_d = 5.3 +/- 0.5 kpc; the parameter errors reflect the range between various decomposition methods. The bulge parameter, n, is rather insensitive to bandpass effects and its value suggests a first rapid formation via mergers followed by secular growth from the disk. The M31 halo has a 2D power-law index ~ -2.5 +/- 0.2. We find that the M31 bulge light is mostly dominant over the range R_minor < 1.2 kpc. The disk takes over in the range 1.2 kpc < R_minor < 9 kpc, whereas the halo dominates at R_minor > 9 kpc. The stellar nucleus, bulge, disk, and halo components each contribute roughly 0.05%, 23%, 73% and 4% of the total light of M31 out to 200 kpc along the minor axis. Nominal errors for the structural parameters of the M31 bulge, disk and halo amount to 20%. If M31 and the Milky Way are at all typical, faint stellar halos should be routinely detected in galaxy surveys reaching below mu_i ~ 27 mag/arcsec^2. We stress that our results rely on this photometric analysis alone. Structural parameters may change when other fundamental constraints, such as those provided by abundance gradients and stellar kinematics, are considered simultaneously (Abridged).

preprint2010arXiv

Integrated Stellar Populations: Confronting Photometry with Spectroscopy

We investigate the ability of spectroscopic techniques to yield realistic star formation histories (SFHs) for the bulges of spiral galaxies based on a comparison with their observed broadband colors. Full spectrum fitting to optical spectra indicates that recent (within ~1 Gyr) star formation activity can contribute significantly to the V-band flux, whilst accounting for only a minor fraction of the stellar mass budget which is made up primarily of old stars. Furthermore, recent implementations of stellar population (SP) models reveal that the inclusion of a more complete treatment of the thermally pulsating asymptotic giant branch (TP-AGB) phase to SP models greatly increases the NIR flux for SPs of ages 0.2-2 Gyr. Comparing the optical--NIR colors predicted from population synthesis fitting, using models which do not include all stages of the TP-AGB phase, to the observed colors reveals that observed optical--NIR colors are too red compared to the model predictions. However, when a 1 Gyr SP from models including a full treatment the TP-AGB phase is used, the observed and predicted colors are in good agreement. This has strong implications for the interpretation of stellar populations, dust content, and SFHs derived from colors alone.

preprint2010arXiv

Stellar Population Trends in S0 Galaxies

We present stellar population age and metallicity trends for a sample of 59 S0 galaxies based on optical SDSS and NIR J & H photometry. When combined with optical g and r passband imaging data from the SDSS archive and stellar population models, we obtain radial age and metallicity trends out to at least 5 effective radii for most of the galaxies in our sample. The sample covers a range in stellar mass and light concentration. We find an average central light-weighted age of ~ 4 Gyr and central metallicity [Z/H] ~ 0.2 dex. Almost all galaxies show a negative metallicity gradient from the center out, with an average value of Delta[Z/H]/Delta(log(r/Re)) = -0.6. An age increase, decrease, and minimal change with radius is observed for 58%, 19%, and 23%, respectively, for a mean age gradient of Delta(age)/Delta(log(r/Re)) = 2.3 Gyr dex^{-1}. For 14 out of 59 galaxies, the light-weighted age of the outer region is greater than 10 Gyr. We find that galaxies with both lower mass and lower concentration have younger light-weighted ages and lower light-weighted metallicities. This mass-metallicity relation extends into the outer regions of our S0 galaxies. Our results are consistent with the formation of S0 galaxies through the transformation of spiral galaxy disks. Determining the structural component that makes up the outer region of galaxies with old outksirts is a necessary step to understand the formation history of S0 galaxies.

preprint2005arXiv

Origin of the Joint Distribution of Structural Parameters in Disc Galaxies

We use a simple model of disc-galaxy formation in the LCDM cosmology to simultaneously reproduce the slopes, zero-points, scatter and uncorrelated residuals in the observed velocity-luminosity (VL) and radius-luminosity (RL) relations. Observed I-band luminosities are converted to stellar masses using the IMF-dependent relation between stellar mass-to-light ratio, M/L, and color. The model treats halo concentration, spin parameter, and disc mass fraction as independent log-normal random variables. Our main conclusion is that the VL and RL zero-points and the uncorrelated residuals can only be reproduced simultaneously if adiabatic contraction is avoided. One or the other could be fixed by appealing to unrealistic values for M/L, halo concentration c, and spin parameter lambda, but not all together. The small VL scatter is naturally determined by the predicted scatter in c and in M/L, quite independent of the large scatter in lambda. However, the RL scatter, driven by the scatter in lambda, can be as low as observed only if sigma_lambda simeq 0.25, about half the value predicted for CDM haloes. This may indicate that discs form in a special subset of haloes. The VL slope is reproduced once star formation occurs only above a threshold surface density, and the RL slope implies that the disc mass fraction is increasing with halo mass, consistent with feedback effects. A model that incorporates the above ingredients provides a simultaneous fit to all the observed features. In particular, the elimination of adiabatic contraction allows 80% disc contribution to the observed rotation velocity of bright discs at 2.2 disc scale lengths. The lack of halo contraction may indicate that disc formation is not as smooth as typically envisioned, but instead may involve clumpy, cold streams.

preprint2002arXiv

Structure of Disk Dominated Galaxies I. Bulge/Disk Parameters, Simulations, and Secular Evolution

(Abridged) A robust analysis of galaxy structural parameters, based on the modeling of bulge and disk brightnesses in the BVRH bandpasses, is presented for 121 face-on and moderately inclined late-type spirals. Each surface brightness (SB) profile is decomposed into a sum of a generalized Sersic bulge and an exponential disk. The reliability and limitations of our bulge-to-disk (B/D) decompositions are tested with extensive simulations of galaxy brightness profiles (1D) and images (2D). Galaxy types are divided into 3 classes according to their SB profile shapes; Freeman Type-I and Type-II, and a third ``Transition'' class for galaxies whose profiles change from Type-II in the optical to Type-I in the infrared. We discuss possible interpretations of Freeman Type-II profiles. The Sersic bulge shape parameter for nearby Type-I late-type spirals shows a range between n=0.1-2 but, on average, the underlying surface density profile for the bulge and disk of these galaxies is adequately described by a double-exponential distribution. We confirm a coupling between the bulge and disk with a scale length ratio r_e/h=0.22+/-0.09, or h_bulge/h_disk=0.13+/-0.06 for late-type spirals, in agreement with recent N-body simulations of disk formation and models of secular evolution. This ratio increases from ~0.20 for late-type spirals to ~0.24 for earlier types. The similar scaling relations for early and late-type spirals suggest comparable formation and/or evolution scenarios for disk galaxies of all Hubble types.

preprint1996arXiv

The Bulge/Disk Connection in Late-type Spirals

Recent ground-based photometric investigations suggest that central regions of late-type spirals are closely coupled to the inner disk and probably formed via secular evolution. Evidence presented in support of this model includes the predominance of exponential bulges, the correlation of bulge and disk scale lengths, blueness of the bulge and small differences between bulge and central disk colors, detection of spiral structure into the core, and rapid rotation. Recent HST observations show that our own bulge and that of M31, M32, and M33 probably harbor both an old and intermediate-age populations in agreement with models of early collapse of the spheroid plus gas transfer from the disk. Secular evolution provides a mechanism to build-up central regions in late-type spirals; mergers or accretion of small satellites could explain the brighter, kinematically distinct bulges of Sa's and SO's.