EDP Sciences
Home All issues Volume 495 / No 2 (February IV 2009) A&A, 495 2 (2009) 379-387 Full HTML
Free Access
Issue
A&A
Volume 495, Number 2, February IV 2009
Page(s) 379 - 387
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361:200810731
Published online 20 November 2008

Galaxy evolution in Abell 85

I. Cluster substructure and environmental effects on the blue galaxy population

H. Bravo-Alfaro1 - C. A. Caretta1 - C. Lobo2,3 - F. Durret4 - T. Scott5


1 - Departamento de Astronomía, Universidad de Guanajuato, Apdo. Postal 144, Guanajuato 36000, Mexico
2 - Centro de Astrofísica da Universidade do Porto, Rua das Estrelas 4150-762 Porto, Portugal
3 - Depto. de Matemática Aplicada, Faculdade de Ciências, Univ. do Porto, R. do Campo Alegre 687, 4169-007 Porto, Portugal
4 - Institut d'Astrophysique de Paris, CNRS, UMR 7095, Université Pierre et Marie Curie, 98bis Bd Arago, 75014 Paris, France
5 - Center for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, AL10 9AB, UK

Received 1 August 2008 / Accepted 25 September 2008

Abstract
In this series of papers we explore the evolution of late-type galaxies in the rich cluster Abell 85. In this first paper we revisit the complex dynamical state of A 85 by using independent methods. First, we analyze the galaxy redshift distribution towards A 85 in the whole range 0-40 000  km s-1, and determine the mean redshifts of the background clusters A 87 and A 89, very close in projection to A 85. Then we search for substructures in A 85 by considering the 2D galaxy distribution of its members (13 000-20 000  km s-1) and by applying the kinematical 3D ${\Delta }$-test to both projected positions and radial velocities. This clearly reveals several substructures: one close to the cluster core and three more projected towards the southeast, along the region where an X-ray filament has been extensively studied. We also analyse the distribution of the brightest blue galaxies across a major fraction of the cluster volume, considering if they are gas-rich or poor. We report a very asymmetric distribution of the blue member galaxies, with most of them to the east and southeast, namely in the region joining the core of A 85 to its farthest substructure in this direction - dubbed the SE clump. By matching our sample of bright blue member galaxies with H  I detections reported in the literature, we identify gas-rich and gas-poor ones. As expected, the H  I-rich blue galaxies follow the same trend as the parent sample, with most of them projected on the eastern side of the cluster as well. Interestingly no blue objects have been detected in H  I up to a projected radius of 2 Mpc in this zone. We finally estimate the ram pressure stripping exerted by the intra-cluster medium as a function of the projected distance from A 85, in order to quantify how important this mechanism might be in sweeping the gas out of the infalling spirals.

Key words: galaxies: evolution - galaxies: clusters: individual: Abell 85, Abell 87

1 Introduction

Important evolution is known to occur between $z \sim 1$ and the present time. Galaxy clusters keep evolving throughout these epochs by accreting galaxies and groups of galaxies, mainly along filamentary structures (e.g. Adami et al. 2005, and references therein), which gives rise to merging episodes in a hierarchical fashion. Many clusters display virialized cores surrounded by infall regions of galaxies and/or groups, which are gravitationally bound but not in equilibrium with the main cluster body. The presence (or absence) of such substructures is relevant since it provides tight constraints on models of structure formation and evolution at large scales. Important changes are also evident for individual galaxies: the cosmic star formation rate may suffer a change of an order of magnitude since z = 1 (e.g. Lilly et al. 1996; Madau et al. 1996). Further major evidence of galaxy evolution is the decreasing number of spirals as one goes from the field to the central regions of clusters - and the corresponding rising population of lenticulars - in the nearby Universe (the morphology-density relation, Oemler 1974; Dressler 1980,2004, and references therein).

Nowadays an interesting debate has arisen about this last topic: while some authors discard major environmental effects on structural parameters such as size and mass concentration (Li et al. 2006), extensive multifrequency surveys over large volumes of galaxy clusters have, on the other hand, provided new insights into the environmental effects on both star formation activity and morphology of individual galaxies (see for example Cayatte et al. 1990; Chung et al. 2007; van Gorkom 2004; Vollmer et al. 2004). Furthermore, the presence of groups of E+A galaxies in very specific regions of clusters like Coma strongly suggests that they have suffered recent - and structural - environmental processing.

In this work we revisit the dynamical state of Abell 85 (hereafter A 85) and the strength of ram pressure stripping to produce the gas deficiency observed to large distances from the cluster center. In a forthcoming paper (Paper II) we will study the possible gravitational effects exerted on individual galaxies infalling from the SE, by considering their gas content, activity features and stellar morphology seen in NIR. The goal of this series of papers is to explore the transformation of late-type galaxies as a function of the cluster environment under different scenarios such as effects produced by the intra cluster medium (ICM) and/or gravitational preprocessing.

The richness class 1 cluster A 85, located at redshift z = 0.055, is a very interesting system that provides the perfect laboratory for studying the effects of environment on the galaxy population and on the properties of the galaxies themselves (gas content, SFR, AGN activity, etc.). This cluster has been extensively studied in the radio continuum, optical and X-rays (e.g. Kempner et al. 2002; Durret et al. 1998b; Giovannini et al. 1999; Durret et al. 2005). Many works have revealed a dynamically young system undergoing several mergings. For example, Kempner et al. (2002) proposed the existence of two subclusters around the inner part of A 85, and studied one of them based on Chandra observations; they report a subcluster much less massive than A 85, merging with the latter from the south. Two other Abell clusters, A 87 and A 89, are very close in projection to A 85 (Abell et al. 1989). Durret et al. (1998b) concluded that A 89 is rather constituted by two sheets of galaxies on the background, while A 87 would not be a separate cluster but rather a series of individual groups moving into the main body of A 85 from the SE. The global cluster orientation, as traced by X-rays and by the optical galaxy distribution, indicates an elongation towards the A 87 region, at a position angle of 160$^{\circ}$, coincident with the large scale structure in which the cluster is embedded. This elongation also coincides with the orientation of the southern X-ray filament seen in different surveys.

