The Color-Magnitude Relation in Coma: Clues to the Age and Metallicity of Cluster Populations

We have observed three fields of the Coma cluster of galaxies with a narrow band (modified Stromgren) filter system. Observed galaxies include 31 in the vicinity of NGC 4889, 48 near NGC 4874, and 60 near NGC 4839 complete to M_5500=-18 in all three subclusters. Spectrophotometric classification finds all three subclusters of Coma to be dominated by red, E type (ellipticals/S0's) galaxies with a mean blue fraction, f_B, of 0.10. The blue fraction increases to fainter luminosities, possible remnants of dwarf starburst population or the effects of dynamical friction removing bright, blue galaxies from the cluster population by mergers. We find the color-magnitude (CM) relation to be well defined and linear over the range of M_5500=-13 to -22. After calibration to multi-metallicity models, bright ellipticals are found to have luminosity weighted mean [Fe/H] values between -0.5 and +0.5, whereas low luminosity ellipticals have [Fe/H] values ranging from -2 to solar. The lack of CM relation in our continuum color suggests that a systematic age effect cancels the metallicity effects in this bandpass. This is confirmed with our age index which finds a weak correlation between luminosity and mean stellar age in ellipticals such that the stellar populations of bright ellipticals are 2 to 3 Gyrs younger than low luminosity ellipticals.Comment: 26 pages AAS LaTeX, 6 figures, accepted for publication in A

Hierarchical galaxy formation and substructure in the Galaxy's stellar halo

We develop an explicit model for the formation of the stellar halo from tidally disrupted, accreted dwarf satellites in the cold dark matter (CDM) framework, focusing on predictions testable with the Sloan Digital Sky Survey (SDSS) and other wide-field surveys. Subhalo accretion and orbital evolution are calculated using a semi-analytic approach within the Press-Schechter formalism. Motivated by our previous work, we assume that low-mass subhalos (v < 30 km/s) can form significant populations of stars only if they accreted a substantial fraction of their mass before the epoch of reionization. With this assumption, the model reproduces the observed velocity function of galactic satellites in the Local Group, solving the ``dwarf satellite problem'' without modifying the popular LCDM cosmology. The disrupted satellites yield a stellar distribution with a total mass and radial density profile consistent with those observed for the Milky Way stellar halo. Most significantly, the model predicts the presence of many large-scale, coherent substructures in the outer halo. These substructures are remnants of individual, tidally disrupted dwarf satellite galaxies. Substructure is more pronounced at large galactocentric radii because of the smaller number density of tidal streams and the longer orbital times. This model provides a natural explanation for the coherent structures in the outer stellar halo found in the SDSS commissioning data, and it predicts that many more such structures should be found as the survey covers more of the sky. The detection (or non-detection) and characterization of such structures could eventually test variants of the CDM scenario, especially those that aim to solve the dwarf satellite problem by enhancing satellite disruption.Comment: 12 pages, 8 figures, Submitted to Ap

On the Structural Differences between Disk and Dwarf Galaxies

Gas-rich dwarf and disk galaxies overlap in numerous physical quantities that make their classification subjective. We report the discovery of a separation between dwarfs and disks into two unique sequences in the mass (luminosity) versus scale length plane. This provides an objective classification scheme for late-type galaxies that only requires optical or near-IR surface photometry of a galaxy. Since the baryonic Tully-Fisher relation for these samples produces a continuous relation between baryonic mass and rotational velocity, we conclude that the difference between dwarfs and disks must be because of their distribution of stellar light such that dwarfs are more diffuse than disk galaxies. This structural separation may be due to a primordial difference between low and high mass galaxies or produced by hierarchical mergers where disks are built up from dwarfs. Structural differences between dwarf and disk galaxies may also be driven by the underlying kinematics where the strong rotation in disks produces an axial symmetric object that undergoes highly efficient star formation in contrast to the lower rotation, more disordered motion of dwarfs that produces a diffuse, triaxial object with a history of inefficient star formation.Comment: 16 pages, 2 figures, AJ in press, AASTeX5.

How special are Brightest Group and Cluster Galaxies?

We use the Sloan Digital Sky Survey to construct a sample of 625 brightest group and cluster galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, redshift, and color. We investigate how the systematic properties of BCGs depend on stellar mass and on their privileged location near the cluster center. The groups and clusters that we study are drawn from the C4 catalogue of Miller et al. (2005) but we have developed improved algorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense environments, we have developed a correction for this effect which can be readily applied to the published catalog data. We find that BCGs are larger and have higher velocity dispersions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger fraction of dark matter. In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs does not vary as a function of galaxy luminosity. Hence BCGs lie on a different fundamental plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber-Jackson relation than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers along preferentially radial orbits. We find tentative evidence that this steepening is stronger in more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs of the same mass, but they have somewhat higher alpha/Fe ratios, indicating that star formation may have occurred over a shorter timescale in the BCGs. Finally, we find that BCGs are more likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are less likely to host an optical AGN. The differences we find are more pronounced for the less massive BCGs.Comment: Replaced with slightly modified version accepted by MNRAS. 28 pages, 25 figures. Version with full resolution figures available at http://www.mpa-garching.mpg.de/~anja/bcgs_avdl.pd

Halo Substructure and the Power Spectrum

In this proceeding, we present the results of a semi-analytic study of CDM substructure as a function of the primordial power spectrum. We apply our method to several tilted models in the LCDM framework with n=0.85-1.1, sigma_8=0.65-1.2 when COBE normalized. We also study a more extreme, warm dark matter-like spectrum that is sharply truncated below a scale of 10^10 h^-1 Msun. We show that the mass fraction of halo substructure is not a strong function of spectral slope, so it likely will be difficult to constrain tilt using flux ratios of gravitationally lensed quasars. On the positive side, all of our CDM-type models yield projected mass fractions in good agreement with strong lensing estimates: f \sim 1.5% at M \sim 10^8 Msun. The truncated model produces a significantly smaller fraction, f \lsim 0.3%, suggesting that warm dark matter-like spectra may be distinguished from CDM spectra using lensing. We also discuss the issue of dwarf satellite abundances, with emphasis on the cosmological dependence of the map between the observed central velocity dispersion of Milky Way satellites and the maximum circular velocities of their host halos. In agreement with earlier work, we find that standard LCDM over-predicts the estimated count of Milky Way satellites at fixed Vmax by an order of magnitude, but tilted models do better because subhalos are less concentrated. Interestingly, under the assumption that dwarfs have isotropic velocity dispersion tensors, models with significantly tilted spectra (n \lsim 0.85, sigma_8 \lsim 0.7) may under-predict the number of large Milky Way satellites with Vmax \gsim 40 km/s.Comment: 5 pages, 2 figures. Poster contribution to the 13th Annual Astrophysics Conference in Maryland, The Emergence of Cosmic Structur

Reionization and the abundance of galactic satellites

One of the main challenges facing standard hierarchical structure formation models is that the predicted abundance of galactic subhalos with circular velocities of 10-30 km/s is an order of magnitude higher than the number of satellites actually observed within the Local Group. Using a simple model for the formation and evolution of dark halos, based on the extended Press-Schechter formalism and tested against N-body results, we show that the theoretical predictions can be reconciled with observations if gas accretion in low-mass halos is suppressed after the epoch of reionization. In this picture, the observed dwarf satellites correspond to the small fraction of halos that accreted substantial amounts of gas before reionization. The photoionization mechanism naturally explains why the discrepancy between predicted halos and observed satellites sets in at about 30 km/s, and for reasonable choices of the reionization redshift (z_re = 5-12) the model can reproduce both the amplitude and shape of the observed velocity function of galactic satellites. If this explanation is correct, then typical bright galaxy halos contain many low-mass dark matter subhalos. These might be detectable through their gravitational lensing effects, through their influence on stellar disks, or as dwarf satellites with very high mass-to-light ratios. This model also predicts a diffuse stellar component produced by large numbers of tidally disrupted dwarfs, perhaps sufficient to account for most of the Milky Way's stellar halo.Comment: 5 pages, 2 figures, Submitted to Ap

The Nature of the Halo Population of NGC 5128 Resolved with NICMOS on the Hubble Space Telescope

We present the first infrared color-magnitude diagram (CMD) for the halo of a giant elliptical galaxy. The CMD for the stars in the halo of NGC 5128 (Centaurus A) was constructed from HST NICMOS observations of the WFPC2 CHIP-3 field of Soria et al. (1996) to a 50% completeness magnitude limit of [F160W]=23.8. This field is located at a distance of 08'50" (~9 kpc) south of the center of the galaxy. The luminosity function (LF) shows a marked discontinuity at [F160W]=20.0. This is 1-2 mag above the tip of the red giant branch (TRGB) expected for an old population (~12 Gyr) at the distance modulus of NGC 5128. We propose that the majority of stars above the TRGB have intermediate ages (~2 Gyr), in agreement with the WFPC2 observations of Soria et al. (1996). Five stars with magnitudes brighter than the LF discontinuity are most probably due to Galactic contamination. The weighted average of the mean giant branch color above our 50% completeness limit is [F110W]-[F160W]=1.22+-0.08 with a dispersion of 0.19 mag. From our artificial-star experiments we determine that the observed spread in color is real, suggesting a real spread in metallicity. We estimate the lower and upper bounds of the stellar metallicity range by comparisons with observations of Galactic star clusters and theoretical isochrones. Assuming an old population, we find that, in the halo field of NGC 5128 we surveyed, stars have metallicities ranging from roughly 1% of solar at the blue end of the color spread to roughly solar at the red end, with a mean of [Fe/H]=-0.76 and a dispersion of 0.44 dex.Comment: Accepted for publication in AJ, 23 pages of text, 13 figures, uses aastex v5.

The Age of Cluster Galaxies from Continuum Colors

We determine the age of 1,104 early-type galaxies in eight rich clusters ($z = 0.0046$ to $0.175$) using a new continuum color technique. We find that galaxies in clusters divide into two populations, an old population with a mean age similar to the age of the Universe (12 Gyrs) and a younger population with a mean age of 9 Gyrs. The older population follows the expected relations for mass and metallicity that imply a classic monolithic collapse origin. Although total galaxy metallicity is correlated with galaxy mass, it is uncorrelated with age. It is impossible, with the current data, to distinguish between a later epoch of star formation, longer duration of star formation or late bursts of star formation to explain the difference between the old and young populations. However, the global properties of this younger population are correlated with cluster environmental factors, which implies secondary processes, post-formation epoch, operate on the internal stellar population of a significant fraction of cluster galaxies. In addition, the mean age of the oldest galaxies in a cluster are correlated with cluster velocity dispersion implying that galaxy formation in massive clusters begins at earlier epochs than less massive clusters.Comment: 35 pages, 10 figures, accepted by Ap

Dark Matter and Stellar Mass in the Luminous Regions of Disk Galaxies

We investigate the correlations among stellar mass (M_*), disk scale length (R_d), and rotation velocity at 2.2 disk scale lengths (V_2.2) for a sample of 81 disk-dominated galaxies (disk/total >= 0.9) selected from the SDSS. We measure V_2.2 from long-slit H-alpha rotation curves and infer M_* from galaxy i-band luminosities (L_i) and g-r colors. We find logarithmic slopes of 2.60+/-0.13 and 3.05+/-0.12 for the L_i-V_2.2 and M_*-V_2.2 relations, somewhat shallower than most previous studies, with intrinsic scatter of 0.13 dex and 0.16 dex. Our direct estimates of the total-to-stellar mass ratio within 2.2R_d, assuming a Kroupa IMF, yield a median ratio of 2.4 for M_*>10^10 Msun and 4.4 for M_*=10^9-10^10 Msun, with large scatter at a given M_* and R_d. The typical ratio of the rotation speed predicted for the stellar disk alone to the observed rotation speed at 2.2R_d is ~0.65. The distribution of R_d at fixed M_* is broad, but we find no correlation between disk size and the residual from the M_*-V_2.2 relation, implying that this relation is an approximately edge-on view of the disk galaxy fundamental plane. Independent of the assumed IMF, this result implies that stellar disks do not, on average, dominate the mass within 2.2R_d. We discuss our results in the context of infall models of disk formation in cold dark matter halos. A model with a disk-to-halo mass ratio m_d=0.05 provides a reasonable match to the R_d-M_* distribution for spin parameters \lambda ranging from ~0.04-0.08, and it yields a reasonable match to the mean M_*-V_2.2 relation. A model with m_d=0.1 predicts overly strong correlations between disk size and M_*-V_2.2 residual. Explaining the wide range of halo-to-disk mass ratios within 2.2R_d requires significant scatter in m_d values, with systematically lower m_d for galaxies with lower $M_*$.Comment: 18 pages, 2 tables, 7 figures, Accepted to ApJ, Table 1 updated, otherwise minor change

The Velocity Function of Galaxies

We present a galaxy circular velocity function, Psi(log v), derived from existing luminosity functions and luminosity-velocity relations. Such a velocity function is desirable for several reasons. First, it enables an objective comparison of luminosity functions obtained in different bands and for different galaxy morphologies, with a statistical correction for dust extinction. In addition, the velocity function simplifies comparison of observations with predictions from high-resolution cosmological N-body simulations. We derive velocity functions from five different data sets and find rough agreement among them, but about a factor of 2 variation in amplitude. These velocity functions are then compared with N-body simulations of a LCDM model (corrected for baryonic infall) in order to demonstrate both the utility and current limitations of this approach. The number density of dark matter halos and the slope of the velocity function near v_*, the circular velocity corresponding to an ~L_* spiral galaxy, are found to be comparable to that of observed galaxies. The primary sources of uncertainty in construction of Psi(log v) from observations and N-body simulations are discussed and explanations are suggected to account for these discrepancies.Comment: Latex. 28 pages, 4 figures. Accepted by Ap