EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 521 (October 2010) A&A, 521 (2010) A78 Full HTML
Free Access
Issue
A&A
Volume 521, October 2010
Article Number A78
Number of page(s) 12
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201014839
Published online 22 October 2010
A&A 521, A78 (2010)

The complex structure of Abell 2345: a galaxy cluster with non-symmetric radio relics[*]

W. Boschin1 - R. Barrena2,3,4 - M. Girardi4,5

1 - Fundación Galileo Galilei - INAF, Rambla José Ana Fernández Perez 7, 38712 Breña Baja (La Palma), Canary Islands, Spain
2 - Instituto de Astrofísica de Canarias, C/vía Láctea s/n, 38205 La Laguna (Tenerife), Canary Islands, Spain
3 - Departamento de Astrofísica, Universidad de La Laguna, Av. del Astrofísico Francisco Sánchez s/n, 38205 La Laguna (Tenerife), Canary Islands, Spain
4 - Dipartimento di Fisica dell'Università degli Studi di Trieste, Sezione di Astronomia, via Tiepolo 11, 34143 Trieste, Italy
5 - INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy

Received 21 April 2010 / Accepted 4 June 2010

Abstract
Context. The connection of cluster mergers with the presence of extended, diffuse radio sources in galaxy clusters is still debated.
Aims. We aim to obtain new insights into the internal dynamics of the cluster Abell 2345. This cluster exhibits two non-symmetric radio relics well studied through recent, deep radio data.
Methods. Our analysis is based on redshift data for 125 galaxies acquired at the Telescopio Nazionale Galileo and on new photometric data acquired at the Isaac Newton Telescope. We also use ROSAT/HRI archival X-ray data. We combine galaxy velocities and positions to select 98 cluster galaxies and analyze the internal dynamics of the cluster.
Results. We estimate a mean redshift $\left<z\right>=0.1789$ and a line-of-sight (LOS) velocity dispersion $\sigma_{V}\sim 1070$ km s-1. The two-dimensional galaxy distribution reveals the presence of three significant peaks within a region of $\sim$1  h70-1 Mpc$\;$(the E, NW, and SW peaks). The spectroscopic catalog confirms the presence of these three clumps. The SW and NW clumps have similar mean velocities, while the E clump has a larger mean velocity ( $\Delta {V}_{\rm rf}\sim 800$ km s-1); this structure causes the presence of the two peaks we find in the cluster velocity distribution. The difficulty in separating the galaxy clumps leads to a very uncertain mass estimate $M\sim 2$  $\times 10^{15}\;h_{70}^{-1}\;M_{\odot}$. Moreover, the E clump well coincides with the main mass peak as recovered from the weak gravitational lensing analysis and is off-set to the east from the BCG by $\sim$1.3$^\prime$. The ROSAT X-ray data also show a very complex structure, mainly elongated in the E-W direction, with two (likely three) peaks in the surface brightness distribution, which, however, are off-set from the position of the peaks in the galaxy density. The observed phenomenology agrees with the hypothesis that we are looking at a complex cluster merger occurring along two directions: a major merger along the $\sim$E-W direction (having a component along the LOS) and a minor merger in the western cluster regions along the $\sim$N-S direction, roughly parallel to the plane of the sky. The eastern radio relic is elongated in the direction perpendicular to that of the major merger, while the peculiar, western radio relic is elongated in the direction perpendicular to the bisecting of the two merger directions.
Conclusions. Our scenario for the internal dynamics of Abell 2345 strongly supports the use of the ``outgoing merger shocks'' model to explain the two radio relics, suggesting a consistent justification for their asymmetry.

Key words: galaxies: clusters: individual: Abell 2345 - galaxies: clusters: general - galaxies: distances and redshifts

1 Introduction

Diffuse radio emission in galaxy systems is a rare phenomenon presently known in only a few tens of (mostly) rich clusters (e.g. Giovannini et al. 1999; see also Giovannini & Feretti 2002; Feretti 2005; Giovannini et al. 2009). These radio sources (named halos and relics) are very interesting for their possible connection with cluster mergers (see Feretti et al. 2002, for a review). They show a typical synchrotron spectrum, a clear sign of the existence of large-scale cluster magnetic fields and of widespread relativistic particles. Cluster mergers have been suggested to provide the large amount of energy necessary for electron reacceleration up to relativistic energies and for magnetic field amplification (Feretti 1999, 2002; Sarazin 2002). Radio relics (``radio gischts'' as referred by Kempner et al. 2003), which are polarized and elongated radio sources located in the cluster peripheral regions, seem to be directly associated with merger shocks (e.g., Ensslin et al. 1998; Roettiger et al. 1999; Ensslin & Gopal-Krishna 2001; Hoeft et al. 2004). Radio halos, unpolarized sources which permeate the cluster volume similarly to the X-ray emitting gas, are more likely to be associated with the turbulence following a cluster merger (Cassano & Brunetti 2005). However, only recently the number of diffuse radio sources discovered in clusters has grown enough to allow their study on the basis of a sufficient statistics and to attempt a classification (e.g. Kempner et al. 2003; Ferrari et al. 2008).

There is growing evidence of the connection between diffuse radio emission and cluster mergers based on X-ray data (see Buote 2002; Feretti 2006, 2008 and references therein). 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 kinematical 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 really complementary to X-ray information since galaxies and intra-cluster medium react on different time scales during a merger (see, e.g., numerical simulations by Roettiger et al. 1997). In this context we are conducting an intensive observational and data analysis program to study the internal dynamics of clusters showing extended radio emission by using member galaxies (DARC, Dynamical Analysis of Radio Clusters; Girardi et al. 2007; http://adlibitum.oat.ts.astro.it/girardi/darc).

Among the clusters hosting diffuse radio emission, only a very small fraction is characterized by double radio relics: Abell 1240 (Giovannini et al. 1999; Bonafede et al. 2009, hereafter B09); Abell 2345 (Giovannini et al. 1999; B09); Abell 3667 (Röttgering et al. 1997); Abell 3376 (Bagchi et al. 2006); RXCJ 1314.4-2515 (Feretti et al. 2005; Venturi et al. 2007). According to van Weeren et al. (2009) cluster ZwCl 2341.1+0000 could host two radio relics, too. However, recent data by Giovannini et al. (2010) suggest a more complex picture, in which the presumed relics (separated in the high-resolution image of van Weeren et al.) are connected to form a large-scale filament of diffuse radio emission. The relics of Abell 1240 and 3667 were explained by the ``outgoing merger shocks'' model (Roettiger et al. 1999; B09). Observations of Abell 3376 agree with both the ``outgoing merger shocks'' and the ``accretion shock'' models (Bagchi et al. 2006).

In the case of Abell 2345 (hereafter A2345), the observations are more difficult to be reconciled with the theoretical scenarios. In fact, the two relics, separated by $\sim$2  h70-1 Mpc, are not perfectly symmetric with respect to the cluster center. Moreover, while the eastern relic (A2345-2) is well explained in the scenario of the ``outgoing merger shocks'' model, the western relic (A2345-1) is quite peculiar (B09). In particular, the spectral index of A2345-1 steepens in the southern and external parts while diffusive shock acceleration models predict a steepening of the radio spectrum towards the cluster center as a consequence of the electron energy losses after shock acceleration (e.g., Ensslin et al. 1998). B09 suggested that the radio properties of A2345-1 could be affected by a possible ongoing merger with another group lying in the external cluster regions, at the NW with respect to the center of A2345 and visible in ROSAT/PSPC X-ray data (the ``X1'' group).

During our observational program we have conducted an intensive study of the cluster A2345 performing spectroscopic and photometric observations with the Telescopio Nazionale Galileo (TNG) and the Isaac Newton Telescope (INT), respectively.

A2345 is a very rich, X-ray luminous, Abell cluster: Abell richness class = 2 (Abell et al. 1989); $L_{\rm X}$(0.1-2.4 keV) =  $4.28\times 10^{44} ~ h_{70}^{-2}$ erg s-1 recovered from ROSAT data (Böhringer et al. 2004). Optically, the cluster is classified as Bautz-Morgan class III (Abell et al. 1989), while it is classified as cD cluster according to Rood-Sastry scheme (Struble & Rood 1987).

Weak gravitational lensing analysis has been performed by Dahle et al. (2002; hereafter D02) and by Cypriano et al. (2004) for the region close to the brightest cluster galaxy (hereafter BCG). D02 found that the highest peak in the mass map is off-set to the east from the central BCG by $\sim$1.5 \ensuremath{{\rm ^\prime}} and noticed that ROSAT/HRI archival data show a large amount of substructure.

Few galaxies having redshift are reported in the field of A2345 (see NED) and the value usually quoted in the literature for the cluster redshift comes from the spectrum of 2124-120, a radio galaxy at z=0.176 (Owen et al. 1995; Owen & Ledlow 1997).

This paper is organized as follows. We present our new optical data and the cluster catalog in Sect. 2. We present our results about the cluster structure in Sect. 3. We discuss our results and give our conclusions in Sect. 4.

Unless otherwise stated, we give errors at the 68% confidence level (hereafter c.l.). Throughout this paper, we use H0=70 km s-1 Mpc-1 in a flat cosmology with $\Omega_0=0.3$ and $\Omega_{\Lambda}=0.7$. In the adopted cosmology, 1 \ensuremath{{\rm ^\prime}} corresponds to $\sim$181  h70-1 kpc $\;$at the cluster redshift.

2 New data and galaxy catalog

Multi-object spectroscopic observations of A2345 were carried out at the TNG telescope in August 2008. We used DOLORES/MOS with the LR-B Grism 1, yielding a dispersion of 187 Å/mm. We used the new E2V CCD, a matrix of $2048\times2048$ pixels with a pixel size of 13.5 $\mu$m. In total we observed four MOS masks for a total of 147 slits. Total exposure times were of 1 h for three masks and 1.5 h for the remaining mask. Wavelength calibration was performed using mercury-neon lamps. Reduction of spectroscopic data was carried out with the IRAF[*] package.

Radial velocities were 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 was 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 $\cal R$, i.e., the parameter given by RVSAO and related to the signal-to-noise ratio of the correlation peak, was chosen. Moreover, all spectra and their best correlation functions were examined visually to verify the redshift determination. In two cases (IDs. 34 and 108; see Table 1) we took the EMSAO redshift as a reliable estimate of the redshift. Our spectroscopic survey in the field of A2345 consists of 132 spectra with a median nominal error on cz of 36 km s-1. The nominal errors as given by the cross-correlation are known to be smaller than the true errors (e.g., Malumuth et al. 1992; Bardelli et al. 1994; Ellingson & Yee 1994; Quintana et al. 2000; Boschin et al. 2004). Double redshift determinations for seven galaxies allowed us to estimate real intrinsic errors in data of comparable quality taken with the same instrument (e.g. Barrena et al. 2007a,b). 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, when using the nominal cross-correlation errors, the small value of $\chi^2$-probability indicates a poor fit, suggesting the errors are underestimated. Only when nominal errors are multiplied by a $\sim$1.4 factor the observed scatter can be explained. Therefore, hereafter we assume that true errors are larger than nominal cross-correlation errors by a factor 1.4. For the seven galaxies with two redshift measurements we used the mean of the two redshift determinations and the corresponding errors.

Finally, our spectroscopic survey in the field of A2345 consists of spectra for 125 galaxies with a median error on cz of 50 km s-1.

As far as photometry is concerned, our observations were carried out with the Wide Field Camera (WFC), mounted at the prime focus of the 2.5 m INT telescope. We observed A2345 with $g^{\prime }$, $r^{\prime }$and $i^{\prime }$ Sloan-Gunn filters on December 2008 in photometric conditions and seeing of $\sim$1.6 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}}( $g^{\prime }$and  $r^{\prime }$images) and $\sim$2.5 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}} ( $i^{\prime }$ image).

