W. Boschin1 - M. Girardi1 - R. Barrena2,3 - A. Biviano4 - L. Feretti5 - M. Ramella4
1 - Dipartimento di Astronomia, Università degli Studi di
Trieste, via Tiepolo 11, 34131 Trieste, Italy
2 - Instituto de Astrofisica
de Canarias, 38200 La Laguna, Tenerife, Spain
3 - INAF - Telescopio
Nazionale Galileo, Roque de Los Muchachos, PO box 565, 38700 Santa
Cruz de La Palma, Spain
4 - INAF -
Osservatorio Astronomico di Trieste, via Tiepolo 11, 34131 Trieste,
Italy
5 - Istituto di Radioastronomia
del C.N.R., via Gobetti 101, 40129 Bologna, Italy
Received 19 September 2003 / Accepted 6 November 2003
Abstract
We present the results of the dynamical analysis of the
rich, hot, and X-ray very luminous galaxy cluster A2219, containing a
powerful diffuse radio-halo. Our analysis is based on new redshift
data for 27 galaxies in the cluster region, measured from spectra
obtained at the TNG, with the addition of other 105 galaxies recovered
from reduction of CFHT archive data in a cluster region of
radius (
0.8 h-1 Mpc
at the cluster distance) centered on
the cD galaxy. The investigation of the dynamical status is also
performed using X-ray data stored in the Chandra archive. Further,
valuable information comes from other bands - optical photometric,
infrared, and radio data - which are analyzed and/or discussed, too.
We find that A2219 appears as a peak in the velocity space at
z=0.225, and select 113 cluster members. We compute a high value for
the line-of-sight velocity dispersion,
km s-1, consistent with the high average X-ray
temperature of 10.3 keV. If dynamical equilibrium is assumed, the
virial theorem leads to
for the global
mass within the virial region. However, further investigation based
on both optical and X-ray data shows significant signs of a young
dynamical status. In fact, we find strong evidence for the elongation
of the cluster in the SE-NW direction coupled with a significant
velocity gradient, as well as for the presence of substructure both in
optical data and X-ray data. Moreover, we point out the presence of
several active galaxies. We discuss the results of our
multi-wavelength investigation suggesting a complex merging scenario
where the main, original structure is subject to an ongoing merger
with a few clumps aligned in a filament in the foreground oriented in
an oblique direction with respect to the line-of-sight. Our
conclusion supports the view of the connection between extended radio
emission and merging phenomena in galaxy clusters.
Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 2219 - galaxies: distances and redshifts - intergalactic medium - cosmology: observations
In the hierarchical scenario for large-scale-structure formation, mergers are an essential ingredient of galaxy cluster evolution (e.g. White 1997; Evrard & Gioia 2002).
One of the most recent aspects of the merging phenomenology is the
possible connection of cluster mergers with the presence of extended,
diffuse radio sources, halos and relics (Feretti 1999; Buote
2002; Giovannini & Feretti 2002). Radio-halos are
located at the cluster center, with typical sizes of 0.5 h-1 Mpc,
regular shape, steep radio spectra and no significant polarization;
relics are located in peripheral regions of the clusters, with
irregular shape and generally highly polarized spectra (e.g. Feretti
& Giovannini 1996; Giovannini & Feretti 2002). The
synchrotron radio emission of halos and relics demonstrates the
existence of large-scale cluster magnetic fields, of the order of
0.1-1
G, and of widespread relativistic particles of energy
density 10-14-10-13 erg cm-3. Although the
properties of halos and relics are not well understood, several
suggestions for the mechanism transferring energy to the relativistic
electrons have been made (e.g. Ensslin 2000).
Cluster mergers were suggested to provide the large amount of energy
necessary for electron re-acceleration and magnetic field
amplification (Feretti 1999; Feretti 2002; Sarazin
2002). However, the precise radio-halo formation scenario is
still debated (e.g. Ensslin & Röttgering 2002). Diffuse
radio sources are quite uncommon and only recently have we become able
to study these phenomena on the basis of sufficient statistics (
30 clusters up to
,
e.g. Giovannini et al. 1999;
see also Giovannini & Feretti 2002).
Presently, growing evidence of connection between diffuse emission and cluster merging is based on X-ray data (e.g. Böhringer & Schuecker 2002; Buote 2002). Studies based on a large number of clusters have found a significant relation between the radio and the X-ray surface brightness (Govoni et al. 2001) and connections between the presence of radio-halos/relics and irregular and bimodal X-ray surface brightness distribution (Schuecker et al. 2001). However, we are far from understanding the whole picture. Unfortunately, to date, optical data is lacking or poorly exploited. Sparse literature concerns some few individual clusters (e.g. Colless & Dunn 1996; Gómez et al. 2000; Barrena et al. 2002; Mercurio et al. 2003).
Optical data are a powerful way to investigate the presence and the dynamics of cluster mergers (e.g. Girardi & Biviano 2002), too. The spatial and kinematic analysis of member galaxies allow us to detect and measure the amount of substructure, to identify and analyze possible pre-merging clumps or merger remnants. This optical information is complementary to X-ray information since galaxies and ICM react on different time scales during a merger (see numerical simulations by Roettiger et al. 1997). Moreover, additional information comes from spectral types of member galaxies, since cluster mergers could stimulate star formation in cluster galaxies; thus the spectral signatures of past activity are useful to determine the relevant time-scales (e.g. Bekki 1999; Terlevich et al. 1999).
To investigate the connection between diffuse radio sources and
cluster mergers, we have performed spectroscopic observations at TNG
and a multi-wavelength analysis of the Abell cluster A2219 (Abell et al. 1989). This cluster shows the presence of a powerful
radio-halo (Giovannini et al. 1999; Bacchi et al. 2003) and is a rich, X-ray luminous, hot cluster at
moderate redshift
- Abell richness =3 (Abell et al. 1989);
(0.1-2.4 keV)
erg s-1 (Ebeling et al. 1996);
keV (Rizza et al. 1998). Possible
evidence of the young dynamical status of A2219 comes from the
ROSAT/HRI X-ray data analysis and strong/weak lensing features (Rizza
et al. 1998; Smail et al. 1995; Dahle et al. 2002).
This paper is organized as follows. We present the new redshift data of A2219 in Sect. 2 and the analysis of optical data in Sect. 3. Our analysis of recent X-ray Chandra archival data is shown in Sect. 4. Active galaxies are analyzed and discussed on the basis of multi-wavelength data using mid-IR and radio data in Sect. 5. Finally, in Sect. 6, we summarize and discuss our results presenting a tentative picture of the dynamical status of A2219.
Unless otherwise stated, we give errors at the 68% confidence level (hereafter c.l.)
Throughout the paper, we assume a flat cosmology with
and
.
For the sake of simplicity in rescaling
we adopt a Hubble constant of 100 km s
h-1 Mpc. In this context, 1
corresponds to
0.152 h-1 Mpc.
Multi-object spectroscopic observations of A2219 were carried out at
the TNG telescope in May 2003 (proposal AOT-7 ID:18). We used
DOLORES/MOS with the LR-B Grism 1, yielding a dispersion of 187 Å mm-1, and the Loral CCD of
pixels (pixel size of 15
m). This combination of grating and detector results in
dispersions of 2.8 Å pix-1. We have taken 1 MOS mask with 39 slits.
We acquired two exposures of 1800.0 s each. Wavelength calibration
was done using arc lamps before or after each exposure
(Helium-Argon).
Reduction of spectroscopic data was carried out with IRAF
. The signal-to-noise ratio per pixel was
computed from the ratio between the signal at 6000 Å and the rms noise within
4000-7000 Å. The resulting average
signal-to-noise ratio of our spectra is
.
Radial velocities are determined using the cross-correlation
technique (Tonry & Davis 1979) implemented in the RVSAO
package (developed at the Smithsonian Astrophysical Observatory
Telescope Data Center). Each spectrum is correlated against six
templates for a variety of galaxy spectral types: E, S0, Sa, Sb, Sc,
Ir (Kennicutt 1992). The template producing the highest value
of
,
i.e., the parameter given by RVSAO and related to the
signal-to-noise of the correlation peak, is chosen. Also, all
spectra and their best correlation functions are examined visually to
verify the redshift determination. In some ambiguous cases, generally
very late-type galaxies, we use the EMSAO package to confirm the
estimated redshift, and in two cases we take the EMSAO redshift as a
more reliable estimate of the redshift.
Out of the 39 spectra of objects, 9 turned out to be objects at
(included alignment stars for masks), while for 3 cases,
with
,
we cannot determine the redshift. Thus, we
obtain spectra for 27 galaxies.
We added to our data observations stored in the CFHT archive (proposal ID: 01AF37, P.I.: J.P. Kneib): 1 MOS mask of one 1800 s -exposure with 155 slits. We reduced the spectra with the same procedure adopted for TNG ones. We obtained spectra for 111 galaxies, 6 of which are in common with TNG galaxies. The average S/N of these CFHT spectra is 8.
To combine our TNG data with those from the CFHT archive, we have
to check for a possible zero-point correction. We use the
six galaxies with double redshift determination. This allows us to
obtain a more rigorous estimate for the redshift errors since the
nominal errors as given by the cross-correlation are known to be
smaller than the true errors (e.g. Malamuth et al. 1992;
Bardelli et al. 1994; Ellingson & Yee 1994; Quintana
et al. 2000). Thus, for the six galaxies with double redshift
determination, we fit the first determination vs. the second one by
using a straight line and considering errors in both coordinates
(e.g. Press et al. 1992). The fitted line agrees with the one
to one relation (
,
), but the small value of
-probability indicates a
poor fit, suggesting that the nominal cross-correlation errors are
underestimated. Only when nominal errors are multiplied by
factor can the observed scatter be explained. Therefore,
hereafter we assume that true errors are larger than nominal
cross-correlation errors by a factor of 2.3. Slightly smaller correction
factors are obtained for nearby clusters (e.g. Malamuth et al. 1992; Bardelli et al. 1994; Quintana et al. 2000).
Our spectroscopic catalogue consists of 132 galaxies sampling a cluster
region of 5
around the dominant galaxy (No. 65), which is a
cD galaxy (e.g. Allen et al. 1992).
We also determine the equivalent widths (EW hereafter) of the
emission line [OII] and the absorption line H
,
which are good
indicators of current and recent star formation, respectively.
Hereafter all the EW are given in angstroms and are positive if the
line is in absorption and negative if the line is in emission. We
estimate the minimum measurable EW as the width of a line spanning 2.8 Å (our dispersion) in wavelength, with an intensity three times the
rms noise in the adjacent continuum.
A variety of spectroscopic-classification schemes are presented in
the literature (see, e.g., Couch et al. 1994;
Dressler et al. 1999 and references therein). We use a conservative approach
leading to a sparse classification (
of the sample, see
Table 1). We define e-type galaxies as those showing strong
active star formation as indicated by the presence of an [OII] line
with an equivalent width of EW([OII])
Å (e.g.
Hammer et al. 1997; Postman et al. 1998). Out of
galaxies having S/N> 10 we define as a-type those with strong
Balmer absorption
Å (e.g. the
"post-starburst'' galaxies of Couch et al. 1994). We find
fourteen "active'' galaxies: four of e-type - measured EW([OII])
in the range
[-39.3,-14.3] Å; ten of a-type - measured
)
in the range [4.0; 12.4] Å. Out of non-active galaxies
with S/N> 10 those (37) with measured
Å are classified as passive galaxies (k-type) and the remaining nine
galaxies, intermediate between a and k are classified as
i-type. Non-active galaxies do not show significant [OII] emission
and only two a-galaxies show a very small [OII] emission (-3.5 Å).
As far as photometry is concerned, we obtained exposures on a
unvignetted field centered on
,
(J2000) through B and R Harris filters
using the CCD camera mounted on the 1m JKT telescope (at Roque de los
Muchachos Observatory, La Palma) in July 2003. We performed several
exposures for total integration times of 3000 and 7200 s in R and B band, respectively, with seeing conditions about 1.4
.
We
carried out the image reduction (bias, flat and fringing corrections)
using the IRAF package and photometry using the AUTOMAG SExtractor
procedure (Bertin & Arnouts 1996). This photometry was
calibrated using Landolt's standard fields at the JKT and TNG
telescopes and transformed into the Johnson filter system. We estimate
that our photometric sample is complete down to R= 19.7 and B=21.4within the observed field.
We assign R magnitudes to 131 out of the 132 galaxies of our
spectroscopic catalogue: the only exception is galaxy No. 52, which is
too close to a very bright star. Out of these 131 galaxies, 129 have
assigned magnitudes in B band, too. We measure redshift for galaxies
down to
mag, but a high level of completeness is
obtained only for bright galaxies (
66% completeness for
R<19.7 mag).
Table 1 lists the velocity catalogue
(see also Fig. 1): identification number of each galaxy,
ID (Col. 1); right ascension and declination,
and
(J2000, Cols. 2-7); R magnitudes (Col. 8); B-R colours (Col. 9);
heliocentric radial velocities,
(in km s-1,
Col. 10) with assumed errors,
,
i.e., the nominal
ones given by cross-correlation technique multiplied by 2.3
(Col. 11); the code for the spectral classification SC, where "I''
and "R'' indicate mid-IR and radio emitting galaxies, respectively
(see Sect. 5), (Col. 12); the telescope used to obtain the spectra
(Col. 13, T and C indicate TNG and CFHT, respectively).
![]() |
Figure 1: R-band image of A2219 (North at top and East to left). Galaxies with successful velocity measurements are labeled as in Table 1. Circles and boxes indicate cluster member and non-member galaxies, respectively. Labels in bold indicate active galaxies (see text). |
| Open with DEXTER | |
According to the analysis of the redshift distribution based on the
one-dimensional adaptive kernel technique (Pisani 1993, see
also Fadda et al. 1996 and Girardi et al. 1996) A2219 appears as a peak at
in redshift
space (Fig. 2). Out of 132 galaxies with redshift,
fourteen are foreground objects and three are background objects. In
particular, the foreground galaxies form two peaks, one of six
galaxies at
and one of eight galaxies at
.
To reject possible interloper galaxies in the main peak we apply the
procedure of the "shifting gapper'' of Fadda et al. (1996),
which combines velocity and position information. This procedure
rejects galaxies that are too far from the main body of galaxies in
rest-frame velocity (
)
within a fixed radial bin shifted along
the distance from the cluster center. According to the Fadda et al.
prescriptions, we use a gap of 1000 km s
and a bin of 0.4 h-1 Mpc, or
large enough to include 15 galaxies. For the cluster center, based
on the two-dimensional kernel analysis, we find that the densest peak
is centered at
RA
and
Dec
,
very close to the cD galaxy
position; thus hereafter we assume the cD position as the cluster
center RA
and
Dec
.
The shifting gapper
procedure rejects two galaxies as non-members. They are indicated by
crosses in Fig. 3, which shows the velocity histogram of the
113 selected cluster members. To check the reliability of the
member selection procedure, we also apply the method by den Hartog &
Katgert (1996). Of the initial 132 galaxies with redshift
we find the same 113 cluster members.
![]() |
Figure 2: Square root of the redshift galaxy density, as provided by the adaptive-kernel reconstruction method. Units on the y axis are arbitrary. |
| Open with DEXTER | |
![]() |
Figure 3:
Left panel: rest-frame velocities
|
| Open with DEXTER | |
The mean cluster velocity
km s
is computed using the biweight estimator (ROSTAT package; Beers et al. 1990). We estimate the line-of-sight velocity dispersion,
,
using the biweight estimator and applying the
cosmological correction and the usual correction for velocity errors
(Danese et al. 1980). We obtain
km s-1, where errors are obtained with the
bootstrap technique. Assuming that the system is in dynamical
equilibrium, the value of
leads to a value of the
radius of the collapsed, quasi-virialized region of
h-1 Mpc
(Carlberg et al. 1997) and a virial mass estimate of
(Girardi et al. 1998; Girardi & Mezzetti 2001).
The high value of
is confirmed if we apply a
tighter member selection (i.e., using a bin of 0.25 h-1 Mpc
in the
shifting gapper): only another three galaxies are rejected -those with
the lowest velocities- resulting in
km s
for the remaining galaxies. Moreover, we compute the
differential and integral mean velocity and velocity dispersion
profiles in Fig. 4: this analysis shows that the global
values of
and
are already reached
within the central cluster region of 0.1-0.2 h-1 Mpc. The conclusion of
the above analyses is that contamination by obvious field interlopers
and/or close cluster companions cannot explain the high value of the
global velocity dispersion. More probably, this value is connected
to the peculiarity of the internal dynamics of the cluster itself,
which will be analyzed below.
![]() |
Figure 4:
Differential (big squares) and integral (small points) mean velocity
and line-of-sight velocity dispersion profiles. For the differential
profiles, the results for four annuli from the cluster center, each of 0.2 h-1 Mpc, are shown. For the integral profiles, the mean and dispersion at
a given (projected) radius from the cluster center is estimated by
considering all galaxies within that radius. The error bands
at the |
| Open with DEXTER | |
As a first step in the study of internal dynamics we analyze the velocity distribution. We consider three estimators of Gaussianity: the kurtosis; the skewness, and the scaled tail index (see, e.g., Beers et al. 1991). We find no significant evidence that the velocity distribution differs from a Gaussian. Moreover, we do not find any evidence of a peculiar velocity of the cD galaxy with respect to the average velocity of cluster members (see Gebhardt & Beers 1991).
To investigate the velocity field we divide galaxies in low and
high-velocity samples using the median cluster velocity and check
the difference between the two distributions of galaxy positions.
Figure 5 shows that low and high-velocity galaxies are
segregated roughly along the E-W direction. The two distributions
are different at the
c.l. according to the two-dimensional
Kolmogorov-Smirnov test (hereafter 2DKS-test; see Fasano &
Franceschini 1987, as implemented by Press et al. 1992). To estimate the direction of the velocity
gradient we perform a multiple linear regression fit to the observed
velocities with respect to the galaxy positions in the plane of the
sky (see also den Hartog & Katgert 1996; Girardi et al. 1996): we find a
degrees (measured
counter-clock-wise from the North, see Fig. 5).
![]() |
Figure 5: Spatial distribution on the sky of 113 cluster members. Open and solid circles indicate low-and high-velocity galaxies: the larger the symbol, the smaller is the radial velocity. The plot is centered on the cluster center (coincident with the cD galaxy). The solid and faint lines indicate the position angle of the cluster velocity gradient and relative errors, respectively. |
| Open with DEXTER | |
Moreover, we compute the mean velocity and LOS velocity dispersion
separately for each of the four samples corresponding to the NE, SE,
SW and NW cluster sectors. Table 2 shows that mean
velocities are larger in the Western than in the Eastern sectors in
agreement with the above analyses. For
,
this is
largest in the SE sector and smallest in the NW sector: the difference
is significant at the
c.l. according to the F-test (Press et al. 1992).
Table 2: Mean galaxy velocities and velocity dispersion (in km s-1) in four cluster sectors (see text). N indicates the number of member galaxies in each sector.
To check for the presence of three-dimensional substructure we
compute the
-statistics devised by Dressler & Schectman
(1988) and establish its significance running 1000 Monte Carlo
simulations, in which we randomly shuffle galaxy velocities. The
signal of substructure is significant at the 95% c.l. In Fig. 6
we plot the distribution on the sky of all galaxies, each marked by a
circle: the larger the circle, the larger is the deviation
of the local parameters from the global cluster
parameters, i.e., there is more evidence for substructure. The most
likely substructure lies close to the cluster center roughly to the SE
of the cD galaxy.
![]() |
Figure 6: Spatial distribution on the sky of the 113 cluster members, each marked by a circle: the larger the circle, there is more evidence for substructure (according to the Dressler & Schectman test, see text). The plot is centered on the cluster center. |
| Open with DEXTER | |
We check for possible colour segregation of galaxies using the
110 member galaxies for which we have B-R colours. We
obtain a correlation between
and colour (at
the
c.l.). We divide the sample in two subsamples with
colours bluer and redder than the median B-Rcolour = 2.12. We obtain that the velocity dispersions of the two
subsamples differ at the
c.l. (
and
km s
for
the B-R< 2.12 and
B-R > 2.12 galaxies, respectively).
![]() |
Figure 7:
B-R vs. R diagram for galaxies with available spectroscopy:
small dots and crosses denote cluster and field members, respectively.
Large symbols denotes spectroscopic classified galaxies: circles,
triangles, squares, and stars indicate k-, i-, a-, and
e-galaxies, respectively. The solid line gives the best-fit CMR;
the dashed lines are drawn at |
| Open with DEXTER | |
To further investigate the above difference we also use the
colour-magnitude relation (hereafter CMR), which indicates the
early-type galaxy locus,
| (1) |
Table 3: Mean velocities and velocity dispersions (in km s-1) for galaxies with different colors (see text). N is the number of galaxies of each type.
We also analyze the two-dimensional spatial distribution of galaxies.
First, we consider bright cluster members, i.e. the sample of 80 galaxies with R<19.7 for which we have a good level of spatial
completeness (see Sect. 2). Figure 8 shows the presence of
two clumps in the cluster center which defines a SE-NW elongated
structure (
): the main peak roughly coincides with
the cD position, and the secondary peak - at the West - traces the
filament of high velocity galaxies well visible in Fig. 1
(from No. 49 to No. 28). The cD also shows an elongation in the SE-NW
direction described by the position angle on the celestial sphere
degrees (Smail et al. 1995).
Then, to study fainter galaxies and to work with much more robust
statistics, we consider the B and I photometric data of Smail et al. (1998) and define likely members using the CMR. Specifically, we select 484 likely cluster members considering
objects (I<22, stellar index <0.9) within 0.25 mag from the
.
The top panel of
Fig. 9 shows a main structure centered on the cD galaxy and
elongated in the SE-NW direction, with a PA of
for
the central region (along the filament described above) and a larger PA for the external region. A similar structure is also traced by the
subsample of very red galaxies - those above the CMR
(Fig. 9 middle panel). Instead, the distribution of not
very red galaxies shows three main clumps aligned in the NE-SW
direction (Fig. 9 bottom panel). In particular, the densest
peak is located at about
South-East of the cD, in the
position of the substructure suggested by the Dressler-Schectman test
(Sect. 3.2).
![]() |
Figure 8:
Spatial distribution on the sky and relative isodensity
contour map of 80 bright cluster (R<19.7 galaxies), obtained with
the adaptive-kernel method (cf. Pisani 1993, 1996).
The plots are centered on the cluster center. The 1
|
| Open with DEXTER | |
![]() |
Figure 9: Same as Fig. 8 but for likely cluster members selected on the basis of the B-I CMR in the photometric sample of Smail et al. (1998). Top panel: the 484 I<22galaxies. Middle and bottom panels: two subsamples corresponding to very red and not very red galaxies (273 and 211 galaxies, respectively). |
| Open with DEXTER | |
The X-ray analysis of A2219 is performed on the archival data of the
Chandra ACIS-S observation 800072 (exposure ID #896, P.I.:
J. Houck). The pointing has an exposure time of 42.8 ks. Data
reduction is performed using the package CIAO
(Chandra Interactive
Analysis of Observations) on chip S3 (field of view
). First, we remove events from the level 2 event list with a status not equal to zero and with grades one, five
and seven. Then, we select all events with energy between 0.3 and 10 keV. In addition, we clean bad offsets and examine the data, filtering
out bad columns and removing times when the count rate exceeds three
standard deviations from the mean count rate per 3.3 s interval. We
then clean S3 chip for flickering pixels, i.e., times where a pixel
has events in two sequential 3.3 s intervals. The resulting exposure
time for the reduced data is 42.3 ks.
In Fig. 10 we plot an R-band image of the cluster with
superimposed the X-ray contour levels of the Chandra image and the
radio contour levels of a low-resolution VLA image by Bacchi et al.
(2003). The elongated shape of the cluster is clearly
visible. Using the CIAO package Sherpa we fit an elliptical Beta
model to the X-ray photon distribution to quantify the departure from
the spherical shape. The model is defined as:
![]() |
(2) |
To detect possible substructures in A2219 we perform a
wavelet analysis by running the task CIAO/Wavdetect on the chip S3. The task was run on different scales to search for
substructures with different sizes. The significance
threshold
was set at 10-6. The
results are shown in Fig. 10. Thick ellipses represent
four significant substructures found by Wavdetect. The principal one,
located at RA
and Dec
,
is centered on the cD
galaxy. There are two more significant peaks in the core of the
cluster located at
RA
and
Dec
,
and at
RA
and
Dec
,
respectively.
Finally, outside the core, we find a fourth structure at
RA
and
Dec
,
about 2
South-East of
the dominant galaxy. The presence of a photon excess in this position
is also evident from Fig. 11, where we plot the residual
map of the elliptical Beta model fit. A peak of positive residuals is
observed at the position of the fourth structure.
For the spectral analysis of the cluster X-ray photons, we first
compute a global estimate of the ICM temperature. The temperature is
computed from the spectrum of the cluster within a circular aperture
of 3
radius around the cluster center. Fixing the absorbing
galactic hydrogen column density at
cm-2,
computed from the HI maps by Dickey & Lockman (1990), we fit
a Raymond-Smith (1977) spectrum using the CIAO package Sherpa
with a
statistics. We find a best fitting temperature of
10.3
keV, in agreement with previous
estimates by Allen (2000) and White (2000).
To detect the possible presence of temperature gradients in the
cluster we divide the cluster in eight sectors, as shown in
Fig. 12, and compute the temperature in each of them.
We find a temperature gradient from the cooler SE sectors to the
hotter NW ones. This pattern is also confirmed by a softness ratio
analysis of the Chandra image. We define the softness ratio as
,
where S is the number count of photons with
soft energies in the range 0.5-2 keV, while H is the number count
of photons with hard energies in the range 2-7 keV. Then we count
soft and hard photons in circles with 50 pixel aperture radius
sliding on a
pixel grid. Background photon counts are
subtracted using an ACIS background event file by
Markevitch. The background-subtracted values of the softness
ratios are plotted in Fig. 13. In this figure different
grey levels identify regions with different
s. High
(lower
temperature) regions are concentrated mainly in sectors 1, 2, 3 and 4,
in agreement with the temperature map shown in Fig. 12.
![]() |
Figure 10: R-band image of the cluster A2219 with, superimposed, the contour levels of the Chandra image (blue) and of a low-resolution VLA image (green) by Bacchi et al. (2003). The red ellipses represent the structures detected by Wavdetect. North is at top and East to left. |
| Open with DEXTER | |
The star formation properties of the galaxies in our sample are
examined on the basis of the spectral classification introduced in
Sect. 2 (see Table 1). Out of 60 classified galaxies, 46 are
member galaxies: 29/37 k-galaxies, 6/9 i-galaxies, 1/4
e-galaxy, and 10/10 a-galaxies. The relative fraction of
cluster a-galaxies is
.
For different EW(H
)
levels we
find
and
(>5 and >3 Å, respectively). For
comparison, Poggianti et al. (1999) found
of k+a/a+kgalaxies with EW(H
)>3 Å and Balogh et al. (1999)
found a fraction
of K+A with
EW(H
) > 5 Å. An open question is whether strong H
galaxies are typical of cluster environments (cf., e.g., Poggianti et al. 1999 with Balogh et al. 1999) since these galaxies
could be produced by a cluster merger (Bekki 1999). Very
interestingly, we find 0/14 a-galaxies in the field compared to 10/46 in cluster environment.
![]() |
Figure 11: Residual map of the elliptical Beta model fit (see text). Solid and dotted contour lines indicate regions with positive and negative residuals, respectively. Ellipses represent the structures detected by Wavdetect. |
| Open with DEXTER | |
![]() |
Figure 12: Temperature map of A2219. The cluster area is divided in eight numbered sectors ( top panel). Sectors are superimposed to the contour levels of the Chandra image. Ellipses show the four structures detected in the X-ray image. The graph in the bottom panel shows the temperature computed in each sector. |
| Open with DEXTER | |
![]() |
Figure 13: Softness ratio map of A2219 (see text). Superimposed to the map are the eight sectors (see Fig. 12). Dashed ellipses are the structures detected in the X-ray image and crosses show the positions of the cluster members. |
| Open with DEXTER | |
We check for possible spectral segregations of cluster galaxies both
in position and in velocity space by comparing the dynamics of passive
k-galaxies with those of active e- plus a-galaxies. We find a
larger velocity dispersion in the case of active galaxies, but the
difference is not significant. Moreover, we verify that the two
subsamples do not differ in the distribution of galaxy positions by
using the 2DKS-test. Figure 14 shows the distribution of
spectroscopically-classified galaxies: 4/25 and 7/21 are classified
as active in the central and external cluster regions (within and
outside 2
). Three central active galaxies are located
in the NW filament of high velocity galaxies in the central region,
corresponding to the NW subclump identified in the two-dimensional
analysis of Sect. 3.3. The fourth galaxy (No. 81) lies at the SE of
the cD galaxy.
For the interpretation of spectra, the only detected e-galaxy is
very blue (the bluest within this sample) as generally found for
starburst models. Nine of ten a-galaxies are classified as blue according to our definition (see Fig. 7). Thus their
spectra can be explained by "post-starburst'' models which reproduce
strong H
and no significant emission lines in a quiescent
phase soon after a starburst (about between a few Myr and 1.5 Gyr;
e.g. Poggianti et al. 1999). Alternatively, strong H
could be obtained during a period of 1-2 Gyr which follows the
end of star formation after a 1 Gyr period of constant star formation
(Morris et al. 1998).
![]() |
Figure 14: Spatial distribution on the sky of 113 cluster members. Large symbols denotes spectroscopic classified galaxies: circles, triangles, squares, and stars indicate k-, i-, a-, and e-galaxies, respectively. |
| Open with DEXTER | |
Galaxy No. 81 shows the strongest H
,
and it is therefore
spectroscopically classified as "a''. However, it also shows evidence
of mild [OII] emission,
Å. It would have been
classified "e(a)'' in the scheme of Poggianti et al. (1999). The galaxy coincides with one of the X-ray peaks,
identified at the smallest wavelet scale. The X-ray emission could
therefore be point-like, and consistent with emission from an AGN. The lack of strong emission lines could be explained if there is
substantial dust obscuration. This interpretation is consistent with
the fact that e(a) galaxies are generally associated with dusty
star-forming galaxies (Poggianti et al. 1999). The AGN
emission should then be re-emitted in the infrared (IR). Although
galaxy No. 81 was not detected at 15
m with ISOCAM by Barvainis
et al. (1999), one can see low-level isocontours of mid-IR
emission centered on this galaxy in Fig. 4 of Barvainis et al. (1999). Given that their rms noise level is
0.11 mJy, and their faintest detected source has a 15
m flux density
of 0.53 mJy, the 15
m flux density of galaxy 81 should be
0.4-0.5 mJy.
Of the 5 mid-IR sources detected by Barvainis et al. (1999) in the
field of A2219, three are in our spectroscopic sample. These are the
galaxies No. 42, 45 and 72 (respectively No. 1, 2, and 3 in Barvainis
et al. Table 1). Only galaxy No. 42 is a cluster member. The other
two are foreground galaxies. Galaxy No. 42 is an a-type galaxy, with
an H
Å. Its spectrum shows evidence of a mild [OII] emission,
Å. It would have been classified
"e(a)'' in the scheme of Poggianti et al. (1999). The flux in the [OII]
line can be used to compute the star formation rate of the galaxy,
using the relation of Kennicutt (1998). We find
yr-1.
We have examined the optical-mid-IR spectral energy distribution (SED) of galaxy No. 42 comparing the observed data (taken from
Barvainis et al. 1999) with the GRASIL models of Silva et al. (1998). The GRASIL model takes into account the effects of
several kinds of dust particles on the reprocessing and obscuration of
the stellar radiation (see:
http://web.pd.astro.it/granato/grasil/grasil.html). We find that the
observed SED of galaxy No. 42 is best-fitted by a model of a young (3 Gyr-old) early-type spiral, with significant ongoing star formation
(see Fig. 15). The best-fit model SED is used to compute
the galaxy total IR luminosity,
.
The galaxy star formation rate, SFR, is derived from the
total IR luminosity using the relation of Kennicutt (1998). We
obtain
yr-1. This
value is about 4 times higher than the estimate we obtain from the [OII] line flux. This is expected because of significant dust
extinction at short wavelengths. From the best-fit model SED we
estimate that dust extinction reduces the flux in the [OII] line by a
factor
3.
We conclude that galaxy No. 42 is a dusty star-forming galaxy, rather
than a post-starburst galaxy. We thus find that both e(a) galaxies of
our sample are dusty, active galaxies, a finding that supports the
interpretation given by Poggianti et al. (1999) for these kind of
galaxies.
![]() |
Figure 15: Optical-mid IR spectral energy distribution of galaxy No. 42. The best-fit GRASIL model is also shown (solid tree). |
| Open with DEXTER | |
With our spectral classification of a subset of A2219 galaxies, we can try to shed light on the nature of the "blue'' galaxies which are responsible for the Butcher-Oemler (BO, hereafter) effect (Butcher & Oemler 1978). We have translated the original Butcher & Oemler definition of "blue'' galaxies to the cluster redshift and the Band R filters, using the tables of Fukugita et al. (1995), and taking into account the Galactic absorption in the direction of A2219 (as given in the NASA Extragalactic Database). The dotted lines in Fig. 7 delimit the region of the colour-magnitude diagram where the BO galaxies are located. One can see that the BO population is of a mixed nature. First, there is a substantial number of field galaxies, of all spectroscopic types. Among BO cluster members, the dominant population is that of a-type galaxies. As we have shown, two of these a-type galaxies are in fact dusty active galaxies, rather than pure post-starburst objects. We thus conclude that the BO population of blue galaxies in the cluster A2219 is mostly composed of galaxies having experienced a recent starburst. In some of these galaxies the star-forming activity is not over yet, although partly hidden by dust.
For radio galaxies in the cluster, Bacchi et al. (2003)
confirm the presence of three radio galaxies at the cluster center of
A2219 (see Owen et al. 1992) identified in our catalogue with
No. 51, 77, 85 (i.e. S3, S4, and S6 respectively in Fig. 9 by Bacchi
et al. 2003). The northernmost galaxy No. 77, which lies in
the low-velocity tail of the cluster velocity distribution
(v=62364 km s-1), shows a tailed structure pointing towards the
NW, with a sharp bend towards the NE (see Owen & Ledlow
1997). The galaxy No. 85 is the second brightest galaxy of the
cluster (
17 R mag) and lies at <
SE separation from
the cD galaxy.
We analyze the internal dynamics of A2219 on the basis of
spectroscopic data for 132 galaxies in a cluster region of
(
0.8 h-1 Mpc
at the cluster distance) centered on the cD
galaxy. Out of 132 spectra, 27 galaxies come from new observations
carried out at the TNG. We find that A2219 appears as a peak in the
velocity space at z=0.225, and select 113 cluster members. The
investigation of the dynamical status is also performed using
X-ray data stored in the Chandra archive. Further valuable
information comes from other bands - optical photometric, infrared,
and radio data - which are analyzed and/or discussed, too.
The value we obtain for the line-of-sight velocity dispersion
km s
is in the high-tail of the distribution
of cluster velocity dispersions (see Fadda et al. 1996; Mazure
et al. 1996; Girardi & Mezzetti 2001). This global
value of
is consistent with the average value of
keV coming from the X-ray analysis
(
,
see also
Fig. 4) and with the value of
erg s-1 (Ebeling et al. 1996), converted to
erg s-1 (see
-
relation by,
e.g., Wu et al. 1999; Girardi & Mezzetti 2001).
Therefore, on the basis of global properties only, we would conclude
that A2219 is not far from the dynamical equilibrium and would trust
the large virial mass estimate M
2.2 h-1 Mpc
.
However, the high values of
,
,
and
could be due to the expected
increase during a cluster-merging phase (e.g. Schindler &
Böhringer 1993; Schindler & Müller 1993). Their
agreement could be a coincidence, and our virial mass estimate an
overestimate of the true value.
We find observational evidence that the cluster is not dynamically relaxed. Our optical and X-ray analyses show that A2219 is elongated in the SE-NW direction. This elongation is seen in the spatial distribution of the colour-selected likely cluster members, the shape of the cD galaxy, the X-ray contour levels of the Chandra image, the gradient in the velocity dispersion and in the X-ray temperature. In particular the ellipticity of X-ray isophotes is significantly higher than 0 and corresponds to an axis-ratio = 0.66, which is much lower than the usual value for galaxy clusters (median axis-ratio with 99% c.l. errors = 0.82+0.04-0.02 from Mohr et al. 1995). Note that Dahle et al. (2002) found that the mass distribution of this cluster recovered from a weak lensing analysis appears even more elongated than the light and galaxy number density distribution.
Very interestingly, the optical data show a somewhat smaller PA than
X-ray data (
and
,
respectively). A
different direction in the elongation of light and hot gas
distribution could be the result of a cluster merger in an advanced
phase as shown by numerical simulations where collisional and
collisionless components react in different way (e.g. Ricker &
Sarazin 2001; Schindler 2002). In fact, the
evidence of an ongoing merger in the cluster core was first suggested
by Smail et al. (1995) to explain the difference they find in
the orientation of the mass and gas distributions (
and 128
degrees), as obtained from gravitational lensing analysis and
ROSAT/HRI X-ray data, respectively.
For the presence of the velocity gradient roughly in the E-W
direction, it suggests that the structures are moving into the plane
of the sky at
(e.g. Quintana et al. 1996;
Roettiger & Flores 2000). Intermediate angles are the most
suitable to allow detection of the cluster elongation and the velocity
gradient during a merging phase.
The hypothesis of an advanced-phase merging is also supported by the lack of a cool center in the cluster (see also Allen & Fabian 1998), as expected since strong cluster mergers disrupt the cool cores (Peres et al. 1998; Roettiger et al. 1996), and the existence in the core of several star-forming or post-starburst galaxies, whose activity could be produced by the cluster merger (e.g. Bekki 1999). In particular, three of these active galaxies are part of the filament NW of the cD, also visible as a cool component in the X-ray softness ratio map (see Figs. 8 and 13), suggesting that this structure might be already involved in a merger event. Moreover, the shape of the radio tail of the radio galaxy No. 77 indicates that the velocity of this galaxy was originally pointing towards the cluster center (cD galaxy), but had suddenly changed direction, probably because of the cluster merger.
The structure we find in the cluster core could be related to the
cluster merger, too. Considering optical data, the Dressler-Schectman
test shows the possible existence of a substructure located at SE,
very close to the cD galaxy. For X-ray data, our wavelet
multi-scale analysis reveals three significant emission peaks located
within 1
from the cluster center. Of these, two peaks are
possibly connected to intrinsic X-ray emission from the cD and the
active galaxy No. 81, but no optical counterpart is found for the
third, intermediate peak. Our wavelet X-ray analysis also detects an
external diffuse substructure located 2
SE of the cD galaxy.
We suggest that this structure is a premerging clump. It coincides
with the cold SE X-ray region and thus could be well explained by the
two-dimensional superposition of a cold clump in the foreground of
the main cluster. The presence of this foreground clump could also
explain why, on average, the velocity to the East is lower than the
West since at least for some galaxies the low velocity is of
cosmological rather than of kinematic origin. These anomalous low
velocity galaxies could also produce a spurious enhancement of the
observed velocity dispersion in the SE region (see Table 2). A
population of very close foreground galaxies could also explain why
very red galaxies strangely differ from not-very red
galaxies both in kinematic and spatial properties (see
Sect. 3.3). While the former are really associated with the main
structure of A2219, the latter could be mixed with the clump members,
i.e. a population of foreground - and thus less red - galaxies.
Unfortunately, we have very poor redshift information in the SE region
to confirm our hypothesis and identify clump members.
On a much larger scale, the cluster is not so well isolated in the
redshift space, and the next possible merging clump might be the one at
(see Fig. 2).
Summarizing the above evidence, we suggest a complex merging scenario for A2219 involving many clumps, possibly in different dynamical states. We suggest that the main, original structure (hosting the cD galaxy) is subject to an ongoing merger with a few clumps aligned in a foreground filament obliquely oriented with respect to the line-of-sight; the filament projected on the sky lies in the SE-NW direction. In this scenario, the merging with one or more clumps is already in a very advanced phase, well after the first core passage, while one or more clumps are still in a pre-merging phase.
In the context of the above scenario, we argue that the presence of a
radio-halo in A2219 is related to the existence of an ongoing, very
advanced phase of merging between the main cluster and one or more
groups. In fact, the elongation of radio map contours (see
Fig. 10) is in the same direction as the X-ray surface
brightness contours. Moreover, the time necessary to re-accelerate
the electrons producing the radio-halo as a consequence of the
cluster merger (
1 Gyr) is comparable to that needed to obtain
post-starburst spectral signatures in galaxies of which A2219 is rich
(see Sect. 5).
In conclusion, our scenario supports the view of the connection between extended radio emission and merging phenomena in galaxy clusters.
Acknowledgements
We thank Luca di Fabrizio for B and R calibration observations at the TNG telescope and the anonymous referee for useful suggestions and comments. We particularly thank Ian Smail for providing us accurate Palomar 5-m telescope photometry. Work 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''), by the Italian Space Agency (ASI), and by INAF (Istituto Nazionale di Astrofisica) through grant D4/03/IS. This publication makes use of data accessed as Guest User at the Canadian Astronomy Data Center, which is operated by the Dominion Astrophysical Observatory for the National Research Council of Canada's Herzberg Institute of Astrophysics (http://cadcwww.dao.nrc.ca/cfht/cfht.html), as well as of data obtained from the Chandra data archive at the NASA Chandra X-ray center (http://asc.harvard.edu/cda/).