We follow a systematic strategy to detect substructures around A 85 based on a thorough analysis applied to a large set of optical photometric and spectroscopic data (described in Sect. 2). Our first approach is based on revisiting the galaxy redshift distribution in the direction of A 85, up to 40 000  km s-1 ( $z \sim 0.133$), and identifying the different systems along the line of sight (Sect. 3.1). We then probe, in Sect. 3.2, the 2D spatial distribution of the member galaxies of A 85 in the same field of view. For a better assessment of three-dimensional sub-structures, we perform a 3D kinematical analysis of the cluster making use of the ${\Delta }$-test (Dressler & Shectman 1988), which is one of the most sensitive algorithms to trace substructures. This is described in Sect. 3.3. As a complementary method for probing substructures, we carry out an analysis of the distribution of the brightest blue galaxies in A 85 (in Sect. 3.4), considering if they are H  I-rich or poor by matching positions with H  I detections reported previously in the literature. We discuss in Sect. 4 the results yielded by the different methods and trace a scenario of galaxy evolution in A 85 where ram pressure seems to have played an important role in the evolution of the bright blue galaxies. This is supported by our estimates of the strength of ram pressure stripping as a function of the projected distance from the cluster center. We summarize our results in Sect. 5.

Throughout this paper we assume $\Omega_{\rm M} = 0.3$, $\Omega_\Lambda = 0.7$, and H0 = 75 km s-1 Mpc-1. In this cosmology, the scale is 10$^{\prime}$ = 0.6 Mpc at the cluster mean redshift.

2 Observational data

We obtained position and optical magnitude data for galaxies in the A 85/87/89 complex from the SuperCOSMOS Science Archive, available from the Royal Observatory of Edinburgh[*]. The area covers a region of 50 arcmin radius around the center of A 85, taken as the position of its cD galaxy, MCG-02-02-086, which also coincides with the peak of the X-ray emission of the cluster. Magnitudes in $b_{\rm J}$ and $r_{\rm F}$ bands came from the digitization of UKST Survey plates (IIIa-J + GG395 and IIIa-F + OG590, respectively) by the SuperCOSMOS machine, using 10 $\mu$m size pixels (0.67'') (Hambly et al. 2001b). These magnitudes were placed on an absolute scale using SDSS-EDR (Stoughton et al. 2002), EIS (Arnouts et al. 2001) and ESO-Sculptor Survey (Arnouts et al. 1997) data, and are good to about 0.1 mag (rms). Note that a new calibration has been applied that supersedes the original SuperCOSMOS Sky Survey. The detection limit of SuperCOSMOS data is estimated to be about $b_{\rm J} = 20.5$ and $r_{\rm F} = 19.5$, at a completeness level above 95%, while the astrometry is accurate to better than 0.3 $^{\prime\prime}$ in both right ascension and declination (Hambly et al. 2001a).

This photometric catalogue was then matched with radial velocity data from the compilation of Andernach et al. (2005)[*], the largest compilation of velocities for Abell cluster member galaxies to date. For the A 85/87/89 complex (the area inside a radius of 50 arcmin around the center of A 85) this compilation contains 1693 entries, including multiple observations for most of the members. After discarding measurements with large deviations we averaged the radial velocities for each galaxy. This resulted in redshifts for 574 galaxies in the area of the A 85/87/89 complex in the range 0-40 000  km s-1. In the range 13 000-20 000  km s-1, corresponding to A 85 (e.g. Durret et al. 1998a), the compilation yields 367 redshifts for member galaxies. We found a match of 98.9% between these galaxies and SuperCOSMOS data.

In order to study the global H  I content in A 85 we consider the positions of twelve H  I-detected galaxies in this cluster recently reported by Bravo-Alfaro et al. (2008). This was obtained from an H  I survey of a region of about 1.5 $\times$ 1.5 degrees centered on A 85, observed with the VLA in its C-configuration, pointing the array at multiple positions and central frequencies. In this fashion the whole complex A 85/87/89 has been homogeneously surveyed with an average noise of rms $\sim $ 0.25 mJy beam-1, within the velocity range 14 580-18 440  km s-1. These observations have, on average, an H  I mass detection threshold of 7 $\times$ 108 $M_{\odot}$, corresponding to a 6$\sigma$ detection level. More details on the observations, on the complex data reduction process, and the gaseous properties of individual galaxies will soon appear (van Gorkom et al., in prep.).

3 Results on the galaxy distribution

3.1 The velocity distribution towards A 85

 \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig1.eps} \end{figure} Figure 1:

Histogram showing the distribution of 574 velocities in the region of the A 85/87/89 complex. The most conspicuous concentrations are marked according to their cluster and/or group association.
Open with DEXTER

 \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig2.eps} \end{figure} Figure 2:

Histogram showing the velocity distribution of galaxy-members of A 85, within the range 13 000-20 000  km s-1.
Open with DEXTER

