C. Ferrari 1,2 - C. Benoist 1 - S. Maurogordato1 - A. Cappi 3 - E. Slezak 1
1 - Laboratoire Cassiopée, CNRS/UMR 6202, Observatoire de la Côte d'Azur, BP 4229,
06304 Nice Cedex 4, France
2 - Institut für Astrophysik,
Technikerstraße 25, 6020 Innsbruck, Austria
3 - INAF,
Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna,
Italy
Received 6 August 2004 / Accepted 31 August 2004
Abstract
We present the analysis and results of a new VRI photometric and
spectroscopic survey of the central
region of the galaxy cluster A3921 (z=0.094). We detect
the presence of two dominant clumps of galaxies with a mass ratio of
5: a main cluster centred on the Brightest Cluster Galaxy
(BCG) (A3921-A), and an NW sub-cluster (A3921-B) hosting the second
brightest cluster galaxy. The distorted morphology of the two
sub-clusters suggests that they are interacting, while the velocity
distribution of 104 confirmed cluster members does not reveal strong
signatures of merging. By applying a two-body dynamical formalism to
the two sub-clusters of A3921, and by comparing our optical results
to the X-ray analysis of A3921 based on XMM observations (Belsole
et al. 2005), we conclude that A3921-B is probably tangentially
traversing the main cluster along the SW/NE direction. The two
sub-clusters are observed in the central phase of their merging
process (
0.3 Gyr), with a collision axis nearly perpendicular
to the line of sight. Based on the spectral features of the galaxies
belonging to A3921 we estimate the star formation properties of the
confirmed cluster members. Substantial fractions of both
emission-line (
13%) and post-star-forming objects (so called
k+a's,
16%) are detected, comparable to those measured at
intermediate redshifts. Our analysis reveals a lack of bright
post-star-forming objects in A3921 with respect to higher redshift
clusters, while the fraction of k+a's increases towards fainter
magnitudes (
). Similar results were obtained in
the Coma cluster by Poggianti et al. (2004), but at still fainter
magnitudes, suggesting that the maximum mass of actively
star-forming galaxies increases with redshift ("downsizing
effect''). The spatial and velocity distributions of k+a galaxies do
not show significant differences to those of the passive population,
and to the whole cluster. Most of these objects show relatively red
colours and moderate Balmer absorption lines, which suggest that
star formation has ceased
1-1.5 Gyr ago. Their presence is
therefore difficult to relate to the on-going merging event. We find
that star-forming galaxies share neither the kinematics nor the
projected distribution of the passive cluster members. Moreover,
most emission-line galaxies are concentrated in A3921-B and in the
region between the two sub-clusters. We therefore suggest that the
ongoing merger may have triggered a star-formation episode in at
least a fraction of the observed emission-line galaxies.
Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 3921 - galaxies: distances and redshifts - cosmology: observations
In the standard cosmological scenario of hierarchical structure formation, bound objects form from the collapse of initial density fluctuations that grow under the influence of gravity through merging of smaller structures that have formed before. As optical and X-ray studies reveal that clusters of galaxies are still forming at the present epoch (e.g. Jones & Forman 1992; West et al. 1995; Donnelly et al. 2001), merging clusters provide a unique tool to test or analyze the physics of structure formation and evolution.
Table 1: Imaging: summary of the observations.
Major cluster-cluster collisions are the most energetic events that
have occurred in the Universe since the Big Bang, as they release
total energies up to
erg (Sarazin 2003), and
their effects on all cluster components are far from being fully
understood. The intra-cluster gas experiences compression, rarefaction
and shock waves (Schindler & Müller 1993), and shows characteristic
features such as cold fronts (Markevitch et al. 2000). Additionally,
as regards galaxy distribution, the velocity dispersion of the cluster
members can increase up to a factor of two during the merger event
(Schindler & Böhringer 1993). While the effects on the
intra-cluster medium and on cluster internal dynamics have been
analyzed in some detail, the effects on the galaxies are still
debated. Caldwell et al. (1993) showed that strong Balmer-line
absorption galaxies are distributed between the main Coma cluster and
the merging SW sub-cluster, suggesting that star formation is
triggered by the merging event. However, various physical mechanisms
have been shown to affect the process of star formation within
clusters. Changes of the tidal gravitational field during cluster
merging could affect star formation in galaxies (e.g. Bekki 1999). The
increasing external pressure following the infall of galaxies into the
dense intra-cluster gas may trigger star formation (Dressler & Gunn
1983; Evrard 1991), while gas stripping in galaxies due to
ram-pressure exerted by the ICM could weaken the starburst phenomenon
during cluster-cluster collisions (Fujita et al. 1999). So far it is
not clear which of these competing effects is the dominant one, since
results in favour of the former (e.g. Abraham et al. 1996; Wang et al. 1997; Moss & Whittle 2000; Flores et al. 2000;
Gavazzi et al. 2003; Poggianti et al. 2004) and of the latter (e.g.
Tomita et al. 1996; Balogh et al. 1997, 1998; Baldi et al.
2001) were obtained.
Combined optical and X-ray studies have been particularly successful in revealing the complex dynamics of merging clusters (e.g. Davis et al. 1995; Lemonon et al. 1997; Roettiger et al. 1998; Durret et al. 1998; Arnaud et al. 2000; Donnelly et al. 2001; Ferrari et al. 2003). The X-ray properties of merging clusters are currently being investigated through a set of XMM observations (XMM guaranteed time, Sauvageot et al. 2001), together with the properties of their galaxy distribution through optical observations by our group, hence providing a unique combined analysis. In this paper we will concentrate on the optical analysis of the merging cluster Abell 3921, and our results will be compared to conclusions obtained from the X-ray analysis of XMM data by Belsole et al. (2005).
Abell 3921 is a R=2, BM II Abell cluster at z=0.094 (Katgert
et al. 1998). Previous ROSAT and Ginga observations revealed the
presence of a main cluster and a substructure with a very perturbed
morphology, interpreted as falling onto the main component (Arnaud
et al. 1996). Following XMM-Newton/EPIC observations (Sauvageot
et al. 2001) confirmed the presence of a two-component structure. This
motivated new optical observations in order to better characterize the
merger scenario of this complex cluster. In this paper we analyze the
spectroscopic and photometric galaxy catalogues of the central
region of the cluster, with new data
of multi-object spectroscopy (239 new spectra) and VRI-band imaging,
obtained with EFOSC2 at the 3.6 m ESO telescope and with WFI at the
2.2 m ESO telescope respectively. Section 2 briefly describes the
observations, the data processing technique and the completeness level
achieved by the new spectroscopic sample. In Sect. 3 cluster members
are identified and in Sect. 4 the optical morphology of the cluster is
studied. In Sect. 5 we perform a kinematical and dynamical analysis of
the cluster. The optical masses of the two main subclusters are
estimated in Sect. 6, and we solve the two-body problem for these two
systems. The photometric and spectral properties of the cluster
members are investigated in Sect. 7, making it possible to define several
subsamples whose spatial and velocity distributions are analyzed. The
main results and their interpretation are summarized in the final
Sect. 8. All numbers are expressed as a function of h75, the
Hubble constant in units of 75 km-1 s-1 Mpc. We have used the
CDM
model with
and
,
thus 1 arcmin
corresponds to
0.097
h75-1 Mpc in the following.
The optical observations of Abell 3921 were carried out using the Wide
Field Imager (WFI) mounted at the Cassegrain focus of the MPG/ESO 2.2 m
telescope at La Silla observatory. WFI is a mosaic camera with
4 2 CCD chips covering a total area of 34
33
arcmin2. To cover the gaps between the eight individual chips of
the camera we adopted a standard dithering sequence. The field centred
on
was observed in
service mode in the V, R and I passbands. The data, including
photometric calibration images, were obtained between August 20th and 23rd, 2001. In Table 1 we summarize the observations.
The data reduction was performed using the package "alambic'' developed
by B. Vandame based on tools available from the multi-resolution visual
model package (MVM) by Bijaoui and collaborators (Bijaoui & Rué
1995; Rué & Bijaoui 1997). The photometric calibration was obtained
using several standard stars from Landolt (1992) over a large range of
airmasses, leading to an accuracy of the photometry of 0.05 mag.
The SExtractor software package (Bertin & Arnouts 1996) was used to identify
all sources in the field as well as to classify them and measure their
magnitudes. The latter were corrected for galactic absorption using
E(B-V)=0.027 as derived from Schlegel et al. (1998), yielding AV=0.09 mag,
AR=0.07 mag and AI=0.05 mag. They were also transformed to the AB system
given by the following relations:
,
,
.
In this paper we present the results of the analysis of new
spectroscopic data obtained through two sessions of observations at
the ESO 3.6 m telescope (2 nights in September 2001, 2 nights in
October 2002). We used the ESO Faint Object and Camera (EFOSC2) with
grism#03 and a punching head of 1.35'', obtaining a spectral
resolution of
Å over the wavelength range
3050-6100 Å. For each frame we made at least two science exposures
in order to eliminate cosmic rays, with an integrated exposure time of
5400 s for brighter objects (
)
and of 7200 s for
fainter ones (
). We made a standard
spectroscopic reduction using our automated package for multi-object
spectroscopy based on the task "apall'' in IRAF
. Spectra were wavelength calibrated using the
helium-argon lamp spectra taken after each science exposure. We
determined redshifts using the cross-correlation technique (Tonry &
Davis 1981) implemented in the task "xcsao'' of the RVSAO
package. Spectra of late-type stars were used as radial velocity
standards.
In Table A (available at the CDS) we
list the results of our spectroscopic observations in the following way:
Col. 1: identification number of each target galaxy; Cols. 2 and 3: right
ascension and declination (J2000.0) of the target galaxy; Cols. 4-6:
the V, R and I band magnitudes in the AB system; Cols. 7 and 8: best estimate
of the heliocentric redshift (expressed as cz) and associated error from the
cross-correlation technique (those values have been set to "-2'' if the
object is a star and to "-1'' if we have no redshift information); Col. 9:
run of observations; Col. 10: a quality flag for the redshift determination:
1 = good determination (R parameter of Tonry & Davis 3), 2 = uncertain
determination, 3 = no determination; Cols. 11 and 12: the equivalent widths
(negative or positive in the case of emission or absorption features
respectively) for [OII]
3727 and
(
Å)
lines, measured only for the cluster members; Col. 13: the spectral
classification of the galaxy members.
A total of 239 new spectra has been obtained, among which 56 are
stars, 100 are galaxies with a very good redshift determination and 83
are spectra with R<3. In the following analysis, we also consider
the redshifts measured in the region of A3921 by other authors (ENACS
collaboration, Katgert et al. 1996; Mazure et al. 1996). We compare
the values of the redshifts of 11 galaxies obtained from the common
sample to those published (35) and we obtain a mean difference of
km s-1, which shows a good consistency between the two
datasets. The final sample includes our 239 redshifts and 24 available redshifts from the ENACS catalogue (Katgert et al. 1996;
Mazure et al. 1996), among which 22 with quality flag = 1 following our
criterion. It covers a field of
,
slightly
larger than the size of our own spectroscopic follow-up
(
).
The R-band magnitude distribution of the whole spectroscopic sample
of galaxies is presented in Fig. 1 (207 galaxies), and
compared to the distribution of galaxies with a good redshift
determination, i.e. quality flag = 1 (122 galaxies), and among them,
those belonging to A3921 (104 galaxies, see Sect. 3).
![]() |
Figure 1: R-band magnitude distributions of the galaxies of our spectroscopic sample (207 objects - dotted line), of all the galaxies with good velocity determination (122 - dashed line), and, among them, of those belonging to A3921 (104 - shaded area). |
Open with DEXTER |
The ratio of the number of galaxies with measured redshift to the
total number of galaxies detected within the central field of
covered by our last observations is
plotted in Fig. 2 as a function of the R-band
magnitude. The spectroscopic catalogue is complete to better than 50%
up to
(
).
![]() |
Figure 2:
Completeness of the spectroscopic sample as a function of the
R-band magnitude in the central
![]() |
Open with DEXTER |
Among the 263 objects of our spectroscopic sample we will consider in the following only the 122 galaxies with redshifts determined with a quality flag = 1 (i.e. our new 100 high quality spectra plus 22 published previously).
![]() |
Figure 3: Apparent radial velocity (cz) histogram of the galaxies in the central field of the cluster A3921, with a binning of 500 km s-1. |
Open with DEXTER |
Figure 3 shows the radial velocity (cz) distribution of
this dataset in bins of 500 km s-1. The bulk of the cluster is
concentrated between 25 400 and 30 400 km s-1. Four objects have czhigher than 75 000, one galaxy appears to be in the foreground, while
13 are located between 31 000 and 50 000 km s-1 and most of them (12) are
also spatially concentrated in the region around the two brightest
objects of the West side of the cluster, BG2 and BG3. If we consider
only the 8 galaxies around 40 000 km s-1, they correspond to a peak
located at
km s-1 (
)
with a velocity
dispersion
km s-1, indicating that they could be a
background group.
In order to eliminate the galaxies not belonging to the cluster, we
apply the standard iterative 3
clipping (Yahil & Vidal 1977).
104 cluster members with cz between 25 400 and 30 400 km s-1 are
selected. Note that the background peak being located
10 000 km s-1behind the mean velocity of the cluster is excluded as a subgroup of
the cluster. Figure 4 shows the location of the 104 identified
members superimposed on a fraction of the R-band image.
![]() |
Figure 4:
Galaxies identified as members of A3921 (squares), superimposed on a
fraction of the R-band image (
![]() |
Open with DEXTER |
![]() |
Figure 5:
Projected galaxy density maps (with
![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Velocity distributions of the main structures of the projected
galaxy density map (same as the right panel of
Fig. 5). The dotted lines in the histograms indicate
the relative velocities of the neighbouring brightest galaxy. At the
level of
![]() |
Open with DEXTER |
![]() |
Figure 7:
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Projected galaxy density maps (using red sequence galaxies and excluding
galaxies known not to be cluster members from spectroscopy, as in
the right panel of Fig. 5) on a
![]() ![]() |
Open with DEXTER |
In order to investigate the projected spatial morphology of A3921
several density maps of the galaxy distribution have been built on the
basis of a multi-scale approach. The adopted algorithm is a 2D generalization of the algorithm presented in Fadda et al. (1998). It
involves a wavelet decomposition of the galaxy catalogue performed on
five successive scales from which the significant structures are
recombined into the final map (following the Eq. (C7) of Fadda et al. 1998). These significant structures are obtained by thresholding
each wavelet plane at a level of three times the variance of the
coefficients of each plane except for the two smallest scales for
which the threshold is increased to four and five times the variance
in order to reduce false detections due to the very low mean density
of the Poisson process at these scales (0.01 for a chosen grid of
).
Such density maps are presented in Fig. 5 for three input
catalogues. In the left panel, all galaxies with
(
)
are used, leading to a map revealing two
dominant clumps in the central part of the field (hereafter A3921-A in
the centre and A3921-B to the NW) in the central
.
In order to avoid possible
projection effects, we isolate galaxies likely to be early types at
the same redshift on the basis of their colour properties. Indeed,
despite the bimodal structure of A3921 one can notice in the two
colour magnitude diagrams of Fig. 7 a well defined red
sequence, the characteristic linear structure defined by the bulk of
early-type galaxies in a cluster. The determination of the slope,
intercept and width of these "red sequences'' is performed using the
technique described in Appendix A. The central panel of
Fig. 5 shows the resulting red sequence density map
(keeping only galaxies at
around the red sequence),
where one can notice that several small-scale clumps in the cluster
core disappear as well as in its surroundings, leaving the
larger-scale structure of the two main clumps unaffected.
In order to go one step further in avoiding projection effects we take
advantage of our spectroscopic follow up. Even though it is complete down
to a level of 50% at
,
most of the clumps seen
visually have been spectroscopically sampled allowing us to discard
possible external groups. For example, as mentioned in the previous
section, the background redshift-peak group may lead to an
over-density unrelated to the cluster in the iso-density
map. Therefore, from the sample of red sequence galaxies we
additionally exclude galaxies known from spectroscopy not to be
cluster members. The result is shown in the right panel of
Fig. 5, which is quite similar to the map in the central
panel except that the location of clump B is shifted towards the South and
its central peak is shifted close to the position of the BG2. If we
can trust at a high confidence level that the sub-structures still
visible are part of A3921, it is still not clear whether or not they
are projection effects within the cluster. In Fig. 6
the relative velocity distributions of these various sub-structures
are presented using the available redshifts. At the level of
these clumps actually have full spectroscopic coverage
except for the eastern group with 5 confirmed cluster members out of 7
galaxies belonging to the red sequence. The strong clustering in
redshift space for all the sub-structures confirms the accuracy of the
picture reflected by this density map.
![]() |
Figure 9:
Projected galaxy density map (
![]() ![]() |
Open with DEXTER |
In Fig. 8 the same maps are presented for three magnitude cuts. The global bi-modal structure of the cluster remains unchanged with luminosity. This is particularly true for the central parts of clumps A and B. However, note that the eastern group, more compact than clump B, is stronger at fainter levels.
The results given above are based on the only red sequence galaxies. In Fig. 9 the galaxies bluer than the red sequence are presented, superimposed on the red sequence iso-density map. These blue galaxies appear much less clustered than the redder ones over the whole field and show only little correlation with the red density peaks except in the case of clump B showing an excess of blue galaxies relative to the other clumps and in particular to the center of clump A. This asymmetry is also present in the distribution of emission line galaxies as will be shown and discussed in Sect. 7.3.1.
To summarize, A3921 is characterized by a) a bimodal morphology; b) the presence of several substructures inside each of the two main clumps; c) an offset of the Brightest Cluster Galaxy (BCG) from the main density peak of clump A, and d) an excess of blue galaxies around the clump B. These results suggest that this system is out of dynamical equilibrium and that it is probably composed of a main cluster interacting with at least two groups, one to the North-East (clump B) and one to the East, the latter being significantly less luminous than the former. In the following section we will analyze the dynamical properties of A3921 in order to understand which phase of the merging process we are witnessing.
Using the biweight estimators for location and scale (Beers
et al. 1990), we find a mean apparent velocity of
km s-1, corresponding to a mean redshift of
,
and a velocity dispersion of
km s-1 (at
significance level, see
Table 2). Figure 10 shows the histogram of
the cosmologically and relativistically corrected velocity offsets
from the mean cluster redshift (
)
of the 104 cluster members. Contrary to the projected morphology on
the sky we do not find any bi-modal structure in the velocity
distribution. In dissipationless systems, gravitational interactions
of cluster galaxies over a relaxation time generate a Gaussian
distribution of their radial velocities; possible deviations from
Gaussianity could provide important indications of ongoing dynamical
processes. In the following, we are therefore interested in measuring
a possible departure of the observed cluster velocity distribution
from a Gaussian.
A velocity distribution with slight tails can indicate the presence of two or more overlapping Gaussian components in the whole velocity histogram, while an asymmetric distribution can be the result of the contemporary presence of subclusters with different numbers of galaxies (Ashman et al. 1994). We use two kinds of shape estimators: the traditional third and fourth moments, i.e. skewness and kurtosis, and the asymmetry and tail indices (Bird & Beers 1993). In Table 2 we report the corresponding values and significance levels, estimated from Table 2 in Bird & Beers (1993) under the null hypothesis of a Gaussian distribution.
The values obtained for Skewness and AI show that the velocity
histogram is quite symmetric (significance level >10%). Values
observed for kurtosis and TI indicate a heavy-tailed distribution,
rejecting the Gaussian hypothesis for the first (i.e. kurtosis) at
10% significance level, and for the second at even better than 1%. While light-tailed distributions indicate multi-modality, heavily
populated tails could be due to contamination by non-cluster
galaxies. We think that the analysis of Sect. 3
excludes this possibility, but we apply an additional test in order to
definitively reject the presence of outliers in our velocity sample of
104 galaxies. Extensive data available for low-z clusters show that
most (
95%) of the galaxies in the central region of the
clusters have radial velocities within
3500 km s-1 of the mean
cluster redshift (Postman et al. 1998). We thus verify that all the
104 galaxies of our final sample have
km s-1 (with
defined at the beginning of this section). In any case,
the shape parameters do not show strong evidence of subcomponents in
velocity space.
7 of the 13 normality tests performed by ROSTAT reject the Gaussian hypothesis at better than 10% significance level (see Table 3). Moreover, we find one significant weighted gap in our dataset (Beers et al. 1990) at cz=27381.9 km s-1 (shown in Fig. 10) and characterized by a normalized size of 2.47 and a significance of 3%. On the contrary, among the 6 tests performed by ROSTAT that do not reject the Gaussian hypothesis, the DIP test (Hartigan & Hartigan 1985) accepts the unimodal hypothesis at better than 99% level, in agreement with the conclusions obtained with the symmetry tests (i.e. Skewness and AI).
Classical tests of Gaussianity therefore give controversial results for the kinematical properties of A3921.
Because of the uncertain results obtained in the previous section we apply the KMM mixture-modeling algorithm of McLachlan & Basford (1988) which, using a maximum-likelihood technique, assigns each galaxy to a possible parent population, and evaluates the improvement in fitting a multiple-component model over a single-one. First of all, we try to allocate the 104 galaxies of our dataset into two possible subclusters, respectively with mean velocity lower and higher than the gap position. Several tests are then performed by changing both the estimated mean velocities and the estimated mixing proportion for each group, but we always obtain the same result: the KMM algorithm tries to fit a 2-group partition from these guesses, obtaining intermediate confidence levels that cannot reject the null hypothesis of unimodal distribution. This result actually agrees both with the DIP test and with what is observed from the shape parameters of the velocity distribution, the absence of light-populated tails and of a significant asymmetry in the histogram excluding the possible presence of several overlapping subunits populated differently.
Classical statistical tests for sub-clustering are then applied to the 104 cluster members in order to look for the existence of correlated substructures in velocity and spatial distribution (Dressler & Shectman 1988; Bird 1994; West & Bothun 1990). The results, obtained by the bootstrap technique and 1000 Monte Carlo models, are presented in Table 4.
Table 2:
Properties of the cz distribution and significance levels of
shape estimators for various subsamples of A3921.
and
are the mean velocity and the velocity dispersion of the
different distributions (biweight estimators for location and scale,
Beers et al. 1990).
Neither
nor
parameters find evidence of
substructures with a high significance level, while the
test
has an intermediate and inconclusive value. Figure 12
shows graphically the results obtained with the
test.
In conclusion, while no bi-modality is shown in the velocity distribution, the various tests of Gaussianity and subclustering do not reach a consensus, but do not show "extreme'' departures from Gaussianity.
The analysis of the mean velocity and particularly of the radial
profile of the velocity dispersion provides a useful tool for
investigating the dynamics of galaxy clusters, as they can reveal
signs of subclustering and ongoing merging (Quintana et al. 1996;
Muriel et al. 2002). Such integrated measurements are therefore
performed with the 104 members of A3921 up to a distance of
2
h75-1 Mpc from the cluster centre (taken to be at the position
of BG1), and are presented in the left panels of Fig. 11. The
mean velocity has a very constant value (
28 000 km s-1) from the
centre of the cluster to its outer edges. Following the
classification of the radial profiles of the velocity dispersions of
den Hartog & Katgert (1996), A3921 presents an "inverted'' shape: the
cumulative velocity dispersion shows an initial increase with radius
up to
4 arcmin away from the cluster centre and then decreases
nearly to its first bin value (
900 km s-1). Finally, the velocity
dispersion becomes flat in the external regions of the cluster
(
1
h75-1 Mpc). This result suggests that this final value is
representative of the total kinetic energy of the cluster members
(Fadda et al. 1996). A similar shape was detected by den Hartog &
Katgert (1996) and Nikogossyan et al. (1999) for the galaxy cluster
A194, and these authors interpreted such a profile as originating in a
nearly relaxed region. Therefore, neither the integrated radial profiles
of the velocity dispersion nor those of the mean velocity reveal the
presence of significant velocity gradients that could be produced by
the presence of internal sub-structures. These results are confirmed
by the right panels of Fig. 11, where mean velocities and
velocity dispersions are estimated in rings containing the same number
of galaxies (13). The shapes of these differential profiles agree with
the integrated ones, again presenting an inverted shape. Its minimum
velocity dispersion corresponds to the ring containing the group of
galaxies located around BG2. Within the errors, the mean velocities are
nearly constant around 28 000 km s-1, as for the integrated profile.
As shown above, the dynamical and kinematical properties of the whole
cluster do not reveal strong signatures of merging. More indications
on the dynamical state of the cluster could come from the analysis of
the velocity distribution of the two main subclusters detected
separately on the iso-density maps. For this, our velocity sample is
divided into two datasets, containing the confirmed cluster members in
the two circles displayed in Fig. 14, chosen as the largest
not intersecting. The radius of the two circles is 0.34 Mpc
and they are centred on the density peaks previously detected
(Sect. 4). We plot in Fig. 13 the ratio of
the number of galaxies with very good redshift determinations (QF = 1)
to the total number of objects detected in the two circles of
Fig. 14 as a function of R-band magnitude; the
spectroscopic sampling of these two regions, and in particular of
A3921-A, appears to be quite good, with completeness levels
50%
for
(
).
Figure 15 shows the velocity distributions of the two
datasets and in Table 2 we quote their velocity means and
dispersions. Both datasets show a mean velocity very close to each
other and to the whole cluster value (the velocity offset
between the mean velocities of A3921-A and A3921-B is only
89 +155
-177 km s-1); clump A is characterized by the highest velocity
dispersion.
Most (8) of the 13 normality tests contained in the ROSTAT package accept the hypothesis of Gaussianity for the velocity distribution of the A3921-B component, and actually none of them rejects the null hypothesis for the radial velocities of the clump A3921-A (Table 3). Shape parameters (Table 2) accept the Gaussian hypothesis at more than 10% significance level, with the exception of the skewness and of the kurtosis, which indicates the possible presence of an asymmetric velocity distribution (significance level lower than 10%) and of heavily populated tails (significance level lower than 5%) in the A3921-B velocity histogram.
Summarizing, Gaussianity tests suggest that the merging event has not
strongly affected the internal dynamics of the two sub-clusters.
![]() |
Figure 10:
Top: stripe density plot of the radial velocity offsets of the
104 A3921 members from ![]() |
Open with DEXTER |
Table 3: 1-D statistical tests performed in the ROSTAT package that exclude the hypothesis of a single Gaussian distribution for the different velocity datasets considered in the paper. In Cols. 1 and 2 we report the name and the value of the statistics, while Col. 3 indicates their significance levels.
Table 4:
3-D substructure indicators for the sample of 104 objects with
quality flag = 1 in the velocity range 25 400 30 400 km s-1.
Determining the respective masses of A3921-A and A3921-B requires idealized assumptions which cannot be strictly valid for an interacting system. However, as seen in Sect. 5.4, the dynamics of the central regions of the two clumps appear to be relatively unaffected by the merging event.
We have therefore assumed that each subcluster is virialized. As in
the case of velocity dispersions, the virial radius of each clump was
estimated selecting only those galaxies within a projected distance of
0.34
h75-1 Mpc from its center (see Fig. 14). In order
to have better statistics and at the same time to minimize biases due
to spectral incompleteness and background contamination, we included
all the galaxies belonging to the red sequence (thus assuming that
early-type galaxies trace the mass profile; see Katgert et al. 2004), except for those with a measured redshift and identified
as outliers on the basis of our previous analysis.
![]() |
Figure 11:
Profile of the velocity means ( top) and dispersions ( bottom) in the
cluster: left: as a cumulative function of galaxies centred on the
position of BG1 - right: differential values as a function of radius
in rings with 13 galaxies; the last point is far from the others due
to the lower degree of completeness reached by spectroscopic data in
the outer regions of the cluster (beyond ![]() |
Open with DEXTER |
The mass was calculated with the classical virial equation:
In Eq. (2) U is the potential energy
of the system and
is the mean harmonic radius, defined as:
where
is the separation between the ith and jth
galaxies, and N is the total number of objects in the system).
We have to apply the above relations to projected separations, so that
we can estimate the projected mean harmonic radius
and obtain
the corresponding projected virial radius
.
Assuming
spherical symmetry,
(see Limber &
Mathews 1960) and
,
where
is
the radial velocity dispersion of the system.
The above method is based on the pairwise estimator of the harmonic
radius: an alternative is the so-called ringwise estimator (see
Carlberg et al. 1996). We have applied both methods, finding no
significant difference (notice that this is not true in general).
![]() |
Figure 12: Projected positions, centred on the coordinates of BG1, of the galaxies with vr=25 400-30 400 km s-1. Each galaxy is represented by a circle, whose dimension is weighted by the estimator of Dressler & Shectman (1988). Filled circles (red in the electronic version of the paper) are centred on the positions of the three BGs. |
Open with DEXTER |
![]() |
Figure 13: Completeness of the spectroscopic sample for the two clumps A3921-A (solid line) and in A3921-B (dashed line). |
Open with DEXTER |
![]() |
Figure 14:
Iso-density contours of the projected distribution of the
red-sequence galaxies (with
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 15: Velocity histogram, with a binning of 200 km s-1, of the galaxies in A3921-A (solid line) and in A3921-B (dotted line), according to the division showed in Fig. 14. The relative Gaussian best-fits to the velocity distributions are superimposed. Arrows show the radial velocities of the three brightest cluster galaxies. |
Open with DEXTER |
The results obtained applying the pairwise method along with the
corresponding 1 errors, are shown in Table 5. The
values of the virial radii are relatively small (only slightly larger
than the window radius): we have to take into account the possibility
that they are underestimated, although they could be appropriate for
subcluster systems.
Given our cosmology with
,
we would
expect
(see e.g. Eke et al.
1996), where
is the mean density
within
and
is the critical density. The ratio
between
these two quantities is related to the virial radius according to the
following equation (always derived for our flat cosmology):
Putting into the above equation the mean cluster redshift and
the estimated values of the three-dimensional virial radii,
,
we find
and
respectively for A3921-A and A3921-B, which
indicates that we are probably underestimating the virial radii. A
further uncertainty is associated to the surface pressure term, which
acts in the opposite direction with respect to the previous correction
for the underestimate of the virial radius (Carlberg et al. 1996, and
references therein). However, while the exact values of the virial
masses might be larger, it is clear that A3921 is made up of two
unequal mass systems, with
,
a value which essentially
reflects the factor 2 ratio between the velocity dispersion of the two
clumps.
Table 5:
Columns 1 and 2: projected virial radii
of A3921-A and A3921-B; Cols. 3 and 4: virial mass estimates for A3921-A and
A3921-B; Col. 5: mass ratio.
Masses are in
units and radiii in h-1 Mpc units.
In this section we apply the two-body dynamical formalism (Gregory
& Thomson 1984; Beers et al. 1992) to the two sub-clusters
of A3921. By assuming that at t=0 the two subclusters were at zero
separation parametric solutions to the equations of motion can be
derived both in the bound:
where
is the separation of the two subclusters at maximum
expansion, M is the total mass of the system,
is the developmental
angle and
is the asymptotic expansion velocity. The relative
velocity V and the spatial separation R between the substructures are
respectively related to their radial (
)
and projected components (
)
through the following relations:
We close the system of equations taking our measured values of
Mpc and
km s-1, and assuming different
values for t0, i.e. the epoch of the last encounter between the two
clumps (Gregory & Thompson 1984; Barrena et al. 2002). This
implies in each case a relationship between the total mass of the
system M and the angle
which is displayed in Fig. 16.
In the first case, we set
t0 = 12.6 Gyr, i.e. the age of the Universe in
our cosmology. This suggests that subclusters are moving apart or together for
the first time. We are then in the pre-merger scenario. In the case of A3921,
the measured total mass of the two sub-clusters is
(Sect. 6.1). The top left
panel of Fig. 16 shows that the possible solutions for this value are
bound. In the incoming case, two solutions are possible, with respectively a
very high and a very low value for the projection angle
,
while in the
outgoing case only a high value of
is allowed (cases (a)
Table 6).
We then suppose we are witnessing the sub-clusters after the collision,
testing the possible solutions for t0 = 0.3 Gyr (case (b)), 0.5 Gyr (case (c)), and 1 Gyr
(cases (d)). In all the cases, the possible solutions are
bound as shown in Fig. 16 (top, right and bottom panels). For both the
smaller values of t0 (0.3 and 0.5 Gyr), there is only one outgoing solution
possible, leading to increasing values of
for higher t0 (see
Table 6). At larger times after merging (t0=1 Gyr) three
values of
are possible, two for incoming solutions and one in the
outgoing configuration (see Table 6). At this stage, apart from
the fact that the sub-clusters are bound with a very high probability, it is
difficult to discriminate between the various merging scenarios from the
dynamical analysis of the galaxies alone.
![]() |
Figure 16:
The total mass of the two clumps as a function of the angle ![]() ![]() |
Open with DEXTER |
In the following we wish to assess the star formation properties of
the galaxies belonging to A3921, and in particular to identify
galaxies with signatures of recent star activity from an older
population, as well as galaxies with emission lines, reflecting
present star formation. This analysis is based on our high S/N spectra
catalog, therefore it is restricted to 83 A3921 members. To reach
this goal, the first step is to classify objects of our spectroscopic
catalog of galaxies. Various methods have been used in the
literature, using generally the presence and strength of the [OII]
(
Å) line and of Balmer lines (typically one or a
combination of
,
,
and
lines). In principle, as shown by Newberry et al. (1990), the most
robust approach to detect post-starburst galaxies would make full use
of all three Balmer lines. However, our limited spectral range does
not allow us to include
in numerous cases, and the S/N
ratio of the
line is generally poor. Therefore, we
use the combination of [OII] and of the
line to
establish our spectral classification, as in Dressler et al. (1999).
We apply two methods for measuring the equivalent widths of these
lines, using respectively the tasks "splot'' and "sbands'' in IRAF
on our continuum normalized spectra. The continuum is fitted by a
spline function of high order (generally 9), and the quality of the
best fit verified for each spectrum. In particular, special attention
is paid to the good quality of the fit in the vicinity of the line to
be measured, which is crucial for a reliable measurement of the EW. While "splot'' fits a
Gaussian for the line, "sbands'' measures
the decrement or increment in signal relative to the fitted continuum
in an interval containing the line to be measured. We check the
consistency of the measures with the two methods, and choose the
Gaussian fitting technique because it gives more accurate measurements
for the following analysis. These measurements are listed in
Table A.1, with equivalent widths of absorption and emission
features defined as positive and negative respectively. We estimate
the minimum measurable EW of each spectrum as in Barrena et al. (2002), taking the width of a line spanning 3.1 Å (our
dispersion) in wavelength, and with an intensity three times the noise
rms in the adjacent continuum. A minimum (maximum) measurable EW of
2.8(-2.8) Å results in the case of absorption (emission)
lines.
The fraction of poststarburst galaxies in galaxy clusters, (the so
called E+A objects; Dressler & Gunn 1983, also referred to more
recently as "k+a'' and "a+k'' types by Franx 1993; Dressler et al. 1999; or
-strong galaxies by Couch & Sharples
1987) has been a quite controversial topic in the last decade. While
the general property of these objects is the presence of a young
stellar population (A type) superimposed on an older one, a strict
consensus does not yet exist on the criteria used for their definition
as shown in Table 3 of Tran et al. (2003). This can lead to
significant differences when estimating and comparing the fraction of
these objects in various samples of clusters. The high frequency of
these objects at high redshift (up to
26% in the MORPHS survey,
Dressler et al. 1999) is in contrast with the very low fraction
obtained at low redshift (around
,
Dressler 1987). These objects
are considered to have undergone recent (younger than 1.5 Gyr) star
formation activity, followed by a quiescent phase. The situation at
intermediate redshifts is still debated, with low fractions obtained
by Balogh et al. (1999) for the CNOC1 survey (
for
), to higher values of
obtained
by Tran et al. (2003). However, these discrepancies could be explained
by differences in the selection criteria and in the correction for
incompleteness. This population of E+A galaxies was suggested to
consist of recently accreted field galaxies, with suppression of the
star formation activity by ram pressure stripping by the ICM. The
evolution of the fraction of such objects with redshift has often been
interpretated as a consequence of both a higher star formation density
and infall rate (Kauffmann 1995) at high redshift. Recent analyses
also show a luminosity/mass effect, with an observed decrease in the
characteristic E+A mass (see Poggianti et al. 2004;
Tran et al. 2003). We will present hereafter the analysis of the E+A galaxies
in A3921.
Table 6: Column 1: name of the scenario in the text - Col. 2: time since last interaction between the two clumps - Col. 3: angle between the plane of the sky and the line connecting the centres of the two clumps - Col. 4: relative velocity between the two clumps - Col. 5: spatial separation between the two systems - Col. 6: state of the systems for the possible solutions.
In a first step, we adopt the criteria defined by Dressler et al. (1999), dividing E+A galaxies in k+a's and a+k's, and defining k+a
and a+k objects from spectra presenting an intermediate and strong
line in absorption respectively, with no detectable
[OII] line in emission (with the criteria
Å,
and 3 Å
Å for k+a,
Å for a+k). We detect 6 k+a
objects in our sample, corresponding to a fraction of
of the
sample, and no a+k galaxies. As we are interested in all objects with
signs of recent star formation activity, we also consider spectra
presenting a clear inversion of the intensities of the K and H calcium
lines, due to the presence of a blend of the H line with the Balmer
line
.
We measure the ratio of the line intensities
,
shown by Rose et al. (1985) to have
a low value for young populations and to reach a plateau at
1.2
for older ones. We select a sample of objects with a ratio of
as candidates of k+a galaxies (coded
in our Table A.1 with k+a?). We detect 11 objects, with four
of them included in the previous sample (
Å), and the remaining ones spreading within the range (2.5 Å
Å). When taking these objects to be candidates
of galaxies having recently undergone star formation, we obtain a
fraction of 15.6% of our sample. A quantitative comparison of the
observed fraction of active galaxies in A3921 to those measured in
other clusters is however difficult as it can be affected by several
factors like the size of the field where spectroscopy was performed,
the spatial sampling of the galaxies or the way active galaxies were
defined. Nevertheless, the measured fraction of E+A galaxies seems to
be relatively large for a cluster at low redshift.
![]() |
Figure 17: The histogram in absolute magnitude in the RAB band of our high S/N spectroscopic sample (83 objects). The various fractions of galaxies of different spectral types in each bin are encoded as follows: passive (white), k+a (less shaded), k+a? (more shaded) and emission-line galaxies (black). The dashed and dotted lines show the magnitude limits corresponding to the MORPHS (Dressler et al. 1999) and "giants/dwarfs'' (Poggianti et al. 2004) cuts respectively. Note that the contribution of active (k+a, k+a? and e = emission line) galaxies creates a second peak at fainter magnitudes. |
Open with DEXTER |
In order to compare our results to those of Dressler et al. (1999), we
convert their spectroscopic limit
to our
band, taking into account the different cosmological models
used. This leads to
(using a mean V-R =0.5as in Poggianti et al. 1999,
and
). Applying this magnitude cutoff, we obtain respectively 1 secure k+a galaxy, and 3 k+a candidates.
As seen from the histogram displayed in Fig. 17, our secure
k+a and "k+a candidates'' (k+a?) are spread in the magnitude range
,
with the majority of
objects lying between
.
Taking
as a boundary in absolute magnitude, we consider
two sub-samples corresponding to faintest and brightest objects, and
analyse the respective fraction of different spectral types. The
results are listed in Table 7. While the k-spectrum
galaxies dominate the "bright'' sample, the fraction of secure k+a
and k+a candidates increases clearly in the "faint'' one. The
luminosity distribution of active galaxies in A3921 appears
intermediate between the MORPHS sample (Dressler et al. 1999) and the
Coma one (Poggianti et al. 2004). The fraction of bright k+a's
(
secure case,
when including candidates) is lower than
in MORPHS but higher than in Poggianti et al. (2004) who found no
object at such magnitudes. In the latter case however, the k+a
distribution is dominated by a population of dwarf galaxies (fainter
than
in our case) which is scarcely tested by
our shallower survey. We clearly detect an evolution of the typical
luminosity of our k+a population, which is brighter than in Coma but
fainter than that detected by MORPHS.
In contrast to Coma, Fig. 7 shows that most of the k+a
galaxies of our spectroscopic sample have red colours, as indicated by
their positions in the A3921 colour-magnitude diagram (circles in
Fig. 7). This can also be seen in Fig. 20, which
displays the
equivalent width as a function of the
colour as in Poggianti et al. (1999). Among the
13 detected k+a and k+a-candidate galaxies, 11 have
colours spanning the same range as the k galaxies (0.35-05), while
only 2 show clearly bluer colours. Moreover, only one k+a galaxy
shows
Å (Fig. 20), most objects
displaying relatively low
,
comparable to those
of the "redder k+a'' Coma subsample of Poggianti et al. (2004) and
significantly lower than the EWs of
lines detected
in their bluer k+a galaxies. We do not detect this young k+a blue
population in A3921, where galaxies with k+a spectra look rather like
the result of an episode of star formation or a star-burst occurring
ago, with redder colours and weaker Balmer
lines.
Table 7:
Number of spectra in different magnitude sub-samples. Bright
and faint galaxies correspond to objects with
and
respectively. Bright objects correspond to the same
sampling in absolute magnitude as in the MORPHS spectroscopic
catalogue (Dressler et al. 1999).
![]() |
Figure 18:
Projected galaxy density map (34![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 19: From left to right, velocity distribution (with a binning of 400 km s-1) of k, k+a and emission line galaxies (according to the definition of Dressler et al. 1999, see text for details). In the central figure, the white component corresponds to secure k+a's and the shaded part to k+a candidates. In the right panel, the shaded component of the histogram corresponds to e(a)'s, the black component to e(b)'s, the white part to e(c) type objects. |
Open with DEXTER |
In Fig. 18 we show the projected positions on the sky of secure k+a's (squares) and k+a candidates (circles). From a simple visual inspection of this figure, most of k+a galaxies seem to be concentrated in projection around the central field of A3921-A. This is particularly clear in the case of k+a candidates, as only one of them is located elsewhere, in the cluster outskirts. Confirmed k+a's are less clustered around the central field of A3921-A. They are however concentrated in the highest density regions of the whole cluster and roughly aligned along an East/West direction. Applying a bi-dimensional Kolmogorov-Smirnov test, we cannot exclude that the spatial distribution of k+a's is significantly different from that of k-type objects with a quite high significance level (87%). An intermediate and therefore inconclusive significance level (26%) is on the contrary obtained by comparing the projected positions of k+a candidates and passive galaxies.
The velocity distribution of the whole k+a sample is shown in the
middle panel of Fig. 19. According to a visual inspection and
the usual statistical tests (i.e. 1D-KS, Rank-Sum, Hoel 1971; and
F-tests, Press et al. 1992), it is not significantly different from
the velocity distribution of the k-types (left panel). k+a candidates
(shaded in Fig. 19) are characterized by a significantly
higher velocity dispersion than passive and confirmed k+a galaxies
(96% and 93% c.l. respectively according to an F-test), very close
to the whole cluster value (1188 km s-1 vs. 1008 km s-1, see
Tables 8 and 2). Mean velocities of the
different sub-samples are largely comparable within 1 errors.
The emission-line population also shows interesting features, both in
its fraction and colour distribution. We find 11 emission line
galaxies, among which there are: a) 2 starbursts (e(b)), with very
strong [OII] emission, b) 4 spirals forming stars at a constant rate
(e(c)), with weak Balmer absorption and weak to moderate [OII] emission, and c) 5 e(a)-type objects, with strong
absorption and measurable [OII] emission, suggested to be dusty
starburst galaxies (Poggianti et al. 1999). The global fraction of
emission line galaxies (13%) is lower than that estimated at high
redshift from MORPHS (26.5%), but comparable to that measured at
intermediate redshifts (
0.18<z<0.55) in the CNOC1 survey
(
16%, Ellingson et al. 2001). Fig. 21 shows the
as a function of
colour. The
sample can be divided into two groups: the bluest ones with
,
and typically large absolute values of
,
(<-10 Å), and the reddest ones with colours typical
of the passive population
,
and lower
absolute values of
,
Å), except for one
object with an exceptionally strong [OII] line (
Å), and a red colour index
,
which is
identified with BG3. The presence of this red population of galaxies
presenting emission lines is quite unusual, and reminiscent of what
was observed in the case of the merging cluster RX J0152.7-1357 at
z=0.837 analysed by Demarco et al. (2004). This suggests that an
on-going episode of star formation or bursting is occurring in an
older population.
![]() |
Figure 20:
Colour-
![]() |
Open with DEXTER |
Table 8: Properties of the cz distribution for the different spectral types. * In the case of the two separate k+a and k+a? subsamples we use the classical variance estimator for velocity dispersion and not the biweight indicator, owing to the low number of objects (Beers et al. 1990).
![]() |
Figure 21: Colour-[OII] diagram for galaxy members of A3921. Crosses indicate k type galaxies, squares k+a type galaxies, circles galaxies with H-K inversion and triangles emission line galaxies. |
Open with DEXTER |
![]() |
Figure 22:
Projected galaxy density map (
![]() ![]() |
Open with DEXTER |
Most (8 out of 11) star-forming galaxies seem to be spatially concentrated in the central region of the whole cluster, i.e. on the NW side of A3921-A, on the SE side and in the centre of A3921-B and between the two sub-clusters (see Fig. 22). A Rank-Sum test confirms this visual impression, as star-forming objects are found to be more concentrated than k-type galaxies around a central position between the two subclusters with a very high c.l. (99.7%). Moreover, according to a bidimensional KS test, the spatial distributions of e-type and k-type galaxies are not drawn from the same parent population (2% significance level).
Figure 19 shows that emission line galaxies have a much more
dispersed velocity distribution than do passive objects. According to
the F-test, the velocity dispersion of emission line galaxies (
1607
+314-219 km s-1) is significantly higher than the corresponding
value for k-type objects (
770 +119-96 km s-1), with a
100% confidence level.
Most star-forming objects are therefore distributed between the two sub-clusters, and their velocity distribution is much dispersed than that of the relaxed k-type galaxies (whose velocity distribution does not differ significantly from a Gaussian), but with a comparable mean velocity. This suggests that, at least for a fraction of them, the star formation episode could be the result of an interaction with the ICM, connected with the merging process of the cluster.
We derive here the star formation rate (SFR) of the galaxies within
A3921 from our measurements of the [OII] equivalent widths of our
galaxy spectra. The uncertainties in the SFR estimated from [OII] are
known to be larger than in that derived through
.
This
is due to the fact that, while flux estimates from
scale with the number of ionizing stars, estimates from [OII] are
hampered by systematics due to large variations in excitation,
implying substantial dispersions in the [OII]/
and in
the [OII]/
+[NII] correlations. In our case, however,
[OII] is the only tracer of the SFR available in our wavelength
domain, and is used with caution. Using the equation derived by
Kennicutt (1992):
![]() |
(11) |
the SFR is directly estimated from the
integrated broad-band B luminosity, the [OII] equivalent width,
assuming a reasonable value for the extinction correction
.
Converting our V-band magnitudes to the B-band using
for the objects belonging to the red sequence and
for bluer objects, and taking a solar luminosity
MB=5.48, as well as the canonical value for extinction (1 mag for
)
used in Kennicutt (1992), we obtain an estimate of the
SFR for our 11 galaxies with [OII] emission lines. From the majority
of spectra within A3921, the SFR ranges from 0.14 to
,
with however one object with a very high SFR (
). This is the galaxy BG3
previously quoted. We estimate a median SFR of
0.58 and a mean
SFR of
(BG3 excluded). The
cumulated SFR for the whole cluster reaches
(BG3 excluded). As expected, the values show that star
formation in A3921 is suppressed compared to field galaxies (Tresse
& Maddox 1998). The highest optical SFR in our data apart from the
BG3 (
), is intermediate between the
values obtained in A114 at z=0.32 by Couch et al. (2001) (
), and that obtained by Duc et al. (2002)
(
)
in Abell 1689 at z=0.18. However,
the SFRs in Couch et al. (2001) were measured from
,
and, as shown in Duc et al. (2002), measurements of the SFR with [OII]
can be seriously underestimated due to dust extinction. Therefore our
values can be considered as a lower limit of the SFR in A3921.
BG3 shows an exceptionally high star formation rate, suggesting a very
active star-bursting object. The spectrum shows strong absorption
lines typical of an old stellar population (Ca K and H, G band, 4000 Å break) as well as emission lines ([OII],
and
). Additional observations are planned to further
investigate the nature of this object and to test if the emission is
due to star formation or to an active galactic nucleus (only the bluer
part of the spectrum was obtained in these observations and the
equivalent width of [OIII] is necessary to perform the test).
Using our new spectroscopic and VRI-band photometric catalogues of
the central
region of A3921, we
detect two main sub-clusters of galaxies: a South-East clump (A3921-A)
that hosts the BCG (BG1), and, at
0.74 Mpc distance, a
North-West system (A3921-B), hosting the second brightest cluster
galaxy (BG2) and the third brightest object (BG3) in its
outskirts. From the analysis of the projected density distribution,
there are clear signs of merging events within this cluster. The
density distribution of galaxies of A3921-A shows an elliptical shape,
while A3921-B is characterized by a very morphologically disturbed
distribution, with an extended low-density tail extending towards
South. This morphology suggests that sub-cluster B is approaching
cluster A from the North or exiting after colliding with it from the
South. The fact that the sub-cluster B is well defined would suggest
that the merging has not yet occurred, but in case of a high impact
parameter, or of a low mass ratio between the secondary and the main
component, the density structure of the sub-clusters may survive. On
the other hand, the presence of the Southern tail of subcluster B
rather suggests a post-merger event. In this case, the tail would be
made of faint galaxies lagging behind the brightest objects during
their tangential crossing of subcluster A.
However, in spite of clear merging signatures in the density
distribution, the kinematical and dynamical properties both of the
whole cluster and of the two subclusters do not show strong signatures
of merging. Moreover, the two sub-clusters show very similar mean
projected velocities. When applying a two-body dynamical model to the
two components, several solutions could explain the observed dynamics
of A3921 allowing both the pre-merging and the post-merging cases. A
comparison with the X-rays properties of the cluster can allow us to
discriminate between them. The analysis of XMM-Newton observations by
Belsole et al. (2005) reveals that the X-ray emission of A3921-A can be
modelled with a 2D-model, leaving a distorted residual
structure toward the NW, coincident with A3921-B (see
Fig. 23). The main cluster detected in X-ray is centred on the
BCG position (BG1), while the X-ray peak of the NW clump is offset
from the brightest galaxy of A3921-B (BG2) (see Belsole
et al. 2005). The temperature map of the cluster shows an extended hot
region oriented parallel to the line joining the centres of the two
sub-clusters. A comparison of this image with numerical simulations by
Ricker & Sarazin (2001) suggests that we are observing the central
phases of an off-axis merger between two unequal mass objects, where
clump A is the most massive component (Belsole et al. 2005), consistent
with optical results.
![]() |
Figure 23:
X-ray residuals after subtraction of a 2D-![]() ![]() |
Open with DEXTER |
By combining the signatures of merging derived both from the optical
iso-density map and from X-ray results, we can now reconsider the
solutions of the two-body dynamical model (Table 6). In
the pre-merger case, the high-angle solutions of cases (a1) and (a2)
would imply a very large real separation for the two sub-clusters
(6-7 Mpc), which is very unlikely taking into account the clear
signs of interaction between the two clumps observed both in the
optical and in X-rays. In the "recent'' post-merger case, we can also
exclude the solution (c), as we expect a higher value of the relative
velocity (
1000 km s-1) between the two clumps for obtaining so
clearly a hot bar in the temperature map. Finally, the comparison of
observed and simulated galaxy density and temperature maps
(e.g. Schindler & Böhringer 1993; Ricker & Sarazin 2001) clearly
exclude an older merger (e.g. t0=1 Gyr), as we would not expect to
observe a clear bimodal morphology in the optical any longer and we
should not detect such obvious structure in the temperature map.
Therefore, only the solutions corresponding to the very central phases
of merging (
Gyr) can explain all our observational results, implying a
collision axis nearly perpendicular to the line of sight. This is
consistent with the absence of strong merging signatures in the
observed projected velocity distribution of cluster members.
The superposition of the X-ray residuals onto the optical iso-density map (Fig. 23) shows that the bulk of X-ray emission in A3921-B is offset towards SW from the main concentration of galaxies. As numerical simulations show that the non-collisional component is much less affected by the collision than the gas distribution (e.g. Roettiger et al. 1993), this offset suggests that A3921-B is tangentially traversing A3921-A in the SW/NE direction, with its galaxies in advance with respect to the gaseous component. The off-axis collision geometry has probably prevented total assimilation of the B group into the main cluster A. This off-axis collision scenario is also consistent with the shape of the feature in the temperature map (Belsole et al. 2005).
We detect very few k+a galaxies (1 secure and 3 candidates) at the
bright absolute magnitude cut of MORPHS (Dressler et al. 1999). Most
of our k+a and k+a? objects have
absolute magnitudes in
the range [-19/-20]. This effect is very similar to the trend in
luminosity detected by Poggianti et al. (2004) in Coma, but the
typical luminosity of our k+a/k+a? objects is at least one magnitude
brighter. However, we do not sample the faint end of the luminosity
function, which may contain part of the undetected population as in
the case of Coma. We therefore confirm the increase of the typical
luminosity of the k+a population with redshift as detected in
Poggianti et al. (2004), comparing Coma to the MORPHS clusters. We
also find a significant fraction of k+a in our low redshift cluster
lying between the fractions obtained at low and high redshift. The
mechanism responsible for this "redshift'' effect in A3921 is not
immediately clear. This can be a consequence only of redshift, due to
some cosmic "downsizing'' effect, as it has been shown that the star
formation activity at high z was more efficient for more massive
galaxies. In this case, this effect could be explained only by the
infall of the galaxies into the cluster. On the other hand, A3921
shows very strong sub-clustering and signatures of merging, which can
be suspected to be responsible, at least in some fraction, for the
enhancement of the fraction of such objects. This will be discussed in
detail in the following.
The typical EW(H)
of the k+a/k+a? galaxies detected in A3921
is moderate and most of the objects have red colours,
indistinguishable from red-sequence objects. We fail to detect the
population of blue k+a with strong Balmer lines as detected in Coma by
Poggianti et al. (2004), which can only be reproduced by a starburst
in the recent past. In contrast, our objects reflect the evolution of
galaxies having undergone starburst or starforming activity which has
been suppressed by some physical process, and now are in the second
half of their lifetime (typically <1.5 Gyr), with redder colours,
and fainter Balmer lines, before reaching a k-type spectrum. This
population can be reproduced by simply halting continuous star
formation, without evoking a strong star burst. The k+a? objects are
probably galaxies that have undergone a still older and fainter last
episode of star formation, as they do not have Balmer lines strong
enough to enter the k+a sample, but show clear signatures of past
activity.
A fraction of emission line galaxies comparable to that found in
intermediate redshift clusters has also been detected. The spatial and
kinematical distributions of k+a and k+a? as well as emission line
objects have been compared to the substructures found both in the
whole galaxy and in the gas distribution. While k+a/k+a? galaxies are
mostly distributed over the main cluster A, the emission line galaxies
mostly lie in the region of the sub-cluster B, and in the region in
between A and B. A similar spatial distribution is shown by blue
galaxies, which are more clustered in the central region of A3921-B
than in the whole field and, in particular, than in the center of the
more massive clump A (see Fig. 9). The comparison of the
observed distributions of blue and emission line objects with the
merging scenario presented previously (0.3 Gyr) suggests that
the interaction with the ICM during the passage of the sub-cluster B
on the edge of cluster A may have triggered star-bursting. This
hypothesis is supported by the significant difference in the radial
velocity dispersions of emission-line and passive galaxies, suggesting
that emission-line objects are a dynamically younger population than
the general cluster members. In contrast, the k+a/k+a? population
shows the signature of older star formation activity which can hardly
be related to the on-going merger, but may be understood either as
previously infalling galaxies, or relics from another older merging
process, but now at rest within the main cluster, as shown from the
velocity distribution.
Ongoing multi-wavelength (optical, X-ray and radio) follow-up will allow us to study in detail the environment and star formation properties of this cluster.
Acknowledgements
We thank the anonymous referee for useful suggestions that helped to improve this paper. We warmly thank all our collaborators working on the X-ray analysis of A3921: Elena Belsole, Jean-Luc Sauvageot, Gabriel W. Pratt and Hervé Bourdin. We thank B. Vandame for providing us with his imaging reduction package "alambic'' as well as his support in using it and G. Mars for his contribution to the reduction. We thank Albert Bijaoui, Andrea Biviano and Philippe Prugniel for fruitful discussions and helpful comments. We thank Christophe Adami for providing spectra of ENACS galaxies. This research was supported in part by Marie Curie individual fellowships MEIF-CT-2003-900773 (C.F.), by the PNC ("Programme National de Cosmologie'') grants in 2002 and 2003, and by the "Conseil Scientifique'' of the "Observatoire de la Côte d'Azur''. C.F. thanks François Mignard and the CERGA department of the "Observatoire de la Côte d'Azur'' for having financially supported part of this work. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology under contract with the National Aeronautics and Space Administration.
The technique adopted to isolate the red sequence of the elliptical galaxies of A3921 and calculate its parameters (slope, intercept and width) is described in this Appendix.
The Colour Magnitude Diagram (CMD) of the central field of A3921
(
)
is shown in Fig. 7. Due to the
high asymmetry of the galaxy distribution and the heavy contamination
of field objects in our CMD, a simple linear regression plus an
iterated 3
clip (Gladders et al. 1998) does not give
satisfactory results. We thus use the median absolute deviation (MAD)
as scale indicator of the object distribution in our CMD, as it is
extremely efficient (
90%) in the case of heavy-tailed
distributions (Beers et al. 1990). In fact, we expect that, in the
CMD, the deviations of the red-sequence galaxies from the optimal
linear fit follow a Gaussian distribution, while deviations of field
objects do not and populate the tails of the distributions instead.
We moreover apply an asymmetric clip (
)
on both sides of the distribution
(
). We then
follow a procedure similar to that described in Gladders et al.,
iterating our linear fitting and clipping until convergence on a final
solution was obtained. The red sequence can be described by the
linear equation
with a width
of
(see Fig. 7).