A&A 439, 997-1011 (2005)
DOI: 10.1051/0004-6361:20047012
T. H. Puzia1,2, - M. Kissler-Patig3 -
D. Thomas4,5 - C. Maraston4,5 - R. P. Saglia4 - R. Bender2,4 - P. Goudfrooij 1 - M. Hempel3
1 - Space Telescope Science Institute, 3700 San Martin Drive,
Baltimore, MD 21218, USA
2 -
Sternwarte der Ludwig-Maximilians-Universität,
Scheinerstr. 1, 81679 München, Germany
3 -
European Southern Observatory, 85749 Garching bei München,
Germany
4 -
Max-Planck-Institut für extraterrestrische Physik,
Giessenbachstrasse, 85748 Garching bei München, Germany
5 -
University of Oxford, Astrophysics, Keble Road, Oxford, OX1 3RH, UK
Received 5 January 2004 / Accepted 13 May 2005
Abstract
An analysis of ages, metallicities, and [/Fe] ratios of globular
cluster systems in early-type galaxies is presented, based on Lick index
measurements summarized in Puzia et al. (2004, A&A, 415, 123, Paper I of this series). In
the light of calibration and measurement uncertainties, age-metallicity
degeneracy, and the relative dynamic range of Lick indices, as well as
systematics introduced by abundance ratio variations (in particular
variations of [
/Fe] ratios), we find that the most reliable age
indicator for our dataset is a combination of the Lick Balmer-line indices
H
,
H
,
and H
.
[MgFe]
is used
as a spectroscopic metallicity indicator which is least affected by
[
/Fe] variations. We introduce an interpolation routine to
simultaneously derive ages, metallicities, and [
/Fe] ratios from
diagnostic grids constructed from Lick indices. From a comparison of
high-quality data with SSP model predictions, we find that
2/3 of
the globular clusters in early-type galaxies are older than 10 Gyr, up to
1/3 have ages in the range
5-10 Gyr, and only a few cluster are
younger than
5 Gyr. Our sample of globular clusters covers
metallicities from [Z/H]
-1.3 up to
0.5 dex. We find that
metal-rich globular clusters show on average a smaller mean age and a
larger age scatter than their metal-poor counterparts. [
/Fe] diagnostic plots show that globular cluster systems in early-type galaxies have super-solar
/Fe abundance ratios with a mean [
/Fe] = 0.47
0.06 dex and a dispersion of
0.3 dex. We
find evidence for a correlation between [
/Fe] and metallicity, in
the sense that more metal-rich clusters exhibit lower
-element
enhancements. A discussion of systematics related to the Lick index system
shows that the method suffers to some extent from uncertainties due to
unknown horizontal branch morphologies at high metallicities. However,
these systematics still allow us to make good qualitative statements. A
detailed investigation of indices as a function of data quality reveals
that the scatter in Balmer index values decreases for higher-quality data.
In particular, extremely low Balmer index values that are lower than any
SSP model prediction tend to disappear. Furthermore, we find that observed
photometric colors are in good agreement with computed SSP colors using
ages and metallicities as derived from the spectroscopic line indices.
Key words: galaxies: star clusters - galaxies: general
Two prominent models of early-type galaxy formation are lively discussed
in the literature. In the monolithic-collapse scenario
(Tinsley 1972; Larson 1975; Silk 1977; Arimoto & Yoshii 1987, etc.) the majority of
stars in early-type galaxies forms early, at redshifts .
Empirical
scaling laws, such as the fundamental plane, the Mg-
relation,
and the [
/Fe]-
relation support such an early formation
epoch (Trager et al. 2000a; Bender 1996; Kuntschner 2000; Bower et al. 1992; Thomas et al. 2005; Trager et al. 2000b; Peebles 2003; Bender et al. 1996; Treu et al. 1999). On the other hand, the hierarchical merging picture (Baugh et al. 1998; Kauffmann et al. 1993; Somerville et al. 2001; White & Rees 1978; Cole et al. 2000) sees early-type galaxies as the result of
multiple merging and accretion events of smaller units over an extended
period of time until the very recent past. In this way a significant
fraction of stars is formed below a redshift of unity. Observed ongoing
mergers and accretion events (e.g. Whitmore & Schweizer 1995; Ibata et al. 1995; Koo et al. 1996)
and the existence of kinematically decoupled cores (Jedrzejewski & Schechter 1988; Bender 1988; Franx & Illingworth 1988; Surma & Bender 1995; Davies et al. 2001) are clearly enforcing arguments for
this scenario.
These models describe two antipodal paradigms of galaxy formation. The quantification of their importance as a function of redshift, environment, and galaxy morphology is necessary to provide detailed insight in galaxy formation and assembly. A major difference between the two pictures are the different star-formation histories of early-type galaxies. While in the hierarchical picture most massive galaxies are thought to experience long assembly time scales, the monolithic collapse scenario predicts very early and short bursts of star formation. Clearly, the predicted star-formation histories stand in marked contrast and are an important piece of evidence to differentiate between these two models.
However, to recover star-formation histories, one ideally has to resolve the underlying stellar populations which build up a galaxy. Only a few early-type systems are close enough so that their halos can be resolved into single stars (e.g. in M 32: Grillmair et al. 1996, NGC 3115: Elson 1997, NGC 5128: Harris & Harris 2002, NGC 3379: Gregg et al. 2004). Most photometric and spectroscopic studies of the unresolved diffuse light are hampered by the mix of ages, metallicities, and abundance ratios in the stellar populations (e.g. Maraston & Thomas 2000). In combination with the well-known age-metallicity degeneracy (Faber 1972; O'Connell 1976; Worthey 1994), it is extremely difficult to disentangle even the different major stellar populations from studies of the diffuse light only, let alone to reconstruct a detailed star formation history.
Globular cluster systems are very useful tools to study star formation histories of galaxies. Several arguments support the hypothesis that globular clusters trace major star-formation events in galaxies. (1) The formation of massive star clusters, which are likely to survive a Hubble time as globular clusters, is observed in interacting/merging and starburst galaxies (e.g. Homeier et al. 2002; Whitmore & Schweizer 1995; Schweizer 1997; Goudfrooij et al. 2001; Johnson et al. 1999). (2) Young massive star clusters are observed in "simmering'' late-type galaxies. Their number is correlated with the star formation rate per unit area in these systems (Larsen & Richtler 2000). (3) Normalizing the number of globular clusters to the total baryonic mass of the host, reveals a surprisingly constant value (McLaughlin 1999) and points towards a tight link between star and globular cluster formation.
In other words, globular clusters are fossil records of major star-formation episodes in galaxies. Their ages, metallicities, and abundance ratios can provide detailed information on the formation histories of their host systems. As simple stellar populations, globular clusters consist of stars characterized by one age and one metallicity. Hence, the interpretation of their observed colours and spectra is less ambiguous than for the diffuse light. Their ubiquity in all galaxy types and their high surface brightness make them easy to observe out to large distances.
Spectroscopy opens an independent way to access ages and metallicities of
globular clusters besides photometry. Low resolution spectroscopy
(R
1000) provides information on the strength of prominent diagnostic
features, such as the Balmer series and some relatively strong iron and
other metal features. Moreover, it can provide clues on the basic
chemistry of globular cluster systems.
The Lick index system standardizes the measurement of spectroscopic line
indices for many strong absorption features (Worthey et al. 1994; Burstein et al. 1984; Worthey & Ottaviani 1997; Trager et al. 1998). The outstanding role of the Lick system is its
provision of index measurements for many stars with a wide range in
,
and [Fe/H]. With this information, so-called fitting
functions
are
computed which are an essential ingredient of theoretical SSP model
predictions. Using index response functions
(Tripicco & Bell 1995),
recent SSP models (Thomas et al. 2003; Trager et al. 2000a; Thomas et al. 2004) provide also
information on fundamental abundance ratios, such as [
/Fe]. The
combination of age-sensitive (e.g. Balmer-line indices) and
metallicity-sensitive indices (
Fe
,
Mgb, [MgFe]
,
etc.) allows in principle to derive accurate spectroscopic ages and
metallicities. This approach is less affected by the age-metallicity
degeneracy than broadband photometry. In addition, combining indices
sensitive to the abundance of
-elements and iron can provide
important clues on star-formation time scales (Tinsley 1979; Greggio 1997; Matteucci 1994; Thomas et al. 1999).
The present paper makes use of the spectroscopic data presented in
Puzia et al. (2004, hereafter Paper I) to derive global ages,
metallicities, and [/Fe] ratios for globular clusters in
early-type galaxies. A more detailed study of these global parameters
as a function of internal galaxy properties, such as morphological
type, environment, galaxy mass, etc. will be presented in subsequent
papers of this series.
We describe the selection of high-quality spectroscopic data in
Sect. 2. The best combination of Lick indices to achieve most
reliable age and metallicity estimates is discussed in
Sect. 3. In Sect. 4 we derive ages and
metallicities of globular clusters in the studied early-type
galaxies. Global [/Fe] ratios and correlations with age and metallicity
for these globular clusters are presented in Sect. 5.
We discuss our results in Sect. 6.
Table 1: Basic information on host galaxies. According to the column the references are: (1) de Vaucouleurs et al. (1991); (2) Tonry et al. (2001); (3) Tully (1988); (4) and (5) Kissler-Patig et al. (1997), Ashman & Zepf (1998).
Our entire sample contains 143 globular clusters in seven early-type
galaxies. About half of our spectra satisfy the ideal S/N standards
(30 per Å) to derive accurate ages, metallicities, and
[
/Fe] ratios. It is possible to achieve an age resolution of
1 Gyr (at ages
Gyr) only for the brightest globular
clusters in our sample. The typical separation between the 15 and 14 Gyr
isochrone in current SSP models is of the order
H
Å, and
0.1 Å for the higher-order
Balmer indices. This separation increases towards younger ages. Hence, the
final sample has to be built from a compromise between age/metallicity
resolution and sample size. We, therefore, set the selection to clusters
with a statistical measurement uncertainty of
0.4 Å and
0.6 Å for H
and higher-order Balmer line indices, respectively. This
selection corresponds to a minimum age resolution
.
As
a metallicity indicator, we use the composite [MgFe]
index (see below). An error cut at 0.2 Å for this index guarantees a
metallicity resolution of
0.25 dex at high and
0.4 dex at low
metallicities. The above selection criteria leave 71 globular cluster
spectra in our sample which correspond to
50% of the initial data.
Henceforth, we consider only these selected globular clusters and combine
the data of all sample galaxies to increase statistical significance in
the discussed relations. Note, that this sample is biased towards the
bright end of the globular cluster luminosity function, i.e. most massive
clusters. Otherwise, it spans a wide range in globular cluster color and
host galaxy properties (see Table 1).
Our original colour selection criteria prevented us from observing
globular clusters with a combination of low metallicities [Z/H]
-1.3and relatively young ages (
Gyr). However, these objects are
essentially non-existent in our sample galaxies (see Paper I).
Figure 1 illustrates the age and metallicity ranges imposed
by our colour selection. All Milky Way globular clusters would be selected
by our colour cuts
,
,
and
(we refer to Paper I for details). The applied colour cuts
also suggest that very metal-rich old globular clusters might be missing in our
final sample. We point out that, within our luminosity cut (
mag),
the fraction of such objects amounts only a few percent. In order to fill
the slit-masks we included cluster candidates redder than our colour cuts.
However, only two genuine globular clusters were found among
30"mask-filler'' objects. Unlike the expected old ages, we find
intermediate ages (5-10 Gyr) for these two globular clusters. The fraction
of old to intermediate-age, very metal-rich globular clusters is likely to
increase in a sample reaching fainter magnitudes. However, given our
luminosity selection, the number of missed objects is negligible.
![]() |
Figure 1: Illustrated are colours which were used to select globular cluster candidates as a function of age and metallicity and taken from model predictions of Maraston (2005). Upper and lower colour cuts are indicated as horizontal lines. Curves in left panels are parametrized by metallicity, curves in right panels are parametrized by age and indicated accordingly. |
Open with DEXTER |
Our dataset does not guarantee a strict consistency of the derived mean
ages, metallicities and [/Fe] ratios among the galaxies of our
sample. These parameters are subject to change, since neither the sampled
fraction of each globular cluster system, nor the sampling of colour
distributions and luminosity functions is identical from galaxy to galaxy.
Given the still relatively small numbers of globular clusters per galaxy,
a peculiar age and/or metallicity distribution in one galaxy can influence
the relation between age, metallicity and [
/Fe] of the complete
sample, although we verified that no galaxy represents a clear outlier
(see future papers in these series).
In the following we determine the best combination of indices as diagnostics for age and metallicity. Taking into account the uncertainty of our line index measurements, the mean uncertainties of the Lick system, and the limits on the predicting power of SSP models, we construct the relatively best combination of diagnostic diagrams from Lick line indices. The found combination maximally reduces the internal uncertainties of the Lick system and the age-metallicity degeneracy of line indices as well as the influence of abundance ratio variations.
![]() |
Figure 2:
The behavior of Lick Balmer indices as a function the horizontal branch morphology (HBR). This parameter is defined in Lee et al. (1994) as HBR (B-R)/(B+V+R), where B and R are the number of stars blue-wards and red wards of the instability strip. V is the number of variable stars inside the instability strip. HBR = 1 indicates an entirely blue and HBR = -1and entirely red horizontal branch. Here we plot data for globular
clusters in the Milky Way ( squares: data from Puzia et al. 2002 and
Trager et al. 1998), M 31 ( invertedtriangles: Puzia et al. 2005,
Rich et al. 2003, and Trager et al. 1998), and the Large Magellanic Cloud
triangles: Beasley et al. 2002b). Filled symbols show globular
clusters with a metallicity [Z/H] -0.6 (note, all have HBR = -1).
NGC 6388 and NGC 6441 are indicated by stars at HBR = -0.7. A small arrow
at each panel's right ordinate indicates the most extreme Balmer index
value for globular clusters with [Z/H] > -0.6. Several clusters show a
tendency for higher Balmer indices at HBR ![]() ![]() ![]() |
Open with DEXTER |
To empirically estimate the effect of a varying horizontal branch
morphology on Balmer indices at high metallicities, which might be
responsible for the increased Balmer indices of metal-rich globular
clusters, we parametrize the HB morphology with the HBR parameter
(Lee et al. 1994). Figure 2 shows the strength of Balmer
line indices as a function of HB morphology (HBR). In general, blue
horizontal branches produce significantly stronger Balmer indices for
globular clusters in the Milky Way (Puzia et al. 2002; Trager et al. 1998),
M 31 (Puzia et al. 2005; Trager et al. 1998), and the LMC (Beasley et al. 2002b). Each
panel shows a sequence of metallicity, where the scatter can be
attributed to the "second parameter''. We use the globular clusters
NGC 6388 and NGC 6441 (which host the bluest horizontal branches among
metal-rich Galactic globular clusters, indicated by stars in
Fig. 2) and clusters at [Z/H]
-0.6 (filled
symbols) which have entirely red HBs (e.g. NGC 6356 and NGC 6637) to
derive a representative "second-parameter'' variation of Balmer line
indices at high metallicities. As this approach is fully empirical and
based on the largest HB morphology fluctuation locally observed, we cannot
rule out that even more extreme HB morphologies for globular clusters at a
given metallicity exist outside the Local Group. We find, in the extreme
case, offsets of 0.4 Å in H
,
3.3 Å in H
,
1.4 Å in H
,
2.0 Å in H
,
and 1.0 Å in H
between metal-rich globular cluster with the weakest Balmer index
(with entirely red HBs) and NGC 6388 and 6441. The HB morphology has
negligible effect on the [MgFe]
index. Consequently, increased
Balmer indices at high metallicities might well be, at least partly, the
result of HB morphology variations.
![]() |
Figure 3:
Influence of [![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In the following we use the Balmer indices as age indicators with
confidence at low metallicity ([Z/H]
-0.6), because HB morphology is
included in our SSP models and under control. At high metallicity
([Z/H]
-0.6), we have the warning in mind that ages could be
degenerate with the presence of unresolved blue HBs. For this reason we
will refer to such young ages as "formal''.
![]() |
Figure 4:
Passband definitions for Balmer-line, Mgb, Fe5270, and Fe5335 Lick indices
with their feature and adjacent continuum passbands. The over-plotted
spectrum is a high-S/N spectrum of the Galactic globular cluster NGC 6284
(Puzia et al. 2002). The resolution is ![]() |
Open with DEXTER |
Recently, Thomas et al. (2003,2004) calculated new theoretical Lick index
predictions which are parametrized for well-defined [/Fe] ratios for a wide range of ages and metallicities. These models take into
account the effects of changing element abundance ratios on Lick indices,
hence give Lick indices not only as a function of age and metallicity, but
also as a function of the [
/Fe] ratio. They are based on the
evolutionary population synthesis code of Maraston (1998). The impact
from element-ratio changes is computed with the help of the
Tripicco & Bell (1995) and Korn et al. (2005) response functions, using an extension
of the method introduced by Trager et al. (2000a). Because of the inclusion of
element-ratio effects, the influence of [
/Fe] on Balmer indices
can be studied, and is illustrated in Fig. 3.
Here we show the influence of [
/Fe] variations on frequently used
age/metallicity diagnostic diagrams. In general, variations of Balmer-line
indices for isochrones with [
/Fe] ratios between solar and +0.5 dex are of the order
0.007-0.25 Å for low, and
0.18-3.63 Å for high metallicities, and introduce a relative age uncertainty in the range
.
These variations are due to different contaminations by metal absorption features inside the Balmer index
passband definitions, as illustrated in Fig. 4.
The figure shows the higher-order Balmer indices generally include more
metal absorption lines in their feature and pseudo-continuum passbands
than the H
index, which is the Balmer index with the narrowest passband
definitions. Hence, among all Lick Balmer line indices, the least
[
/Fe]-sensitive index is H
,
followed by H
,
H
,
H
,
H
(see
also Thomas et al. 2003,2004). This exercise demonstrates that at high
metallicities the impact of [
/Fe] variations on age/metallicity
determinations (in particular, those of early-type galaxies) can
significantly alter the results and/or introduce spurious correlations.
The Balmer line series provides the best spectroscopic age indicator
among the set of Lick line indices. The Lick system defines five indices (H,
H
,
H
,
H
,
and H
)
for three Balmer lines (Worthey et al. 1994; Worthey & Ottaviani 1997). Figure 4 shows the passband
definitions for all Balmer indices. In combination with a
metallicity diagnostic, these higher-order Balmer line indices are
widely used to determine (luminosity-weighted) ages and metallicities
of galaxies (e.g. Thomas et al. 2003; Trager et al. 2000a; Kuntschner et al. 2002a; Poggianti et al. 2001; Kuntschner 2000; Trager et al. 2000b,1998).
Table 2: Summary of the coefficients relevant to Eq. (1). The coefficients in Cols. 2-8 are given in units of Å. Columns 9-14 are given in dex/Gyr. Note that for the higher-order Balmer indices Col. 5 gives the arithmetic mean of the uncertainties for cool and warm stars (see Tables 3 and 4 in Worthey & Ottaviani 1997).
However, different types of diagnostic plots based on different
Balmer-line indices are used throughout the literature. Although the age
predicting power of an arbitrarily chosen diagnostic plot (most common
versions include the H
and Mg2 or
Fe
indices
might yield accurate-enough results for a specific scientific goal (e.g.
the mean age difference between two different galaxy samples), the choice
of a specific diagnostic plot is still subject to observational
constraints and individual assessment, and makes comparisons between
studies difficult. As a consequence, most authors use several diagnostic
plots with different Balmer line indices and assign equal importance to
the results derived from each of those.
In the following we provide a recipe to define a quantity from which the relatively best Balmer-line age indicator can be determined. This quantity takes into account the quality of a given dataset and the diagnostic power of theoretical predictions from which one intends to derive the age and metallicity. It does not take into account systematic uncertainties in the fitting functions for a particular index.
In particular, the age sensitivity of an index is a function of the following parameters:
The highest
indicates the best age indicator with least
age-metallicity and age-
/Fe degeneracy. In
Table 2 we provide values for
at two different metallicities and for
at [
/Fe] = 0.3 and four age-metallicity combinations for each Balmer line index. Since SSP models do not provide continuous but discrete
predictions the partial derivatives are substituted by difference ratios,
e.g.
.
The quotients
are determined by linear interpolation of SSP models.
We determine the relatively best age indicator from the set of five Lick
Balmer indices by combining the mean dynamic range
,
the mean age-metallicity and age-
/Fe degeneracy
parameter
,
and the total index
uncertainty which is the denominator in Eq. (1). The final
mean
is documented in the last column of
Table 2. We find that the relatively best age
diagnostic for our data is the H
index followed by the
indices H
and H
.
H
and H
have the
smallest
values and are not considered to be reliable age indicators.
It is instructive to see that despite the relatively large age-metallicity
degeneracy of the H
index, it still provides the most
accurate age predictions. This fact is primarily due to the large dynamic
range of H
compared to its mean measurement uncertainty.
H
,
on the other hand, has a relatively large total uncertainty and
the measurements will therefore be more scattered over the diagnostic
plot's parameter range. However, H
shows by far the highest
value (Col. 13 in
Table 2; see also Fig. 3) and
is therefore least sensitive to the cumulative effects of age-metallicity
and age-
/Fe degeneracy. In general, the higher-order Balmer lines
require less S/N to guarantee a similar total index accuracy as H
,
which is primarily due to the narrow passband definition of H
(see
Fig. 4). If our dataset would be infinitely
accurate (i.e.
in Eq. (1)), the order of
from
the best to worst Balmer index would remain unchanged. The value of
is predominantly governed by uncertainties in the fitting functions of
the respective index. To vary this order, the mean measurement
uncertainties have to be very discrepant and the SSP-model predictions
have to deviate significantly from the model used here. It is expected
that the relative accuracy of Balmer index measurements is comparable
between different datasets as they are usually derived from one optical
spectrum. The relative age scale of SSP models appears to be quite
stable against the choice of different stellar evolutionary tracks for
ages >1 Gyr (Trager et al. 2000a; Charlot et al. 1996; Maraston 2003). This scale is
used in our above prescription. It can therefore be expected that no large
fluctuation in
will arise from the use of different SSP model
predictions
.
The index with the highest metallicity sensitivity and minimal
age-sensitivity could in principle be found in a comparable way as it
was done for the relatively best age diagnostic. The major impact of
typical metallicity tracers, such as
,
Mg2,
and Mgb, on the absolute metallicity scale is expected to arise from
changing abundance ratios. To reduce the influence of [
/Fe] variations on age and metallicity determinations, Thomas et al. (2003) modify the [MgFe] index
to obtain an [
/Fe]-insensitive metallicity index,
It is known that 40-60% of early-type galaxies show indications of
emission in their absorption spectra (Phillips et al. 1986; Caldwell 1984; González 1993; Goudfrooij et al. 1994). In a narrow-band imaging survey,
Macchetto et al. (1996) find ionized gas in
80% of early-type galaxies
including flocculent H
[N II] emission in NGC 3379, NGC 5846,
and NGC 7192, well within
,
but no significant
emission in NGC 3115. This gas is located in the central parts and
distributed in a rather regular way, suggestive of a disk. If dominant,
all Balmer indices would be affected, along with potential contamination
of Fe5015 by [O III] (
5007 Å) and of Mgb by [N I] (
5199 Å) (Goudfrooij & Emsellem 1996). In that case, our measurements
would indicate too old ages. This effect rapidly decreases from H
towards H
and H
(Osterbrook 1989), i.e. higher-order Balmer indices are less affected and should be
preferentially used for age determinations in the presence of ionized gas.
In order to exclude a major effect of ionized gas on line-strength measurements in our globular cluster data, we performed several tests/estimates.
Since line emission is concentrated in the central parts of galaxies
(Macchetto et al. 1996), we expect a correlation of Balmer indices with
galacto-centric radius if line-emission contamination is significant. We
find no evidence that Balmer indices are correlated with galacto-centric
distance. A more detailed analysis of background spectra shows that most
clusters are located within
and that the flux
level of the diffuse galaxy light is well below the object flux. In
particular, we find no correlation inside one effective radius, where line
emission is expected to be strongest. Furthermore, we find no correlations
of Balmer indices for globular clusters with Balmer indices measured on
corresponding background spectra.
![]() |
Figure 5:
Computed vs. measured photometric colors for our sample
globular clusters. Computed colors were derived from spectroscopic ages
and metallicities using the SSP models of Maraston (2005). The circles
show all selected globular clusters for which optical/near-infrared colors
were obtained (see Paper I for details). Solid circles show the
high-quality sample, which was selected by a more constrained error
cut (see text for details). Errorbars indicate 1-![]() |
Open with DEXTER |
Visual re-inspection of the background subtraction process for some
low-H
outliers (see Fig. 6) underline the good
quality of background modeling and subtraction. However, problems with
accurate background subtraction might occur in very few cases when line
emission has a very filamentary structure. For instance, the worst case
scenario would be when a globular cluster overlaps with a filament of
ionized gas while the slit is aligned perpendicular to such a filament.
However, a filamentary emission pattern is not found in the
Macchetto et al. (1996) study for our host galaxies. We conclude that line
emission has no measurable effect on the scatter in the age/metallicity
and [
/Fe] diagnostic plots.
Having determined the relatively best metallicity diagnostic, we determine the best combination of Balmer-line indices as our prime age indicator. In deciding whether a Balmer index will be chosen as part of this most reliable age proxy, we inspect the following points for each Balmer index and assign priorities in descending order:
Since the SSP model predictions allow full control over [/Fe] variations within the diagnostic grids, we use an iterative approach in
combination with a second diagnostic grid from which we derive the
[
/Fe] ratios (Mg2 vs.
Fe
,
see
Sect. 5). As a first step the
-enhancement for each
individual globular cluster is derived. This value is used to interpolate
the age/metallicity diagnostic grid for the correct
-element
enhancement. From so adjusted grids ages and metallicities are computed
using linear interpolation of the model predictions
employing a least-square technique. These ages and metallicities are then
used to adjust the [
/Fe] diagnostic grid, which itself is slightly
dependent on age. This procedure is iterated until the age, metallicity,
and [
/Fe] values converge. Extensive tests show that this goal is
reached within a few iteration steps. Along with these values we also
determine the statistical errors due to index measurement uncertainties.
Typically the errors are asymmetric which is a result of the skewness of
diagnostic grids (see Fig. 6).
For objects which fall outside the diagnostic grid, we assign the most
likely extreme grid value and do not extrapolate the model
predictions, i.e. we project the data onto the closest grid point along
its error vector. Note that due to the increasing influence of hot
blue horizontal-branch stars, the iso-age tracks for old (t>8 Gyr)
stellar populations with metallicities below -1.0 dex tend to overlap.
This introduces an ambiguity in assigning ages and metallicities to
individual globular clusters and artificially broadens the age
distributions at old ages. However, only six objects are affected by this
ambiguity which has no effect on the following results. For these six objects the routine randomly assigns either young or old ages and the corresponding metallicities.
The comparison of derived formal ages and metallicities of individual
globular clusters shows that different diagnostic diagrams make age
predictions which are somewhat inconsistent with one another. The average
relative age uncertainty varies in the range
which we attribute to systematic errors in the age scale of model
predictions and/or systematics in the calibration of the data. By
combining age predictions from different diagnostic plots we smooth out
these systematic effects which influence age determinations derived from
single Balmer indices.
This age inconsistency is partly surprising as the Maraston et al. (2003) models are well calibrated for Milky Way globular clusters and their ages derived from different Balmer indices are consistent with each other.
![]() |
Figure 6:
Age-metallicity diagnostic plots (panels a)- e)) constructed from
Balmer indices H
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
On the other hand, we find very good agreement of metallicity
predictions from all diagnostic plots within an average uncertainty of 0.2 dex.
Note also, that a striking feature of all age/metallicity diagnostic plots in Fig. 6 is that a significant fraction of globular clusters lie below the oldest iso-age track. This is observed in other samples of extragalactic globular clusters in the literature as well (e.g. Larsen et al. 2003; Cohen et al. 1998; Kissler-Patig et al. 1998; Puzia et al. 2000). Systematics in the data reduction have been deemed extremely unlikely given the good Lick standard calibration (see Paper I) and are also unlikely given the general problem in the literature. We could also exclude emission filling as source for the Balmer index inconsistencies which is presented in Sect. 3.5.
To test whether the spectroscopic age and metallicity predictions for these
outliers are consistent with their observed photometric colors (see Paper I for details), we use SSP models of Maraston (2005) to compute
photometric colors from the derived spectroscopic ages and metallicities.
Figure 5 shows the comparison of computed and observed
data. We find good agreement for the ,
,
and
colors. We find no systematics in color residuals as a function of Balmer index with respect to the one-to-one relation. However, we find significantly smaller
residuals towards more metal-rich objects.
We also investigate if higher-quality data alleviates these
inconsistencies and select a sub-sample with statistical errors
H
Å and
H
Å for
higher-order Balmer lines, which is shown as solid dots in
Fig. 5. For this sub-sample, we indeed find a better
agreement between computed and observed photometric colors. We take this
as an indication that the age and metallicity predictions for the
higher-quality dataset are more robust. Based on our photometric
cross-check, we only consider the ages and metallicities derived for our
high-quality sub-sample as trustworthy. However, to illustrate the
difference between the high-quality and the rest of the data we plot the
remaining dataset in the following diagnostic plots as open symbols in
Fig. 6 and show their distributions as hatched
histograms in Fig. 7.
In the following, we discuss ages and metallicities for the high-quality
sub-sample derived from diagnostic plots using the Balmer-line indices
H,
H
,
and H
(see
Sect. 4.2). The corresponding diagnostic diagrams are shown
in Fig. 6 along with the other two diagrams
constructed from the Balmer indices H
,
and H
.
Since for seven globular clusters one or more of the three Balmer indices H
,
H
,
and H
are not available, we exclude these objects from the subsequent analysis to avoid systematics, which leaves 17 globular clusters in our high-quality
sub-sample. Their age and metallicity distributions are shown in Fig. 7 as solid histograms.
Our method reveals that a significant fraction of clusters in the
high-quality sample (6/17 or 35%, see solid histogram in
Fig. 7) has formal ages between 5 and 10 Gyr. Only one object
has a formal age below 5 Gyr (3.8 Gyr) and there are no globular
clusters with extremely young ages (
1 Gyr). We split the sample at
[Z/H] = -0.6, corresponding to the dip in the Milky Way globular cluster
metallicity distribution (Harris 1996), into metal-poor and metal-rich
globular clusters. With this distinction we find an increase in age spread
from metal-poor to metal-rich clusters. We derive a dispersion of 0.5 Gyr
for metal-poor and a dispersion of 3.5 Gyr for metal-rich globular
clusters. A weighted linear least-square fit to our high-quality dataset
reveals a weak age-metallicity relation, in the sense that more metal-rich
globular clusters appear on average younger. The significance of this
relation needs to be tested with larger samples.
Our high-quality globular cluster sample covers metallicities in the range
-1.3
[Z/H]
+0.5, with a mean -0.36
0.13 and dispersion
dex. Four out of 17 globular clusters have formally
super-solar metallicities. All clusters that are younger than 10 Gyr have
metallicities [Z/H]
-0.4, which excludes their formation from
primordial gas clouds.
![]() |
Figure 7:
Histograms for ages, metallicities, and [![]() |
Open with DEXTER |
We note that no super metal-rich, intermediate-age counterpart globular cluster population is found in the Milky Way (e.g. Harris 2001). However, globular clusters with intermediate ages are found in M 31 (e.g. Puzia et al. 2005; Burstein et al. 2004; Beasley et al. 2004; Barmby et al. 2000), whose bulge-to-disk ratio is larger than that of the Milky Way. An interesting idea is that this mode of globular cluster formation may not have manifested itself in our Galaxy due to Milky Way's smaller bulge (Goudfrooij et al. 2003).
In summary, the conclusions that can be drawn from our previous analysis of globular clusters in early-type galaxies are the following:
In this Section we derive [/Fe] ratios for our high-quality sample globular
clusters using a diagnostic diagram which is least sensitive to
age/metallicity variations. Such a diagram can be constructed from the
indices
Fe
and Mg2, which primarily trace the
abundances of iron and the
-element magnesium (Tripicco & Bell 1995).
We note that among the three Mg-sensitive indices, Mg1, Mg2, and Mgb, theoretical index predictions show a relatively large spread in
Fe
and Mg2 for
[
/Fe] ratios between solar and
0.5 dex, at high mean
metallicities. Given the quality of our data we can expect a good
discrimination between enhanced and solar-type [
/Fe] ratio for
individual globular clusters at metallicities [Z/H]
-0.8 and a good
estimate of the mean [
/Fe] ratio at lower metallicities.
Figure 6 (panel f) shows that the Mg2 vs. Fe
diagnostic diagram is not entirely free from the age/metallicity degeneracy. Iso-[
/Fe] tracks for three different ratios (0.0, 0.3, and 0.5 dex) are plotted for two ages (3 and 13 Gyr, indicated by dotted and solid lines, respectively). It is obvious that age information is needed to choose the correct set of tracks for a reliable [
/Fe] determination. The diagnostic grid is interpolated in our
iterative fitting routine using information derived from age/metallicity
diagnostic grids (see Sect. 4.1). In other words, [
/Fe] ratios are determined simultaneously with ages and metallicities. The averaging is identical to the one for ages and
metallicities in Sect. 4.1.
A histogram of [/Fe] ratios for our high-quality sample is shown
in Fig. 7. All globular clusters are consistent with
super-solar [
/Fe] ratios. The mean [
/Fe] of the sample is
0.47
0.06 dex, with a dispersion of 0.26 dex.
Using a weighted linear least-square fit, we find evidence for a
[/Fe]-metallicity relation in the sense that more metal-rich
globular clusters have lower [
/Fe] ratios. However, due to the
reduced [
/Fe] resolution of Lick indices at low metallicities and
our modest sample size (especially at the metal-poor end), this trend
needs to be confirmed with more data. Consistent with the previously found
age-metallicity and [
/Fe]-metallicity relation, we find no
evidence for a [
/Fe]-age correlation. Globular clusters at all
ages appear to have on average super-solar [
/Fe] ratios.
Recent measurements of [/Fe] ratios in globular cluster systems in
other early-type galaxies reveal very similar results. For instance, the
data of Kuntschner et al. (2002b) for globular clusters in NGC 3115 show that
most clusters are consistent with [
/Fe]
0.3 over the
entire range of sampled metallicities. Larsen et al. (2002b) find a super-solar
mean [
/Fe] ratio of +0.4 dex for globular clusters of all
metallicities in NGC 4594 (Sombrero). Using SSP models with non-constant
[
/Fe] ratios, Forbes et al. (2001) argue that at least some globular
clusters in NGC 1399 exhibit super-solar [
/Fe] ratios. However,
the found [
/Fe]-metallicity correlation is seen in this sample for
the first time. This is mainly due to the higher quality of our data and a
more reliable age/metallicity determination compared to previous studies.
Note that super-solar [/Fe] ratios have also been reported for the
diffuse light of early-type galaxies (e.g. Kuntschner et al. 2002a; Thomas et al. 2005; Trager et al. 2000b; Davies et al. 2001,1993). In subsequent papers of
this series we will compare the index measurements of globular clusters
with indices derived from our long-slit spectra of the hosts' diffuse
stellar light.
In the following we summarize the major findings of this section:
![]() |
Figure 8: Top row: comparison of globular cluster ages and metallicities in our sample galaxies with those in a simulated host galaxy with a similar luminosity (mass). All globular cluster in our sample are shown as hatched histograms, outliers in diagnostic plots are indicated by double-hatched histograms (see Sect. 4.2), and high-quality data is shown as solid histograms. Observed and simulated globular cluster samples were split into metal-poor and metal-rich globular clusters at [Z/H] = -0.9 and the shading is accordingly transfered to the age distribution plot. Bottom row: average ages and metallicities of simulated globular cluster systems in faint and bright host galaxies. All predictions were taken from Beasley et al. (2002a). |
Open with DEXTER |
If the dispersion in Balmer-[MgFe]
diagnostic plots is entirely
driven by age
, we find
indications for a significant fraction of relatively young globular
clusters in early-type galaxies. About
1/3 of our sample globular
clusters have formal ages younger than 10 Gyr, implying formation
redshifts
(in a
CDM concordance universe).
In Fig. 8 we compare the derived ages and metallicities
of our high-quality sample with the corresponding distributions of
globular cluster systems simulated according to a hierarchical merging
scenario (Beasley et al. 2002a). As the average luminosity of our sample host
galaxies is
0.6 (see
Table 1), we use the globular cluster system of a typical
simulated galaxy with a similar luminosity
MB=-20.42 (see also
Fig. 6 in Beasley et al. 2002a) for comparison. We split both cluster samples
(simulated and observed) in metal-poor and metal-rich clusters at
[Z/H] = -0.9. The top right panel of Fig. 8 shows that
our sample is biased towards metal-rich globular clusters, and the
comparison of metal-poor globular clusters is limited to the statement
that both observed and simulated clusters are all consistently old. For
metal-rich globular clusters there is good agreement in the age range
covered by data and simulations.
The somewhat naive prediction of more extended formation timescales
of more massive structures in the hierarchical picture holds only if the
fraction of gas-poor to gas-rich mergers, the so-called dry to mixed
merger ratio, is constant with redshift. Khochfar & Burkert (2003) predict that
this ratio depends on galaxy mass in the sense that most massive
ellipticals formed early in rather dissipationless (dry) mergers of
bulge-dominated precursors (see also Kauffmann & Haehnelt 2000). This would
imply that on average low-luminosity ellipticals would harbor a higher
fraction of young globular clusters, while most massive galaxies would
preferentially host old globular cluster systems. This is reflected in the
simulations of Beasley et al. (2002a). The lower panels of
Fig. 8 illustrate the mean ages and metallicities of
globular cluster systems in faint and bright galaxies with a cut at
MB=-21.0. The plots show that while host galaxy luminosity affects
the mean metallicity of globular cluster systems only marginally, a
significant fraction (25%) of low-luminosity galaxies is predicted
to host globular cluster systems with an average age younger than
10 Gyr, compared to a tiny fraction (
2%) in bright galaxies. Unfortunately, the transition between dry and mixed merger-dominated evolution is predicted to occur for early-type galaxies in the luminosity
range between
and -21 mag, where all of our sample galaxies reside.
To first order, previous spectroscopic studies find evidence for mostly old globular cluster systems in the massive Fornax galaxy NGC 1399 (Kissler-Patig et al. 1998), and the two Virgo galaxies M 87 (Cohen et al. 1998) and M 49 (Cohen et al. 2003), which appears to fit into the framework of dissipationless merging with an early assembly. We will get back to this point in a subsequent paper in which globular cluster systems will be investigated as a function of the host galaxy property.
Stellar populations with super-solar [/Fe] ratios, as they are
observed in massive elliptical galaxies, are generally interpreted as the
result of very short formation time scales. However, in the case of
dissipative merging, hierarchical merging models predict frequent merging
events, which is expected to result in relatively low [
/Fe] ratios
(e.g. Thomas et al. 1999). Although somewhat on the high side, our formal
mean [
/Fe] ratio of
0.47 dex for globular clusters in
early-type galaxies is in line with values measured for the diffuse light.
Hence, the formation time scales of field stars and globular clusters
appear to be similar in elliptical galaxies and not immediately compatible
with hierarchical merging models.
The found evidence for a [/Fe]-metallicity relation in early-type
globular cluster systems is consistent with a generic chemical enrichment
scenario in which more metal-rich stellar populations have lower
-element enhancements. Two types of supernovae contribute to the
build-up of
-peak (type II) and iron-peak elements (type Ia) in
the interstellar medium on different time scales
(e.g. Matteucci 1994; Thomas et al. 1999), because of the evolutionary delay
of
1 Gyr of type-Ia supernova progenitors (Greggio 1997). Since
globular clusters are not massive enough to support significant
self-enrichment, their [
/Fe] ratios reflect the large-scale
chemical conditions during their formation. However, to allow a detectable
increase in [
/Fe], for high-metallicity stellar populations it
requires a larger number of type II supernovae to outweigh any previous
metal-enrichment by type Ia supernovae compared to lower metallicities.
This implies that high [
/Fe] ratios at high metallicities are
likely the result of shorter formation timescales and/or higher
star-formation rates than similar [
/Fe] ratios at lower
metallicities. Careful modeling of the chemical evolution of a globular
cluster system is needed to quantify these formation timescales. The
relation between [
/Fe] and metallicity found here provides the
first important constraint for this exercise.
It is also important to test whether a steeper [/Fe]-age relation
can be found for globular cluster systems in more massive galaxies, which
would be expected if these systems formed on more extended timescales, as
predicted by hierarchical merging models.
We have conducted a study of ages, metallicities, and [/Fe] ratios
of extragalactic globular clusters in early-type galaxies, based on the
Lick index system. We find that up to
1/3 of our sample of globular
clusters have ages formally younger than 10 Gyr. This result is not
biased by one single globular cluster system in our galaxy sample and
appears representative for early-type galaxies in general. We cannot state
with confidence whether the younger ages (i.e. t < 10 Gyr) are real or
due to an unexpected blue horizontal branch morphology at high
metallicities. If the younger ages are real, the found fraction should be
taken as an upper limit, since our data only probe the bright end of the
globular cluster luminosity function where relatively young globular
clusters, if present, are expected to reside. In general, less than
1/3 of the brightest
10% of globular cluster systems in early-type galaxies which we sample with our study could have formed at redshifts
.
We find that the formal age scatter increases and the mean age decreases from metal-poor to metal-rich globular clusters, resulting in an age-metallicity relation. For our high-quality sample globular clusters with metallicities [Z/H] < -0.6 we find an age dispersion of 0.5 Gyr. The metal-rich sub-sample has a dispersion of 3.5 Gyr.
Our high-quality sample spans a wide range in metallicity between
[Z/H]
-1.3 and
+0.5 dex. We find evidence for an
[
/Fe]-metallicity relation in the sense that more metal-rich
globular clusters have lower
-element enhancements, which needs to
be confirmed with more data that better sample the metal-poor regime.
However, there is no indication for an [
/Fe]-age correlation.
[/Fe] ratios are found to be on average super-solar with a mean
0.47
0.06 dex and a dispersion of
0.3 dex, which indicates
formation timescales shorter than
1 Gyr. In other words, the
progenitor clouds of globular clusters were predominantly enriched by
type-II SNe.
Acknowledgements
We are grateful to Michael Rich for providing colour magnitude diagrams of M 31 globular clusters prior to publication. Many thanks go to Michael Beasley for sending electronic tables of his hierarchical-clustering simulations. We also thank Scott Trager for a very constructive referee report. T.H.P. gratefully acknowledges the support by the German Deutsche Forschungsgemeinschaft, DFG project number Be 1091/10-2, and the support in form of an ESA Research Fellowship.