We revisit the velocity distribution in the direction of A 85 by analyzing the redshift catalogue compiled in Sect. 2. The distribution of these radial velocities is shown in Fig. 1. As expected there are many concentrations of galaxies along this line of sight. The most remarkable one is the A 85 spike, within the velocity range 13 000-20 000  km s-1. This peak corresponds to 64% of the galaxies in the cone. We find a biweight mean radial velocity of 16 583 (+60, -52)  km s-1 and a velocity dispersion of 1122 (+54, -54)  km s-1 for the 367 member galaxies in this system, using a robust statistical ( ROSTAT) estimation (Beers et al. 1990). Their velocity distribution, enlarged in Fig. 2, is not far from a Gaussian. It shows a small positive skewness of 0.11, probably because of the small sub-peak between 19 000 and 20 000  km s-1. Most of the galaxies in this velocity sub-peak have positions coincident with a subcluster located in the western region of A85, that may be infalling into the main cluster. The kurtosis of the distribution is also positive (0.49), revealing a slightly peaked shape.

The next two prominent peaks along the line of sight to A 85, at 23 500 and 28 500  km s-1 respectively, whose centers are projected in the direction of A 89 (especially the former one), refer to the near and remote sheets identified by Durret et al. (1998b). Our sample consists of 77 and 48 galaxies for these two structures, respectively. The last peak ($\sim $39 000  km s-1) is centered on the A 87 region. This is most probably the system of galaxies associated with the A 87 cluster identified by Abell et al. (1989). A 87 is a distance class 5 cluster, which, following Abell's estimation, corresponds to about 39 000  km s-1 on average ( $z \sim 0.131$). With 32 member galaxies, the ROSTAT calculation gives a biweight mean of 39 014 (+113, -143)  km s-1 and a velocity dispersion of 737 (+259, -110)  km s-1 for this system of galaxies.

In general, our values of central velocity and velocity dispersion are in good agreement with recent works, for instance that of Christlein & Zabludoff (2003). In spite of a slight difference in the velocity range we consider (13 000-20 000  km s-1 against 13 423-19 737  km s-1 in their work), our systemic velocity (16 583  km s-1) is consistent with theirs (16 607  km s-1). Instead, a small difference is seen between our reported velocity dispersion (1122  km s-1) and Christlein & Zabludoff's value (993  km s-1), likely due to our larger velocity range and sample of galaxies (367 in this work against 280 in theirs).

3.2 The 2D galaxy distribution in A 85

 \begin{figure} \par\includegraphics[width=8.5cm,clip]{0731Fig3.eps} \end{figure} Figure 3:

Iso-density contours drawing the distribution of galaxies in the direction of A 85; red thin contours in the background correspond to all galaxies with magnitudes $b_{\rm J} < 20.5$. Black thick contours indicate the distribution of the 367 members of A 85 (13 000-20 000  km s-1). Main clusters and substructures of A 85 are indicated.
Open with DEXTER

To better understand the structure of A 85 we constructed an isodensity contour map for the cluster members, i.e., those 367 galaxies with radial velocities in the range 13 000-20 000  km s-1. This is shown in Fig. 3 with black contours. The galaxies were counted within square cells of 4.1$^{\prime}$ side, each cell displaced only half of its size from the previous one in order to smooth the contours. This provides a resolution of about 0.125 Mpc at the distance of A 85. The plotted contours correspond to 3, 5, 7, 9 and 11 galaxies per cell (the lowest level corresponds to 0.2 galaxies/arcmin2, while the peak reaches 0.8 galaxies/arcmin2). These contours are superimposed on the isodensity map of all galaxies with magnitudes $b_{\rm J} < 20.5$ in the A 85/87/89 complex area (red contours). The A 89 concentration appears only in red contours, whereas the A 87 region presents a superposition of two structures. These structures are the A 87 cluster (at 39 000  km s-1, seen in red contours) and, slightly displaced to the south, an infalling subcluster of A 85 (seen in black contours and labeled SE in Fig. 3), the latter coincident with one of the clumps found by Durret et al. (1998b) in the filament region.

The black isodensity contours constitute the first galaxy density map of A 85 built exclusively with galaxy members. It reveals three significant peaks: the first represents the core of the main cluster (indicated with A85 in the figure), nearly coincident with the cD galaxy; a secondary peak is coincident with the Southern Blob (SB as in Kempner et al. 2002); and a third, and less conspicuous one, appears some 35$^{\prime}$ (about 2.1 Mpc) to the south-east of the cluster center (indicated as SE in the figure), which corresponds to the subcluster superposed on A 87 already mentioned above. We analyze in more detail this SE subcluster in the next section.

3.3 The 3D ${\Delta }$-test

 \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig4.eps} \end{figure} Figure 4:

Substructures detected around A 85 following the kinematical analysis of the $\delta $ parameter. Larger black circles indicate higher probabilities of the presence of substructures. The most conspicuous substructures detected are indicated. The red cross marks the cluster center, coincident with the cD.
Open with DEXTER

As pointed above the galaxy density map of A 85 members with velocities between 13 000 and 20 000  km s-1 already shows a non-uniform distribution. Now we move a step forward in quantifying the presence of substructures by applying a more powerful statistical test that considers both positions and velocities. One of the most sensitive algorithms proposed in the literature (e.g. Pinkney et al. 1996) is the ${\Delta }$-test proposed by Dressler & Shectman (1988), which searches for deviations of the local velocity mean and dispersion from the global values. These statistics produce a deviation measure ($\delta $ parameter) for each galaxy, which can be used to find subclumps of ``deviant'' galaxies:

\begin{displaymath}% \delta^2 = \frac{(N_{\rm neighb} + 1)}{\sigma^2} [(v_{\rm local} - v)^2 + (\sigma_{\rm local} - \sigma)^2] \end{displaymath} (1)