The WFC consists of a four-CCD mosaic covering a 33 \ensuremath{{\rm ^\prime}} $\times $ 33 \ensuremath{{\rm ^\prime}} field of view, with only a 20% marginally vignetted area. We took 9 exposures of 720 s in $g^{\prime }$ band, 13 of 360 s in $r^{\prime }$ band and 10 of 600 s in $i^{\prime }$band. So, the total exposure times were 6480 s, 4680 s and 6000 s in  $g^{\prime }$, $r^{\prime }$ and $i^{\prime }$ band, respectively. By developing a dithering pattern we were able to build a ``supersky'' frame that was used to correct our images for fringing patterns (Gullixson 1992). In addition, the dithering helped us to clean cosmic rays and avoid gaps between the CCDs in the final images. The complete reduction process (including flat fielding, bias subtraction and bad-column elimination) yielded a final coadded image where the variation of the sky was lower than 0.8% in the whole frame. Another effect associated with the wide field frames is the distortion of the field. In order to match the photometry of several filters, a good astrometric solution is needed to take into account these distortions. Using the imcoords IRAF tasks and taking as a reference the USNO B1.0 catalogue, we were able to find an accurate astrometric solution (rms $\sim$0.4 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}}) across the full frame. The photometric calibration was performed observing the SA92 Landolt standard field, that is also calibrated in the Sloan-Gunn system (Smith et al. 2002).

We finally identified galaxies in our $g^{\prime }$, $r^{\prime }$ and  $i^{\prime }$ images and measured their magnitudes with the SExtractor package (Bertin & Arnouts 1996) and AUTOMAG procedure. In a few cases (e.g. close companion galaxies, galaxies close to defects of the CCD) the standard SExtractor photometric procedure failed. In these cases we computed magnitudes by hand. This method consists in assuming a galaxy profile of a typical elliptical galaxy and scaling it to the maximum observed value. The integration of this profile gives us an estimate of the magnitude. This method is similar to PSF photometry, but assumes a E-type galaxy profile, more appropriate in this case.

\begin{figure} \par\includegraphics[width=16cm,clip]{14839f01.eps} \end{figure} Figure 1:

INT $r^{\prime }$-band image of the cluster A2345 (North at the top and East to the left) with, superimposed, the contour levels of the (smoothed) ROSAT/HRI archival image (pointing ID: US800769H, thick contours) and the contour levels of a VLA radio image at 1.4 GHz (see B09; beam 50 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}} $\times $ 38 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}}, thin contours starting at 0.24 mJy/beam and then spaced by a factor 2). Arrows highlight the contour levels of the two radio relics.

Open with DEXTER

As a final step, we estimated and corrected the galactic extinction from Burstein & Heiles's (1982) reddening maps. The values of the extinction coefficients are $A_{g^{\prime}} \sim0.15$, $A_{r^{\prime}} \sim0.12$and $A_{i^{\prime}} \sim0.09$ for the $g^{\prime }$, $r^{\prime }$ and $i^{\prime }$ Sloan-Gunn filters, respectively.

We estimated that our photometric sample is complete down to $g^{\prime}=23.3$ (24.0), $r^{\prime}=22.7$ (23.3) and $i^{\prime}=21.4$ (22.5) for S/N=5 (3) within the observed field.

We measured redshifts for galaxies down to magnitude $r^{\prime}\sim 21$, but a high level of completeness is reached only for galaxies with magnitude $r^{\prime }<20$ (>60% completeness).

Table 1 lists the velocity catalog (see also Fig. 2): identification number of each galaxy, ID (Col. 1); right ascension and declination, $\alpha$ and $\delta$(J2000, Col. 2); $g^{\prime }$, $r^{\prime }$ and $i^{\prime }$magnitudes (Cols. 3-5); heliocentric radial velocities, ${v}=cz_{\hbox{$\odot$ }}$ (Col. 6) with errors, $\Delta {v}$ (Col. 7); member assignment (Col. 8; 1:A2345, 2:background/foreground).

2.1 Radio galaxies in the field of A2345

Our spectroscopic catalog lists some radio galaxies observed in the sky region of A2345. Contour levels in Fig. 1 show a pointlike radio source $\sim$40 \ensuremath{\ensuremath{{\rm ^\prime}}\hskip -0.1em\ensuremath{{\rm ^\prime}}\;} W of the BCG. Indeed, a high resolution radio image by B09 (see their Fig. 1) reveals that this source is a blend of two head-tail galaxies (ID. 70 and ID. 76 in our catalog). At $\sim$5 \ensuremath{{\rm ^\prime}} SW of the BCG, Fig. 1 shows another pointlike source. B09 image reveals that this is the radio source 2124-120 cited above, another head-tail galaxy. Its optical counterpart in our catalog is galaxy ID. 35. Finally, $\sim$5 \ensuremath{{\rm ^\prime}} NE of the BCG there is a third pointlike source (NVSS J212725-120623). Contrary to galaxies IDs. 70, 76 and 35, which are obvious cluster members (see Sect. 3.1 for details on member selection), this radio galaxy has no counterpart in our optical image and is a likely background object.

3 Analysis and results

3.1 Member selection and the question of the cluster center

