In this section we use the single-age single-metallicity stellar
population models from Vazdekis et al. (1996)
to infer luminosity weighted mean ages, metallicities, and abundance ratios
from the data.
We use the models with a bi-modal IMF with high mass slope
.
It is not surprising that
from a single observational quantity, such as (B-r),
it is not possible to determine both age and metallicity
(e.g. Worthey 1994).
For example, a color of
,
which is typical for
E and S0 galaxies,
can be matched by both
a low age and high metallicity, say 3 Gyr and Z = 0.05,
an intermediate age and solar metallicity, say 8 Gyr and Z = 0.02,
and a high age and low metallicity, say 13-17 Gyr and Z = 0.008.
The numbers quoted are from the Vazdekis et al. models.
Given two observational quantities we could hope to determine both the age and metallicity. Unfortunately, for a number of such color-color, color-index, and index-index diagrams the effects of age and metallicity are nearly degenerate (e.g. Worthey 1994, Faber et al. 1995).
It turns out, that in the two-dimensional
-
and
-
diagrams
the effects of age and metallicity
are not degenerate.
This was mentioned by Faber et al. (1995).
We assume homology and are then able to calculate
from the data
(Eq.
; cf. Sect.
, p.
).
Our data are plotted in the
-
and
-
diagrams
in Fig.
.
Overplotted are the predictions from the Vazdekis et al. models.
The first thing to note is that the models span the data quite well,
the measurement errors taken into account.
We are free to shift the values of
up and down since they depend on the two unknown quantities
H0 and the fraction of dark matter.
The used values of
and
give a good match to the data,
and we do not apply any offset.
The Vazdekis et al. models predict
,
,
and
for a grid of 45 (age,[M/H]) points,
with the 15 age values ranging from 1.00 to 17.38 Gyr, and the three
[M/H] values being -0.4, 0.0, and 0.4.
Recall that
,
with
.
Note, that not all the age values are shown on Fig.
.
The inverse problem, i.e. given observed values of
and
(or
and
)
determine the corresponding age and metallicity,
is a question of interpolation in an irregular grid.
That the grid is irregular is apparent in Fig.
.
We did this irregular interpolation in three steps.
First, a Delaunay triangulation of the irregular grid points
was established.
Second, a large
(say 1000
1000) regular grid of interpolated or extrapolated
values of both
and [M/H] was calculated,
with the grid being chosen to
span the data. The interpolation and extrapolation was done using the
Akima's quintic polynomials.
Third, a standard bilinear interpolation was used to get the final
values of
and [M/H].
These calculations were done using IDL (Interactive Data Language).
For details, see the help pages for triangulate and trigrid.
Estimates of uncertainties on
and [M/H]
were obtained by in turn keeping
and
(or
)
constant
while varying the other by plus/minus the observational error, calculating
and [M/H] for those four points, and taking half the
min-max variation as the estimate of uncertainty.
For the interpolation in the
-
diagram,
we omitted the three galaxies with an uncertainty on
larger than
0.065.
To discuss the results from the interpolation in
the
-
and
-
diagrams,
we introduce the following notation
The above notation indicates our first order assumptions:
we assume that the metallicity inferred from
gives the magnesium
abundance [Mg/H], and that
the metallicity inferred from
gives the iron
abundance [Fe/H].
This is despite the fact that
the Vazdekis et al. models have solar abundance ratios,
including [Mg/Fe] = 0.
That these are reasonable approximations is supported by the
work of Tripicco & Bell (1995), and Weiss et al. (1995),
as described in
Sect.
(p.
).
In summary,
Tripicco & Bell (1995) found
that the
index depended strongly on the magnesium abundance,
and that the
index depended strongly on the iron abundance,
although is was just as sensitive to changes in the total metallicity.
Weiss et al. (1995) found from isochrones with [Mg/Fe]
0
that the effect on e.g. the luminosity
of changing the abundance ratios
while keeping the total metallicity constant was small.
If the models provided an adequate description of the data,
the metallicity difference
and the age difference
should be zero within their errors caused by the observational errors.
However, this is not the case, and
we already suspected this discrepancy from the
-
diagram, Fig.
(p.
).
The metallicity difference
is large compared with the two metallicities, typically 70%,
while the age difference
is small compared with the two ages, typically only 5%.
We regard the differences in ages as an indication of the
limitations in the method, and use
our
as an estimate of the true
.
Note,
that
has to be non-zero
since
is non-zero, because
appears in both diagrams.
The Vazdekis et al. models for ages > 5 Gyr can be
approximated by
(J97; Eq.
).
Applying this to the two diagrams, we get
,
or
.
A fit to the combined HydraI+Coma sample (N=83) gives
,
in agreement with this.
The scatter is small,
.
Note, that
since
,
,
and the ages all depend on
,
and since we can only determine
to within an offset (cf. above),
these four quantities are also only determined to within an offset.
For
,
on the other hand, the effect of this unknown offset cancels out,
at least to some extend.
In Figure
we show histograms over the distribution of
,
,
,
,
,
and
.
is included just to show that also the absolute range in
in this quantity is rather small.
We tested the HydraI and the Coma data against each other by means of
Kolmogorov-Smirnov tests.
This test gives the probability
that the two samples are drawn from the
same underlying distribution.
For the HydraI sample (N=41) and
the Coma sample with
data (N=42)
we find
in the range 26%-98%
for the above six quantities.
In other words, we do not find any significant differences between
the HydraI and the Coma samples.
For the HydraI sample and
the full Coma sample (N=111)
we find
% for
and 41% for
;
again we find no significant differences between
the HydraI and the Coma samples.
Properties of E and S0 Galaxies in the Clusters HydraI and Coma
Master's Thesis, University of Copenhagen, July 1997
Bo Milvang-Jensen (milvang@astro.ku.dk)