where $N_{\rm neighb}$ is the number of neighbors used to estimate the local mean velocity ( $v_{\rm local}$) and velocity dispersion ( $\sigma_{\rm local}$), while v and $\sigma$ are the global mean velocity and velocity dispersion, respectively.

A correlated spatial and kinematic variation, i.e., the presence of a region (RA, Dec) containing several galaxies that deviate from the global kinematic values, indicates a high probability of substructure in this zone. In the following discussion we indicate these clumps of galaxies as showing high kinematical deviations from the global behavior. The cumulative value of $\delta $ for the entire cluster (${\Delta }$) can be compared to Monte Carlo simulations to estimate the confidence level of the general deviation of the cluster.

We applied the ${\Delta }$ statistics to the 367 members of A 85, using $N_{\rm neighb} = 10$, and the corresponding results are shown in Fig. 4 (we obtained similar solutions with values of $N_{\rm neighb}$ = 15, 19 $\simeq$ $\sqrt N$, and 37 $\simeq$ N/10). In Fig. 4 each galaxy is marked by a circle with a size proportional to $\delta $ (instead of $e^{\delta}$, as usual, since the signal in A 85 is very high). We take the arbitrary value $\delta $ = 2.0 as a convenient threshold to discriminate objects with kinematics that deviate significantly from the cluster one. Therefore, galaxies with $\delta > 2.0$ are indicated with black circles in Fig. 4, while green circles represent galaxies with $\delta < 2.0$. The cumulative ${\Delta }$ parameter for the cluster is larger than the one obtained in 99.3% of 1000 Monte Carlo random simulations.

We note five prominent regions displaying enhanced variation in kinematics with respect to the mean; these are spots of larger black circles in Fig. 4: the first one is close to the center, a few arcmin to the west, and we identify it as C2. Another substructure is found towards the south, coincident with the Southern Blob (SB), and two more appear in the southeast region of A 85: one along the X-ray filament detected in previous works, namely with ROSAT and XMM-Newton data (F, see Durret et al. 2003, and references therein), and the other along an extension of the latter (indicated with SE). We detect a last substructure to the west (W) of the cluster.

3.4 The distribution of the brightest blue galaxies in A 85

Another effective method to trace substructures in clusters consists of analysing the distribution (both in projected position and velocity) of their brightest blue members in a first step, and then verifying their gas content. This method has successfully revealed substructures at different merging stages in other systems (see, for example, van Gorkom et al. 2003; Dickey 1997; Bravo-Alfaro et al. 2000). We expect that spirals already undergoing strong environmental effects present different degrees of H  I deficiency, while the newcomers to the cluster neighborhood display normal gas content. In order to apply this method to A 85 we built a catalogue of the brightest blue member galaxies[*] based on the following strategy. (1) Magnitude and color: taking superCOSMOS data we produce a color-magnitude diagram to exclude the red sequence. In this fashion we select all galaxies brighter than $b_{\rm J} = 18.0$ and bluer than $b_{\rm J} - r_{\rm F} = 1.2$. (2) Position: we take only the objects within the reported region surveyed in H  I, as shown in Fig. 5. (3) Velocity: we only consider galaxies with known redshift, and within the effective H  I velocity coverage 14 600-18 400  km s-1; this represents a velocity interval roughly centered on the systemic value for A 85 (16 600  km s-1) and encompassing roughly three times the cluster velocity dispersion.

Table 1:   The sample of brightest blue member galaxies of A 85. Objects with an asterisk are H  I detections. ``#'' indicates strong H${\alpha }$ emission, and ``##'' indicates high probability of an AGN (see text).

There are 42 galaxies matching these criteria which constitute our sample of the brightest blue members of A 85; they are listed in Table 1. After the running number (first column), the next three columns list galaxy names and positions. Galaxies with an asterisk were detected in H  I and have reported gas masses of over 109 $M_{\odot}$. The symbol ``#'' indicates galaxies with strong H${\alpha }$ emission (as detected by Boué et al. 2008), and ``##'' indicates a high probability of harboring an AGN, according to coincident positions with X-ray point-like detections by Sivakoff et al. (2008). A particular case is the blue H  I-rich galaxy J0043017-094721 which, due to its projected position very close to A85[DFL98]425, had remained previously uncatalogued as an independent source. These two objects are roughly coincident with the position that Abell et al. (1989) reported for the cluster A 87, but are separated in velocity by more than 10 000  km s-1. The position of J0043017-094721 was determined in this work from a MEGACAM image taken at the 3.6-m CFHT (Boué et al. 2008). In Col. 5 we give the optical radial velocity of the galaxies, obtained from the compilation described in Sect. 2. The respective estimated error (rms) in velocity is given in Col. 6. Column 7 gives the $b_{\rm J}$ SuperCOSMOS magnitude of the galaxies. Column 8 gives the $b_{\rm J}-r_{\rm F}$ color. Column 9 gives the projected distance to the cluster center in Mpc.

 \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig5.eps} \end{figure} Figure 5:

Distribution of the brightest blue galaxies through A85, indicating the position of the H  I rich (encircled crosses) and the H  I poor galaxies (crosses). The center of the X-ray emission (red contours) as seen by ROSAT coincides with the cD galaxy. The ellipse roughly draws the position of the X-ray filament mentioned in the text. Large dashed circles correspond to the 1.0 Mpc and 2.4 Mpc radius. The blue external polygon indicates the VLA-H  I surveyed area.
Open with DEXTER