To select cluster members out of 125 galaxies having redshifts, we follow a two steps procedure. First, we perform the 1D adaptive-kernel method (hereafter DEDICA, Pisani 1993 and 1996; see also Fadda et al. 1996; Girardi et al. 1996). We search for significant peaks in the velocity distribution at >99% c.l. This procedure detects A2345 as two close peaks at $z\sim0.1775$ and 0.1806 populated by 55 and 46 galaxies, but largely overlapped since many galaxies (66/101) have a non-negligible probability of belonging to both peaks (see Fig. 3). Out of 24 non members, 12 and 12 galaxies are foreground and background galaxies, respectively.

\begin{figure} \par\includegraphics[height=18.5cm,clip]{14839f02.eps} \vspace*{3mm} \end{figure} Figure 2:

INT $r^{\prime }$-band image of the cluster A2345 (West at the top and North to the left). Circles and boxes indicate member and non-member galaxies, respectively (see Table 1). Labels indicate the IDs of cluster galaxies cited in the text.

Open with DEXTER

All the galaxies assigned to the cluster are analyzed in the second step which uses the combination of position and velocity information: the ``shifting gapper'' method by Fadda et al. (1996). This procedure rejects galaxies that are too far in velocity from the main body of galaxies and within a fixed bin that shifts along the distance from the cluster center. The procedure is iterated until the number of cluster members converges to a stable value. Following Fadda et al. (1996) we use a gap of 1000 km s$^{-1} \;$- in the cluster rest-frame - and a bin of 0.6  h70-1 Mpc, or large enough to include 15 galaxies.

The choice of the center of A2345 is not obvious. No evident dominant galaxy is present in our catalog, e.g. the BCG (ID. 82) is only 0.4 mag brighter in the $r^{\prime }$ band than the second brightest galaxy (ID. 2). However, ID. 2 is located in the very north-western cluster region. The eastern central region is someway dominated by the BCG, since the other luminous galaxy (ID. 70) is $\sim$1 mag fainter. Probably, that is why this cluster was classified as a cD cluster in the past. Other two luminous galaxies, ID. 35 (the radio galaxy 2124-120 discussed above) and ID. 33, lie in the south-western region. From the X-ray point of view, the center is not obvious, too. Böhringer et al. (2004) list a center very close ($\sim$40 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}}) to the position of the BCG, but Ebeling et al. (1996) list a $\sim$4 \ensuremath{{\rm ^\prime}} NE center. Indeed, the X-ray map is quite substructured as shown by the the analysis of ROSAT data (B09; D02; see also our following analysis). In view of these difficulties we decide to adopt from the operative and graphical point of view the location of the BCG, which is well defined, as the cluster center [RA =  $21^{{\rm h}}27^{{\rm m}}13~\rlap{\hbox{${\rm ^s}$ }}{\hbox{$.$ }}~70$, Dec =  $-12\ensuremath{{\rm ^\circ}}09\ensuremath{{\rm ^\prime}} 46.1\ensuremath{\ensuremath{{\rm ^\prime}}\hskip -0.1em\ensuremath{{\rm ^\prime}}\;} $ (J2000.0)]. After the ``shifting gapper'' procedure we obtain a sample of 98 fiducial cluster members (see Fig. 4, top panel).

\begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f03.eps}} \end{figure} Figure 3:

Redshift galaxy distribution. The solid line histogram refers to the 101 galaxies assigned to A2345. The redshift galaxy density, as provided by the adaptive kernel reconstruction method DEDICA, is showed in the insect panel, with an arbitrary normalization.

Open with DEXTER
\begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f04.eps}} \end{figure} Figure 4:

Top panel: rest-frame velocity vs. projected clustercentric distance for the 101 galaxies in the two peaks forming the cluster (see Fig. 3). Crosses indicate galaxies rejected as interlopers by our ``shifting gapper'' procedure. Large squares indicate the five brightest galaxies. Middle and bottom panels: integral profiles of mean velocity and LOS velocity dispersion, respectively. The mean and dispersion at a given (projected) radius from the cluster-center are estimated by considering all galaxies within that radius - the first values computed on the five galaxies closest to the center. The error bands at the $68\%$ c.l. are also shown.

Open with DEXTER
\begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f05.eps}} \end{figure} Figure 5:

The 98 galaxies assigned to the cluster. Upper panel: velocity distribution for all galaxies and for galaxies with $r^{\prime }<20$ (solid and dashed histograms, respectively). The big arrows indicate the position of the two peaks found by the 1D DEDICA algorithm. The Gaussians corresponding to the three galaxy clumps detected in the $r^{\prime }<20$ sample through the 3D DEDICA algorithm are shown (NW, SW and E clump from left to right, respectively). Notice that the separation of these three groups is much more evident in the 2D galaxy distribution (see Fig. 9 in the following). The small arrows show the velocities of the five brightest galaxies indicating the correspondent parent clump. Lower panel: stripe density plot where the arrows indicate the positions of the significant gaps.

Open with DEXTER

3.2 Global cluster properties

By applying the biweight estimator to the 98 cluster members (Beers et al. 1990, ROSTAT software), we compute a mean cluster redshift of $\left<z\right>=0.1789\pm 0.0004$, i.e. $\left< {v}\right>=\left<cz\right>=(53645\pm$108) km s-1. We estimate the LOS velocity dispersion, $\sigma_{V}$, by using the biweight estimator and applying the cosmological correction and the standard correction for velocity errors (Danese et al. 1980). We obtain $\sigma_{V}=1069_{-62}^{+75}$ km s-1, where errors are estimated through 1000 bootstrap simulations performed with the ROSTAT routine by Beers et al. (1990), which uses bias-corrected percentile intervals.

To evaluate the robustness of the $\sigma_{V}$ estimate we analyze the velocity dispersion profile (Fig. 4, bottom panel). The integral profile rises out to $\sim$0.3  h70-1 Mpc$\;$and then flattens suggesting that a robust value of $\sigma_{V}$ is asymptotically reached in the external cluster regions, as found for most nearby clusters (e.g., Fadda et al. 1996; Girardi et al. 1996).

3.3 Velocity distribution

We analyze the velocity distribution to look for possible deviations from Gaussianity that might provide important signatures of complex dynamics. For the following tests the null hypothesis is that the velocity distribution is a single Gaussian.

We estimate three shape estimators, i.e. the kurtosis, the skewness, and the scaled tail index (see, e.g., Bird & Beers 1993). We find no evidence that the velocity distribution departs from Gaussianity.

Then we investigate the presence of gaps in the velocity distribution. We follow the weighted gap analysis presented by Beers et al. (1991, 1992; ROSTAT software). We look for normalized gaps larger than 2.25 since in random draws of a Gaussian distribution they arise at most in about $3\%$ of the cases, independent of the sample size (Wainer & Schacht 1978). We detect two significant gaps (at the $99.4\%$ c.l.) which divide the cluster in three groups of 48, 5 and 45 galaxies from low to high velocities (see Fig. 5, lower panel).

We also use here the Kaye's mixture model (KMM) test to find a possible group partition of the velocity distribution (as implemented by Ashman et al. 1994). The KMM algorithm fits a user-specified number of Gaussian distributions to a dataset and assesses the improvement of that fit over a single Gaussian. In addition, it provides the maximum-likelihood estimate of the unknown n-mode Gaussians and an assignment of objects into groups. This algorithm is usually used to analyze the velocity distribution where theoretical and/or empirical arguments indicate that the Gaussian model is reasonable. Here the one-dimensional KMM test fails in confirming both the two groups suggested by the DEDICA analysis and the three groups detected by the weighted gap method.

3.4 2D cluster structure

To analyze the two dimensional galaxy distribution we first use our photometric data sample which covers a larger spatial region.

In our photometric catalog we select likely members on the base of both ( $r^{\prime }$- $i^{\prime }$) and ( $g^{\prime }$- $r^{\prime }$) colours. Goto et al. (2002) showed that there is a small tilt in the color-magnitude relations ( $r^{\prime }$- $i^{\prime }$) vs. $r^{\prime }$ and ( $g^{\prime }$- $r^{\prime }$) vs. $r^{\prime }$ and that the scatter in the latter relation is roughly the double than the scatter in the first one. Out of our photometric catalog we consider as likely cluster members those objects with a SExtractor stellar index $\le $0.9 having ( $r^{\prime }$- $i^{\prime }$) and ( $g^{\prime }$- $r^{\prime }$) lying within $\pm $0.15 and $\pm $0.3 from the median values of $r^{\prime }$- $i^{\prime}=0.47$ and $g^{\prime }$- $r^{\prime}=1.18$ colours of the spectroscopically cluster members (see Fig. 6).

