A&A 408, 57-65 (2003)
DOI: 10.1051/0004-6361:20030971
A. Mercurio1 - M. Massarotti 2 - P. Merluzzi 2 - M. Girardi 1 - F. La Barbera 2 - G. Busarello2
1 - Dipartimento di Astronomia, Università degli Studi di Trieste,
via Tiepolo 11, 34100 Trieste, Italy
2 - INAF - Osservatorio Astronomico di Capodimonte,
via Moiariello 16, 80131 Napoli, Italy
Received 27 February 2003 / Accepted 19 June 2003
Abstract
We derive the luminosity functions in three bands (BVR) for
the rich galaxy cluster ABCG 209 at z=0.21. The data cover an area of 78 arcmin2 in the B and R bands, while a mosaic of three
pointings was obtained in the V band, covering an area of
approximately 160 arcmin2. The galaxy sample is complete to B = 22.8 (
), V = 22.5 (
)
and R = 22.0 (
). Although the fit of a
single Schechter function cannot be rejected in any band, the
luminosity functions are better described by a sum of two Schechter
functions for bright and faint galaxies, respectively. There is an
indication for a presence of a dip in the luminosity functions in the
range
V = 20.5-21.5 and
R = 20.0-21.0. We find a marked luminosity
segregation, in the sense that the number ratio of bright-to-faint
galaxies decreases by a factor 4 from the center to outer regions. Our
analysis supports the idea that ABCG 209 is an evolved cluster,
resulting from the merger of two or more sub-clusters.
Key words: galaxies: clusters: individual: Abell 209 - galaxies: photometry - galaxies: luminosity function, mass function
Galaxy luminosity function (LF) is a powerful tool to constrain galaxy
formation and evolution, since it is directly related to the galaxy
mass function and hence to the spectrum of initial perturbations.
Hierarchical clustering models predict a mass distribution
characterised by a cut-off above a given mass M* and well
described by a power law at low masses (Press & Schechter
1974). Starting from these results Schechter (1976)
analysed the luminosity distribution of 14 galaxy clusters observed by
Oelmer (1974), by introducing an analytical description in the
form:
Although investigated in several works, the universality of the LF
faint-end slope is still controversial. The value of the faint-end
slope turns out to be
for field galaxies (e.g.,
Efstathiou et al. 1988; Loveday et al. 1992), while
clusters and groups seem to have steeper slopes,
(e.g., De Propris et al. 1995; Lumsden et al. 1997;
Valotto et al. 1997), suggesting the presence of a larger
number of dwarf galaxies (but see also e.g., Lugger 1986;
Colless 1989; Trentham 1998). Changes in the slope of
the faint-end of the LF in clusters can be related to environmental
effects. An increase of the steepening of the LF faint-end in the
cluster outer regions was actually observed (Andreon 2001) and
explained taking into account that the various dynamical processes
which can destroy dwarf galaxies act preferentially in the
higher-density cores.
Lopez-Cruz et al. (1997) showed that clusters with a flat LF (
)
are a homogeneous class of rich clusters with a
single dominant galaxy, symmetrical single peaked X-ray emission and
high gas masses. Irregular clusters have a steeper faint-end, in
particular, the LFs of ABCG 1569 and Coma which present
substructures, can be suitably fitted with the sum of two Schechter
functions with
and
(Lopez-Cruz et al.
1997). Trentham (1997) also suggested that the
faint-end slope of the LF flattens as clusters evolve because of the
destruction of dwarf galaxies by merging with giants galaxies.
The density of the environment seems to affect the distribution of
galaxy luminosity in the sense that Schechter fits are poor for
galaxies in dense environment, where there is an indication of a dip
(e.g., Driver et al. 1994; Biviano et al. 1995; Wilson
et al. 1997; Molinari et al. 1998; Garilli et al. 1999; N
slund et al. 2000; Yagi
et al. 2002). Trentham & Hodgking (2002) identified
two types of galaxy LF, one for dynamically evolved regions
(i.e. region with a high elliptical galaxy fraction, a high galaxy
density, and a short crossing time), such as Virgo cluster and Coma
cluster, and one for unevolved regions, such as the Ursa Major cluster
and the Local Group. A dip is present in the LF of Virgo and Coma and
is absent in LFs of Ursa Major and Local Group.
The differences in shape of the LFs from cluster to cluster could be explained assuming that the total LF is the sum of type specific luminosity functions (hereafter TSLFs), each with its universal shape for a specific type of galaxies (Binggeli et al. 1988). The total LF then assumes a final shape which can be different from cluster to cluster according to the mixture of different galaxy types. Therefore, different mixtures of galaxies, induced by cluster-related processes, may be at the origin of the presence and of the different shape of dips seen in cluster LFs and may be responsible for the differences seen in the total LFs among field, groups and galaxy clusters.
Indeed, dips are found in several clusters, occurring roughly at
the same absolute magnitude (
or
M* + 2.5), within a range of about one magnitude, suggesting that
clusters have comparable galaxy population. However the dips may have
different shapes, and also depend on the cluster region. This could be
related to the relative abundances of galaxy types, which depend on
the global properties of each cluster and on the local density (Durret
et al. 1999).
In this work we study the Abell galaxy cluster ABCG 209 at z = 0.21 (Kristian et al. 1978; Wilkinson & Oke 1978; Fetisova
1981; Mercurio et al. 2003, hereafter Paper I) which
is a rich, X-ray luminous, and massive cluster (richness class R = 3,
Abell et al. 1989;
,
Ebeling et al. 1996;
,
Rizza et al. 1998; M(<
,
Paper I).
The cluster shows an elongation and asymmetry in the X-ray emission
with two main clumps (Rizza et al. 1998), but no strong
cooling flow is detected. The dynamical analysis presented in Paper I
showed that ABCG 209 is characterized by a very high value of
the line of sight velocity dispersion:
km s-1 and by a preferential
SE-NW direction as indicated by: a) the presence of a velocity
gradient in the velocity field; b) the elongation in the spatial
distribution of colour-selected cluster members; c) the elongation of
the X-ray contour levels in the Chandra image; d) the elongation of
the cD galaxy. There is significant evidence of substructure, as shown
by the Dressler & Schectman test. The two-dimensional distribution
of the colour-selected members shows a strong luminosity
segregation. Furthermore, the young dynamical state is also
indicated by the possible presence of a radio halo (Giovannini et al. 1999), possibly a remnant of a recent cluster merger
(Feretti 2002).
This observational scenario suggests that ABCG 209 is presently undergoing strong dynamical evolution with the merging of two or more sub-clumps along the SE-NW direction, but could not allow us to discriminate between two alternative pictures (Paper I). The merging might be either in a very early dynamical state, where the clumps are still in a pre-merging phase, or in a more advanced state, where luminous galaxies trace the remnant of the core-halo structure of a pre-merging clump hosting the cD galaxy.
In order to further investigate the cluster dynamical state and to
discriminate between the previous pictures, we derived the LFs by
using new photometric data for ABCG 209 based on ESO-NTT imaging in
the B, V and R wavebands. The new photometric data are presented in
Sect. 2. In Sect. 3 we describe the data
reduction, and the photometric calibrations. The aperture photometry
is presented in Sect. 4, whereas Sect. 5 deals with
the LFs and Sect. 6 with the spatial distribution of
galaxies of different luminosity. Section 7 is dedicated to
the summary and the discussion of the results. In this work we assume
H0 = 70 km s Mpc,
,
.
According to this cosmology, 1 arcmin corresponds to 0.205 Mpc at z = 0.209.
New observations of the galaxy cluster ABCG 209 were carried out at the ESO New Technology Telescope (NTT) with the EMMI instrument in October 2001. The data include B-, V- and R-band imaging, plus multi-slit spectroscopy (EMMI-NTT) for 112 cluster members. The spectroscopic data are presented in Paper I.
A field of
(
Mpc2), was observed in the B, V and R bands around
the cluster center. In order to sample the cluster at large distance
from the center, we observed other two fields in V band. The total
observed area in V band, accounting for overlapping, is
160 arcmin2 and is shown in Fig. 1. The relevant
information on the photometry are summarized in Table 1. Standard stars from Landolt (1992) were also
observed before and after the scientific exposures, and were used for
the photometric calibration.
![]() |
Figure 1:
Area covered by the photometry in the region of
ABCG 209. All the galaxies brighter than V = 21.5 are marked with
crosses. The central field (solid contour) was observed in B, V and R bands, while the two adjacent fields (dashed contours) were
observed only in V band. Circles mark the regions analysed in
Sect. 6.1. North is up and East is on the left. The
origin of the coordinates coincides with the cluster center
(
![]() ![]() |
Open with DEXTER |
Standard procedures were employed for bias subtraction, flat-field
correction, and cosmic ray rejection using the IRAF package. For each waveband the flat-field was obtained
by combining twilight sky exposures. After bias subtraction and
flat-field correction the images were combined using the IRAF task
IMCOMBINE with the CRREJECT algorithm. Residual cosmic rays and hot
pixels were interpolated applying the IRAF task COSMICRAYS. The
resulting images show a uniform background with typical rms of 2.3%, 1.4%, and 1.1% for the B, V, and R bands, respectively.
Table 1: Information on the photometric observations.
The photometric calibration was performed into the
Johnson-Kron-Cousins photometric system by using the Landolt standard
fields. The instrumental magnitudes of the stars were measured in a
fixed aperture by using the IRAF packages APPHOT and DAOPHOT. The
aperture size was chosen in order to i) enclose the total flux, ii) obtain the maximum signal-to-noise ratio. By comparing the
magnitudes of the stars in different apertures, we found that a
reasonable compromise is achieved with an aperture of radius 10 pixels
(cf. Howell 1989). For each band, we adopted the following
calibration relation:
As a test we compare the (B-R, V-R) diagram of the Landolt standard stars with that of the stars in the cluster field. We measured the magnitude of these stars by using the software SExtractor (Bertin & Arnouts 1996). We verified that the observed distribution of stars in our images matches that of the Landolt stars in the (B-R, B-V) plane, proving the accuracy of the photometric calibration.
Table 2: Results of the photometric calibration of BVR data.
For each image, a photometric catalog was derived by using the
software SExtractor (Bertin & Arnouts 1996). We measured
magnitudes within a fixed aperture of 5.0'', corresponding to 17 kpc at z = 0.209, and Kron magnitudes (Kron 1980),
for which we used an adaptive aperture with diameter
,
where
is the Kron radius and a is a
constant. We chose a = 2.5, yielding
94% of the total source
flux within the adaptive aperture (Bertin & Arnouts
1996). The measured magnitudes were corrected for galactic
extinction following Schlegel et al. (1998). The uncertainties
on the magnitudes were obtained by adding in quadrature both the
uncertainties estimated by SExtractor and the uncertainties on the
photometric calibrations.
The completeness magnitudes were derived following the method of
Garilli et al. (1999), as shown in Fig. 2. We
estimated the completeness magnitude as the magnitude where galaxies
start to be lost because they are fainter than the brightness
threshold in the detection cell. The completeness magnitudes are
(
),
(
), and
(
).
In order to derive the LF we considered the objects brighter than the
completeness limit and adopted the Kron magnitude, since this is the
best estimate of the total magnitude. The star/galaxy classification
was based on the SExtractor stellar index (SG), defining as stars the
sources with
.
The catalogs for the central region of the cluster were cross-correlated by using the IRAF task XYXYMATCH. For each object, colours were derived within the fixed aperture.
![]() |
Figure 2:
Completeness magnitudes of the B-, V-, and R-band
images are estimated by comparing magnitudes in the fixed aperture (
![]() ![]() ![]() ![]() |
Open with DEXTER |
In order to measure the cluster LF in each band we used all the galaxy
photometric data up to the completeness magnitude, and removed the
interlopers by statistically subtracting the background
contamination. We used galaxy counts in B, V and R bands from the
ESO-Sculptor Survey (Arnouts et al. 1997; de Lapparent et al. 2003) kindly provided to us by V. de Lapparent. These data
cover an area of 729 arcmin2 observed with EMMI instrument. The data reduction was performed by using the procedures
similar to those adopted in the present work. In particular, the
photometric catalog was obtained by SExtractor and the total
magnitudes were estimated from the Kron magnitudes defined by adopting
the same aperture.
We assumed Poisson statistics for the background and cluster field
galaxy counts. The errors on the cluster LFs were computed by adding
in quadrature Poissonian fluctuations. Since the field of view of our
observations is 10 times smaller than the area covered by the
background counts, the errors on the cluster LFs are dominated by
Poissonian errors of cluster counts.
We derived the LFs for the central field (see Fig. 1) in
the B, V and R bands by fitting the galaxies counts with a single
Schechter function. In V and R bands the fits were also computed with
the sum of two Schechter functions in order to describe the bright and
the faint populations (Sect. 5.1). We compared the LFs with
counts obtained selecting red sequence galaxies
(Sect. 5.2). All the fit parameters and the statistics, are listed in Table 3.
Table 3:
Fits to the Luminosity Functions. Errors on the
M* and
parameters can be obtained by the
confidence contours shown in Fig. 3.
Figure 3 (solid lines) shows the LFs in the B, V and R bands for the central cluster region, obtained by a weighted
parametric fit of the Schechter function to the statistically
background-subtracted galaxy counts (filled circles). The parameters
of the fit are:
B* = 20.06,
,
V* =
18.29,
,
R* = 17.78,
.
We evaluated the quality of the fits by means of the
statistics (see Table 3). The single Schechter function
gives a fair representation of the global distribution of the data,
that is the single Schechter fit cannot be rejected in any band even
at the 10% c.l.
On the other hand, there is indication of a dip in the
distribution at
and
.
According to the
fitted single Schechter function, there should be 121 and 125 galaxies
in the range
V = 20.5-21.5 and
R = 20.0-21.0 respectively, whereas
in our counts we find
and
galaxies. If we define the dip amplitude as:
The presence of a dip was discussed in the literature by considering
the properties of the TSLFs (e.g., Binggeli et al. 1988) in
nearby clusters. In a study of the Virgo cluster, Sandage et al. (1985) showed the presence of two distinct classes of
galaxies (normal and dwarfs), with different dynamical and luminosity
evolution. Since the bright-end of the LF is well studied, we can use
a priori the information about the LF shape for bright galaxies
to fit our data with two different Schechter functions, representing
bright and faint galaxies. We assumed a Schechter model with
and
(from N
slund et al. 2000 results, scaled
according to the adopted cosmology). Using the V-R colour term as
derived from Eq. (4), we also fix
and
.
Figure 4 shows the cluster LFs in V and R bands modelled with
the two Schechter functions. The fit procedures yields
,
.
For
the slope of the faint galaxies function the V- and R-band data
provide poor constraints, so that we have only an indication of a
steep faint-end. According to the
statistics (see Table 3), combining two Schechter functions for bright and faint
galaxies the quality of the fit increases both in V and R bands.
![]() |
Figure 3:
Luminosity function in the B, V and R bands in the central
field of
![]() ![]() ![]() |
Open with DEXTER |
We obtained the Colour-Magnitude (CM) relation by fitting the
photometric data of the spectroscopically confirmed cluster members
(see Paper I) with a biweight algorithm (Beers et al. 1990):
By selecting galaxies on the CMR (within the observed scatter), the background contamination is expected to be negligible. So we can directly compare galaxy counts (open circles in Fig. 3) with those derived in Sect. 5.1 (filled circles in Fig. 3). In both V and R bands the counts obtained with the two different approaches are very similar.
Figure 3 (dashed lines) shows the LFs in the V and R bands for the central cluster region, obtained by a weighted
parametric fit of the Schechter function to the red galaxy counts
(open circles). The parameters of the fit are:
V* = 17.87,
,
R* = 17.50,
.
According to the
statistics (see Table 3)
we can reject the fit with a single Schechter function in the V band
at 98% c.l. At the same time, the dip amplitude in the V band,
%, is unchanged respect to the case of the global LF (Sect. 5.1). It turns out that, independently from the fitted
Schechter functions, the ratio (
60%) of observed counts in
the bins inside the dip and in the bins adjacent to the dip regions is
the same for both galaxy samples, while the overall distribution of
red galaxies cannot be described with a single Schechter
function. According to the
statistics (see Table 3), combining two Schechter functions the quality of the
fit increases both in V and R bands, and becomes acceptable also in
the V band.
![]() |
Figure 4:
Luminosity function in the V and R bands in the central
field of
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5: V-R vs. R CM diagram for all the galaxies within the completeness limit in the central field of ABCG 209. Galaxies of the red sequence (see solid line) are plotted as filled circles. The solid line defines the CM sequence for spectroscopically confirmed cluster members. |
Open with DEXTER |
The data in the V band, covering an area of 160 arcmin2,
corresponding to a circular region with equivalent radius 0.6
(
Mpc; see
Paper I), allow to study the environmental dependence of LF and the
spatial distribution of galaxies as a function of the clustercentric
distance.
Figure 6 (upper panel) shows the V-band LF in the
whole observed area, modelled by using a weighted parametric fit to a
single Schechter function, with best fit values
V* = 18.45(
M*V = -22.03) and
.
The LF shape is
very similar to that obtained in the central field (Fig. 3)
and also in this case the Schechter function overestimates the
observed counts in the range
V = 20.5-21.5. The dip amplitude is
%.
![]() |
Figure 6:
Luminosity function in the V band. In the upper panel are
shown the LF in the whole observed area (160 arcmin2) and the 1, 2 and 3![]() ![]() |
Open with DEXTER |
Table 4: Best fit value of the LFs measured in three regions around the center of the cluster.
In order to compare the abundances of galaxies of various luminosities
in different regions we computed the LFs in three areas at different
distances from the center. First we considered the cluster
counts in an area of 0.5
h-270 Mpc2 around the center of
the cluster (region A), and then in two concentric circular rings
around the first central region, each in an area of 1.0
h-270 Mpc2, respectively at
2
and
4.7
from the center (region B and C). We fitted the counts using a single
Schechter function with
fixed at the best fit value obtained
from the LF computed over the whole observed field. Table 4
reports the relative fit values in the different regions and
Fig. 6 shows the fitted functions. We also measured the
luminosity weighted ratio of the number of objects brighter and
fainter than the V = 21.0 (Table 4), that is the magnitude
where the dip occurs. This luminosity ratio increases from region B to the center by a factor
,
indicating a significant
luminosity segregation, whereas the ratio does not vary from the
region B to the region C.
In Fig. 7 (left panel) the number density radial profile is shown, i.e. the number of galaxies measured in concentric rings around the center. It has been computed for the whole galaxy population (open circles) and for galaxies brighter (filled circles) and fainter (triangles) than V = 21.0, the magnitude where the luminosity distribution shows the distinctive upturn. The number of faint galaxies is determined up to the completeness limit. The densities are statistically background-subtracted by using background galaxy density measured in the ESO-Sculptor Survey (Arnouts et al. 1997). Errors are assumed to be Poissonian.
The radial profile of all galaxies is centrally peaked, showing the
highest gradient for bright galaxies (Fig. 7, left
panel). Indeed, the bright galaxy density decreases by a factor 5.4 between the first and the second bin, to be compared with a
decrease of a factor
1.4 for faint galaxies and of a factor
2.6 when all galaxies are considered. For faint galaxies (V >
21.0), the distribution decreases slowly with clustercentric distance.
Figure 7 (right panel) shows the bright-to-faint ratio (BFR) as a function of the clustercentric distance. The BFR shows a maximum at the center, where the number of bright galaxies is two times the number of faint galaxies, and decreases rapidly between the first and the second bin. Outside the central region of the cluster (r > 1 arcmin) the BFR distribution is flat and there are roughly two faint galaxies for each bright galaxy.
The analysis performed on ABCG 209 can be summarized by the following results:
![]() |
Figure 7:
Background-subtracted radial profile of ABCG 209 is
shown in the left panel. Open circles represent the distribution
of all galaxies in the cluster, filled circles are bright galaxies
(
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Dips are found at comparable absolute magnitudes, in different clusters, suggesting that these clusters may have similar galaxy population. However, the dips have not necessarily the same shape: the dips found in the the LF of Shapley 8 and Coma are broader than those found in Virgo, ABCG 963, and ABCG 85. As already mentioned in the introduction, this can be explained with the hypothesis of an universal TSLF (Binggeli et al. 1988), whereas the different relative abundances of galaxy types, induced by cluster-related processes, are at the origin of the different dip shapes (see also Andreon 1998).
Moreover, as shown in Table 5, there is no correlation between
the dip position and the redshift, suggesting little or no
evolutionary effect in the dip, in agreement with N
slund et al. (2000), which directly
compare the LFs of Coma, ABCG 963 and MS 2255. They detected no
qualitative difference between nearby and distant clusters.
Studying the R-band LFs for a photometric sample of 10 clusters at different redshifts and with different richness classes, Yagi et al. (2002) demonstrated that the dips seen in the LFs are almost entirely due to r1/4-like galaxies. This evidence is in agreement with those shown in Fig. 3 (open circles). In fact by comparing the counts of red sequence galaxies with those obtained with a statistical background subtraction, the dip in V band appears more pronounced.
Environmental effects can be at the origin of the different positions of the dips seen in different clusters, since the transition from the bright-dominated to the faint-dominated parts of the LF can occur at different magnitudes in different environment. The cluster morphological mix and the morphology-density relation (Dressler 1980) should give rise to LFs with different shape when subdividing a sample galaxies respect to the cluster environments.
The effect of the environment can be seen in ABCG 209 by comparing
the shape of the LF in V band in different regions around the cluster
center (see Fig. 6 and Sect. 6.1). The fitted M* value is shifted toward fainter magnitudes going from the inner (A) to the outer (C) region and the luminosity ratio of
bright-to-faint galaxies decreases by a factor 2, indicating
a luminosity segregation.
The bright galaxies are markedly segregated in the inner 0.2 h-170, around the cD galaxy. This suggest that bright galaxies could trace the remnant of the core-halo structure of a pre-merging clump.
Table 5: Compilation of dip positions in the restframe R band in cluster LFs from literature. Data were transformed to R-band absolute magnitudes by using k- and evolutionary correction, and average colours for early-type galaxies of Poggianti (1997).
Although ABCG 209 is a cD-like cluster, with cD galaxy located in
the center of a main X-ray peak, it shows an elongation and asymmetry
in the galaxy distribution (Paper I). Moreover, the faint-end slope
turns out to be
at more than 3
c.l. in both V and R bands, thus reconciling the asymmetric properties of X-ray
emission with the non flat-LF shape of irregular systems as found by
Lopez-Cruz et al. (1997).
These results allow to discriminate between the two possible formation scenarios suggested by the dynamical analysis (Paper I). We conclude that ABCG 209 is an evolved cluster, resulting from the merger of two or more sub-clusters, while the elongation and asymmetry of the galaxy distribution (of the X-ray emission) and the shape of the LFs show that ABCG 209 is not yet a fully relaxed system.
Our analysis is in agreement with the existence of i) an universal LF for bright galaxies, which is well described by a flat Schechter function, and ii) a steep Schechter function for faint galaxies. A definitive conclusion regarding the faint-end slope of the LF needs deeper photometry able to sample the luminosity distribution of dwarf galaxies.
Acknowledgements
We thank Valerie de Lapparent who provided us with the galaxy counts used to derive the cluster LFs. A. M. thanks Massimo Capaccioli for the hospitality at the Osservatorio Astronomico di Capodimonte, and Francesca Matteucci for support during this work. This work has been partially supported by the Italian Ministry of Education, University, and Research (MIUR) grant COFIN2001028932: clusters and groups of galaxies, the interplay of dark and baryonic matter, and by the Italian Space Agency (ASI).