One of the most important results of this paper is shown in Fig. 5, where a very striking asymmetry in the blue-galaxy distribution through the whole cluster volume is seen. Crosses indicate the 42 brightest blue (cluster-member) galaxies plotted over a composite view of A 85 seen in optical and X-rays. As expected, the 12 H  I detections (encircled crosses) are within this sample, and the remaining 30 constitute the category of gas deficient (H  I-non detected) blue galaxies. To check if this asymmetric distribution is really a property of the blue members of A 85, we examined the distribution of non-emission line cluster galaxies taken from a sample of 232 objects with SDSS spectra, which should be representative of red members. As expected, these display a very regular distribution, concentrated around the cluster center (Caretta et al. 2008).

The spatial distribution of the H  I detections follows the same asymmetric trend as blue galaxies: all but two are projected east of the center of A 85 where they are roughly distributed in two regions. The first group of five H  I detections is projected at an average distance of $\sim $30$^{\prime}$ NE from the center of A 85. An interesting feature to notice is that these five H  I-rich blue galaxies (A85[DFL98]374/461/491/502 and [SDG98]3114) display an intriguingly low velocity dispersion of just 120  km s-1, but are spread (in projection) across a large area of $\sim $1.2 Mpc.

The second group of H  I detections is projected SE of A 85, around the position of our detected SE substructure (see Fig. 4). Five galaxies are detected in H  I in this region with at least four of them having higher probabilities of belonging to a real subcluster when compared to those in the NE (see below). In addition to these H  I detections (A85[DFL98]323, 486, 347 and J004301-094721), many other (gas poor) blue objects appear in the SE quadrant. The H  I-rich galaxies are projected, on average, some 35$^{\prime}$ SE of the center of A 85, along the extension of the X-ray filament (indicated by an ellipse in the figure). Similar to the five galaxies detected NE of A 85, these are also spread over a large area (more than one Mpc, Bravo-Alfaro et al. 2008), and they form a grouping with low velocity dispersion (see Table 1).

The asymmetry in the distribution of blue galaxies is better quantified if we refer to the quadrants indicated in the figure. SW $\to$ NW $\to$ NE $\to$ SE gives the sequence from the poorest to the richest in blue galaxies. Relative to the other quadrants, the SW region is almost empty of blue galaxies, containing only four of these objects (all H  I-poor). In comparison, six blue objects are seen in the NW and nine in the NE, where two and five H  I detections were found respectively. By far the richest region in blue galaxies is the SE, where 23 blue objects are distributed from the highest density zone of the cluster ICM, going through the X-ray filament area (indicated by an ellipse), and then spread around the position of the SE substructure. Only five of those are detected in H  I. Most of the blue galaxies projected on this region display velocities in the range 15 000-16 300  km s-1 (Table 1), and several display wide H${\alpha }$ equivalent widths (larger than 25 Å). The presence of this indicator of star formation (Boué et al. 2008, and Paper II in prep.) strongly supports the hypothesis of a merging between a subcluster of galaxies and A 85.

4 Discussion

4.1 Substructures across A 85

The results of the previous section clearly reveal the following substructures.


The core: the central region of the cluster appears to be dynamically very active and presents strong signs of past activity (Kempner et al. 2002; Durret et al. 2005). We confirm this by detecting the most prominent substructure (C2) projected $\sim $5$^{\prime}$ west of the cluster center and consisting of two blobs. The ridge defined by these two blobs is oriented in a SE-NW axis, and the position angle of this axis is roughly the same as the orientation of the X-ray filament and its extension to the SE substructure. Most of the galaxies in C2 present velocities around 15 800  km s-1, well under the systemic value of A 85 (16 600  km s-1).

The southern part of C2 has been seen in X-rays as a diffuse source with Chandra (see Fig. 1 of Kempner et al. 2002, the structure which they call ``the southwest subcluster''). Here we report that the clump of galaxies associated with this structure further extends some 10$^{\prime}$ to the N (see Fig. 4), with another clump centered on the second brightest galaxy of A 85 (KAZ 364). This is a very peculiar spiral reported as an individual X-ray source by Sivakoff et al. (2008), and classed an AGN. In spite of the large area covered by this survey, the authors only report individual X-ray sources; no diffuse emission is detected around KAZ 364, equivalent to that seen by Kempner et al. (2002) in the southern tail of this grouping. This object also has the lowest radial velocity in the whole A 85 sample and, being outside the VLA velocity coverage, it is not included in Table 1. The $\delta $ signals associated with C2 decreases if we remove this galaxy from the sample, but the difference is not significant. The whole C2 structure, with diffuse X-ray emission on its southern tail, strongly suggests a subcluster undergoing an advanced merging process with the main cluster core, in agreement with the general scenario proposed by Durret et al. (2005).


The southern blob and filament: region SB, projected some 10$^{\prime}$ south of the cD, probably consists of an almost disrupted group of galaxies merging from the south (Kempner et al. 2002). Intriguingly enough, the galaxies within this structure display velocity values around 17 000  km s-1, i.e. larger than the substructures C2, F and SE. Clearly identified in Fig. 3, this substructure has the most accentuated density enhancement, even comparable to the core of A 85, while its kinematics do not deviate much from the global pattern - this is indicated with small green circles in Fig. 4 at the same position. The opposite occurs with clump F, a substructure projected midway along the X-ray filament described by Durret et al. (2003), about 20$^{\prime}$ away from the cluster center. This is not detected as a significant enhancement in the isodensity plot (Fig. 3) but deviates very much from the global kinematics (see the large black circles in Fig. 4). This structure is between the subgroup 3 in Durret et al. (1998b) and the southern blob.