\begin{figure} \par\includegraphics[width=8cm,clip]{14839f06.eps} \end{figure} Figure 6:

Upper panel: $r^{\prime }$- $i^{\prime }$ vs. $r^{\prime }$ diagram for galaxies with available spectroscopy. Circles and crosses indicate cluster members and non members. The solid line gives the median colour determined on member galaxies; the dashed lines are drawn at $\pm $0.15 mag from this value. Lower panel: $g^{\prime }$- $r^{\prime }$ vs. $r^{\prime }$ diagram for galaxies with available spectroscopy. The dashed lines are drawn at $\pm $0.3 mag from the median value determined on member galaxies.

Open with DEXTER

Table 2:   Substructure from the photometric sample detected with DEDICA.

\begin{figure} \par\includegraphics[width=7.8cm,clip]{14839f07.eps} \vspace*{2mm} \end{figure} Figure 7:

Spatial distribution on the sky and relative isodensity contour map of likely cluster members with $r^{\prime }\le 21$, obtained with the 2D DEDICA method (external regions are not shown for the sake of clarity). The BCG is taken as the cluster center. The location of the five brightest galaxies are indicated by squares. Crosses indicate the peaks found in the ROSAT/HRI image. The circles indicate the position of the main (the largest circle) and the secondary peaks in the projected mass distribution as obtained from gravitational lensing analysis (see D02, notice that they analyzed only a small region around the BCG). The two relics are indicated in a schematic way by the ``+'' symbols. Solid lines highlight the directions of the two mergers. The dashed line is the bisecting of the merger directions (see text).

Open with DEXTER

Figure 7 shows the contour map for the 579 likely cluster members having $r^{\prime }\le 21$. The 2D DEDICA procedure finds for significant peaks. In order of relative density, they are the E peak, the NW peak, the SW peak and the N peak. For each of these subclumps, Table 2 lists the number of assigned members, $N_{\rm S}$; the peak position; the density (relative to the densest peak), $\rho_{\rm S}$; the $\chi^2_{\rm S}$ value. Ramella et al. (2007) tested the 2D DEDICA procedure with Monte-Carlo simulations reproducing galaxy clusters. They show that the physical significance associated to the subclusters is based on the statistical significance of the subcluster (recovered from the $\chi^2$ value) and the $r_{\rm CS}=N_{\rm C}/N_{\rm S}$ parameter, where $N_{\rm S}$is the number of members of the substructure and $N_{\rm C}+N_{\rm S}$is the total number of cluster members. The peaks listed in Table 2 have $r_{\rm CS}$ in the range of 4-6 and, as expected for real comparable substructures, $\chi^2\ga 25$ (see Fig. 2 of Ramella et al. 2007). We can estimate that these four subclumps have a small probability to be false detections (<$5\%$). Similar results are found for samples based on more conservative magnitude limits (e.g. $r^{\prime}\le 20.5$ and $r^{\prime}\le 20$).

To furtherly probe the robustness of these detections, we also apply the Voronoi Tessellation and Percolation (VTP) technique (e.g. Ramella et al. 2001; Barrena et al. 2005). This technique is non-parametric and does not smooth the data. As a consequence, it identifies galaxy structures irrespective of their shapes. For our purposes we run VTP on the same sample of 579 likely members as above. The result of the application of VTP is shown in Fig. 8. VTP is run three times adopting three detection thresholds: galaxies identified as belonging to structures at 95%, 98% and 99% c.ls. are shown as open squares, asterisks and solid circles respectively. VTP confirms the existence of the E, NW and SW clumps found with DEDICA, but doesn't find any significant structure in the position of the N clump. Since the existence of this clump is not confirmed we ignore it in the following. Instead, $\sim$9 \ensuremath{{\rm ^\prime}} south of the BCG, VTP detects a small galaxy peak, possibly a poor group infalling in the potential well of A2345.

\begin{figure} \par\includegraphics[width=8.5cm,clip]{14839f08.eps} \end{figure} Figure 8:

Galaxies belonging to structures as detected by the Voronoi Tessellation and Percolation technique. The algorithm is run on the sample of likely members with $r^{\prime }<21$ (see text). Open squares, asterisks and solid circles indicate galaxies in structures at the 95%, 98% and 99% c. ls., respectively. The E, NW and SW subclumps and the S group are indicated. As a comparison, big crosses indicate the positions of the subclumps found with DEDICA.

Open with DEXTER

We also use a ROSAT/HRI archival image (pointing US800769H) to study the 2D structure of the cluster. The image has an exposure time of $\sim$12.4 ks. As already noticed by D02, these data show a very large amount of substructure. X-ray contours in Fig. 1 show as the cluster structure is elongated in the E-W direction and the likely presence of three subclumps. The wavelet multiscale analysis performed on this image confirms the significance of the two most eastern subclumps, centered at RA =  $21^{{\rm h}}27^{{\rm m}}13~\rlap{\hbox{${\rm ^s}$ }}{\hbox{$.$ }}~99$, Dec =  $-12\ensuremath{{\rm ^\circ}}09\ensuremath{{\rm ^\prime}} 53.3\ensuremath{\ensuremath{{\rm ^\prime}}\hskip -0.1em\ensuremath{{\rm ^\prime}}\;} $ (J2000.0; close to the BCG position) and RA =  $21^{{\rm h}}27^{{\rm m}}05~\rlap{\hbox{${\rm ^s}$ }}{\hbox{$.$ }}~98$, Dec =  $-12\ensuremath{{\rm ^\circ}} 10\ensuremath{{\rm ^\prime}} 44.3\ensuremath{\ensuremath{{\rm ^\prime}}\hskip -0.1em\ensuremath{{\rm ^\prime}}} $ (J2000.0), respectively. The centers of these two wavelet structures, as well as a rough center indicator for the most western X-ray peak, are compared to the position of galaxy subclumps in Fig. 7.

3.5 3D-analysis

The existence of correlations between positions and velocities of cluster galaxies is a footprint of real substructures. Here we use several different approaches to analyze the structure of A2345 combining position and velocity information of the spectroscopic sample.

We analyze the presence of a velocity gradient performing a multiple linear regression fit to the observed velocities with respect to the galaxy positions in the plane of the sky and perform 1000 Monte Carlo simulations to assess the significance of this velocity gradient (e.g., Boschin et al. 2004 and refs. therein). We find no significant velocity gradient.

We also combine galaxy velocity and position information to compute the $\Delta$-statistics devised by Dressler & Schectman (1988; see also Boschin et al. 2006, for a recent application) which is very useful to detect small, compact substructures. We find no significant indication of substructure.

The spectroscopic sample covers only the regions of the E, SW, and NW peaks detected in the galaxy density distribution (see the above section). To study these regions we select galaxies within 0.3  h70-1 Mpc$\;$from each peak, where the limiting radius is chosen to avoid any member superposition, and obtain three subclumps containing 18, 8 and 9 galaxies respectively. According to the standard means-test (e.g., Press et al. 1992) the E and NW clumps have different mean velocities at the $93\%$ c.l., with the E clump having a larger velocity than the NW clump ( $\left<{v_{\rm E}}\right>=(54573\pm 232$) km s-1; $\left<{v_{\rm NW}}\right>=(53487\pm 407$) km s-1). We obtain a similar result comparing the E clump with the combined sample of SW+NW clumps ( $\left<{v_{\rm NW+SW}}\right>=(53~713\pm300$) km s-1).

The existence of a velocity difference between subclumps detected in two-dimensions prompted us to use the 3D KMM analysis to properly assign galaxies to the three subclumps, although the use of a Gaussian for the galaxy distribution is very approximate. We apply the 3D KMM algorithm both on the whole cluster sample and the subsample of 88 galaxies having $r^{\prime }<20$, i.e. with a relevant level of magnitude completeness. In both cases the algorithm fits a three-group partition, at the $99.4\%$ c.l. (at the $98.4\%$ c.l.) according to the likelihood ratio test, leading to three groups of 52, 17 and 29 galaxies (45, 13 and 30 galaxies), corresponding to the E, SW, and NW clumps. The results for the three KMM groups coming from both the analyses are shown in Table 3 (see also Figs. 5 and  9 for the $r^{\prime }<20$ sample).

Table 3:   Kinematical properties of cluster and subclumps.

\begin{figure} \par\resizebox{8.5cm}{!}{\includegraphics{14839f09.eps}} \end{figure} Figure 9:

Spatial distribution on the sky of the cluster galaxies (with $r^{\prime }<20$) showing the three groups recovered by the 3D KMM analysis. Solid circles, open circles and squares indicate the galaxies of E, SW and NW group, respectively. The BCG is taken as the cluster center. Large squares indicate the five brightest cluster members.

