gas cooling plays a key role in galaxy formation . according to the hierarchical structure formation model , galaxies form in gravitationally collapsed dark matter halos .
the dissipative baryonic matter cools by emitting radiation and condenses into halos to form galaxies @xcite .
subsequently stars form through the radiative cooling of interstellar medium ( ism ; e.g. , * ? ? ?
the cooling rate of gas varies with its chemical abundance , because different atoms and molecules have different cooling rates .
the chemical abundance of cosmic gas evolves as a function of time and environment . in the early universe , the primordial gas consists of mainly hydrogen ( mass fraction of @xmath3 ) and helium ( mass fraction of @xmath4 ) with a negligible amount of light metals ( e.g. * ? ? ?
* ) . throughout this paper ,
all the elements heavier than helium are collectively called ` metals ' .
most metals are produced inside stars and spread out to the space by galactic winds driven by supernova ( sn ) explosions . in the metal - enriched gas , the total cooling rate
is significantly enhanced relative to that of the primordial gas by the atomic emission lines from recombination of ionised metals .
metals may influence galaxy formation in two different ways .
firstly , the sf efficiency is increased owing to the shorter gas cooling time .
star formation generally takes place in dusty metal - enriched ism , therefore , the enhanced cooling rate by metals would boost up the sf efficiency .
however , the energy / momentum feedback by sne may suppress the subsequent star formation after the initial starburst by heating up the ambient gas .
this complicates the situation , and it is not clear whether the net effect of metal cooling and sn feedback would be negative or positive .
the feedback process is highly nonlinear , therefore , a direct numerical simulation would be a useful tool to explore the effects of feedback and metal cooling .
secondly , the metals dispersed into the igm by sn feedback also enhance the cooling of igm .
this can lead to the increase of igm accretion onto galaxies , because colder gas can sink into galaxies more easily . both of the above two effects may affect the galaxy growth and cosmic sf history significantly .
therefore it is essential to include the effects of metal cooling and chemical enrichment by feedback in the studies of galaxy formation and evolution , and capture the complex two - way interactions between galaxies and igm .
in addition , metal enrichment also affects the equation of state by altering the mean molecular weight of gas , which influences the hydrodynamic calculation of gas thermal state .
cosmological hydrodynamic simulations are widely used in the studies of galaxy formation and cosmic sf history ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
most of these works included the treatments of radiative cooling by h & he , star formation , and sn feedback .
however , some of the simulations did not include the effects of metal cooling and/or chemical enrichment by galactic wind .
furthermore , the effects of metal cooling on galaxy growth and cosmic sf history has not been explored in detail using cosmological hydrodynamic simulations and presented in the literature , as it is costly to run large cosmological simulations with and without the effect of metal cooling .    in this paper
, we investigate the effects of metal cooling and metal enrichment on galaxy growth and cosmic sf history using a series of cosmological hydrodynamic simulations with and without metal cooling .
the aim of this paper is to single out the effects of metal cooling among our simulations .
this work can be regarded as our initial step towards the long - term goal of developing more complete cosmological hydrodynamic code with physically motivated models of star formation and feedback .
the paper is organised as follows . in
[ sec : method ] , we describe our simulation method focusing on how we implement the metal cooling . in   [ sec : global ] , we study the metal cooling effects on the global properties , such as cosmic sfr , phase space ( @xmath5 vs. @xmath6 ) distribution of gas , and the evolution of four phases ( hot , warm - hot , diffuse , and condensed ) of baryons .
we then study the galaxy mass functions (   [ sec : mf ] ) and gas mass fractions (   [ sec : gasfrac ] ) .
finally , we discuss and summarise our findings in   [ sec : summary ] .
we use the updated version of the tree - particle - mesh ( treepm ) smoothed particle hydrodynamics ( sph ) code gadget-2 @xcite for our cosmological simulations .
the gravitational dynamics is computed by a treepm algorithm , which uses a particle - mesh method @xcite for the long - range gravitational force and a tree method for the short - range gravitational force @xcite .
this hybrid algorithm makes the gravitational force calculation faster than a tree method and allows better force resolution than a pm method in dense regions .
the gas dynamics is computed by an sph method .
the sph is particularly useful if the simulation needs to resolve large dynamical range , which is an inevitable requirement for the study of galaxy formation in a cosmological context .
therefore a treepm - sph simulation can provide a fast and high resolution calculation for both gravitational dynamics and hydrodynamics .
the gadget-2 code adopts the entropy - conservative formulation @xcite , which alleviates the overcooling problem that previous sph codes suffered from .
our basic simulations include radiative cooling and heating processes for hydrogen and helium using a method similar to @xcite .
an external uv background radiation is treated as a spatially uniform photoionising radiation @xcite , and modified to match the ly@xmath7 forest observations @xcite .
implementing star formation and feedback from first principles is not feasible in current cosmological simulations , because the spatial & mass - scales of molecular clouds are not resolved .
star formation and sn feedback are represented by the subgrid multiphase ism model developed by @xcite . in the multiphase scheme ,
a single gas particle represents both hot and cold gas .
the stars are formed in the cold portion when the density exceeds a given threshold , @xmath8 , which is derived self - consistently within the multiphase ism model .
the related parameters , such as the normalisation of gas consumption time and evaporation efficiency of cold gas , are set to satisfy the empirical kennicutt - schmidt law @xcite .
the sn feedback returns some fraction of the cold gas to the hot phase , and increases the thermal energy of the hot gas .    as an extension to the multiphase ism model ,
the simulation includes a phenomenological model for sn - driven galactic wind @xcite .
the galactic wind is particularly important for distributing the metals produced by sne into the igm .
we use the strong kinematic wind with a velocity of 484kms@xmath9 .
it has been shown that this model produces favourable results for the luminosity function of lyman - break galaxies at the bright - end @xcite and the hi column density distribution function @xcite at @xmath1 , when compared to the runs without the wind .
however , @xcite pointed out the problems of this galactic wind model by comparing with the observations of civ absorption lines in quasar spectra , and suggested that the momentum - driven wind model @xcite is a more viable model , which can carry more metals with lower wind velocities . in this paper
we choose not to modify our galactic wind model in order to single out the effects of metal cooling , and to allow direct comparisons to the previous works @xcite .
here we focus on the effects of metal cooling on galaxy growth , while @xcite focused on the civ statistics of the igm .
the metallicity of gas particles are also tracked by the code , assuming a closed box model for each gas particle .
the yield ( @xmath10 ) of 0.02 is assumed . in principle
, there could be a time delay between sn explosions and chemical enrichment of the ambient gas .
unfortunately , current cosmological simulations do not have sufficient resolution to track the detailed mixing process of metals .
in our simulations , we ignore this time delay and assume an instantaneous mixing within each gas particle .
this assumption is reasonable , because the overall time step of the simulations is longer than the mixing time - scale on small scales .
ideally we would like to track individual metals and compute the cooling rate of each element at each time - step of the simulation .
this approach has been used to study the evolution of individual galaxies @xcite . for cosmological simulations
, it requires a large memory to track individual metals and the computation becomes exceedingly expensive ( but see * ? ? ?
* ; * ? ? ?
* for such an attempts ) .
our simulations collectively track the total mass fraction of metals ( @xmath11 ) as the measure of metallicity .
we obtain the chemical abundance and cooling rates for a given @xmath12 and temperature by interpolating the table given in ( * ? ? ?
* hereafter sd93 ) for the standard collisional ionisation equilibrium model , and create a lookup table .
metal cooling is implemented in the full range of metallicity in the table from sd93 , from [ fe / h]=-3 to [ fe / h]=1 .
in addition , the metal cooling outside of this range is computed by extrapolating the table .
we require that the extrapolation does not go lower than the primordial cooling rate .
the primordial abundance pattern is used where [ fe / h ] @xmath13 @xcite , and the solar abundance pattern is used where [ fe / h ] @xmath14 @xcite . for @xmath15 [ fe / h ] @xmath16 ,
the abundance pattern is computed by interpolating between the primordial and solar abundance patterns .
note that , [ fe / h ] = @xmath17 - @xmath18 , where @xmath19 is the number density .
figure [ fig : xyz_table ] shows the relationships between the metal parameters ( [ fe / h ] , @xmath20 , @xmath21 , and @xmath22 ) and the total metallicity @xmath12 . because the abundance pattern is fixed for a given @xmath12 , [ fe / h ] is a monotonically increasing function of @xmath12 .
the main difference between the solar and the primordial abundance pattern is the oxygen or @xmath7-particle enhancement . in the early phase of metal enrichment ,
type ii sne are the main source of metals and they tend to generate more @xmath7-particles than sn ia .
therefore , it is expected that a low metallicity gas to have more @xmath7-particles than the high metallicity one .
as mentioned above , we use two different chemical abundance pattern ( at low- and high - end of the @xmath12 values ) to model this evolution in chemical composition .
we compute the oxygen abundance based on the table 4 of sd93 . in our code
, we incorporate the metal cooling effect by adding the additional cooling contribution from metals on top of the primordial cooling rate , which is computed in a similar fashion as described in @xcite .
metal enrichment changes the mean molecular weight @xmath23 , which is needed to compute the temperature and internal energy of gas . in order to compute @xmath23
, we need to know the electron number density @xmath24 , which varies with metallicity .
we obtain the values of @xmath24 from sd93 in the same way as we obtained the metal cooling rates .
our electron number density estimate ignores the electrons from photoionised metals by the uv background radiation , but we expect that the number of electrons from ionised metals will be small .    varying the mean molecular weight @xmath23 also influences the sf threshold density @xmath8 in the multiphase ism model of @xcite , as well as the cold and hot gas fractions of multiphase gas particles . in this paper , we consider two different treatments of @xmath8 . in the ` constant @xmath8 ' model ,
we take the value of @xmath8 as in the original formulation by @xcite and do not modulate @xmath8 by the varying @xmath23 . in the ` varying @xmath8 ' model , we take into account of the modulation of @xmath8 by the change in @xmath23 .
we use three series of simulations with varying box size and resolution .
the effects of metal enrichment and metal cooling are tested by comparing the runs listed in table  [ tbl : simulation ] .
the effects of galactic wind are also tested by turning on and off the wind .
the adopted cosmological parameters of all simulations are @xmath25 , which mostly concur with the results of the wilkinson microwave anisotropy probe ( wmap1 ) @xcite .    [ cols="<,^,^,^,^,^,^,<,^",options="header " , ]     our naming conventions for the simulations are as follows .
the first part of the run name denotes the particle number and box size used for the simulation series .
the run with the extension of ` mc ' ( for ` metal cooling ' ) adopts the ` constant @xmath8 ' model for metal cooling , the run with the extension ` mv ' ( for ` metal , varying ' ) adopts the ` varying @xmath8 ' model , the run with the extension ` nw ' has no galactic wind , and the run without any extension uses the cooling rates for the primordial chemical composition ( i.e. , only h and he ) .
the ` n216l10 ' series is the fiducial simulation set in our study .
we implement one simulation without galactic wind ( ` n216l10nw ' run ) .
the ` n216l10 ' run has the same physical models as the ` q4 ' run used in @xcite and @xcite with no metal cooling .    to examine the resolution effect ,
we implement lower resolution simulations , the ` n144l10 ' series .
owing to the comparably small box size and missing long wavelength perturbations , we evolve the n144l10 and n216l10 series only down to @xmath26 . to evolve the simulation toward lower redshifts ,
we need a simulation with a larger box size .
therefore we implement the ` n288l34 ' series , which has a larger box size than other series , but a lower resolution .
this series is evolved down to @xmath0 .
the runs within the same series use the identical initial condition , which makes the comparison more robust and enables halo - by - halo comparison if needed .      in cosmological simulations ,
dark matter , gas , and stars are represented by particles , and these particles are monte carlo representations of matter distribution . in order to study galaxy formation using a cosmological simulation ,
we need a working definition for simulated galaxies . in this paper
the simulated galaxies are defined as isolated groups of star and gas particles , identified by a simplified variant of the algorithm developed by @xcite . in more detail
, the code first computes a smoothed baryonic density field to identify candidate galaxies with high density peaks .
the full extent of these groups are found by adding gas and star particles to the groups in the order of declining density . if all @xmath27 is the minimum number of gas and star particles that constitute one isolated group . in this paper
we set @xmath28 .
] nearest neighbour particles have lower densities , this group of particles is considered as a new group .
if there is a denser neighbour , the particle is attached to the group to which its nearest denser neighbour already belongs to .
if two nearest neighbours belong to different groups and one of them has less than @xmath27 particles , these the two groups are merged .
if two nearest neighbours belong to different groups and both of them has more than @xmath27 particles , the particles are attached to the larger group , leaving the other group intact .
in addition , the gas particles in groups should be denser than @xmath29 . in this paper
we are not concerned with groups without star particles .
the fraction of such groups are very small , therefore they do not affect our conclusions .
% by metal cooling throughout the entire redshift range . ]    , but for the n288l34 series , which goes down to @xmath0 .
the n216l10mc run is also shown for a resolution comparison .
the green dot - dashed line ( hs03 ) is for the analytic model of ( * ? ? ?
* eq.(2 ) ) . ]
figure  [ fig : sfr_n216l10 ] shows the effect of metal cooling and galactic wind on the cosmic sf history for the n216l10 series .
the comparison of n216l10mc and n216l10 runs clearly shows that metal cooling enhances the cosmic sfr throughout all redshifts by @xmath30% .
previously , ( * ? ? ?
13 ) argued that the sfr density is enhanced by metal cooling mostly at lower redshifts and hardly at high-@xmath31 , however , our simulations suggest otherwise .
( see   [ sec : summary ] for more discussion on this point . )
metallicity needs to be higher than @xmath32 to noticeably increase the gas cooling rate , because the metal cooling rate is much smaller than that by h and he at @xmath33 .
the metallicity of diffuse igm is generally lower than @xmath34 at @xmath35 , and the metal cooling effect on igm is small .
however the local sfr could be enhanced significantly even at high-@xmath31 when the ism is enriched to @xmath36 , which can be easily achieved after a few events of sn ii .
our results suggest that star formation at high-@xmath31 can also be enhanced by metal cooling , and not just at low-@xmath31 .
we also find that the peak of the cosmic sf history hardly shifts by the introduction of metal cooling , consistently with the estimate by @xcite .
the n216l10mv run has a lower sfr at @xmath37 and a higher sfr at @xmath38 than the n216l10mc run .
this is because at high-@xmath31 , the value of @xmath8 is higher in the ` mv ' run than in the ` mc ' run due to lower metallicity and mean molecular weight , leading to a larger fraction of gas not being able to satisfy the sf criteria in the n216l10mv run than in the n216l10mc .
the opposite is true at @xmath39 .
figure  [ fig : sfr_n216l10 ] also shows that the sfr in the n216l10nw run is higher than other runs by an order of magnitude .
this is because the gas can cool very efficiently without being ejected by the galactic wind .
although the current galactic wind models in cosmological simulations are not well established yet and need to be improved , both theoretical and observational studies evidently suggest the important role of galactic wind feedback @xcite .
as expected , the n144l10mc run has lower sfr densities at @xmath40 compared to the n216l10 series .
this is because a lower resolution run is unable to resolve low - mass halos that would otherwise forms stars in them .
when the star formation at high-@xmath31 is underestimated , some of the gas remain unconverted into stars until lower redshift , resulting in a higher sfr at @xmath41 in the n144l10mc run than in the n216l10mc run .
figure  [ fig : sfr_n288l34 ] shows the cosmic sf history down to @xmath0 for the n288l34 series .
again , the sfr density in the n288l34mc run is systematically higher than that in the n288l34 run at all redshifts .
moreover , the difference in sfr density between the two runs increases at @xmath39 .
by @xmath42 the igm metallicity becomes high enough to noticeably increase the cooling rate , and the accretion of igm onto galaxies becomes more efficient , which later fuels the star formation .
our results suggest that the cosmic sfr density is enhanced by metal cooling at @xmath43 by both of the following two effects : 1 ) more efficient igm accretion onto galaxies due to higher metallicity of igm , and 2 ) increased sf efficiency due to higher cooling rate with the contribution from metal cooling .    in figure
[ fig : sfr_n288l34 ] the analytic model of ( * ? ? ?
* eq.(2 ) , hereafter h&s model ) is also shown in the green dot - dashed line .
this model is a result of many cosmological sph simulations with increasing resolution , and the authors gave a theoretical interpretation to the empirical fitting formula for the cosmic sfr density . since our runs adopt a lower value of @xmath44 than the original value used to derive the h&s model ( @xmath45 ) , we modify the h&s model by rescaling the @xmath46 parameter in their eq.(2 ) and eq.(47 ) .
this reduces the sfr density at @xmath47 by @xmath48% , but not much at @xmath49 .
the results of the n288l34 and n216l10mc runs agree well with the h&s model at @xmath41 and @xmath50 , respectively .
this is expected , because the n288l34 run adopts the same physical models as in the original simulations that were used to derive the h&s model . at high-@xmath31
, both runs fall short compared to the h&s model owing to insufficient resolution .
the metal enrichment and metal cooling also change the phase space distribution ( @xmath51 diagram ) of cosmic gas .
figure  [ fig : phase.n216l10 ] shows the state of cosmic gas at @xmath1 for the n216l10 series . as shown in panel ( @xmath52 ) , the baryons in the universe can be broadly categorised into four different phases according to their overdensity and temperature : ` hot ' , ` warm - hot ' , ` diffuse ' , and ` condensed ' @xcite .
the first two phases are mostly the shock - heated gas in clusters and groups of galaxies .
the ` diffuse ' phase is mostly the photoionised igm with lower temperature , which can be observed as the ly@xmath7 forest .
the tight power - law relationship between @xmath5 and @xmath6 at @xmath53 and @xmath54k is governed by the ionisation equilibrium @xcite , where @xmath55 is the mean density of baryons at @xmath56 .
the ` condensed ' phase is the high density , cold gas in galaxies . because the primordial cooling curve has a sharp cutoff at @xmath57k , the temperature of condensed phase decreases only down to this temperature in figure  [ fig : phase.n216l10]a , b .
the characteristic finger - like feature at @xmath58 is due to the pressurisation of the star - forming gas by sn feedback in the multiphase ism model of @xcite .
a fitting formula for this effective equation of state was derived by @xcite , which is shown by the blue dashed lines . above the sf threshold density @xmath8
, the gas particles in the simulation go into the multiphase mode and are allowed to form stars .
the most noticeable and interesting change in the @xmath51 plot by the metal cooling is in the distribution of ` condensed ' gas .
star formation and sn explosions occur in the multiphase star - forming gas , therefore the immediate effect of metal cooling appears in the condensed phase . the presence of gas at @xmath59k with metal cooling results from the effective temperature calculation for the multiphase medium , and not from the direct gas cooling down to @xmath59k . according to the multiphase scheme in @xcite , the star forming gas consists of two phases : hot gas ( @xmath60k @xmath61 @xmath62k ) and cold gas ( @xmath63k ) .
the effective temperature for the multiphase gas is computed from the internal energy of hot and cold gases , which are derived using the cold gas fraction .
the cold gas fraction is computed from the cooling rate of hot gas and the gas density .
the original gadget imposes the effective temperature of @xmath64k at the sf threshold density . in our metal cooling scheme ,
however , the cooling rate of hot gas increases , and the cold gas fraction increases .
therefore , the multiphase gas above the sf threshold density with higher metallicity could have lower effective temperature of @xmath59k in our implementation .    in the n216l10mc run ,
the value of @xmath8 was kept fixed as the original value , which causes a sharp edge at @xmath65 for the multiphase gas . in the n216l10mv run , the value of @xmath8 was varied according to the change in mean molecular weight , therefore the condensed gas has a spread in both temperature and density near @xmath65 .
galactic wind returns the metal - enriched , star - forming gas to the hot and low density igm , as evidenced in the @xmath66 diagrams .
this metal - enriched , low - density ( @xmath67 ) gas is absent in the n216l10nw run , which clearly shows the necessity of galactic wind to explain the metallicity observed in ly@xmath7 forest .
the phase distribution of igm in the n216l10 , n216l10mc , and n216l10mv runs are very similar .
this is presumably because the effects of metal cooling is not so significant yet in the igm at @xmath68 . as we will discuss in the next section ,
this is not be the case at @xmath2 .
the connection between the four phases of cosmic gas is important in galaxy formation , and it is beneficial for us to study their evolution for a better understanding of metal cooling effects .
a conventional view is that the diffuse gas is shock - heated to hot or warm - hot phase , which later cools down to the condensed phase @xcite .
@xcite also suggested that some diffuse gas can migrate to the condensed phase without being shock - heated , which they called the ` cold mode ' accretion .
figure  [ fig : evolv4 ] shows the evolution of mass fractions of different phases .
all of our simulations show that the diffuse phase continuously decreases from high-@xmath31 to low-@xmath31 , while the ` warm - hot + hot ' and ` condensed ' phases continuously increase , consistently with the general expectations .
the igm ( i.e. , hot + warm - hot + diffuse ) is accreted onto galaxies , to become the condensed phase .
the diffuse component ( panel @xmath52 ) decreases from @xmath69% at @xmath70 to @xmath71% at @xmath0 . the ` warm - hot + hot ' component ( panel @xmath72 ) increases from almost zero at @xmath70 to @xmath73% at @xmath0 .
the mass fraction of the total igm ( panel @xmath74 ) decreases from @xmath75% at @xmath70 to @xmath76% at @xmath1 .
this reflects the increase of the condensed phase ( panel @xmath77 ) from a few percent at @xmath70 to @xmath78% at @xmath1 .
these mass fractions are in general agreement with the results of @xcite .
the comparison between n288l34 and n288l34mc runs shows interesting differences . at @xmath68 ,
the mass fractions of all four phases are similar in the two runs . however , at @xmath2 , the igm mass fraction continues to decrease in the n288l10mc run , while it becomes almost constant in the n288l10 run .
a similar behaviour is seen for the condensed phase , and the n288l10mc continues to increase its condensed mass .
this suggests that the metal cooling enhances the igm accretion onto galaxies mostly at @xmath39 , which supplies the fuel for star formation .
the cooling of warm+hot phase into diffuse phase is also enhanced , leading to a larger mass fraction of diffuse component in the n288l34mc run than in the n288l34 run .
( the same is true for n216l10mc and n216l10 runs . )
the n216l10nw run shows stark differences from all other runs .
it has a significantly lower igm mass fraction , and a very large mass fraction in the condensed phase .
this shows that the igm is mostly enriched by galactic wind , and that most of the masses are returned into the ` warm - hot + hot ' phase , but not to the diffuse phase .
galaxy stellar mass function ( gsmf ) tells us how the stellar mass is distributed in different galaxies , and its redshift evolution reflects the growth of structure in the hierarchical universe .
figure  [ fig : mf_star ] compares the gsmf in different simulations .
we only show the range of @xmath80 , which corresponds to the limiting mass of 32 star particles in n216l10 runs .
we chose 32 star particles as our resolution limits because all mass function start to turn over around this mass range .
panel  ( @xmath52 ) shows that the run with no galactic wind ( n216l10nw ) completely overestimates the gsmf , and the run with lower resolution ( n144l10 ) underestimates the number of lower mass galaxies at @xmath81 relative to the n216l10 run .
panel  ( @xmath72 ) shows that the metal cooling increases the masses of galaxies by @xmath82% at @xmath1 . as a result ,
both the n216l10mc and n216l10mv runs show a slight increase in the number of galaxies at the massive end , with a slightly stronger enhancement in the n216l10mv run .    figures  [ fig : mf_star]c , d demonstrate that the metal cooling also changes the redshift evolution of gsmf
. in the n216l10mc run , the peak of mass function is shifting more toward the massive side than in the n216l10 run at @xmath1 . from @xmath83 to @xmath1 , the low - mass galaxies merge to form more massive galaxies , increasing the number of galaxies with @xmath84 furthermore .
it also seems that the formation of low - mass galaxies ( @xmath81 ) is enhanced by metal cooling at early times ( @xmath85 ) .
however we find that the behaviour of mf at the low - mass end is somewhat dependent on the threshold density for the grouping . as mentioned in
[ sec : groups ] , our grouping code imposes a threshold gas density of @xmath29 for a gas particle to be part of a galaxy . when we lower this threshold density , we find that the two mfs for the n216l10 and n216l10mc converged , as more particles are incorporated into galaxies at the outskirts of galaxies .
metal cooling increases the gas density in galaxies , enabling more low - mass galaxies to satisfy the gas density threshold and to be identified as simulated galaxies .
therefore the behaviour of mf at @xmath86 is somewhat dependent on the choice of the outer threshold density of galaxies .
figure  [ fig : mf_comp ] compares the gas and baryonic ( star + gas ) mfs for the n216l10 series at @xmath1 ( panels @xmath52 & @xmath72 ) , and the redshift evolution of the baryonic mf from @xmath87 to @xmath1 in the n216l10 ( panel @xmath74 ) and n216l10mc ( panel @xmath77 ) runs .
here we show only the range of @xmath88 , which corresponds to the limiting mass of 32 gas particles .
panel ( @xmath52 ) shows that more gas is converted into stars in the runs with metal cooling ( n216l10mc and n216l10mv ) compared to the n216l10 run , while panel ( @xmath72 ) shows that the baryonic mf is similar in all the runs . as we discussed in
[ sec : evolv4 ] , the igm accretion onto galaxies is not so much enhanced yet before @xmath89 due to relatively low igm metallicity , which explains the similarity in the total baryonic mf . in the runs with metal cooling ,
the peak of gsmf is shifted toward higher mass , whereas the peak of gas mf is shifted toward lower mass .
this suggests that the enhanced gsmf at @xmath1 is due to the increased sf efficiency by metal cooling within the galaxies .
again , the apparent enhancement in the number of low - mass galaxies in the n216l10mc run compared to the n216l10 run is somewhat dependent on the threshold density of grouping , therefore it should be interpreted with caution .          as the accretion of igm onto galaxies
become more efficient at @xmath2 due to the increased igm cooling by the metals , we expect that the total baryonic mass of galaxies would be more enhanced at @xmath0 than at @xmath1 .
figure  [ fig : mf_comp_n288l34 ] shows the mfs at @xmath1 and @xmath0 in the n288l34 series , and it clearly demonstrates that this expectation is true .
here we show only the range of @xmath91 , which corresponds to the limiting mass of 32 gas particles .
the general trend in the three ( star , gas , and total baryon ) mfs is similar to that we saw in figure  [ fig : mf_comp ] , although the difference at @xmath1 between the n288l34 and n288l34mc runs is slightly smaller than in the n216l10 series due to poorer resolution .
figure  [ fig : mf_comp_n288l34]f shows most prominently the enhancement of the total baryonic mf at @xmath0 in the n288l34mc run owing to the metal cooling . in panels ( @xmath92 ) & ( @xmath93 ) , the n288l34 run actually has a longer tail at the most massive - end than the n288l34mc run .
the reason for this feature is not fully clear , but it may be related to the balance between igm accretion and feedback . owing to metal cooling , igm accretion rate increases in the mc run , and
the sfr is also enhanced , leading to a stronger feedback .
the amount of mass loss is greater in low mass galaxies , but the net heating of igm is more significant for massive galaxies .
the mc run has stronger feedback , therefore its feedback heating may become more significant than igm accretion for very massive galaxies .
the significant igm heating suppresses the growth of massive - end of mass function from @xmath1 to @xmath0 and results in shorter tail for the mc run at @xmath0 as shown in figure  [ fig : mf_comp_n288l34]e , f .
for a fixed galaxy baryonic mass , the enhancement of star formation by metal cooling would decrease the gas mass fraction , @xmath94 .
figure  [ fig : gfrac_evol_n216l10 ] shows @xmath95 as a function of galaxy stellar mass for the n216l10 series .
panel ( @xmath52 ) shows the data only from n216l10mc run at @xmath1 , and each data point corresponds to a simulated galaxy . to characterise the distribution ,
we compute the following two quantities in each logarithmic stellar mass bin : ` average ' and ` median ' .
the ` average ' is the ratio of total gas mass to total baryonic mass for all the galaxies in each mass bin , i.e. , @xmath96 .
the ` median ' case is simply the median of @xmath95 values in each mass bin .
both quantities show a similar trend , however , there is a sharper drop - off at @xmath97 for the ` median ' case in figure  [ fig : gfrac_evol_n216l10]a .
this mass - scale corresponds to 32 star particles in the n216l10 series .
we find that there are many galaxies with @xmath98 above this mass - scale , which causes the sharp drop - off in the ` median ' line . below this limiting mass ,
galaxies are not resolved well , which results in an underestimate of star formation and an overestimate of @xmath95 .
if we had a higher resolution simulation with finer particle masses , this limiting mass - scale would shift to a lower mass .
therefore the location of this sharp drop - off is currently determined by the resolution of our simulation . however
, dark matter halos would stop forming stars at some lower limiting halo mass , if we had an infinitely high - resolution simulation , this lower limit to the galaxy mass is presumably determined by the photoevaporation of gas by the uv background radiation @xcite .
recent works suggest that star formation could be suppressed by the uv background in halos with @xmath99 at @xmath89 .
galaxies with @xmath100 would reside in halos with @xmath101 .
the n216l10 series would resolve such a halo with @xmath102 dark matter particles , and its mass resolution is actually close to the astrophysical limit for dwarf galaxy formation at @xmath89 .
therefore the sharp drop - off in @xmath95 at @xmath100 may not be so far from the true answer .
we find that @xmath103 increases with decreasing @xmath104 at all redshifts , regardless of metal cooling and wind effects .
this trend is qualitatively consistent with current observations . @xcite estimate the gas fraction as a function of stellar mass using the rest - frame uv - selected star - forming galaxies at @xmath105 , and show that the gas fraction decreases with increasing stellar mass . and using the mean gas and stellar mass , they find the average @xmath106 , which agrees with the predicted gas fraction for massive galaxies in our simulation .
in addition , @xcite reported that the average neutral gas fraction is @xmath107 for the local dwarf galaxies selected from the sloan digital sky survey . in our n216l10 series with metal cooling , the ` average ' @xmath95 reaches 0.6 for galaxies with @xmath108 at @xmath1 & 4 .
figures[fig : gfrac_evol_n216l10]b , c , d show the redshift evolution of @xmath95 for the n216l10 series . in the runs with metal cooling ( n216l10mc and n216l10mv ) , @xmath95 is lower than in the n216l10 run by @xmath30% at all redshifts for galaxies with @xmath109 .
this result suggests that the metal cooling reduces @xmath95 owing to more efficient star formation .
the values of @xmath95 seem to be more convergent at the massive - end ( @xmath110 ) .
in addition , we find that galaxies in the n216l10nw run are the most gas - rich , because almost no galactic gas is returned to the igm .    figures  [ fig : gfrac_evol_n216l10]b , c , d also show that @xmath95 is higher at higher redshifts .
for example , the well - resolved galaxies with @xmath111 have @xmath112 at @xmath87 , but @xmath113 at @xmath1 .
as we saw in figure  [ fig : mf_star]d , there is a considerable growth in gsmf from @xmath87 to @xmath1 .
together with the decrease in the gas fraction , these results suggest that most of the gas accreted during @xmath114 has been converted into stars
.    figure  [ fig : mgfrac_resolution]a highlights the redshift evolution of @xmath95 from @xmath87 to @xmath0 in the n288l34mc run .
the gas mass fraction clearly decreases with decreasing redshift as the gas is converted into stars .
the rate of decrease is greater at @xmath115 than at @xmath116 , with @xmath117 & 0.2 for @xmath118 , & 1 , respectively , for galaxies with @xmath119 .
figure  [ fig : mgfrac_resolution]b shows the effect of resolution on @xmath95 .
the location of the drop - off in @xmath95 shifts to lower masses as the resolution is increased from n144l10mc to n216l10mc run .
the location of the drop - off is at around our resolution limit  the vertical lines in _ panel ( b)_. therefore the values of @xmath95 are overestimated in unresolved galaxies .
using cosmological hydrodynamic simulations with metal enrichment and metal cooling , we studied their effects on galaxy growth and cosmic sfr .
owing to metal cooling , the sfr density increases about 20% at @xmath1 and about 50% at @xmath0 .
our results suggest that metal cooling enhances the star formation through two different processes : 1 ) more efficient conversion of local ism into stars ( i.e. , increase of local sf efficiency ) , and 2 ) the increase of igm accretion onto galaxies .
the former process is in effect essentially at all times at @xmath120 , because the local ism can be instantaneously enriched by sn explosions .
this process enhances the sfr , but does not noticeably change the total baryonic mass of galaxies . the latter process ,
on the other hand , can increase the total baryonic mass of galaxies , as well as enhancing the overall sfr density by supplying more gas for star formation .
this process becomes effective only at lower redshifts ( @xmath39 ) , because it takes some time to enrich the igm up to @xmath121 .    in this paper
, we used two different prescriptions for star formation in the runs with metal cooling : 1 ) the ` constant @xmath8 ' scheme and 2 ) the ` varying @xmath8 ' scheme .
the value of @xmath8 was fixed in the former scheme , while it was modulated according to the metallicity in the latter scheme .
both schemes show similar increase of cosmic sfr density ( see figure  [ fig : sfr_n216l10 ] ) , which implies that the overall sf enhancement by metal cooling is not so sensitive to the choice of @xmath8 ( but with some differences as we discussed in
[ sec : cos_sfr ] ) .
the multiphase ism model for star formation by @xcite contains two free parameters : @xmath8 and the normalisation of gas consumption time - scale , @xmath122 .
the values of these two parameters were originally chosen to match the empirical kennicutt - schmidt law @xcite using simulations of isolated disk galaxies . since we did not change the value of @xmath122 in the runs with metal cooling ( as well as @xmath8 in the ` mc ' runs ) ,
it is possible that our simulations may now violate the kennicutt law .
however , as we show in figure  [ fig : kenni ] , the plots of @xmath123 vs. @xmath124 for the n216l10 and n216l10mc runs are not so different from each other .
this can be understood as follows . as figure  [ fig : evolv4]c showed
, metal cooling increases the density and lowers the temperature of star - forming gas . the fundamental scaling relationship between sfr and cold gas density
does not need to change when the metal cooling is introduced ; the star - forming gas particles would simply slide upward along the kennicutt law , @xmath125 .    in the case of n216l10mv run , we varied @xmath8 according to the gas metallicity , primarily scattering it to lower values because higher metallicity increases the cooling rate and lowers @xmath8
. therefore more gas particles become eligible to form stars and @xmath124 is scattered upward above the kennicutt law for a given value of @xmath123 , resulting in a broader deviation from the power - law relationship of the kennicutt law near @xmath126 . at higher values of @xmath123 , @xmath124 is enhanced , but still within the range of the kennicutt law .
this difference can be alleviated by adjusting the sf timescale ( see * ? ? ?
* ) , however , we try to emphasise the effect of metal cooing on sfr by keeping this parameter the same in all the runs .
further investigations of the hi  aspect of our simulations is beyond the scope of the present paper , and we will present the results elsewhere .    upon evaluating the effect of metal cooling , @xcite considered two different gas phases where the cooling become important .
one is the diffuse gas in galactic halos , which must radiate its thermal energy in order to collapse onto the high - density ism , and the other is the multiphase star - forming gas .
they argued in their
5.2 that the cosmic sfr density at @xmath127 would not be affected by the metal cooling very much , because at high-@xmath31 cooling is so efficient that the gas in diffuse atmospheres of halos cools nearly instantly , even without any metal cooling .
they expected their model results to be largely independent of metal enrichment , because in this regime the evolution of sfr density is driven by the fast gravitational growth of the halo mass function .
for the star - forming multiphase ism , they argued that the parameters of their model ( e.g. cold gas evaporation efficiency and gas consumption time - scale ) can be adjusted such that the normalisation of the kennicutt law can be maintained , yielding to first order unaltered model predictions .
however , our simulations , in particular the comparison between the n216l10 and n216l10mc run , suggest that the metal cooling can enhance the sfr density even at @xmath127 by @xmath128% .
this is because the gas density becomes higher and temperature becomes lower ( see figure  [ fig : evolv4]c ) by responding to the enhanced cooling rate by metals .    in summary
, we consider that our study serves to demonstrate the importance of metal enrichment and metal cooling in galaxy formation , beyond our original motivation to perform a simple numerical study .
our study highlights the role of metals at low-@xmath31 and high-@xmath31 to enhance the star formation in the universe , and demonstrates that the metals increase the sf efficiency in the star - forming ism , as well as enabling more igm to accrete onto galaxies and fuel star formation .
we expect that metal cooling would also alter some of the galaxy properties such as the mass - metallicity relationship ( e.g. , * ? ? ?
* ) and specific sfr ( e.g. , * ? ? ?
we plan to investigate these issues in the future , as well as alternative models for star formation and feedback .
we thank v. springel for allowing us to use the updated version of gadget-2 code for our study , and for useful comments on the manuscript .
kn is grateful for the hospitality of institute for the physics and mathematics of the universe , university of tokyo , where part of this work was done .
we are also grateful to the anonymous referee for constructive comments which improved this paper .
this research was supported in part by the national aeronautics and space administration under grant / cooperative agreement no .
nnx08ae57a issued by the nevada nasa epscor program , by the national science foundation through teragrid resources provided by the san diego supercomputer center ( sdsc ) , and by the president s infrastructure award at unlv .
the simulations were performed at the unlv cosmology computing cluster and the datastar at sdsc .
d.  n. , verde l. , peiris h.  v. , komatsu e. , nolta m.  r. , bennett c.  l. , halpern m. , hinshaw g. , jarosik n. , kogut a. , limon m. , meyer s.  s. , page l. , tucker g.  s. , weiland j.  l. , wollack e. , wright e.  l. , 2003 , apjs , 148 , 175


early - type galaxies in clusters exhibit a linear color - magnitude ( cm ) relation indicating that bright galaxies are systematically redder than their faint cluster companions ( visvanathan & sandage 1977 ) .
this remarkable relation shows very small scatter ( @xmath3 magnitude ) in high precision photometry of local clusters such as coma and virgo ( bower , lucey & ellis 1992a , 1992b , hereafter ble92 ) and can be extended to clusters at medium - to - high redshift ( @xmath4 ) ( ellis et al 1997 , stanford , eisenhardt & dickinson 1998 ) . a first attempt at explaining the universality of the cm relation involves using the age of each galaxy as the main determinant of its color .
ageing stellar populations redden progressively as stars with decreasing initial mass evolve off the main sequence .
therefore , if the colors of cluster galaxies are purely controlled by age , the small scatter about the cm relation implies a nearly synchronous star formation process for all galaxies of a given mass , while the slope of the cm relation implies systematically older ages for more massive galaxies .
as shown most recently by kodama & arimoto ( 1997 ) , such a picture is highly unlikely because it does not preserve the slope nor the magnitude range of the cm relation in time .
another important factor that affects the colors of stellar populations is metallicity . at fixed age
, a more metal - rich stellar population will appear redder and fainter than a more metal - poor one ( e.g. , worthey 1994 ) .
hence , increasing metallicity at fixed age has a similar effect on colors as increasing age at fixed metallicity .
this is usually referred to as the _ age - metallicity degeneracy _ ( worthey 1994 ) .
several studies have shown that cm relation of cluster elliptical galaxies could be primarily driven by metallicity effects ( larson 1974 ; matteuci & tornamb 1987 ; arimoto & yoshii 1987 ; bressan , chiosi & tantalo 1996 ; kodama & arimoto 1997 ) .
the physical mechanism usually involved is that of a galactic wind : supernovae - driven winds are expected to be more efficient in ejecting enriched gas , and hence in preventing more metal - rich stars from forming , in low - mass galaxies than in massive galaxies with deeper potential wells .
although age is generally assumed to be the same for all galaxies in these studies , this has not been proven to be an essential requirement .
in fact , scenarios in which e / s0 galaxies progressively form by the merging of disk galaxies ( schweizer & seitzer 1992 ) in a universe where structure is built via hierarchical clustering also predict that the cm relation is driven primarily by metallicity effects ( kauffmann & charlot 1998 )
. moreover , age effects could be important if , for example , there is sufficiently strong feedback from early galaxy formation to bias the luminous mass distribution of subsequent generations of galaxies by the heating of intergalactic gas .
in this paper we present a new , more model - independent approach for evaluating the full range of ages and metallicities allowed by the spectro - photometric properties of early - type galaxies in clusters .
the method is based on the construction of age - metallicity diagrams constrained by the colors of early - type galaxies in the nearby coma cluster and in 17 clusters observed with the _ hubble space telescope _ ( _ hst _ ) at redshifts up to @xmath5 ( stanford et al .
such an analysis has hitherto been hindered because of the lack of both accurate stellar libraries for different metallicities and reliable morphological information on cluster galaxies at medium - to - high redshifts .
our results can subsequently be reframed into specific theories of galaxy formation , since they will be indispensable for any model that seeks to produce galaxies resembling those actually observed . in
2 we present the spectral evolution models used in this paper . the cluster sample is described in  3 . in
4 we construct the age - metallicity diagrams allowed by the observations , and in  5 we compute the corresponding ranges in mass - to - light ratio and in several commonly used spectral indices .
we discuss our main conclusions in  6 .
we compute the spectral evolution of early - type galaxies using the latest version of the bruzual & charlot ( 1998 ) models of stellar population synthesis .
these span the range of metallicities @xmath6 and include all phases of stellar evolution , from the zero - age main sequence to supernova explosions for progenitors more massive than @xmath7 , or to the end of the white dwarf cooling sequence for less massive progenitors .
in addition , the models predict the strengths of 21 stellar absorption features computed using the worthey et al .
( 1994 ) analytic fitting functions for index strength as a function of stellar temperature , gravity and metallicity .
this constitutes the standard `` lick / ids '' system that is often used as a basis for spectral diagnostics in early - type galaxies .
the resulting model spectra computed for stellar populations of various ages and metallicities have been checked against observed spectra of star clusters and galaxies ( bruzual et al .
1997 ; bruzual & charlot 1998 ) .    the uncertainties in the models are discussed in charlot , worthey & bressan ( 1996 ) .
these can reach up to 0.05  mag in rest - frame @xmath8 , 0.25  mag in rest - frame @xmath9 and a 25% dispersion in the @xmath10-band mass - to - light ratio . with these uncertainties in mind
, we will concentrate more on understanding the trends seen in the observations than on inferring absolute age and metallicity values .
it is worth noting that the most massive elliptical galaxies exhibit [ mg / fe ] ratios in excess of that found in the most metal - rich stars in the solar neighborhood ( by @xmath11 dex ; see worthey , faber , & gonzalez 1992 ) . while this may limit the accuracy of the predicted mg@xmath12 indices of bright elliptical galaxies , the recent models of bressan et al .
( 1998 ) convincingly show that an enhancement in light elements at fixed total metallicity has virtually no effect on the other spectrophotometric properties of model stellar populations .
we approximate model early - type galaxies by instantaneous - burst stellar populations .
the reason for this is that we aim at constraining the age and metallicity ranges of stars dominating the light of early - type galaxies , whose photometric properties are well represented by instantaneous - burst populations .
in fact , this is true even if the galaxies underwent subsequent small amounts of star formation or if the epoch of major star formation was extended over several billion years ( e.g. , fig .  1 of charlot & silk 1994 ) .
the predicted colors of our models at fixed age and metallicity agree well with the results of more refined calculations including the effects of infall and galactic winds for corresponding values of the mean age and metallicity ( kodama & arimoto 1997 , bressan , chiosi & tantalo 1996 ) .
for example , adopting metallicities matching the luminosity  weighted metallicities @xmath13 in table  2 of kodama & arimoto yields @xmath14 and @xmath9 colors that agree to better than 0.05  mag with the results from these authors at an age of 15  gyr .
such a discrepancy is well within the errors of current population synthesis models ( charlot et al .
1996 ) . in the remainder of the present paper , the initial mass function ( imf ) is taken from scalo ( 1986 ) and is truncated at 0.1 and 100@xmath15 .
we use the above models to compute the locations in the age
metallicity diagram of stellar populations satisfying specified spectro - photometric properties .
figure  1 shows four such age
metallicity diagrams corresponding to imposed values of the @xmath14 and @xmath9 colors and mg@xmath16 and h@xmath17 spectral indices .
these quantities are chosen here because they can be constrained by many observations of early - type galaxies (  3 and  5 ) . in each panel , the models satisfying the same value of the spectro - photometric property of interest are related by a continuous line , different lines corresponding to different imposed values . with this definition , the slope of a line in the age
metallicity diagrams indicates the relative sensitivity of the color or index under consideration to age and metallicity .
vertical lines would correspond to a sensitivity purely to age , and horizontal lines to a sensitivity purely to metallicity .
figure  1 then shows immediately that the @xmath14 and @xmath9 colors and mg@xmath16 index depend more strongly on metallicity than on age , while the h@xmath17 index depends more strongly on age than on metallicity .
we will return to this point in  4 and  5 .
the relative dependence of the spectro - photometric properties of instantaneous - burst populations on age and metallicity has been previously investigated by worthey ( 1994 ) .
he used the parameter @xmath18 at fixed color or index to represent the ratio of the change @xmath19 in age needed to counterbalance a change @xmath20 in metallicity in order to keep that color or index unchanged .
the difference between worthey s and our approach is that he computed a single effective value of @xmath18 for each spectro - photometric property , while the different lines in figure  1 show the behavior of the @xmath18 slope for different values of the color or index under consideration . for comparison ,
the arrow in each panel of figure  1 indicates the @xmath21 vector obtained by worthey ( 1994 ) . in each case
, the general agreement with the mean slope of the lines is good .
table  1 gives a more quantitative comparison between worthey s ( 1994 ) and our results .
we computed linear fits to all lines in figure  1 and then a linear fit between the derived slopes and their corresponding color or index value .
the slopes and zero points of these relations for each spectro - photometric property are listed in columns ( 2 ) and ( 3 ) of table  1 .
we then evaluated @xmath18 for four values of the @xmath14 and @xmath9 colors ( and corresponding model predictions for the mg@xmath16 and h@xmath17 indices ) matching the properties of early - type galaxies at four magnitudes @xmath22 , @xmath23 , @xmath24 and @xmath25 along the local cm relation ( ble92 ; see  3 ) .
the agreement with worthey s ( 1994 ) predictions is seen to be of the order of @xmath26 % .
it is worth noting that our generalized fits deviate significantly from worthey s `` 3/2 rule '' which takes the @xmath27 ratio to be 1.5 for any color .
c|cc|cccc|c property & slope & zero point & @xmath28= 17.5 & @xmath28= 19.0 & @xmath28= 20.5 & @xmath28= 22.0 & worthey + [ @xmath14 ] & 0.444 & 0.653 & 1.19 & 1.24 & 1.30 & 1.35 & 1.5 [ @xmath9 ] & 0.886 & @xmath311.004 & 1.59 & 1.69 & 1.79 & 1.89 & 1.9 [ mg@xmath16 ] & 0.118 & 0.976 & 1.33 & 1.36 & 1.40 & 1.45 & 1.7 [ h@xmath17 ] & @xmath310.092 & 0.645 & 0.47 & 0.48 & 0.49 & 0.50 & 0.6
observational constraints on the photometric properties of early - type galaxies are taken from the recent sample of stanford , eisenhardt & dickinson ( 1998 ) .
the sample consists of 19 clusters in the redshift range @xmath34 .
these were extracted on the basis of available _ hst _ imaging from a larger , heterogeneous sample of 46 clusters drawn from a variety of optical , x - ray and radio - selected samples .
the 19 clusters studied by stanford et al .
( 1998 ) were imaged in the near infrared @xmath35 , @xmath36 and @xmath37 passbands .
exposure times in all passbands were chosen to achieve 5-@xmath38 detection of objects with the spectral energy distribution of an unevolved present - day elliptical galaxy down to 2  mag fainter than @xmath39 ( corresponding to an apparent limiting magnitude @xmath40 and 20.0 at @xmath41 and 0.895 , respectively ) . for 17 out of these 19 clusters , photometry is also available in two optical passbands , referred to as _
blue _ and _ red _ , that were tuned as a function of redshift to span the 4000    break in the galaxy rest frame spectra .
this subsample of 17 clusters is of considerable interest to us because it sets the most useful constraints on the spectral properties of galaxies .
stanford et al .
morphologically classified the galaxies in these clusters on the basis of _ hst _ wide field and planetary camera 2 ( wfpc2 ) images .
they did not attempt to distinguish between e and s0 galaxies , and for the purpose of the present analysis we also consider these galaxy types together as a single early - type class .
table  2 lists the name , redshift , _ blue _ and _ red _ passbands and number of e / s0 galaxies for each of the 17 clusters in the sample .    on the average
, cluster membership is expected to be secure for over 85% of the e / s0 galaxies in the sample , as estimated either from morphologically - dependent number counts or from statistical field corrections ( see stanford et al .
1998 for a thorough analysis ) .
furthermore , all clusters exhibit a tight cm relation with a slope showing no significant change out to @xmath5 .
this has been taken as evidence that all galaxies shared a common history of star formation ( stanford et al .
1998 ; kodama et al . 1998 ; but see kauffmann & charlot 1998 ) .
we note that there is an open question about the membership of some morphologically - selected early - type galaxies with faint magnitudes and colors far bluer than the cm relation in the clusters ( e.g. , kodama et al .
this can be appreciated most readily from deep _ hst _ wfpc2 imaging by ellis et al .
( 1997 ) of three of the clusters in table  2 , f1557.19tc , cl  0016 + 16 and j1888.16cl .
although these faint ( @xmath42 ) galaxies may be field contaminants , for completeness we should not abandon the possibility that they are blue cluster early - type galaxies that violate the cm relation .
thus , they are also included in our analysis .
it is worth pointing out that stanford et al .
( 1998 ) find no significant difference between the photometric properties of early - type galaxies in clusters of different richness or x - ray luminosity at a similar redshift .
while this finding does not necessarily imply that all early - type galaxies were assembled early ( kauffmann & charlot 1998 ) , it strongly supports analyses like the one presented in  4 in which constraints on photometric evolution are drawn from comparisons of the properties of galaxies in clusters at different redshifts . also , to avoid any bias near the bright end of the cm relation , we have checked that the brightest cluster galaxies in the sample of table  2 match the @xmath37-band absolute luminosity versus redshift relation published by aragn - salamanca et al .
( 1998 ) .
finally , to tighten the models at @xmath43 , we use the cm relation derived by ble92 for e / s0 galaxies in the coma cluster .
this takes the forms @xmath44 and @xmath45 , with standard deviations of 0.055  mag and 0.065  mag in @xmath14 and @xmath9 colors , respectively . to define the absolute magnitude scale of the models , the recession velocity of the coma cluster
is taken to be @xmath46 km  s@xmath47 ( han & mould 1992 ) .
we adopt @xmath48 km  sec@xmath47mpc@xmath47 and @xmath49 , except when otherwise indicated .
the distance modulus to the coma cluster is therefore 35.41  mag .
the small scatter of the cm relation at fixed luminosity has been used by ble92 to constrain the age range of early - type galaxies in clusters ( see also ellis et al .
they computed the maximum spread in star - formation epoch allowed by the observed scatter in @xmath14 color according to the rate of color change @xmath50 predicted by population synthesis models .
the basis of this argument can be understood from figure  2 .
the top panel shows the @xmath14 color evolution of a stellar population in the time interval @xmath51 $ ] as a function of age @xmath52 for several metallicities , while the bottom panel shows the inferred @xmath50 evolution as a function of @xmath52 .
since @xmath50 decreases with model age @xmath52 , ble92 found that the older the galaxies , the larger the allowed spread in star - formation epoch .
on the basis of this analysis , ble92 favored a major star - formation epoch at @xmath53 for early - type galaxies in clusters , spread over a period of roughly 1  gyr .
it is crucial to realize that ble92 s argument utilizing @xmath50 is based on the _ a priori _ requirement that at fixed luminosity , all early - type galaxies in the cm relation have the same metallicity .
in fact , the analysis was conducted using population synthesis models for uniquely solar metallicity .
as the present paper demonstrates , however , at fixed luminosity the photometric constraints on cluster galaxies allow a wide range of ages and metallicities (  4 ) .
hence , the @xmath54 argument already includes restrictive hypotheses on the ages and metallicities of early - type galaxies with respect to the full allowed ranges .
it is not surprising , therefore , that ble92 s conclusions on the epoch of major star formation for early - type galaxies in clusters are unnecessarily constraining when compared to the results of our more complete analysis below .
lccccc cluster & @xmath55 & _ blue _ @xmath31 _ red _ & n(e / s0 ) & & & ac 118 & 0.308 & @xmath56 & 38 ac 103 & 0.311 & @xmath56 & 32 ms 2137.3 - 234 & 0.313 & @xmath56 & 21 cl 2244 - 02 & 0.330 & @xmath56 & 24 cl 0024 + 16 & 0.391 & @xmath56 & 39 gho 0303 + 1706 & 0.418 & @xmath56 & 38 3c 295 & 0.461 & @xmath57 & 25 f1557.19tc & 0.510 & @xmath57 & 29 gho 1601 + 4253 & 0.539 & @xmath57 & 42 ms 0451.6 - 0306 & 0.539 & @xmath57 & 51 cl 0016 + 16 & 0.545 & @xmath57 & 65 j1888.16cl & 0.560 & @xmath57 & 38 3c 220.1 & 0.620 & @xmath57 & 22 3c 34 & 0.689 & @xmath57 & 19 gho 1322 + 3027 & 0.751 & @xmath58 & 23 ms 1054.5 - 032 & 0.828 & @xmath58 & 71 gho 1603 + 4313 & 0.895 & @xmath59 & 23
we now use the models described in  2 and the observational constraints outlined in  3 to compute the regions allowed in age  metallicity space for early - type galaxies in clusters .
we proceed as follows .
we consider each cluster in table  2 , with redshift @xmath55 , and search for all models which can match the photometric properties of galaxies in that cluster at @xmath55 and those of coma galaxies at @xmath43 . in practice , we start by considering models that span a full range of metallicities and with ages at @xmath60 that range between the age of the universe and the lookback time to redshift @xmath55 . out of these models , we select all those matching @xmath14 and @xmath9 colors of galaxies in coma within the scatter of the observed cm relation (  3 ) . for a set of @xmath14 and @xmath9 colors
this defines a range of possible @xmath10 magnitudes .
we then compute the predicted apparent @xmath37 magnitudes and _ blue_@xmath61 , _ red_@xmath61 , @xmath62 and @xmath63 colors of the selected models at the cluster redshift @xmath55 .    to decide whether a model is acceptable , we compare the predicted photometric properties at @xmath55 with the observations . for this purpose , the observed cm relation of each cluster in table  2
was divided into four apparent @xmath37 magnitude bins , roughly 1  mag wide , allowing good sampling from the brightest to the faintest galaxies ( see  3 ) . for each bin ,
the standard deviation around the relation was computed in _
blue_@xmath61 , _
red_@xmath61 , @xmath62 and @xmath63 colors .
this procedure allows us to account for the increased scatter of the observed cm relations towards fainter magnitudes .
a model is retained if it simultaneously falls within @xmath64 of the cm relation in all four colors .
the main reason for adopting a @xmath64 criterion is that we must allow for known uncertainties in the spectral evolution models .
uncertainties in the predicted optical / infrared color evolution over lookback times of a few billion years can reach a few tenths of a magnitude in current population synthesis models ( charlot et al .
this is several times larger than the typical scatter around the cm relation for the clusters of table  2 ( see fig .
5_a _ of stanford et al .
1998 and table  1 of kodama et al . 1998 ) .
figures  3_a _ , 3_b _ and 3_c _ show the results of our analysis for the clusters in the redshift ranges @xmath65 , @xmath66 and @xmath67 , respectively .
metallicity space is parameterized in terms of the formation redshift @xmath68 and the iron abundance computed as [ fe / h]@xmath69 , with @xmath70 and @xmath71 .
we separate galaxies into four luminosity bins at @xmath43 by defining apparent magnitude bins in coma centered on @xmath72 , 14.25 , 13.5 and 12.75 .
these correspond roughly to absolute @xmath10 luminosities @xmath73 , @xmath74 , @xmath75 and @xmath76 , respectively .
the four panels in figures  3_a_3_c _ separately show the results for each of the four luminosity bins . at fixed luminosity ,
figure  3 indicates that the observations of all 17 clusters roughly constrain similar regions in @xmath77})$ ] space .
this is a consequence of the similarity of the cm relation among all clusters in the sample .
an important discriminant , however , is that the condition @xmath78 allows more recent formation epochs and higher metallicities for galaxies in low - redshift clusters than for those in high - redshift clusters .
as expected , different luminosities imply different allowed values of @xmath68 and [ fe / h ] .
since faint galaxies on the cm relation are bluer than bright galaxies , they can be modelled on average by more metal - poor , and to some extent younger , stellar populations .
we also note that the metallicity range of model galaxies in figure  3 , @xmath79 } \simlt + 0.2 $ ] , compares well with the range observed in nearby e / s0 galaxies ( e.g. , worthey , faber & gonzalez 1992 ) .
the upper panel of figure  4 shows the areas corresponding to all four luminosity bins superimposed on a similar diagram for the cluster cl0016 + 16 .
this cluster , at @xmath80 , has one of the best - defined cm relations in the sample ( 65 e / s0 galaxies ; see table  2 ) .
figure  4 readily shows that a pure age sequence can not account for the photometric properties of early - type galaxies because a single horizontal line can not cross all four areas simultaneously in @xmath81})$ ] space .
alternatively , different combinations of age and metallicities can accommodate the data .
the simplest assumption of a pure metallicity sequence ( vertical line ) requires that the bulk of stars in galaxies formed at redshifts @xmath53 , otherwise one can not account for the photometric properties of the brightest galaxies .
this confirms earlier conclusions by ellis et al .
( 1997 ) based on similar data for cl0016 + 16 . as figure  4 demonstrates , however , there is no requirement that the bulk of stars in all early - type galaxies in this cluster formed at the same epoch .
if the most metal - rich galaxies form most recently , as might be expected from simple chemical enrichment arguments , then figure  4 implies that the dominant stellar populations of all early - type galaxies in cl  0016 + 16 were already in place at redshifts @xmath82 . on the other hand , a wide range of scenarios
are also allowed in which faint metal - poor galaxies form as recently as @xmath83 , i.e. , more recently than bright metal - rich galaxies .
the lower panel of figure  4 shows that the alternative cosmology @xmath84 km  sec@xmath47mpc@xmath47 and @xmath85 would lead on average to slightly larger @xmath68 values for model galaxies because of the lower age of the universe implied at redshift 1 .
we emphasize that the constraints on the formation redshift and metal abundance of the stellar populations of early - type galaxies derived from figure  4 apply to cl  0016 + 16 only . as figure
3_a _ shows , constraints on clusters at lower redshifts are consistent with formation redshifts as small as 0.6 for even the brightest , reddest galaxies .
hence , the apparent lack of evolution of the slope and scatter of the cm relation in clusters does not imply by itself a common epoch of star formation for all early - type galaxies .
instead , what is strictly implied is that new galaxies joining the cm relation at low redshifts must be more metal - rich than their older cluster companions of similar luminosity . the required age
metallicity evolution inferred from figure  3 corresponds roughly to a 25% increase in [ fe / h ] from @xmath86 to @xmath87 , i.e. , over a period of 2.5  gyr for the assumed @xmath88 km  sec@xmath89mpc@xmath47 and @xmath49 .
we note that adopting longer timescales of star formation  instead of an instantaneous burst  for early - type galaxies in clusters would increase the formation redshift of the first stars in model galaxies with respect to the results of figure  3 .
whatever the adopted history of star formation , however , the values of @xmath68 in figure  3 are robust limits on the redshift of the last major event of star formation in early - type galaxies ( see  2 ) .    finally , we have explored the possibility that the faint morphologically - selected early - type galaxies with very blue colors in cl  0016 + 16 were actual cluster members ( see  3 ) .
the inferred @xmath68 and [ fe / h ] ranges for these objects is indicated by the heavy shaded region in figure  4 . as expected , the blue colors imply that the galaxies must be both young and metal - poor .
we may speculate that they are young objects which `` violate '' the main sequence cm relation , but which could eventually evolve towards it as they age and undergo chemical enrichment .
however , spectroscopic confirmation as well as a deeper search for faint blue outliers in clusters at higher redshifts are needed before we can draw any definitive conclusions as to their true nature .
the above analysis has enabled us to determine the areas allowed in @xmath90})$ ] space for the stellar populations of e / s0 galaxies in clusters .
while we emphasize that real galaxies may occupy only part of these allowed areas , it is interesting to compute the implied evolution with redshift of other observable properties of early - type galaxies such as mass - to - light ratios and spectral indices . since we do not know which subareas real galaxies might occupy in figure  3 , in the following we investigate properties averaged over the entire allowed areas .
we first investigate the mass - to - light ratios of the model galaxies selected in figure  3 . for each cluster , we compute the logarithmic mean of @xmath91 in each luminosity bin at the cluster redshift by averaging over the @xmath81})$ ] area constrained by the observations . since absolute values of @xmath91
depend sensitively on the assumed low - mass end of the imf , we adopt for all clusters the arbitrary normalization @xmath92 at @xmath93 and focus our investigation on the dependence of the mass - to - light ratio on luminosity .
the mean values of @xmath94 computed in this way are shown as triangles in figure  5 . since the mass - to-@xmath10 light ratio of a stellar population increases at increasing age and metallicity , brighter galaxies in figure  5 generally tend to have larger @xmath94 than fainter galaxies .
this is not always true , however , as can be understood from the dispersion in the results of figure  3 .
also shown in figure  5 are the results from two recent observational studies of the fundamental plane for cluster e / s0 galaxies .
jrgensen , franx , & kjrgaard ( 1996 ) have parameterized the mass - to - gunn  @xmath95 luminosity ratio , @xmath96 , in terms of a combination of half  light radius and central velocity dispersion from observations of a large sample of 226 e / s0 galaxies in 10 nearby clusters .
they conclude that the half
light radius has a negligible effect on determinations of @xmath96 and show that the data follow the simple relation @xmath97 with a scatter of only 25% .
the central velocity dispersion @xmath38 is tightly correlated with luminosity via the faber - jackson relation ( faber & jackson 1976 ) .
we use ble92 s calibration of this relation based on observations of early - type galaxies in the coma cluster in order to relate @xmath98 to @xmath99 in figure  5 .
we then adopt the mean @xmath100 colors of galaxies along the cm relation to convert jrgensen et al.s ( 1996 ) result into an expression involving the mass - to-@xmath10 luminosity ratio .
this yields @xmath101 , corresponding to an increase of less than 3% in logarithmic slope with respect to the relation derived in gunn  @xmath95 .
the horizontal shading in figure  5 then indicates the range of slopes allowed after accounting for the scatter of jrgensen et al.s ( 1996 ) data around the mean relation .
an alternative constraint on the observed range of mass - to-@xmath10 luminosity ratios can be obtained from the relation recently derived by graham & colless ( 1997 ) between half  light radius , central velocity dispersion and mean surface brightness of early - type galaxies .
their analysis is based on accurate fitting of the @xmath10-light profiles of 26 e / s0 galaxies in the virgo cluster using both homologous ( @xmath102 ) and non - homologous ( @xmath103 , with @xmath104 a free parameter ) radial dependence laws .
the results indicate a slight but systematic breaking of the generally assumed homology in the sense that @xmath104 increases with increasing half - light radius .
we can use the virial theorem and faber - jackson relation to reexpress the results of graham & colless ( 1997 ) in terms of a relation involving @xmath91 and @xmath98 .
the outcome is shown in the upper and lower panels of figure  5 for non - homologous and homologous light profiles , respectively . in each case
, the heavy solid line indicates the mean relation and the slanted shading represents the allowed range of slopes when the uncertainties quoted by graham & colless are included .
figure  5 shows that the mass - to - light ratios of model galaxies are all within the range constrained by current observations .
the models appear to be mostly compatible with the graham & colless ( 1997 ) results for non - homologous light profiles .
however , we point out that the model values correspond to purely stellar mass - to - light ratios , whereas there is some observational evidence for the presence of dark halos around early - type galaxies .
this is inferred via the studies of hi kinematics ( e.g. franx et al .
1994 ) , x - ray emission ( e.g. forman et al .
1994 ) , radial velocities of planetary nebulae and globular clusters ( e.g. mould et al .
1990 , hui et al .
1995 ) , gravitational lensing ( e.g. maoz & rix 1993 ) and measurements of the shape of the stellar line - of sight velocity distribution ( e.g. carollo et al . 1995 ) . more recently , rix et al .
( 1997 ) have analyzed the velocity profile of the elliptical galaxy ngc  2434 and show that roughly half the mass within an effective radius is dark .
hence , one must remain cautious in interpreting trends in the mass - to - light ratios of early - type galaxies on the basis of pure population synthesis models .      for completeness
, we also compute the strengths of several commonly used spectral indices implied by the above photometric study for early - type galaxies in clusters .
most observations in this domain have concentrated on measurements of h - balmer , magnesium and iron - dominated indices of the lick system (  2 ) such as h@xmath17 , mg@xmath12 , mg@xmath16 and mgfe ( worthey et al . 1992 ; bender , burstein & faber 1993 ; gonzalez 1993 ; davies , sadler & peletier 1993 ; ziegler & bender 1997 ) .
one of the main motivations for studies of this type is to identify a pair of indices including one mostly sensitive to age and one mostly sensitive to metallicity in order to break the age
metallicity degeneracy from photometric colors .
sensitive h@xmath17 index , however , can be significantly contaminated by nuclear emission in e / s0 galaxies , for which the corrections are uncertain ( e.g. , gonzalez 1993 ) .
it is therefore deemed preferable to study higher - order balmer indices such as h@xmath105 ( e.g. , worthey & ottaviani 1997 ) . also , the apparent systematic increase of [ mg / fe ] from faint to bright e / s0 galaxies challenges analyses based on models with solar - scaled abundance ratios ( e.g. , worthey et al .
1992 ; tantalo , bressan & chiosi 1997 and references therein ) .
worthey ( 1995 ) shows that the metal - sensitive index c4668 of the lick system appears to be significantly less enhanced with respect to iron than mg in bright e / s0 galaxies .
thus , c4668 represents a good alternative to mg - dominated indices for reliable model investigations of early - type galaxy spectra ( e.g. , kuntschner & davies 1998 ) .    by analogy with our approach in  5.1 ,
we compute mean spectral indices for each of the four luminosity bins in each cluster of our sample .
the results are shown in figure  6 for the mg@xmath12 , h@xmath17 , mg@xmath16 , mgfe , c4668 and h@xmath106 indices of the lick system .
as expected , there is a general correlation between index strength and luminosity , which is more pronounced for metal - sensitive ( mg@xmath12 , mg@xmath16 , mgfe and c4668 ) than for age - sensitive ( h@xmath17 , h@xmath106 ) indices .
the reason for this is that luminous galaxies in figure  3 are on average found to be significantly more metal - rich and slightly older than faint galaxies .
also , figure  1 shows that the sensitivity to age of h - balmer indices such as h@xmath17 is significantly weakened after @xmath107  gyr , when a and b stars have evolved off the main sequence . for comparison , the solid and dashed lines in figure  6 show the locations of passively evolving model galaxies with @xmath108 and the metallicities @xmath109 and @xmath110 , respectively .
the mean index strengths of cluster galaxies appear to evolve roughly along these lines , supporting the consistency of the stellar populations with apparently passive evolution .
in fact , ziegler & bender ( 1997 ) find that the evolution from @xmath43 to @xmath111 of the correlation existing between the mg@xmath16 index strength and central velocity dispersion of e / s0 galaxies is consistent with passive evolution .
it should be noted that the model indices computed in figure  6 are global indices averaged over the emission from all stars in a galaxy .
observations of nearby e / s0 galaxies , however , reveal systematic variations of the strengths of spectral indices between the central and outer regions ( e.g. , worthey et al .
we have shown that the tight photometric constraints on early - type galaxies in clusters allow relatively wide ranges of ages and metallicities for the dominant stellar populations .
in particular , the small scatter of the cm relation out to redshifts @xmath2 does not necessarily imply a common epoch of star formation for all early - type galaxies .
it requires , however , that galaxies assembling more recently be on average more metal - rich than older galaxies of similar luminosity . in this context
it is interesting to mention that , based on the spectral indices of nearby e / s0 galaxies , worthey , trager & faber ( 1996 ) favor younger ages for more metal - rich galaxies than for metal - poor ones at fixed velocity dispersion .
the results of our unbiased analysis therefore define the boundaries in age and metallicity that must be satisfied by theoretical studies aimed at explaining the formation and evolution of early - type galaxies in clusters .    the constraints obtained here on the age and metallicity ranges of e / s0 galaxies are consistent with conventional models in which the galaxies all form monolithically in a single giant burst of star formation at high redshift ( e.g. , kodama et al .
1998 , and references therein ) .
in fact , this implies that regardless of the true ages and metallicities of early - type galaxies within the allowed range , their photometric properties will always be consistent with _ apparently _ passive evolution of the stellar populations .
as figure  6 shows , this consistency even extends to spectral index strengths .
our results are also consistent with scenarios in which e / s0 galaxies are formed by the merging of disk galaxies ( schweizer & seitzer 1992 ) in a universe where structure is built through hierarchical clustering ( kauffmann 1996 ; baugh , cole & frenk 1996 ; kauffmann & charlot 1998 ) . for such scenarios ,
figure  3 constrains the metallicity and epoch of the last major event of star formation in e / s0 galaxies and their progenitors ( see  2 and  4 ) . the ages and metallicities of cluster e / s0 galaxies predicted by hierarchical models are found to be consistent with these constraints ( kauffmann & charlot 1998 ) .    to better assess the origin of e / s0 galaxies in clusters one therefore needs to appeal to observational constraints other than their spectro - photometric properties .
for example , conventional models of e / s0 galaxy formation are being challenged by the paucity of red galaxies found at high redshifts in deep surveys ( kauffmann , charlot , & white 1997 ; zepf 1997 ) . also , morphological distinction between e and s0 galaxies and the evolution of the morphology - density relation out to moderate redshifts appear to point to different formation epochs for e and s0 galaxies ( dressler et al .
1997 ) .
the tightness of the cm relation is proof of a stable process in the assembly of cluster early - type galaxies .
however , as we move towards greater redshifts , a drastic change is expected at lookback times that approach the formation of the first e / s0 galaxies .
this change can arise as a systematic blueing , an increased scatter or a slope flattening in the cm relation ( e.g. , aragn - salamanca et al . 1993 ; charlot & silk 1994 ; kauffmann & charlot 1998 ) .
an interesting question is raised by the presence of morphologically - selected early - type galaxies with very blue colors in clusters at moderate redshifts (  3 and  4 ) .
if these objects are true cluster members , our analysis shows that they could be young metal - poor galaxies that will later join the cm relation .
hence , we need to probe deeper down the galaxy luminosity function in distant clusters in order to assess whether these objects can have any fundamental bearing on the origin of early - type galaxies .
we thank adam stanford for sending us a machine readable list of the cluster photometry used in this paper .
thanks the cnrs and the institut dastrophysique de paris for hospitality and financial support , and also acknowledges a ph.d .
scholarship from the `` gobierno de cantabria . ''
j.s . acknowledges support from nasa , nsf and the blaise - pascal chair at the institut dastrophysique de paris .
aragn - salamanca , a. , ellis , r.s . , couch , w.j . ,
& carter , d. 1993 , mnras , 262 , 764 aragn - salamanca , a. , baugh , c. m. & kauffmann , g. 1998 , mnras , in press ( astro - ph/9801277 ) arimoto , n. & yoshii , y. 1987 , , 173 , 23 baugh , c.m . , cole , s. & frenk , c.s . , 1996 , mnras , 283 , 1361 bender , r. , burstein , d. & faber , s. m. 1993 , , 411 , 153 bower , r. g. , lucey , j. r. & ellis , r. s. 1992 , , 254 , 589 bower , r. g. , lucey , j. r. & ellis , r. s. 1992 , , 254 , 601 bressan , a. , chiosi , c. & tantalo , r. 1996 , 311 , 425 bressan , a. et al . 1998 , in preparation bruzual , g. , barbuy , b. , ortolani , s. , bica , e. , cuisinier , f. , lejeune , t. , & schiavon , r.p . 1997 ,
aj , 114 , 1531 bruzual , g. & charlot , s. 1998 , in preparation carollo , c.m .
, de zeeuw , p.t . ,
van der marel , r.p . ,
danziger , i.j .
& qian , e.e .
, 1995 , apj , 441 , l25 charlot , s. , & silk , j. 1994 , apj , 432 , 453 charlot , s. , worthey , g. & bressan , a. 1996 , , 457 , 625 davies , r.l . ,
sadler , e.m .
, & peletier , r.f .
1993 , mnras , 262 , 650 dressler , a. , et al .
1997 , apj , 470 , 577 ellis , r. s. , smail , i. , dressler , a. , couch , w. j. , oemler , a. , butcher , h. & sharples , r. m. 1997 , , 483 , 582 faber , s.m . , & jackson , r. 1976 , apj , 204 , 668 forman , w. , jones , c. & tucker , w. , 1994 , apj , 429 , 77 gonzalez , j.j .
1993 , phd thesis , univ . of santa - cruz graham ,
a. , colless , m. 1997 , , 287 , 221 han , m. , & mould , j.r .
1992 , apj , 396 , 453 hui , x. , ford , h.c .
, freeman , k.c . & dopita , m.a . , 1995 , apj , 449 , 592 jrgensen , i. , franx , m. & kjrgaard , p. 1996
, , 280 , 167 kauffmann , g. 1996 , mnras , 281 , 487 kauffmann , g. , charlot , s. , & white , s.d.m .
1997 , mnras , 283 , 117l kauffmann , g. & charlot , s. 1998 , mnras , in press kodama , t. & arimoto , n. 1997 , 320 , 41 kodama , t. , arimoto , n. , barger , a. j. & aragn - salamanca , a. 1998 , astro - ph/9802245 kuntschner , h. , & davies , r.l .
1998 , mnras , in press ( astro - ph/9710176 ) larson , r. b. 1974 , mnras , 166 , 585 lilly , s. j. , tresse , l. , hammer , f. , crampton , d. & le f ` evre , o. 1995 , , 455 , 108 maoz , d. & rix , h .- w . , 1993 , apj , 416 , 215 matteuci , f. & tornamb , f. 1987 , , 185 , 51 mould , j.r . , oke , j.b .
, de zeeuw , p.t . &
nemec , j.m .
, 1990 , aj , 99 , 1823 rix , h .- w . , de zeeuw , p.t . ,
cretton , n. , van der marel , m. , roeland , p. , & carollo , m. 1997 , apj , 488 , 702 scalo , j.n .
1986 , fundamentals of cosmic physics , 11 , 1 stanford , s. a. , eisenhardt , p. r. & dickinson , m. 1998 , , 492 , 461 schweizer , f. , & seitzer , p. 1992 , aj , 104 , 1039 tantalo , r. , bressan , a. , & chiosi c. 1997 , a&a , submitted ( astro - ph/9710101 ) visvanathan , n. , & sandage , a. 1997 , apj , 216 , 214 worthey , g. 1994 , , 95 , 107 worthey , g. 1995 , in from stars to galaxies , ed . c.  leitherer , u.  fritze - von alvensleben , & j. huchra ( asp conf . series 98 ) , 167 worthey , g. , faber , s.m .
, & gonzalez , j.j .
1992 , apj , 398 , 69 worthey , g. , trager , s.c . , &
faber , s.m .
1996 , in fresh views of elliptical galaxies , ed .
a.  buzzoni , a.  renzini , & a.  serrano ( asp conf .
series 86 ) , 203 worthey , g. , & ottaviani , d.l .
1997 , apjs , 111 , 377 zepf , s.e .
1997 , nature , 390 , 377 ziegler , b.l .
, & bender , r. 1997 , mnras , 291 , 527