The SE region: as outlined in the previous sections, the SE is a very interesting region. The substructure marked SE in Fig. 4 is close to the projected position of the A 87 cluster, roughly in the middle of the X-ray filament as detected in the ROSAT PSPC field (Durret et al. 1998b), or equivalently, at the end of the filament portion detected in the XMM-Newton field (Durret et al. 2003). This is a clump with very prominent kinematic deviations (large black circles) but is barely detected as a galaxy density enhancement in Fig. 3. The most likely explanation is that we are dealing with a relatively poor subcluster undergoing a high velocity infall towards A 85. Many of the galaxies within this substructure display velocities around 15 800  km s-1, the same as substructures F and C2. Most of the SE galaxies have blue colors, indicative of late type morphologies, and the galaxies detected in H  I southeast of A 85, being within this velocity range, are certainly associated with this infalling group. In Paper II we analyse the spectral properties and NIR images of selected galaxies in this area.


The W group: another region with several large $\delta $ values is seen about 35$^{\prime}$ west of the cluster center (marked with W). It is roughly at the same projected distance from the cluster center as the SE clump. These objects are outside the useful VLA-FOV and we have no information on their gas content; furthermore it is outside all the reported X-ray surveys. This group of galaxies displays a velocity around 18 000  km s-1. As this group is not projected along a major large scale filament, it is impossible to say if a dynamical relation exists with A 85.


The NE region: the H  I-detected galaxies in the NE quadrant of A 85 (encircled crosses in Fig. 3), in spite of displaying a tight velocity dispersion, are neither detected in the isodensity plot (Fig. 3) nor as a kinematical structure (Fig. 4). We applied larger values of $N_{\rm neighb}$ unsuccessfully trying to separate these objects from the global kinematics. This is likely due to the very large area over which they are spread and to the small number of galaxies being part of this possible structure.

4.2 H  I-rich vs. H  I-poor galaxies: the radial distribution

Striking observational evidence found in many of the clusters surveyed in H  I is the presence of galaxies and groups, projected at large radius from the cluster center, that have been stripped of a major fraction of their gas component (e.g. Solanes et al. 2001, and references therein). Their large distance from the densest zones of the ICM raises the question of whether the typical ram pressure stripping approach could account for such effects. The best known examples have been reported in Coma (Gavazzi 1987; Bravo-Alfaro et al. 2000) and in Perseus I (Levy et al. 2007). Based on Arecibo data, Gavazzi (1987) reported a number of H  I deficient galaxies at large radii from Coma, projected on the intermediate zone between this cluster and A 1367. In addition to this, Bravo-Alfaro et al. (2000) reported one group of late-type galaxies projected 0.9 Mpc east of the Coma center, with strong H  I deficiency. More recently, Levy et al. (2007) found a very interesting picture in Pegasus I, a system which lacks a harsh ICM and where we do not expect ram pressure to strip large amounts of gas. These authors report a number of spirals, near the cluster core but also in foreground and background groups, ranging from moderate H  I deficiency to truncated gas disks. These cases can be added to other stripped objects in low density environments, such as NGC 4522 in Virgo (Kenney et al. 2004, and references therein) and NGC 2276 in the NGC 2300 group (Davis et al. 1997). With the same trend we have A 85, with many gas-poor blue galaxies lying at 2-2.5 Mpc from the cluster center in projection.

A very important step forward on this issue was made by Tonnesen et al. (2007), who carried out a cosmological simulation of cluster evolution. They found that gas stripping may occur out to the cluster virial radius, as observed in A 85 and the clusters mentioned above. Very interestingly, their model predicts that ram pressure stripping remains the most important mechanism responsible for the loss of inter-stellar gas, and they found that a fraction of galaxies is affected by this process in the intermediate cluster zone (1.0-2.4 Mpc) and up to the periphery (2.4-5.0 Mpc). This is in close agreement with our observational results. With the help of Table 1 and Fig. 5 we see that, among the 42 brightest blue member galaxies of A 85, 30 are deeply gas deficient and constitute our sample of blue gas-poor galaxies. The radial distribution of these H  I-deficient objects is as follows: only 8 are within the inner radius (r < 1.0 Mpc) while 22 are in the transition zone (1.0-2.4 Mpc). These gas-poor blue galaxies along the transition zone of A 85 would match the stripped population predicted by Tonnesen et al. (2007).

 \begin{figure} \par\includegraphics[angle=-90,width=8.7cm,clip]{0731Fig6.eps} \end{figure} Figure 6:

Histogram of the radial distribution of the two populations: H  I-rich (gray columns) vs. H  I-poor (empty columns) galaxies. This clearly shows that, on average, the former are found at larger distances from the cluster center.
Open with DEXTER

Figure 6 clearly shows that projected positions of H  I-rich galaxies (filled columns), measured from the cluster center, follow a different distribution when compared with the H  I-non detections (empty columns, see also Table 1). The histogram displays bins of 0.5 Mpc, where the X-axis corresponds to the radius measured from the cD galaxy. Average positions are 2.07 Mpc and 1.42 Mpc, for H  I-rich and H  I-deficient galaxies, with standard deviations of 0.53 and 0.59 respectively. In order to examine the difference between both distributions we applied the Kolmogorov-Smirnov test, which is sensitive to global shape differences. We find they are different with a significance level of 0.99. We complement this by applying the t distribution to explore the difference of two means (case n1 + n2 - 2 = 40 degrees of freedom). The hypothesis $H_0: \mu_1 = \mu_2$ is discarded with a confidence level of 0.99.

4.3 Ram pressure stripping in A 85

In order to quantify the strength of the ICM-ISM interaction exerted on our sample of blue galaxies around A 85 we compute, as a first approach, the cluster gas density distribution and extrapolate it to radii of 2.0-2.5 Mpc. We follow the hydrostatic-isothermal $\beta$-model of Cavaliere & Fusco-Femiano (1976), given by:

\begin{displaymath}% \rho(r) = \rho_o [1+(r/r_{\rm c})^2]^{-3\beta/2}. \end{displaymath} (2)

We take the parameters of the gas profile of A 85 from Chen et al. (2007): core radius $r_{\rm c}$ = 82 kpc, central density $\rho_o$ = 0.0257 cm-3 (and assuming that this is an electrically neutral gas), and $\beta$ = 0.532. We then compute the ram pressure stripping as a function of radius from the cluster center following Gunn & Gott's (1972) equation: $P_{\rm ram} = \rho_{\rm ICM}$ $\times$ $v_{\rm rel}^2$, where $\rho_{\rm ICM}$ is the density of the ICM at the galaxy position, and  $v_{\rm rel}$ is the galaxy velocity relative to the ICM. Ideally, the latter should be the component of the real velocity perpendicular to the galaxy disk (as defined e.g. in Cayatte et al. 1994), but not having information on the inclination angle nor on the velocity of the galaxy in the plane of the sky we will thus consider radial velocities, relative to the cluster systemic movement, as a first approach. Normally, since galaxy tangential velocity components are expected, this approximation underestimates the real  $v_{\rm rel}$, making our value of $P_{\rm ram}$ a lower limit of the real one.

For the estimation of $P_{\rm ram}(r)$ we consider three values of  $v_{\rm rel}$: first we take 1000  km s-1, roughly matching the (radial) velocity dispersion of the cluster, which is also a lower limit of the (radial) relative velocity between the SE clump ($\sim $15 800  km s-1) and A 85 ($\sim $16 600  km s-1). We include two other velocity values as reference, 500  km s-1 and 1500  km s-1. Figure 7 displays the output of $P_{\rm ram}$ for values between 1.0 and 2.6 Mpc.

In order to quantify the effects of ram pressure stripping on individual galaxies we compare these values with the restoring force, defined as $F_{\rm r}$ = $2\pi~G~\sigma_{\rm tot}~\sigma_{\rm gas}$. This is a measure of the galaxy autogravitational binding force per cm-2 (having units of dyn cm-2, it can be directly compared with the ram pressure). For this comparison we take the typical value of $F_{\rm r}$ = 2 $\times$ $10^{-12}~{\rm dyn~cm}^{-2}$, obtained by Cayatte et al. (1994) for spiral galaxies in Virgo. If we inspect Fig. 7 we notice that, in spite of the low ICM density values in the transition zone (where the SE clump is found), $P_{\rm ram}$ is larger than the restoring force up to 2.0 Mpc from the cluster center for $v_{\rm rel}$ $\ge$ 1000  km s-1. Even the lowest velocity in our estimation ( $v_{\rm rel}$ = 500  km s-1) produces values of ram pressure of the same order of magnitude as $F_{\rm r}$ within the inner 1 Mpc. This unexpected result confirms that we should not discard ram pressure stripping as a mechanism playing a major role in driving galaxy evolution when dealing with high infall velocities relative to a cluster with a dense ICM such as that of A 85.

5 Summary

Our results can be summarized as follows:

1.
We first analysed the distribution of galaxies with measured redshift in the line of sight towards the complex A 85/87/89, within a region of about 1.5 degrees side around the center of A 85. For A 85 we find 367 member galaxies in the velocity range 13 000-20 000  km s-1. We confirm the presence of two background structures near the position of A 89. The farthest cluster in our sample is A 87, at $\sim $39 000  km s-1.

2.
We revisited the 2D (RA, Dec) distribution of galaxies through the whole A 85/87/89 complex and trace for the first time the isodensity contour map only with members of A 85. This highlights the substructures of A 85 from the background clusters. Three major density peaks are revealed: two in the inner part of the cluster and one at $\sim $2 Mpc in the SE direction, superimposed on the projected position of A 87.

3.
We ran the 3D kinematical ${\Delta }$-test (RA, Dec, vel) for the member galaxies of A 85. We detect different clumps of galaxies going from the center of A 85 and extending towards the SE. The first substructure reveals the inner structure of the core of A 85, with a group of galaxies in an advanced merging stage with the main cluster body. We also report a structure coincident with the X-ray Southern Blob, roughly 0.7 Mpc south of the center of A 85. Another substructure appears along the X-ray filament, more than 1 Mpc SE of the cluster center. Finally, the SE clump is projected at more than 2 Mpc from the center of A 85, at $\sim $15 800  km s-1, coincident with the infalling sub-cluster revealed by the 2D analysis.

 \begin{figure} \par\includegraphics[angle=-90,width=8.7cm]{0731Fig7.eps} \end{figure} Figure 7:

The ram pressure stripping, in units of 10-12 dyn cm-2, as a function of distance from the center of A 85, estimated for radii up to $\sim $2.5 Mpc. Three different velocities, relative to the ICM, were applied: 500, 1000 and 1500  km s-1. As a comparison, the solid line indicates the value of the restoring force $F_{\rm r}$ (see text).
Open with DEXTER

4.
We built a sample of the brightest blue galaxies, with 42 members of A 85 in the velocity range 14 600-18 400  km s-1, and we trace their spatial distribution throughout the cluster. By taking into account the region surveyed in H  I and the velocity range covered by Bravo-Alfaro et al. (2008), we obtain information on the gas content of this sample, with just twelve objects (H  I-detections) showing rather normal gas content or mild deficiency; the remaining 30 blue objects represent deficient galaxies. The spatial distribution of the whole 42 blue brightest members is totally unexpected, with about half of them projected on the SE quadrant of the cluster. This strongly supports the presence of sub-clusters in an early stage of merging with A 85, infalling with a relative velocity greater than 1000  km s-1.