Open with DEXTER

4 Discussion and conclusions

Our estimate of the cluster redshift is $\left<z\right>=0.1789\pm 0.0004$ (cf. with z=0.176 reported in the literature as estimated by only one galaxy).

For the first time the internal dynamics of A2345 is analyzed on the basis of member galaxies.

The global value of the velocity dispersion $\sigma_{V}=1069_{-62}^{+75}$ km s$^{-1} \;$corresponds to that of a very rich cluster (e.g. the Coma cluster, Colless et al. 1996; Girardi et al. 1998) although it is not particularly high within the DARC sample. The position of A2345 in the $L_{\rm X, bol}-\sigma_{V}$ plane is well consistent with that of other clusters (see Fig. 5 of Ortiz-Gil et al. 2004, taking into account the different cosmologies).

4.1 Cluster structure

Our analysis of photometric data shows the presence of three significant peaks in the galaxy distribution, i.e. the E, SW, and NW peaks in order of importance. The E-peak is off-set to the east from the BCG by $\sim$1.3 \ensuremath{{\rm ^\prime}}($\ga$0.2  h70-1 Mpc).

The presence of E, SW, and NW galaxy groups is also confirmed including our spectroscopic data in the analysis. We use different techniques to assign galaxies to the corresponding three subgroups. The corresponding kinematical properties are shown in Table 3. The E group is found to have the largest mean velocity $\langle{v_{\rm E}}\rangle\ga 54~000$ km s-1, while the SW and NW groups have similar mean velocities with a rest frame difference of $\Delta {V}_{\rm rf}\la 800$ km s-1 with respect to the E group. The combination of these groups produces the double-peak structure found by our initial analysis of the spectroscopic sample (cf. the three Gaussians with the position of the two 1D DEDICA peaks in Fig. 5).

The complexity of A2345 structure and the small velocity differences among subgroups are the likely cause of the failure of global, standard 3D techniques when searching for substructures (e.g. the velocity gradient and the Dressler-Schectman statistics). Our results suggest that the SW and NW clumps define a direction parallel to the plane of sky, while the relative direction between the E clump and the SW+NW complex have a component along the LOS.

4.2 Cluster mass

The complicated structure of A2345 and the small velocity differences among subgroups also make us difficult to estimate the velocity dispersions for the individual subclumps. Table 3 lists what we nominally obtain in our analysis but, while the values of  $\sigma_{V}$ computed in the 1D DEDICA analysis are obvious underestimates of the true values due to the abrupt division of the velocity distribution, the values computed using other techniques are likely overestimates due to the possible contamination from the other, close groups. The 3D KMM method suggests that the group with the largest $\sigma_{V}$, i.e. the most massive one, is the E group. This also agrees with the fact that the E group coincides with a strong mass peak (see D02) and hosts the BCG and the main X-ray peak. However, the precise value of $\sigma_{V,\rm E}$ is not obvious: Fig. 4 shows as the galaxy group around the BCG has a $\sigma_{V}\sim 800$ km s$^{-1} \;$increasing out to $\sigma_{V}\ga 1100$ km s$^{-1} \;$when considering more and more distant (and so likely more contaminated) regions. This range of $\sigma_{V,\rm E}$ values overlaps with that of the $\sigma_{\rm WL}$ values obtained by the weak gravitational lensing analysis for a singular isothermal sphere model or other approaches (870-965 km s-1, D02; Cypriano et al. 2004). Hereafter, we will assume that $\sigma_{ V,\rm E}\sim 900$ km s$^{-1} \;$for the E group.

The uncertainties are even larger for the other two groups. The relative density in the 2D galaxy distribution, the amount of the population obtained through the 3D DEDICA method and the mass map (from the gravitational lensing, see the section below) suggest that the NW group is more massive than the SW group, but the values of  $\sigma_{V}$ are not conclusive. We decide to assume for both the groups $\sigma_{V,\rm SW}=\sigma_{V,\rm NW}\sim 700$ km s-1, a value somewhat smaller than that for the E group.

Finally, we attempt to estimate the mass of the whole A2345 complex. Following the methodology for the virial mass estimate already outlined in previous papers for other clusters of DARC sample (e.g., Girardi et al. 1998; Girardi & Mezzetti 2001; Boschin et al. 2009), the assumption of the dynamical equilibrium for a cluster having $\sigma_{V}=1069$ km s-1 would lead to a radius of the quasi-virialized region ${R}_{\rm vir}=0.17\times \sigma_{V}/H(z) = 2.4$  h70-1 Mpc$\;$and a mass of $M({<}{R}\sim2.4 \;h_{70}^{-1}~{\rm Mpc})\sim1.8$  $\times 10^{15}\;h_{70}^{-1}\;M_{\odot}$. Considering the combination of the three above subclumps, each assumed to be in equilibrium, we instead obtain $M\sim 1.1+0.5+0.5\sim 2.1$  $\times 10^{15}\;h_{70}^{-1}\;M_{\odot}$.

4.3 Merging phase

To investigate the phase of the cluster merger among the detected galaxy clumps we must resort to the comparison with results from other wavelengths. The complex cluster structure of A2345 is also revealed using X-ray and gravitational lensing analysis. The analysis of ROSAT/HRI data shows the presence of two, likely three, peaks in the X-ray image, with two X-ray peaks lying between E and SW peaks (of which one coincident with the BCG) and one X-ray peak lying between SW and NW peaks. No peak is shown by the HRI image in the northern region reinforcing the global E-W (or better ENE-WSW) elongation of the X-ray surface brightness. Two more peaks are shown by ROSAT/PSPC archival data in the most external NW region, in particular the peak ``X1'' (see B09) is close to the position of the second brightest galaxy in our spectroscopic catalog (ID. 2).

As for the weak gravitational lensing approach, D02 analyze a small cluster region of 6$^\prime$ $\times $ 6$^\prime$, i.e. in practice the region around our E group. The mass map shows a main peak off-set to the east from the BCG by $\sim$1.5$^\prime$, but well coincident with the position of the E peak as found in the 2D galaxy distribution. Unfortunately, both SW and NW peaks lie outside the region analyzed by D02. However, very interestingly, the mass map shows a strong elongation toward the NW region in the direction of the NW peak.

The displacement between X-ray and optical (or mass) peaks suggests that A2345 is in an ``after cross core'' phase of merger. The observational phenomenology indicate two merging directions. The X-ray brightness distribution is mainly elongated towards the E-W (slightly ENE-WSW) direction, the same direction defined by the E peak and the other two galaxy peaks (SW+NW peaks). This direction is roughly perpendicular to the elongation of the (eastern) A2345-2 radio relic. The second direction is a S-N (slightly SSE-NNW) direction as indicated by the SW and NW peaks, and the intermediate X-ray peak. The bisecting line of this and the above merging direction is roughly perpendicular to the direction of the elongation of the (western) A2345-1 radio relic (see Fig. 7). The same bisecting direction is also roughly indicated by the tail of the radio galaxy ID. 35 (Owen & Ledlow 1997, B09).

In the framework of the ``outgoing merger shocks'' model, the E-W merger is the natural origin of the formation of the eastern radio relic (A2345-2), while the formation of the expected western radio relic was disturbed by the presence of the additional S-N merger in the western cluster regions resulting in the formation of the north-western A2345-1 radio relic. This scenario would explain both the asymmetry of the two relics and, qualitatively, the peculiarity of A2345-1. In fact, the morphology of A2345-1 is different from the usual arc-like morphology with the presence of a kind of tail departing from the South of the relic and bended towards the external regions (see left-upper panel of Fig. 2 of B09). Moreover, the southern and external regions of A2345-1 show a higher spectral index than the northern region suggesting that the northern region has received an additional input of energy.

A difficulty in our above scenario is the absence of a southern relic, which we would expect due to the merger of the SW group with the NW group. However, this could be connected with the fact that the southern group is likely the less important galaxy clump in A2345.

4.4 Toy model for the complex merger

In the framework of the observational picture and the above scenario we attempt to quantify the internal cluster dynamics of A2345 through a simple analytical approach.

In the case of a two-groups merger, the ``outgoing merger shocks'' model is successfully supported by the kinematical data on subclumps through the use of the simple analytical bimodal model (e.g. Barrena et al. 2009, for Abell 1240).

As for A2345, we assume a toy model, obviously very approximate, where the cluster is first interested by the merger between the E group and a W group (formed by the SW+NW groups). Then, the western part of the cluster is interested by a second merger between the SW group and the NW group. In practice, according to this toy model, the internal dynamics of A2345 is analyzed by the combination of two successive bimodal models. For both these mergers we apply the two-body model (Beers et al. 1992; Thompson 1982) following the methodology outlined for A1240 (Barrena et al. 2009).

\begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f10.eps}} \vspace*{-4mm} \end{figure} Figure 10:

System mass vs. projection angle for bound and unbound solutions (thick solid and thick dashed curves, respectively) of the two-body model applied to the E and SN+SW subsystems. Labels BI$_{\rm a}$ and BI$_{\rm b}$ indicate the bound and incoming, i.e., collapsing solutions (thick solid curve). Label BO indicates the bound outgoing, i.e., expanding solutions (thick solid curve). Label UO indicates the unbound outgoing solutions (thick dashed curve). The horizontal lines give the range of observational values of the mass system with a 50% error. The thin dashed curve separates bound and unbound regions according to the Newtonian criterion (above and below the thin dashed curve, respectively).

Open with DEXTER

The two-body model assumes radial orbits for the clumps with no shear or net rotation of the system. According to the boundary conditions usually considered, the clumps are assumed to begin their evolution at time t0=0 with a separation d0=0, and are now moving apart or coming together for the first time in their history. In the case of a collision, we assume that the time t0=0 with separation d0=0 is the time of their core crossing and that we are looking at the system a time t after.

As for the first collision, the values of relevant parameters for the two-clump system (the E group and the SW+NW complex) are the relative LOS velocity in the rest-frame, $V_{\rm rf}\sim 800$ km s-1, and the projected linear distance between the two clumps, $D\sim 1$  h70-1 Mpc. To obtain an estimate of t, we use the Mach number of the shock ${\cal M}\sim 2.8$ as inferred by B09 from the radio spectral index of A2345-2. The Mach number is defined to be ${\cal M}={v}_{\rm s}/c_{\rm s}$, where ${v}_{\rm s}$ is the velocity of the shock and $c_{\rm s}$ is the sound speed in the pre-shock gas (see e.g., Sarazin 2002 for a review). The value of $c_{\rm s}$, obtained from our estimate of $\sigma_{ V,\rm E}\sim 900$ km s-1, leads to a value of ${v}_{\rm s}\sim 2.5\times 10^3$ km s-1. Assuming the shock velocity to be a constant, the shock covered a $\sim$0.9  h70-1 Mpc$\;$scale (i.e., the distance of the relic from the cluster center) in a time of $\sim$0.35 Gyr. We assume this time as our estimate of t. Although the velocity of the shock is not constant, studies based on numerical simulations show how the variation in  ${v}_{\rm s}$ is much smaller than the variation in the relative velocity of the subclumps identified with their dark matter components (see Fig. 4 of Springel & Farrar 2007, and Fig. 14 of Mastropietro & Burkert 2008), thus our rough estimate of t is acceptable as a first order approximation.

\begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f11.eps}} \end{figure} Figure 11:

The same that in Fig. 10 but for the NW and SW subsystems. Thick and thin lines give the results for $V_{\rm rf,LOS}\sim 100$ and 50 km s-1, respectively.

Open with DEXTER

The bimodal model solution gives the total system mass $M_{\rm sys}$, i.e. the sum of the masses of the E+SW+NW groups, as a function of $\alpha$, where $\alpha$ is the projection angle between the plane of the sky and the line connecting the centers of the two clumps (e.g., Gregory & Thompson 1984). Figure 10 compares the bimodal-model solutions with the observed mass of the system considering a 50% uncertainty band. The present solutions span the bound outgoing solutions (i.e., expanding), BO; the bound incoming solutions (i.e., collapsing), BI$_{\rm a}$ and BI$_{\rm b}$; and the unbound outgoing solutions, UO. For the incoming case, there are two solutions because of the ambiguity in the projection angle $\alpha$. The BO solution is the only one to be consistent with the observed mass range leading to a $\alpha\sim 40$-50 degrees. The BO solution means that the E group is moving towards East going in the opposite direction with respect to the observer, while the SW+NW complex is moving toward West toward the observer. The angle estimate means that the true spatial distance between the two subclumps is $\sim$1.3-1.5  h70-1 Mpc$\;$and that the real, i.e. deprojected, velocity difference is $V_{\rm rf} \sim 1000{-}1250$ km s-1. We note that the present relative velocity between galaxy clumps is smaller than the shock velocity, i.e., the regime is not stationary, but this is expected when comparing shock and collisionless components in numerical simulations (Springel & Farrar 2007; Mastropietro & Burkert 2008). The deprojected velocity difference of the two clumps at the cross core time is $V _{\rm rf} \sim 2000{-}2800$ km s-1, i.e. comparable to that of the shock as expected.

As for the second collision, the values of relevant parameters for the two-clump system (the SW and the NW groups) are $V_{\rm rf}\sim 0$ km s$^{-1} \;$(since we do not see any significant velocity difference), $D\sim 0.3$  h70-1 Mpc$\;$and $M_{\rm sys} \sim 1.1+0.5\sim 1.6$  $\times 10^{15}\;h_{70}^{-1}\;M_{\odot}$. In our scenario this merger is more recent than the above one and thus we assume a time $t\sim 0.2$ Gyr.

Figure 11 shows the results for $V_{\rm rf}=50$ km s-1 and $V_{\rm rf}=100$ km s-1. The BO solution is well acceptable with $\alpha\sim 2{-}20$ \ensuremath{{\rm ^\circ}}. This means that the SW group is moving toward South and the NW group towards North and, as someway expected from the observations, the NW and SW groups define a direction almost parallel to the plane of the sky. The true spatial distance is similar to the projected one. The real velocity difference has a huge uncertainty due to its dependence from the very uncertain $\alpha$ value. Notice that this merger should be considered of minor importance with respect to the main one described above, thus here the bimodal model is likely a worse representation of the reality than the above case.

In conclusion, we recover from our toy model quantitative results consistent with the global scenario. Thus, our analysis of the internal cluster dynamics well supports the ``outgoing merger shocks'' model proposed by B09. We also present an explanation for the relics asymmetry and the peculiarity of A2345-1. Our explanation is based on the internal structure of A2345, rather than on a possible merger with an external, close group (e.g. the group ``X1'' detected by B09). Before applying more reliable approaches to the study of A2345 internal dynamics, an improvement of the present observational picture would be very useful, e.g. by obtaining redshift information for more galaxies and deeper X-ray data (e.g. with Chandra and/or XMM telescopes) to compute gas temperature maps.

Acknowledgements
We are in debt with Annalisa Bonafede for the VLA radio image she kindly provided us. We thank the anonymous referee for his/her useful comments and suggestions. This publication is based on observations made on the island of La Palma with the Italian Telescopio Nazionale Galileo (TNG) and the Isaac Newton Telescope (INT). The TNG is operated by the Fundación Galileo Galilei - INAF (Istituto Nazionale di Astrofisica). The INT is operated by the Isaac Newton Group. Both telescopes are located in the Spanish Observatorio of the Roque de Los Muchachos of the Instituto de Astrofisica de Canarias.
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.

References

  1. Abell, G. O., Corwin, H. G. Jr., & Olowin, R. P. 1989, ApJS, 70, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Ashman, K. M., Bird, C. M., & Zepf, S. E. 1994, AJ, 108, 2348 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bagchi, J., Durret, F., Lima Neto, G. B., & Paul, S. 2006, Science, 314, 791 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Bardelli, S., Zucca, E., Vettolani, G., et al. 1994, MNRAS, 267, 665 [Google Scholar]
  5. Barrena, R., Ramella, M., Boschin, W., et al. 2005, A&A, 444, 685 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Barrena, R., Boschin, W., Girardi, M., & Spolaor, M. 2007a, A&A, 467, 37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Barrena, R., Boschin, W., Girardi, M., & Spolaor, M. 2007b, A&A, 469, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Barrena, R., Girardi, M., Boschin, W., & Dasí, M. 2009, A&A, 503, 357 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Beers, T. C., Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 [NASA ADS] [CrossRef] [Google Scholar]
  10. Beers, T. C., Forman, W., Huchra, J. P., Jones, C., & Gebhardt, K. 1991, AJ, 102, 1581 [NASA ADS] [CrossRef] [Google Scholar]
  11. Beers, T. C., Gebhardt, K., Huchra, J. P., et al. 1992, ApJ, 400, 410 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bird, C. M., & Beers, T. C. 1993, AJ, 105, 1596 [NASA ADS] [CrossRef] [Google Scholar]
  14. Böhringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bonafede, A., Giovannini, G., Feretti, L., Govoni, F., & Murgia M. 2009, A&A, 494, 429 (B09) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Boschin, W., Girardi, M., Barrena, R., et al. 2004, A&A, 416, 839 [Google Scholar]
  17. Boschin, W., Girardi, M., Spolaor, M., & Barrena, R. 2006, A&A, 449, 461 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Boschin, W., Barrena, R., & Girardi, M. 2009, A&A, 495, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Buote, D. A. 2002, in Merging Processes in Galaxy Clusters, ed. L. Feretti, I. M. Gioia, & G. Giovannini (The Netherlands: Kluwer Ac. Pub.): Optical Analysis of Cluster Mergers [Google Scholar]
  20. Burstein, D., & Heiles, C. 1982, AJ, 87, 1165 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cassano, R., & Brunetti, G. 2005, MNRAS, 357, 1313 [NASA ADS] [CrossRef] [Google Scholar]
  22. Cypriano, E. S., Sodré, L. Jr., Kneib, J.-P., & Campusano, L. E. 2004, ApJ, 613, 95 [NASA ADS] [CrossRef] [Google Scholar]
  23. Colless, M., & Dunn, A. 1996, ApJ, 458, 435 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dahle, H., Kaiser, N., Irgens, R. J., Lilje, P. B., & Maddox, S. J. 2002, ApJS, 139, 313 (D02) [NASA ADS] [CrossRef] [Google Scholar]
  25. Danese, L., De Zotti, C., & di Tullio, G. 1980, A&A, 82, 322 [NASA ADS] [Google Scholar]
  26. Dressler, A., & Shectman, S. A. 1988, AJ, 95, 985 [NASA ADS] [CrossRef] [Google Scholar]
  27. Ebeling, H., Voges, W., Böhringer, H., et al. 1996, MNRAS, 281, 799 [NASA ADS] [CrossRef] [Google Scholar]
  28. Ellingson, E., & Yee, H. K. C. 1994, ApJS, 92, 33 [NASA ADS] [CrossRef] [Google Scholar]
  29. Ensslin, T. A., & Gopal-Krishna 2001, A&A, 366, 26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Ensslin, T. A., Biermann, P. L., Klein, U., & Kohle, S. 1998, A&A, 332, 395 [NASA ADS] [Google Scholar]
  31. Fadda, D., Girardi, M., Giuricin, G., Mardirossian, F., & Mezzetti, M. 1996, ApJ, 473, 670 [NASA ADS] [CrossRef] [Google Scholar]
  32. Feretti, L. 1999, MPE Rep., 271 [Google Scholar]
  33. Feretti, L. 2002, The Universe at Low Radio Frequencies, Proc. IAU Symp., 199, held 30 Nov. - 4 Dec. 1999, Pune, India, ed. A. Pramesh Rao, G. Swarup, & Gopal-Krishna, 133 [Google Scholar]
  34. Feretti, L. 2005, X-Ray and Radio Connections ed. L. O. Sjouwerman, & K. K. Dyer, Published electronically by NRAO http://www.aoc.nrao.edu/events/xraydio, Held 3-6 February 2004 in Santa Fe, New Mexico, USA [Google Scholar]
  35. Feretti, L. 2006, Proc. XLIst Rencontres de Moriond, XXVIth Astrophysics Moriond Meeting: From dark halos to light, ed. L. Tresse, S. Maurogordato, & J. Tran Thanh Van, [arXiv:astro-ph/0612185] [Google Scholar]
  36. Feretti, L. 2008, Mem. SAIt, 79, 176 [NASA ADS] [Google Scholar]
  37. Feretti, L., Gioia, I. M., & Giovannini, G. 2002, AS&S Library, 272, Merging Processes in Galaxy Clusters (The Netherlands: Kluwer Academic Publisher) [Google Scholar]
  38. Feretti, L., Schuecker, P., Böhringer, H., Govoni, F., & Giovannini, G. 2005, A&A, 444, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Ferrari, C., Govoni, F., Schindler, S., Bykov, A. M., & Rephaeli, Y. 2008, Space Sci. Rev., 134, 93 [NASA ADS] [CrossRef] [Google Scholar]
  40. Giovannini, G., & Feretti, L. 2002, in Merging Processes in Galaxy Clusters, ed. L. Feretti, I. M. Gioia, & G. Giovannini (Kluwer Ac. Pub.: The Netherlands), Diffuse Radio Sources and Cluster Mergers [Google Scholar]
  41. Giovannini, G., Tordi, M., & Feretti, L. 1999, New Astron., 4, 141 [NASA ADS] [CrossRef] [Google Scholar]
  42. Giovannini, G., Bonafede, A., Feretti, L., et al. 2009, A&A, 507, 1257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Giovannini, G., Bonafede, A., Feretti, L., Govoni, F., & Murgia, M. 2010, A&A, 511, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Girardi, M., & Biviano, A. 2002, in Merging Processes in Galaxy Clusters, ed. L. Feretti, I. M. Gioia, & G. Giovannini (The Netherlands, Kluwer Ac. Pub.): optical Analysis of Cluster Mergers [Google Scholar]
  45. Girardi, M., & Mezzetti, M. 2001, ApJ, 548, 79 [NASA ADS] [CrossRef] [Google Scholar]
  46. Girardi, M., Fadda, D., Giuricin, G., et al. 1996, ApJ, 457, 61 [NASA ADS] [CrossRef] [Google Scholar]
  47. Girardi, M., Giuricin, G., Mardirossian, F., Mezzetti, M., & Boschin, W. 1998, ApJ, 505, 74 [NASA ADS] [CrossRef] [Google Scholar]
  48. Girardi, M., Barrena, R., & Boschin, W. 2007, Contribution to Tracing Cosmic Evolution with Clusters of Galaxies: Six Years Later conference, http://www.si.inaf.it/sesto2007/contributions/Girardi.pdf [Google Scholar]
  49. Goto, T., Sekiguchi, M., Nichol, R. C., et al. 2002, AJ, 123, 1807 [NASA ADS] [CrossRef] [Google Scholar]
  50. Gregory, S. A., & Thompson, L. A. 1984, ApJ, 286, 422 [NASA ADS] [CrossRef] [Google Scholar]
  51. Gullixson, C. A. 1992, in Astronomical CCD Observing and Reduction techniques, ed. S. B. Howell, ASP Conf. Ser., 23, 130 [Google Scholar]
  52. Hoeft, M., Brüggen, M., & Yepes, G. 2004, MNRAS, 347, 389 [NASA ADS] [CrossRef] [Google Scholar]
  53. Kempner, J. C., Blanton, E. L., Clarke, T. E. et al. 2003, Proc. The Riddle of Cooling Flows in Galaxies and Clusters of Galaxies, ed. T. H. Reiprich, J. C. Kempner, & N. Soker, [arXiv:astro-ph/0310263] [Google Scholar]
  54. Kennicutt, R. C. 1992, ApJS, 79, 225 [Google Scholar]
  55. Malumuth, E. M., Kriss, G. A., Dixon, W. Van Dyke, Ferguson, H. C., & Ritchie, C. 1992, AJ, 104, 495 [NASA ADS] [CrossRef] [Google Scholar]
  56. Mastropietro, C., & Burkert, A. 2008, MNRAS, 389, 967 [NASA ADS] [CrossRef] [Google Scholar]
  57. Ortiz-Gil, A., Guzzo, L., Schuecker, P., Böringer, H., & Collins, C. A. 2004, MNRAS, 348, 325 [NASA ADS] [CrossRef] [Google Scholar]
  58. Owen, F., & Ledlow, M. J. 1997, ApJ, 108, 41 [Google Scholar]
  59. Owen, F., Ledlow, M. J., & Keel, W. C. 1995, AJ, 109, 14 [NASA ADS] [CrossRef] [Google Scholar]
  60. Peres, C. B., Fabian, A. C., Edge, A. C., et al. 1998, MNRAS, 298, 416 [NASA ADS] [CrossRef] [Google Scholar]
  61. Pisani, A. 1993, MNRAS, 265, 706 [NASA ADS] [Google Scholar]
  62. Pisani, A. 1996, MNRAS, 278, 697 [NASA ADS] [CrossRef] [Google Scholar]
  63. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, in Numerical Recipes, Second Ed. (Cambridge University Press) [Google Scholar]
  64. Quintana, H., Carrasco, E. R., & Reisenegger, A. 2000, AJ, 120, 511 [NASA ADS] [CrossRef] [Google Scholar]
  65. Ramella, M., Boschin, W., Fadda, D., & Nonino, M. 2001, A&A, 368, 776 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Ramella, M., Biviano, A., Pisani, A., et al. 2007, A&A, 470, 39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Roettiger, K., Loken, C., & Burns, J. O. 1997, ApJS, 109, 307 [NASA ADS] [CrossRef] [Google Scholar]
  68. Roettiger, K., Burns, J. O., & Stone, J. M. 1999, ApJ, 518, 603 [NASA ADS] [CrossRef] [Google Scholar]
  69. Röttgering, H. J. A., Wieringa, M. H., Hunstead, R. W., & Ekers, R. D. 1997, MNRAS, 290, 577 [NASA ADS] [CrossRef] [Google Scholar]
  70. Sarazin, C. L. 2002, in Merging Processes in Galaxy Clusters, ed. L. Feretti, I. M. Gioia, & G. Giovannini (The Netherlands: Kluwer Ac. Pub.), The Physics of Cluster Mergers [Google Scholar]
  71. Smith, J. A., Tucker, D. L., Kent, S., et al. 2002, AJ, 123, 2121 [Google Scholar]
  72. Springel, V., & Farrar, G. R. 2007, MNRAS, 380, 911 [NASA ADS] [CrossRef] [Google Scholar]
  73. Struble, M. F., & Rood, H. J. 1987, ApJS, 63, 555 [NASA ADS] [CrossRef] [Google Scholar]
  74. Thompson, L. A. 1982, in Early Evolution of the Universe and the Present Structure, ed. G. O. Abell, & G. Chincarini (Dordrecht: Reidel), IAU Symp., 104 [Google Scholar]
  75. Tonry, J., & Davis, M. 1979, ApJ, 84, 1511 [Google Scholar]
  76. van Weeren, R. J., Röttgering, H. J. A., Bagchi, J., et al. 2009, A&A, 506, 1083 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Venturi, T., Giacintucci, S., Brunetti, G., et al. 2007, A&A, 463, 937 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Wainer, H., & Schacht, S. 1978, Psychometrika, 43, 203 [CrossRef] [Google Scholar]