5.
The distribution of the blue sample is also remarkable when measured as a function of the projected distance to A 85. From the sample of 42 blue members, 32 are projected on the transition zone (1.0-2.4 Mpc), and 22 of these are very gas deficient. We compare the positions of the two populations as a function of projected radius from the cluster center. We show that H  I-rich galaxies are projected, on average, much further than H  I-deficient ones: the mean projected distances are 2.1 Mpc and 1.4 Mpc respectively.

6.
We estimated ram pressure stripping in A 85 as a function of radius in order to quantify the environment's ability to decrease the gas fraction in spirals, in particular in the gas-poor blue galaxies reported in the SE region, $\sim $2 Mpc away from the center of A 85. Considering gas density parameters of A 85 ($\beta$ model and X-ray parameters) and different values of velocity relative to the ICM, we found that ram pressure is several times larger than the restoring force for typical spirals, when velocities relative to A 85 are over 1000  km s-1. This would explain the presence of the impressive sample of gas-poor blue-galaxies projected at the edge of the transition zone, i.e. at radius $\ga$2 Mpc. The observational evidence shown here matches theoretical cluster models such as that of Tonnesen et al. (2007) who predict active gas stripping in that transition zone.
These results reveal an even more complex structure of A 85 than previously found, with groups of galaxies at different merging stages, especially in the SE region but in the cluster core as well. The observational evidence shown in this work confirms that the cluster environment is playing an important role in driving galaxy evolution. We show that ram pressure stripping might be effective up to unusually large radial distances from a cluster center, when infalling velocities are large enough. In a forthcoming paper (Paper II) we will quantify the effects of gravitational mechanisms in disturbing stellar disks, in particular in the SE region.

Acknowledgements
We thank the anonymous referee for her/his very useful comments which helped us to strenghten and improve this paper. The authors are also grateful to H. Andernach and E. Tago for providing radial velocity data from their compilation. H.B.A. acknowledges CONACyT grant 50794, and financial support from Institut d'Astrophysique de Paris and Centro de Astrofísica da Universidade do Porto, allocated to carry out working stays at these institutes. C.A.C. also aknowledges support from CONACyT, grant 50921, and the University of Guanajuato (DINPO) for the grant 0102/06. C. Lobo aknowledges finantial support from project PTDC/CTE-AST/66147/2006 (FCT, Portugal).

References

Footnotes

... Edinburgh[*]
http://surveys.roe.ac.uk/ssa/index.html
...Andernach et al. (2005)[*]
Most of the data for A 85 in this compilation comes from Durret et al. (1998a); Malumuth et al. (1992); Smith et al. (2004); Beers et al. (1991); Christlein & Zabludoff (2003) and SDSS (Stoughton et al. 2002).
... galaxies[*]
In the absence of morphological classification, we expect that bright blue objects will very likely be late-type galaxies.

All Tables

Table 1:   The sample of brightest blue member galaxies of A 85. Objects with an asterisk are H  I detections. ``#'' indicates strong H${\alpha }$ emission, and ``##'' indicates high probability of an AGN (see text).

All Figures

  \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig1.eps} \end{figure} Figure 1:

Histogram showing the distribution of 574 velocities in the region of the A 85/87/89 complex. The most conspicuous concentrations are marked according to their cluster and/or group association.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig2.eps} \end{figure} Figure 2:

Histogram showing the velocity distribution of galaxy-members of A 85, within the range 13 000-20 000  km s-1.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8.5cm,clip]{0731Fig3.eps} \end{figure} Figure 3:

Iso-density contours drawing the distribution of galaxies in the direction of A 85; red thin contours in the background correspond to all galaxies with magnitudes $b_{\rm J} < 20.5$. Black thick contours indicate the distribution of the 367 members of A 85 (13 000-20 000  km s-1). Main clusters and substructures of A 85 are indicated.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig4.eps} \end{figure} Figure 4:

Substructures detected around A 85 following the kinematical analysis of the $\delta $ parameter. Larger black circles indicate higher probabilities of the presence of substructures. The most conspicuous substructures detected are indicated. The red cross marks the cluster center, coincident with the cD.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8cm,clip]{0731Fig5.eps} \end{figure} Figure 5:

Distribution of the brightest blue galaxies through A85, indicating the position of the H  I rich (encircled crosses) and the H  I poor galaxies (crosses). The center of the X-ray emission (red contours) as seen by ROSAT coincides with the cD galaxy. The ellipse roughly draws the position of the X-ray filament mentioned in the text. Large dashed circles correspond to the 1.0 Mpc and 2.4 Mpc radius. The blue external polygon indicates the VLA-H  I surveyed area.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[angle=-90,width=8.7cm,clip]{0731Fig6.eps} \end{figure} Figure 6:

Histogram of the radial distribution of the two populations: H  I-rich (gray columns) vs. H  I-poor (empty columns) galaxies. This clearly shows that, on average, the former are found at larger distances from the cluster center.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[angle=-90,width=8.7cm]{0731Fig7.eps} \end{figure} Figure 7:

The ram pressure stripping, in units of 10-12 dyn cm-2, as a function of distance from the center of A 85, estimated for radii up to $\sim $2.5 Mpc. Three different velocities, relative to the ICM, were applied: 500, 1000 and 1500  km s-1. As a comparison, the solid line indicates the value of the restoring force $F_{\rm r}$ (see text).
Open with DEXTER
In the text

Copyright ESO 2009

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.