Online Material

Table 1:   Velocity catalog of 125 spectroscopically measured galaxies in the field of the cluster A2345. ID. 82, in boldface, highlights the BCG.


Footnotes

... relics[*]
Table 1 is only available in electronic form at http://www.aanda.org
... IRAF[*]
IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

All Tables

Table 2:   Substructure from the photometric sample detected with DEDICA.

Table 3:   Kinematical properties of cluster and subclumps.

Table 1:   Velocity catalog of 125 spectroscopically measured galaxies in the field of the cluster A2345. ID. 82, in boldface, highlights the BCG.

All Figures

  \begin{figure} \par\includegraphics[width=16cm,clip]{14839f01.eps} \end{figure} Figure 1:

INT $r^{\prime }$-band image of the cluster A2345 (North at the top and East to the left) with, superimposed, the contour levels of the (smoothed) ROSAT/HRI archival image (pointing ID: US800769H, thick contours) and the contour levels of a VLA radio image at 1.4 GHz (see B09; beam 50 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}} $\times $ 38 \ensuremath {\ensuremath {{\rm ^\prime }}\hskip -0.1em\ensuremath {{\rm ^\prime }}}, thin contours starting at 0.24 mJy/beam and then spaced by a factor 2). Arrows highlight the contour levels of the two radio relics.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[height=18.5cm,clip]{14839f02.eps} \vspace*{3mm} \end{figure} Figure 2:

INT $r^{\prime }$-band image of the cluster A2345 (West at the top and North to the left). Circles and boxes indicate member and non-member galaxies, respectively (see Table 1). Labels indicate the IDs of cluster galaxies cited in the text.

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f03.eps}} \end{figure} Figure 3:

Redshift galaxy distribution. The solid line histogram refers to the 101 galaxies assigned to A2345. The redshift galaxy density, as provided by the adaptive kernel reconstruction method DEDICA, is showed in the insect panel, with an arbitrary normalization.

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f04.eps}} \end{figure} Figure 4:

Top panel: rest-frame velocity vs. projected clustercentric distance for the 101 galaxies in the two peaks forming the cluster (see Fig. 3). Crosses indicate galaxies rejected as interlopers by our ``shifting gapper'' procedure. Large squares indicate the five brightest galaxies. Middle and bottom panels: integral profiles of mean velocity and LOS velocity dispersion, respectively. The mean and dispersion at a given (projected) radius from the cluster-center are estimated by considering all galaxies within that radius - the first values computed on the five galaxies closest to the center. The error bands at the $68\%$ c.l. are also shown.

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f05.eps}} \end{figure} Figure 5:

The 98 galaxies assigned to the cluster. Upper panel: velocity distribution for all galaxies and for galaxies with $r^{\prime }<20$ (solid and dashed histograms, respectively). The big arrows indicate the position of the two peaks found by the 1D DEDICA algorithm. The Gaussians corresponding to the three galaxy clumps detected in the $r^{\prime }<20$ sample through the 3D DEDICA algorithm are shown (NW, SW and E clump from left to right, respectively). Notice that the separation of these three groups is much more evident in the 2D galaxy distribution (see Fig. 9 in the following). The small arrows show the velocities of the five brightest galaxies indicating the correspondent parent clump. Lower panel: stripe density plot where the arrows indicate the positions of the significant gaps.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8cm,clip]{14839f06.eps} \end{figure} Figure 6:

Upper panel: $r^{\prime }$- $i^{\prime }$ vs. $r^{\prime }$ diagram for galaxies with available spectroscopy. Circles and crosses indicate cluster members and non members. The solid line gives the median colour determined on member galaxies; the dashed lines are drawn at $\pm $0.15 mag from this value. Lower panel: $g^{\prime }$- $r^{\prime }$ vs. $r^{\prime }$ diagram for galaxies with available spectroscopy. The dashed lines are drawn at $\pm $0.3 mag from the median value determined on member galaxies.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=7.8cm,clip]{14839f07.eps} \vspace*{2mm} \end{figure} Figure 7:

Spatial distribution on the sky and relative isodensity contour map of likely cluster members with $r^{\prime }\le 21$, obtained with the 2D DEDICA method (external regions are not shown for the sake of clarity). The BCG is taken as the cluster center. The location of the five brightest galaxies are indicated by squares. Crosses indicate the peaks found in the ROSAT/HRI image. The circles indicate the position of the main (the largest circle) and the secondary peaks in the projected mass distribution as obtained from gravitational lensing analysis (see D02, notice that they analyzed only a small region around the BCG). The two relics are indicated in a schematic way by the ``+'' symbols. Solid lines highlight the directions of the two mergers. The dashed line is the bisecting of the merger directions (see text).

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8.5cm,clip]{14839f08.eps} \end{figure} Figure 8:

Galaxies belonging to structures as detected by the Voronoi Tessellation and Percolation technique. The algorithm is run on the sample of likely members with $r^{\prime }<21$ (see text). Open squares, asterisks and solid circles indicate galaxies in structures at the 95%, 98% and 99% c. ls., respectively. The E, NW and SW subclumps and the S group are indicated. As a comparison, big crosses indicate the positions of the subclumps found with DEDICA.

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{8.5cm}{!}{\includegraphics{14839f09.eps}} \end{figure} Figure 9:

Spatial distribution on the sky of the cluster galaxies (with $r^{\prime }<20$) showing the three groups recovered by the 3D KMM analysis. Solid circles, open circles and squares indicate the galaxies of E, SW and NW group, respectively. The BCG is taken as the cluster center. Large squares indicate the five brightest cluster members.

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f10.eps}} \vspace*{-4mm} \end{figure} Figure 10:

System mass vs. projection angle for bound and unbound solutions (thick solid and thick dashed curves, respectively) of the two-body model applied to the E and SN+SW subsystems. Labels BI$_{\rm a}$ and BI$_{\rm b}$ indicate the bound and incoming, i.e., collapsing solutions (thick solid curve). Label BO indicates the bound outgoing, i.e., expanding solutions (thick solid curve). Label UO indicates the unbound outgoing solutions (thick dashed curve). The horizontal lines give the range of observational values of the mass system with a 50% error. The thin dashed curve separates bound and unbound regions according to the Newtonian criterion (above and below the thin dashed curve, respectively).

Open with DEXTER
In the text

  \begin{figure} \par\resizebox{9cm}{!}{\includegraphics{14839f11.eps}} \end{figure} Figure 11:

The same that in Fig. 10 but for the NW and SW subsystems. Thick and thin lines give the results for $V_{\rm rf,LOS}\sim 100$ and 50 km s-1, respectively.

Open with DEXTER
In the text

Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.