BUDHIES III

BUDHIES III: The fate of Hi and the quenching of galaxies in evolving environments

Abstract

In a hierarchical Universe clusters grow via the accretion of galaxies from the field, groups and even other clusters. As this happens, galaxies can lose their gas reservoirs via different mechanisms, eventually quenching their star formation. We explore the diverse environmental histories of galaxies through a multiwavelength study of the combined effect of ram-pressure stripping and group ‘processing’ in Abell 963, a massive growing cluster at from the Blind Ultra Deep Hi Environmental Survey (BUDHIES). We incorporate hundreds of new optical redshifts (giving a total of 566 cluster members), as well as Subaru and XMM-Newton data from LoCuSS, to identify substructures and evaluate galaxy morphology, star formation activity, and Hi content (via Hi deficiencies and stacking) out to . We find that Abell 963 is being fed by at least seven groups, that contribute to the large number of passive galaxies outside the cluster core. More massive groups have a higher fraction of passive and Hi -poor galaxies, while low-mass groups host younger (often interacting) galaxies. For cluster galaxies not associated with groups we corroborate our previous finding that Hi gas (if any) is significantly stripped via ram-pressure during their first passage through the intracluster medium, and find mild evidence for a starburst associated with this event. In addition, we find an overabundance of morphologically peculiar and/or star-forming galaxies near the cluster core. We speculate that these arise as groups pass through the cluster (post-processing). Our study highlights the importance of environmental quenching and the complexity added by evolving environments.

keywords:
galaxies: clusters: general – galaxies: clusters: individual: Abell 963 – galaxies: evolution – galaxies: interactions – galaxies: peculiar.

1 Introduction

It has been established by numerous observational efforts that star formation in galaxies has steeply declined since (e.g. Bell et al., 2005; Daddi et al., 2007; Madau & Dickinson, 2014). The declining star formation activity is accompanied by a decrease in the fraction of late-type galaxies and a comparative increase of the early-type galaxy population (Desai et al., 2007; Poggianti et al., 2009) that is particularly visible in clusters (Dressler et al., 1997; Fasano et al., 2000). In fact, studies of large samples of nearby galaxies show a strong correlation between the passive galaxy fraction and both stellar mass and environment (Baldry et al., 2006; Haines et al., 2007; Peng et al., 2010). Such findings have suggested that, the star formation activity of galaxies, as well as their morphologies, are strongly correlated with stellar mass, and that the quenching of star formation is accelerated in dense environments such as galaxy clusters. All of the aforementioned observations are consistent with a scenario in which galaxies are quenched as clusters grow hierarchically with cosmic time (e.g. De Lucia et al., 2012), and many galaxies transition from the field into groups and clusters. What remains unclear is which are the physical mechanisms driving galaxy evolution in different environments. The most likely mechanisms at play in groups and clusters (see Boselli & Gavazzi, 2006, for a review) can be divided in three broad categories.

  • Gravitational interactions among galaxies. Simulations have shown that a major merger between two disc galaxies can result in a massive elliptical (Toomre & Toomre, 1972; Barnes & Hernquist, 1996), while the accretion of small satellites by a disc galaxy can create an S0 (Walker, Mihos & Hernquist, 1996). Galaxy-galaxy interactions however, are expected to be more frequent in low-mass groups, or the field, where the relative velocities are low enough for the collisions to occur. In clusters, galaxies can instead experience fast encounters with other cluster members. The accumulation of such encounters, known as harassment (Moore et al., 1999), can cause disc thickening and gas fuelling of the central region. Harassment can also strip Hi from galaxies (Duc & Bournaud, 2008), although it is most effective in low-surface-brightness galaxies, and towards the core of galaxy clusters (Moore et al., 1999; Smith et al., 2015).

  • Galaxy-Cluster interaction. The gravitational potential of a cluster can affect the cluster galaxies, by, for example, compressing the gas, inducing inflow, forming bars, and centrally concentrating the star formation (see e.g. Miller, 1986; Byrd & Valtonen, 1990). However, Boselli & Gavazzi (2006) showed that gas is hardly directly removed in such interactions.

  • Interaction between galaxies and the intracluster medium (ICM). Gunn & Gott (1972) first predicted (analytically) that ram-pressure can effectively remove a galaxy’s interstellar medium (ISM) as it passes through a dense ICM, thereby quenching its star formation. The efficiency of the stripping is proportional to the density of the ICM () and to the square of the galaxy’s velocity (). When the ram-pressure is greater than the anchoring force of the galaxy, the diffuse ISM will be stripped (Abadi, Moore & Bower, 1999; Quilis, Moore & Bower, 2000; Vollmer et al., 2001; Jaffé et al., 2015). Some models predict that Hi can be fully depleted after a single transit through the ICM’s core (Roediger & Brüggen, 2007). It is also possible that only the diffuse hot halo gas can be ram-pressure stripped or thermally evaporated (Bekki, Couch & Shioya, 2002). This process, known as starvation or strangulation, halts the funnelling of gas into the disc, provoking a slow decline of the galaxies’ star formation activity. The quenching time-scale of starvation ( Gyr) is much longer than that of ram-pressure stripping (RPS) of the disc gas (see Peng, Maiolino & Cochrane, 2015, and references therein).

It is possible that at least some of the aforementioned mechanisms act at the same time, increasing the difficulty of separating them observationally. One complication is that clusters continuously grow via the accretion of galaxies from the field, filaments and/or groups. In fact, it has been proposed that galaxies can start to be quenched in low-mass groups, that “pre-process” them (via galaxy-galaxy interactions, galaxy-group tidal interactions, and/or RPS) before falling into a more massive cluster (e.g. Zabludoff & Mulchaey, 1998; Balogh et al., 1999; Verdes-Montenegro et al., 2001; Lewis et al., 2002; McGee et al., 2009; De Lucia et al., 2012; Jaffé et al., 2012). Although the relative importance of group pre-processing in the quenching of cluster galaxies is still debated, it has been estimated from cosmological simulations that a cluster at has accreted of its galaxies () from groups more massive than (McGee et al., 2009). Moreover, recent observations show that around half of galaxies in groups, and a significant fraction of galaxies in filaments show signs of quenched star formation activity (Cybulski et al., 2014).

In addition, the simulations of Vijayaraghavan & Ricker (2013) showed that during a group’s pericentric passage within a cluster, galaxy-galaxy collisions are enhanced. Shock waves produced by the merger also enhance the ram-pressure intensity on group and cluster galaxies, while the increased local density during the merger leads to greater galactic tidal truncation. The combination of these effects has been referred to as group ‘post-processing’. Observations suggest that % of clusters at are undergoing mergers with other clusters (Katayama et al., 2003; Sanderson, Edge & Smith, 2009; Hudson et al., 2010), and this fraction is expected to increase with redshift (Mann & Ebeling, 2012). Most studies find that cluster mergers trigger star formation, although there has been some debate (Stroe et al., 2015, and references therein).

Neutral hydrogen (Hi ) is an ideal tool to study environmental processes as it is a very fragile component of galaxies (Haynes & Giovanelli, 1984; Giovanelli & Haynes, 1985). In fact, there is unimpeachable evidence that Hi can be disturbed and stripped in late-type galaxies as they pass through the dense ICM of massive galaxy clusters (e.g. Cayatte et al., 1990; Bravo-Alfaro et al., 2000, 2001; Poggianti & van Gorkom, 2001; Kenney, van Gorkom & Vollmer, 2004; Crowl et al., 2005; Chung et al., 2007, 2009; Abramson et al., 2011; Scott et al., 2010, 2012; Gavazzi et al., 2013) and even groups (e.g. Catinella et al., 2013; Hess & Wilcots, 2013). Simulations suggest that this happens via RPS (see Roediger, 2009, for a review) and/or gravitational interactions (e.g. Vollmer et al., 2001; Vollmer, 2003; Kapferer et al., 2009; Tonnesen & Bryan, 2009). Furthermore, Hi is also known to be sensitive to galaxy-galaxy interactions (see e.g. Yun, Ho & Lo, 1994), common in group environment. Therefore, Hi is the perfect tool to study the effect of environment on galaxies across cosmic time.

Due to technical limitations Hi has not been studied beyond the low-redshift Universe. The Blind Ultra Deep Hi Environmental Survey (BUDHIES; Verheijen et al., 2007, 2010; Jaffé et al., 2013), observed, for the first time, Hi in and around two clusters at with the Westerbork Synthesis Radio Telescope (WSRT) detecting 160 galaxies. Recently, this redshift regime has started to be explored by other Hi surveys. The HIGHz Arecibo survey observed 39 galaxies at in pointed observations (Catinella & Cortese, 2015). CHILES, the COSMOS Hi Large Extragalactic Survey is currently observing a a arcmin field with the Very Large Array (VLA) covering the range . A pilot for this survey covering detected 33 galaxies (Fernández et al., 2013).

In Jaffé et al. (2012) we studied the Hi content in Abell 2192 (A2192_1 hereafter, following the notation of Jaffé et al., 2013), a cluster with mass, and a significant amount of substructure. We found that this intermediate-size cluster is in the process of assembling, and that the incidence of Hi -detections significantly correlates with its substructure. At large clustercentric radii (2-3 virial radii), where the infalling substructures were found, many galaxies are Hi -rich, while at the core of the forming cluster, none of the galaxies are Hi -detected. The reduction of Hi starts to become significant in low-mass groups (with masses as low as ), so that by the time the group galaxies fall into the cluster their Hi reservoir is significantly reduced.

In Jaffé et al. (2015) we concentrated on the more massive cluster Abell 963 (A963_1 hereafter, following previous papers) to study the effect of RPS by the ICM. We compared the location of Hi -detected cluster galaxies in velocity versus position phase-space diagrams with the distribution of dark matter haloes in zoom-in cosmological simulations, and found strong evidence for Hi gas stripping (down to the detection limit of BUDHIES) on first infall into the cluster. In particular, we defined regions in phase-space where galaxies are expected to be experiencing Hi stripping via ram-pressure. We also found however that ram-pressure cannot be the only mechanism at play, since at the outskirts (infall region) of the cluster there is a significant fraction of passive Hi -poor galaxies. We speculated that these galaxies had to have been quenched prior to entering the cluster.

Figure 1: Coverage of the different observations around the central galaxy of A963_1. The field-of-view of the WSRT Hi observations is highlighted with a thick blue box and labelled accordingly, along with the multiwavelength photometry: INT ( and ; red square), Spitzer MIPS (cyan) and IRAC (green regions), XMM-Newton (X-rays; cyan tilted square), Galex ( and ; magenta circle), Subaru ( and ; thin blue square), and UKIRT (dashed black square). In addition, the black dots show the distribution of all the spectroscopic members of A963_1 (from WHT and MMT) and the cyan contours the extended X-ray emission detected by XMM. The horizontal black line at the bottom-left corner of the plot indicates the scale.

In this paper, we combine BUDHIES data around A963_1 with hundreds of new redshifts, as well as optical and X-ray imaging data from the Local Cluster Substructure Survey (LoCuSS) to pursue an in-depth investigation of the combined effect of RPS by the ICM with group pre- and post-processing in the removal of Hi and subsequent star formation enhancement/quenching in cluster galaxies.

The paper is organized as follows: in Section 2 we present the multiwavelength data available in the field of A963_1. In Section 3.1 we study the degree of substructures within the cluster and identify infalling groups. In Sections 3.2 and 3.3 we explore the properties of the cluster members, such as galaxy colours and morphologies as a function of environment (i.e. position within the cluster, groups, and phase-space distribution). In Section 3.4 we study the Hi content of galaxies as a function of environment through Hi deficiencies and Hi stacking. We complement our results in Section 3.5 by performing stacking of optical spectra in different regions of phase-space. Finally, in Section 4 we discuss and summarize our results.

Throughout this paper we use AB magnitudes, and assume a concordance CDM cosmology with 0.3, 0.7, and H70 km s Mpc.

2 Data

2.1 Budhies

This study is based on data from BUDHIES, an ultra-deep Hi survey of galaxies in and around two clusters at . With effective volume depth of 328 Mpc and a coverage on the sky of 1212 Mpc for each cluster, BUDHIES covers a large range of environments. The ultra-deep Hi observations, taken with the WSRT (Verheijen et al., 2007, 2010), yielded 126 Hi detections in the volume of A963_1 and 36 in A2192_1 (Verheijen et al., in preparation). The synthesized beam has a size of arcsec (corresponding to kpc at ), and the velocity resolution is 19 km s, allowing the Hi emission lines to be well resolved. Four Hi detection criteria were applied to data cubes of different resolution, namely, that the Hi signal was above the rms noise in a number of consecutive channels: ; ; ; or , where is the velocity resolution element. The Hi mass () detection limit was typically , although it varies with distance from the pointing centre (due to primary beam attenuation) and frequency.

BUDHIES also has - and images from the Isaac Newton Telescope (INT), GALEX and , Spitzer IRAC and MIPS, UKIRT in the near-infrared, and Herschel SPIRE and PACS for a limited region. SDSS photometry is also available for the galaxies in the field, as well as spectroscopy, although this is scarce at . To alleviate the lack of redshifts in the surveyed volumes, we pursued a spectroscopic campaign with WIYN and the AF2 on the 4.2m William Herschel Telescope (WHT), from which we have assessed cluster memberships, measured the equivalent width of the [O] line when present, and characterized the environment in the surveyed volumes (Jaffé et al., 2013). The spectroscopic campaign had a spectral coverage of 3900–6900Å and limited completeness (% down to in A963_1, the least complete cluster in BUDHIES). To increase the spectroscopic completeness in and around A963_1, in this paper we have included hundreds of new redshifts from LoCuSS (private communication) and Hwang et al. (2014). As part of the LoCuSS cluster sample A963_1 was targeted with Hectospec, a multi-object spectrograph on the 6.5m MMT telescope. The observations provide a wide wavelength range (3650–9200Å) at 6.2Å resolution. Targets were selected from UKIRT/WFCAM -band images down to , with colours consistent with known cluster members over a field of view (Haines et al., 2009), and additional preference was given to those targets detected at 24m (Haines et al., 2013). Hwang et al. (2014) also obtained more redshifts in the volume of A963_1 using MMT/Hectospec. They targeted galaxies brighter than in the field of the cluster, irrespective of colour. The inclusion of all new redshifts have increased our spectroscopic completeness significantly (to above at ), as shown in Fig. 19 of Appendix A.

In addition to the spectroscopic observations, LoCuSS also observed the inner () region of A963_1 with XMM-Newton (Martino et al., 2014) and Subaru/Suprime-Cam imaging in and (with a seeing of  arcsec, corresponding to 2.5 kpc at the cluster redshift; see Okabe et al., 2010). We make use of these data sets in Section 2.2 to study cluster substructure, and in Section 3.3 to characterize the morphologies of the cluster galaxies.

Fig. 1 shows the coverage of the multiwavelength observations in the field of A963_1 used in this paper, along with the spatial distribution of the cluster members.

Figure 2: CMD of all the cluster members of A963_1. The horizontal line divides galaxies into “passive” ( or undetected in , plotted at at the top of the plot for reference; red diamonds), and blue ( and ; light and dark blue circles respectively). The vertical dashed line shows the spectroscopic completeness limit ). We use different symbol shapes, colours and sizes in this plot to easily identify each sample (e.g. passive versus star-forming or bright versus faint) in Fig. 8. The number of galaxies in each CMD region are labelled, and Hi -detections are marked with open grey squares.
Figure 3: SFR (left) and sSFR(right) versus stellar mass for all cluster galaxies, colour-coded by colour. Galaxies not detected in NUV were plotted with the reddest colour of the scale. In both panels, Hi -detected galaxies are highlighted with open cyan squares. The star-forming “sequence” at from Speagle et al. (2014) is plotted for reference with a blue dashed line (and dotted lines for scatter). The solid line (in both panels) corresponds to sSFR, and roughly separates blue Hi -detected galaxies from redder undetected ones.

2.2 A963_1

This paper focuses on A963_1 (RA10h 27m 36.97s; DEC-9d 18m 11.34s), the most massive cluster of BUDHIES. The improved completeness, and the additional 426 galaxies with a redshift allowed us to better characterize the mass of the cluster and the substructure within it (Section 3.1). From these new estimates we also re-assigned cluster memberships considering both their position and velocity, as can be seen in the right-hand side of Fig. 7. The figure shows a projected clustercentric distance versus redshift phase-space diagram in an extended field around the BUDHIES volume. The cluster members are coloured in red, and were chosen to follow the “trumpet” shape characterizing the cluster. The revised mean redshift for the cluster is .

From a sample of 140 spectroscopic cluster members (), in Jaffé et al. (2013) we estimated the cluster’s mass to be ( km s). With the increased spectroscopic sample we now have 566 cluster members. Using the same magnitude-limit as in Jaffé et al. (2013) we re-computed the mass of the cluster to be: ( km s). This is still somewhat higher than the most recent X-ray (XMM-Newton) estimate, 1, but closer to the latest weak lensing mass estimate ( Okabe & Smith, 2015). We note however that our mass estimate (inferred from ) assumes that the cluster is relaxed and in equilibrium. The discrepancy between the measurements could thus be caused by the presence of substructure.

A963_1 has been classified as regular by previous X-ray studies and lensing analysis (Smith et al., 2005; Chon, Böhringer & Smith, 2012). However, in Jaffé et al. (2013) we saw weak evidence for substructure, based on a kinematic analysis. In Section 3.1 we will indeed show that there are several distinct groups feeding the cluster. If we exclude galaxies belonging to subgroups, the dynamical mass drops down to ( km s ), which coincides with the weak lensing mass estimate.

2.3 Colours, stellar masses and star formation rates

Fig. 2 shows the versus colour-magnitude diagram (CMD) of the spectroscopically confirmed members of A963_1. Because traces the ratio of young to old stars, red colour can be associated with passive galaxies. Galaxies not detected in are considered “red” and are plotted in Fig. 2 with arbitrarily red colour (=7.8, redder than the reddest cluster galaxy) for reference. The vast majority of the undetected galaxies lie in the versus red-sequence and do not show emission in their optical spectra. We therefore refer to all galaxies in the red () bin of the CMD, including galaxies not-detected in as “passive” hereafter. In the figure, Hi -detections are also highlighted (squares) and the number of galaxies in the different samples is indicated. Most of the Hi detections are in the bluest part of the diagram, as expected. In particular, it is visible how the Hi detections drop with colour so that almost no galaxy redder than is detected in HiȮur findings are consistent with the work of Cortese et al. (2011), Catinella et al. (2013) and Brown et al. (2015), that show a gradual decrease in Hi content of galaxies with colour. The only exception is one passive galaxy detected in HI, which is an elliptical galaxy with a small nearby galaxy (in projection) that could be responsible for the Hi emission. The CMD of Fig. 2 shows a prominent population of dwarf blue cluster galaxies that extend 1-2 mag fainter than red dwarfs. This reflects the BUDHIES spectroscopic target selection criteria that allocated fibers to galaxies brighter than with the exception of Hi -detected galaxies that could be as faint as (see Fig. 2 in Jaffé et al., 2013). Also, for the faint population, spectroscopic redshifts were easier to obtain when emission-lines were present. For however, our sample is fairly unbiased and complete (see Figure 19).

It is known that the m emission reflects the obscured star formation (Kennicutt et al., 2007), while the observed UV emission traces the un-obscured component (Calzetti et al., 2007). We thus computed total star formation rates (SFR) from a combination of GALEX FUV luminosity and infrared luminosity () derived from MIPS. First, we derive from the MIPS flux and redshift using the calibration of Rieke et al. (2009). For the un-obscured component, we calculate after applying a -correction to the FUV flux. Then, for galaxies detected in both FUV and MIPS , we calculate total SFR following Equation 9 of Murphy et al. (2011). For galaxies detected only in FUV, we calculate SFR using Equation 3 of Murphy et al. (2011), which accounts for only the un-obscured component. And for any galaxies detected in MIPS but not in FUV, we use Equation 4 of Murphy et al. (2011). We note that the computed SFRs are dominated by the infrared component, and thus they trace the obscured component of star formation activity well. Specific star formation rates (sSFRs) were computed dividing SFR by stellar mass (). We measured from the available and INT photometry, using the method described in Zibetti, Charlot & Rix (2009), and a Kroupa initial mass function.

Fig. 3 shows the SFR (left) and sSFR (right) as a function of stellar mass for the cluster members, colour-coded by their colours. The black solid line indicates the sSFR boundary between most Hi -detections (that also roughly corresponds to galaxies with ) and non-detections. Such boundary corresponds to , as shown in the right-hand panel. As evidenced by Fig. 3, is sensitive to recent star formation and it can thus be used as a proxy for specific star formation rate (sSFR, see Salim et al., 2005, 2007; Schiminovich et al., 2007).

3 Results

3.1 Substructure within the cluster

X-rays

Previous Chandra and HST observations of the central few hundred kpc region of the cluster reveal a regular and concentrated X-ray emission, around a cD galaxy (Fig. 6 of Smith et al., 2005). Recent analysis of XMM-Newton observations by Haines et al. (in preparation) has revealed additional substructures surrounding the centre, as shown by the cyan contours of Fig. 1 (see also Fig. B.6. in Zhang et al., 2007). Three of them are at the redshift of the cluster, at a distance of from the cluster centre. These cluster substructures, labelled ‘A’, ‘B’ and ‘C’ and listed in Table 1, have velocity distributions consistent with that of groups, as shown in the top panel of Fig. 5. The most massive group detected in X-rays is group C, that could even be considered a small cluster from its mass. In fact, it has a distinctively bright galaxy at its centre. This group is however over an order of magnitude less massive than A963_1. Groups B and A do not display a distinctive bright central galaxy. Group B in fact hosts many merging galaxies, and the low-mass group A shows a prominent population of blue, unperturbed, late-type galaxies.

To assign group membership we selected galaxies within a radius equal to the of the groups and line-of-sight velocities consistent with the bulk of the other group galaxies, as shown in Fig. 5. was computed for each group from , estimated using the X-ray luminosity- relation of Leauthaud et al. (2010, see Haines et al. in preparation). The of the X-ray groups are plotted as open circles in Fig. 4, and labelled ‘A’, ‘B’ and ‘C’.

The Dressler-Shectman test

An alternative way to identify substructures within a galaxy cluster is through the Dressler-Shectman (DS; Dressler & Shectman, 1988) test, that compares the local (10 nearest neighbour) velocity and velocity dispersion for each galaxy with the global (cluster) values.

In short, the method consists of measuring individual galaxy deviations, (Eq. 6 in Jaffé et al., 2013) for the cluster members, to then quantify the degree of substructure in two ways:

  1. The ‘critical value’ method evaluates .

  2. The probabilities () method instead computes -values by comparing the -value to ‘shuffled’ -values is then calculated as , so that a small value indicates the presence of substructure.

The details of the DS test can be found in Jaffé et al. (2013), where we carried out this analysis in all the BUDHIES clusters with the limited spectroscopic sample available at the time, finding marginal signs of substructure in A963_1. In particular, we found weak evidence for the presence of 4 groups within the cluster (A963_1a, A963_1b A963_1c and A963_1d). It is known however that the mean -deviation substructure index is very sensitive to completeness (see Fig. 3 in Ragone-Figueroa & Plionis, 2007). For this reason, we re-assess here the presence of substructures using the increased spectroscopic sample.

The results of the DS test on the new data reveals a significant degree of substructures (/ and ). In fact, several substructures surrounding the cluster core can be seen in Fig. 4, that highlights with darker colours the galaxies with the largest deviations (i.e. the highest values). Interestingly, the X-ray groups (enclosed by open circles, and labelled ‘A’, ‘B’, and ‘C’) do not show the strongest deviations from the local and global velocity distribution of the cluster. The strongest deviations () instead define several additional substructures or galaxy overdensities that are not identified from the X-rays. The four most prominent “DS groups” are labelled with the Greek letters , , and for distinction from the X-ray groups, and listed in Table 1.

Group (km s) () No. members Method
A 383 10 (7) X
B 234 8 (6) X
C 753 34 (16) X
338 17 (9) DS
1129 22 (12) DS
108 6 (5) DS
197 9 (4) DS
Table 1: Properties of the subgroups found in A963_1. The columns are: Group given name, velocity dispersion (), mass (), number of members (number above the spectroscopic completeness limit), and method of identification: (X for X-rays or DS for dynamical). The observed velocity dispersion () was estimated for galaxies above the spectroscopic completeness limit. was estimated from the X-ray data for the X-ray groups A, B, and C (top). For the DS groups (bottom) it was derived from the measured . We note that groups and have very wide and narrow velocity distributions respectively, so their dynamical mass estimates were not attempted.

The DS groups and correspond to the groups A963_1d and A963_1c, identified in Jaffé et al. (2013, despite the spectroscopic incompleteness). Our revised analysis confirms their reality and allows us to assign secure group memberships to galaxies in these two groups. We also found that the group A963_1b in Jaffé et al. (2013, see Fig. 10) instead, corresponds to the X-ray group A, while group A963_1a is not picked up as a strong DS substructure in the new DS analysis, so we discard it as a real substructure.

With the exception of group , all the DS groups have a narrow velocity distribution as shown in the bottom panel of Fig. 5, confirming they are real associations. None of the DS groups are detected in X-rays, but a plausible explanation exists. The DS groups , and are located outside of the region covered by the XMM-Newton observations (see Fig. 1). Only the group is well inside the XMM-Newton region. It is possible however, that this group is being disrupted by the main cluster due to its vicinity to the cluster core and it is thus not detected in X-rays. This is consistent with its broad velocity distribution.

It has been shown by Cohn (2012) that it is easier to detect a substructure via the DS test for large subgroups, and that the likelihood of detection also depends on the line-of-sight direction. This could mean that the four substructures we found via the DS test in A963_1 are the larger subgroups in A963_1, but there could be more undetected ones. In fact, from Fig. 4 we can see signs of weaker DS structures that could be individual groups or part of disrupted groups, as will be discussed below.

Figure 4: The spatial distribution of all the spectroscopically confirmed members of A963_1 is plotted in solid circles, colour-coded by their deviations (see text for details). The darker the symbols the more deviated the galaxies are from the mean velocity distribution of the cluster. Open circles enclose the galaxy groups identified form the X-rays (Groups A, B, and C; pink), and the DS test (Groups , , , and ; violet), and the dashed circle corresponds to the cluster’s .
Figure 5: Velocity distribution of the galaxies in the main cluster (solid grey histogram), and those in the X-ray groups A, B, and C (left) and the DS groups , , , and (right), as labelled.
Figure 6: Projected phase-space diagram for all the cluster members (small grey crosses). Galaxies in the X-ray or DS groups are highlighted with bigger circles and diamonds respectively, and coloured as in Fig. 5. The ‘virialized’ and ‘recent infall’ zones defined in Jaffé et al. (2015) are indicated (see also description in Section 3.1).

Accretion history of the cluster

Figure 7: Left: Galaxies in the redshift range 0.184-0.224 in a  deg area around A963_1. Galaxies are colour-coded according to their velocities (centred on the cluster’s mean value), as indicated. The large open circle has a radius of ( at ), and roughly corresponds to the coverage of our spectroscopic campaign. Outside of this circle galaxies with a redshift from SDSS are plotted. Inside of the circle, smaller open circles demark the location of the identified X-ray and DS galaxy groups. The figure shows that A963_1 is being fed by at least two filaments (NE and SE), and possibly a third one (W). Galaxies in the filaments are enclosed by open diamonds, triangles and squares. Right: distance from the centre of A963_1 (in Mpc) versus redshift for all the galaxies with a redshift between 0.190 and 0.220 (small open circles) The red filled circles are the galaxies that are considered as cluster members. Their distribution in this diagram fill the typical “trumpet” shape of fairly virialized clusters in phase-space (see also bottom panel of Fig. 8). The bigger symbols correspond to the galaxies in the filaments identified in the left-hand panel, as labelled.

The presence of sub-groups within A963_1 is consistent with cosmological simulations (e.g. Cohn, 2012), which predict that sub-groups host a significant fraction of cluster galaxies, with larger sub-groups more likely to be in more massive clusters. The number of galaxies in groups however, drops significantly with time since infall into the cluster. This prediction suggests that the groups that we are able to detect inside of A963_1 have not been in the cluster for long. Fig. 6 supports this assumption. The horizontal axis of the plot shows the projected distance from the cluster centre, , normalized by (as computed from X-rays), while the vertical axis displays the peculiar line-of-sight velocity of each galaxy with respect to the cluster recessional velocity (), normalized by the velocity dispersion of the cluster ( km s; computed excluding group members). The figure also shows the ‘recent infall’ and ‘virialized’ regions, defined in Jaffé et al. (2015). In short, the ‘recent infall’ region represents the cluster outskirts, where galaxies that were most recently accreted into the cluster are more likely to be. In the ‘virialized’ region instead, % of the galaxies have been in the cluster for over Gyrs (roughly a cluster crossing time). Fig. 6 thus shows how most group galaxies clearly avoid the ‘virialized’ region of the cluster in phase-space, preferring the ‘recent infall’ region.

Although we assume a spherical morphology for the groups (circles in Fig. 4), many of them show a non-spherical spatial distribution, such as elongations or preferred directions that can extend beyond the drawn boundaries. It is possible that some galaxies in the vicinity of these groups were previously bound to them, but are now being ejected, as the group falls deeper into the cluster’s potential well. For example, there are a few galaxies with large values south of group C, west of group and south-west of group (see Fig. 4) that could have well been part of these groups in the past. Simulations indeed show that galaxy groups can be stretched out along the direction of infall before falling into the cluster (Vijayaraghavan & Ricker, 2013). This could be a consequence of the group falling in along a cosmological filament or the effect of the cluster’s tidal field, or a combination of both. The elongations of some of the groups (together with their phase-space location) thus suggest that (at least some of) the groups are on first infall and could be part of a larger filamentary structure.

Figure 8: Spatial (top) and phase-space (bottom) distribution of all the cluster galaxies, colour-coded by and sized by luminosity as in Fig. 2. The open grey squares correspond to the Hi -detected galaxies. Open circles in the top panel mark the groups defined in Fig. 4 and the dashed circle the of A963_1. Shaded areas in the bottom panel indicate several regions of interest such as the ‘recent infall’, ram-pressure ‘stripping’, and ‘virialized’ regions, discussed in the text.

In an attempt to map the large-scale structure (LSS) around A963_1, in the left-hand panel of Fig. 7 we plotted galaxies with SDSS redshifts in the surroundings of the cluster, and in the redshift range . Although SDSS spectroscopy is limited at this redshift, the plot shows evidence of at least two filamentary structures feeding the cluster from the North-East (NE) and the South-East (SE). A possible third filament can also be seen from the West (W) direction. These structures seem to be separated by void or empty regions, and in phase-space (right-hand side of Fig. 7) they show narrow velocity () dispersions at all , which supports the idea that they are indeed cosmic filaments. In fact their phase-space distribution is similar to the large-scale filaments found around Virgo (see Figure 3 in Lee, Kim & Rey, 2015) and Coma (see Figure 4 in Falco et al., 2014). It is possible that the groups (plotted as grey open circles in Fig. 7) were accreted into the cluster from such filamentary structure.

In addition to the groups currently being accreted into the cluster, the X-ray temperature map in the inner region of the cluster (made by AF, not shown here) further reveals a past merger/interaction. These events typically leave inhomogeneities in the temperature map at the cluster core, which are visible in A963_1.

3.2 Distribution of red, blue and Hi -rich galaxies

In the top panel of Fig. 8 we inspect the spatial distribution of all the cluster galaxies, highlighting their colour and Hi content in the same manner as in the CMD of Figure 2. It is evident that the cluster core, as well as groups host a significant population of passive, gas-poor galaxies. Passive galaxies also extend beyond the imposed boundaries () of the groups, most notably in Groups , C-B and , possibly due to group elongation caused by the cluster potential well, as discussed in Section 3.1.

To study more closely the connection between the passive galaxies and environment, in the top panel of Fig. 9 we evaluate the incidence of passive galaxies as function of distance from the cluster centre. The plot clearly shows how passive galaxies (which are also typically the most massive ones) dominate the cluster core, and become more scarce with increasing , as expected. The bottom panel of Fig. 9 further shows, for galaxies in groups with robust mass estimates, how the fraction of passive group galaxies increases with increasing group mass, at the expense of Hi -detections. While the fraction of passive galaxies in the lower mass groups coincides with that of the infall region at the largest clustercentric distances (), at higher group mass, it jumps to . The observed jump in the passive fraction (as well as a comparable decrease of Hi -detections) from low to intermediate-mass groups confirm the existence of a strong connection between the gas content, star formation activity and host halo mass, and suggest that environmental star formation quenching starts being significant in groups more massive than . We note that for the plot in the bottom panel of Fig. 9 we did not make any magnitude cut due to the low number of galaxies. We note however that the stellar mass distribution of the group galaxies are very similar, ranging from to and with mean value .

Figure 9: Top: The fraction of passive cluster galaxies brighter than as a function of distance from the cluster centre. Bottom: fraction of passive (red circles) and Hi -detected ( cyan squares) galaxies in the groups A, B and C, as a function of group mass (groups of similar mass were added to increase the number of galaxies and make the Poisson errors smaller). In this case, no magnitude cut was used.

The bottom panel of Fig. 8 shows the distribution colour and Hi in projected position versus velocity phase-space. A similar plot was already shown in Jaffé et al. (2015) but it is presented here again with the much larger spectroscopic sample. A few areas, defined in Jaffé et al. (2015), are shown. The ‘recent infall’ and ‘virialized’ areas reflect the orbital histories of cluster galaxies, as already discussed in Section 3.1. The ‘stripping’ zone, also defined in Jaffé et al. (2015) is the outcome of embedding Gunn & Gott (1972)’s prescription of RPS, together with a model of A963_1’s ICM, into dark matter cosmological simulations. Our model predicted gas-rich galaxies to be located in the ‘recent infall’ region of phase-space, avoiding the ‘stripping’ and ‘virialized’ zones. The observed distribution of Hi -detections is in striking agreement with the model (see Figure 6 in Jaffé et al., 2015). The idea is that as galaxies fall into the cluster from high clustercentric distances (‘recent infall’ region in phase-space), they can gain velocity as they approach the cluster core, which leads them into the ‘stripping’ zone. After their first pass, they will oscillate in phase-space slowly settling in the potential well of the cluster, building up the ‘virialized’ region.

Fig. 8 further shows that passive galaxies dominate both the central part of the cluster, and the ‘virialized’ part of phase-space, suggesting that many of them have been in the cluster for a long time ( Gyr). However, not all passive galaxies are in the virialized region. In fact, a significant population of passive galaxies can be seen in the infall region of the cluster. It is possible that some of these (especially the most massive ones) have been ‘mass quenched’ (by internal processes such as AGN feedback; see Peng et al., 2010), but additional ‘environmental quenching’ might also be at play, especially for lower mass galaxies. As discussed in Jaffé et al. (2015), the abundance of passive galaxies in the cluster outskirts cannot easily be explained by RPS alone (or backsplashing galaxies). However, a more plausible explanation for this effect is that many of these galaxies started to be quenched in a previous environment. One likely possibility, given the results shown in Fig. 9, is group ‘pre-processing’. In fact, % of the passive galaxies in the infall region of the cluster belong to the groups identified in Section 3.1. This is a lower limit however, as we can only identify the most obvious group members and many group galaxies have probably mixed with the cluster population already. In addition, the limited area of the XMM observations does not allow the identification of X-ray groups at large clustercentric distances.

Figure 10: Example Subaru -band images (arcsecarcsec) of the peculiar galaxies identified in A963_1. Most of the peculiar galaxies show signs of interactions with another galaxy, but a few could be stripped by the ICM. The three most convincing ‘jellyfish’ candidates are shown in the bottom row.

3.3 Morphological peculiarities

In the quest to detect signs of the different physical processes driving galaxy evolution in this actively accreting cluster, two of the authors of this paper (YJ and MMC) visually inspected the morphologies of the galaxies in the central  arcmin region of the cluster using the Subaru and images. We found many galaxies with peculiar morphologies. Most of them show signs of gravitational interactions with another galaxy, such as ‘bridges’ of stars, multiple nuclei, or irregular morphologies, although the origin of the perturbations is not always clear. Some examples are shown in Fig. 10. A few of the peculiar galaxies exhibit one-sided ‘tails’ or ‘jellyfish’ morphologies, typically associated with galaxy-ICM interactions. The most convincing cases are shown in the bottom row of Fig. 10. In our analysis, we treat all peculiar galaxies as ‘interacting’, not attempting to identify the cause of their disturbances.

The top panel of Fig. 11 shows the distribution of the cluster galaxies in the sky. Peculiar galaxies (solid symbols) clearly trace the structures of A963_1: they cluster near the centre (extending to the east), and they are also common in and around the groups. When selecting only group galaxies, we find that peculiar galaxies are more common in the lower-mass groups A, B, and (% on average, for galaxies with ) in comparison with Group C (%). Moreover, when inspecting the colours of the peculiar group galaxies, we find that in grup C, they are all passive, while in the lower-mass groups they are a mixed population.

Figure 11: The spatial (top), and SFR versus (bottom) distribution of all cluster members (no magnitude cut; open circles) are compared with galaxies with peculiar morphology (filled red diamonds). Note that morphologies are only available in the area covered by Subaru (plotted area in top panel). The bigger open circles in the top panel delimit the different groups as in Fig. 8. The dashed circle in the centre corresponds to the of the cluster. The lines in the bottom panel are as in Fig. 3.

The SFR-stellar mass diagram of the bottom panel of Fig. 11, further shows that the vast majority () of peculiar galaxies are actively forming stars.

The abundance of peculiar star-forming galaxies in groups is not a surprise, since most peculiar galaxies in our sample are clear mergers, a phenomena more likely to occur in galaxy systems of low velocity dispersion. The large number of peculiar galaxies in the central part of the cluster however, is harder to explain.

Using cosmological simulations Vijayaraghavan & Ricker (2013) show that during a group’s pericentric passage within a cluster, galaxy-galaxy collisions are enhanced. In addition, merger shock enhances the ram-pressure on group and cluster galaxies and an increase in local density during the merger leads to greater galactic tidal truncation. They call the combination of these effects ‘group post-processing’.

It is possible that the high fraction of peculiar (and/or star-forming) galaxies near the cluster core arise via post-processing by the numerous groups falling into A963_1 and the very high galaxy density. This is supported by the high degree of substructure in this cluster, and the presence of not only mergers, but also ram-pressure stripped galaxies.

Figure 12: Hi -deficiencies as a function of distance from the cluster centre for blue galaxies in A963_1 (blue circles) and the Coma cluster (red triangles; Boselli & Gavazzi, 2006). For A963_1 galaxies, Hi deficiencies for Hi -detections were computed using the average of all morphological types, with error bars corresponding to the Hi -deficiencies of spirals from Sa to types later than Sc. For non detections (open circles) we plot lower limits on the Hi deficiency. All points were corrected for primary beam attenuation. The solid and dashed lines trace the mean Hi deficiency values for the galaxies in A963_1 and Coma respectively, and the shaded regions correspond to uncertainties.
Figure 13: Hi stacking of passive (top) and blue (bottom) cluster galaxies (excluding group members) inside (left) and outside (right) . These stacks exclude group galaxies. In each panel, the thin grey lines correspond to the eight reference offset stacks. The number of galaxies that went into each stack () and number of direct Hi -detections () are indicated in the top-left corner of each panel. The horizontal bar underneath the spectra indicates the rest-frame velocity range used in the flux measurement.

3.4 The Hi content of galaxies as a function of environment

Hi deficiencies

We first assess the Hi content of galaxies in A963_1 by computing Hi deficiencies. From a sample of 324 (mostly spiral) isolated galaxies Haynes & Giovanelli (1984) defined Hi deficiencies as the logarithmic difference between the observed Hi mass and the expected value in isolated objects of similar morphological type and linear size.

Unfortunately, the visual classification of A963_1 galaxies from available ground-based images (see Section 3.3) could only allow the identification of interacting galaxies but not a detail assessment of T-type morphologies. For this reason, we assumed all blue galaxies are spirals and all passives ones are elliptical or S0, and computed Hi deficiencies for all the blue () galaxies in A963_1. Galaxy sizes were determined from the 50% light radii in B-band and the angular distance to A963_1. As a comparison sample for the computation of the deficiencies we used the local field galaxies of Haynes & Giovanelli (1984). For galaxies not directly detected in Hi we compute the lower limit of Hi -deficiency correcting for primary beam attenuation.

Fig. 12 shows the computed Hi -deficiencies as a function of clustercentric distance. The plotted error bars correspond to the Hi -deficiencies of Sa spirals to types later than Sc. In addition, we plot Coma galaxies from Boselli & Gavazzi (2006) as local reference. Note that A963_1 and Coma are similarly massive () clusters. Despite the scatter, the figure clearly shows how galaxies in both clusters become increasingly more deficient as they get closer to the cluster centre. The radial decay of Hi deficiencies in A963_1 agrees within errors with the trend observed in Coma, although close to the cluster core, the z=0.2 galaxies are a factor of 2 less deficient than the local ones.

Hi stacking in and around the cluster

As in Verheijen et al. (2007), we also probe deeper into the gas content of the cluster galaxies by stacking the Hi spectra of galaxies with known optical redshifts in bins of colour and environment. To do this, we extract the spectra around the expected location of the 21 cm line in each individual galaxy using its Hi position and redshift when available, or the optical information otherwise. We do this over an area of roughly the size of the FWHM of the synthesized beam, and 401 channels around the galaxy, corresponding to a window of km s in velocity. This generous selection allows us to quantify the noise around the 21 cm line robustly. We then add all the spectra in a given sample of galaxies to compute an average flux, that is later converted into a assuming all galaxies are at the mean distance of the cluster (): , where  Mpc (the distance to A963_1), = , the flux density in channel , and km s, the rest-frame channel width.

To estimate the errors on the Hi fluxes (and consequently ), for each stack, eight reference stacks at different offset positions were made. We then calculated the variance in the eight reference fluxes, integrated over the same velocity windows as the on-target stack. The values derived from the stacks presented in this section, and their associate errors are listed in Table 3.

Figure 14: Hi -stacking of red (top row) and blue (bottom) galaxies in groups, separated in three bins of group mass: low (left, Groups A and ), intermediate (middle, groups B and ), and high mass (right, group C). No magnitude cut was applied. Lines and labels are as in Fig. 13.
Figure 15: Hi -stacking of passive (top) and blue (bottom) galaxies in the ‘virialized’ (left), “ram-pressure” (middle) and ‘recent infall’ (right) regions of phase-space, as defined in the bottom panel of Fig. 8. Lines and labels are as in Fig. 13. Only galaxies above the spectroscopic completeness limit () were considered.
Figure 16: Results from the Hi stacking in different environments. (A) The average Hi mass of passive and blue galaxies inside and outside is compared for all cluster members with , excluding galaxies in groups. We plot the mean value in each bin against the average . (B) The average in groups of similar mass (see Fig. 14) is plotted against average group mass for passive and blue galaxies. (C) Average in different regions of phase-space for passive and blue () galaxies. In all panels solid lines and solid circles are used when there was a measured whereas dashed lines and inverse triangles are used when only upper limits could be measured.

We demonstrate the power of stacking in Appendix B. The top panel of Fig. 21 shows the results of the Hi stacking in different bins of colour and magnitude. It is noticeable how the 21 cm emission decreases dramatically with redder colour and fainter magnitudes, as expected. To further test whether the direct Hi detections dominate the stacks, we performed the same stacks in colour and magnitude bins, but excluding the direct Hi detections (bottom panel of Fig. 21). The trends seen before stand, although naturally, the total Hi masses in this case decrease. However, the bulk of the galaxies not directly detected in Hi do have significant amounts of gas, in many cases not much below the threshold for the detected ones. The results of this exercise are summarized in Fig. 20.

In the following, we present several Hi stacks as a function of environment, considering Hi -detected and non-detected galaxies combined. We use bins of colour, which is a proxy for star formation activity, and only consider galaxies above our spectroscopic completeness limit () unless otherwise stated. We prefer to bin in observable quantities, to avoid contamination from poor modelling.

Hi stacking in the substructures

In Fig. 13 we compare the Hi content of cluster galaxies in the innermost region of the cluster () with those in the outskirts of the cluster (). We do this using the spectroscopically complete sample in the colour bins defined in Fig. 2 (i.e. blue and red/passive), and excluding galaxies in groups. The stacks show a mild detection of Hi emission ()in the passive galaxies residing in the outskirts of the cluster, and no sign of emission () in the innermost part. In contrast, when stacking blue cluster galaxies a prominent emission is seen. At , , while in the cluster core is reduced by . The average fluxes and of the stacks are listed in Table 3 and the results are summarized in panel A of Fig. 16, where the solid blue line shows the reduction of Hi towards the cluster centre.

In Section 3.1 we identified seven groups falling into the cluster, that show an increasing fraction of passive galaxies with increasing group mass at the expense of Hi -detections (Fig. 9). To further investigate the Hi content of the groups, in Fig. 14, we stack the Hi emission of group galaxies in bins of colour and group mass. We only considered groups with reliable mass estimates, and merged together galaxies in groups of similar mass: A and (average mass of ); B and (), and C (). Due to the reduced number of galaxies in these stacks we did not limit the group samples by magnitude.

Passive group galaxies (top row in Fig. 14) barely host any HiȮnly passive galaxies in the lowest mass groups have a detection (). Blue group galaxies instead have larger Hi gas reservoirs. This can be seen in panel B of Fig. 16. The Hi gas masses of blue galaxies in the two lowest mass groups () are consistent with those of non-group blue galaxies in the infall region of the cluster (left-most blue point in panel C. The Hi masses of blue group galaxies progressively declines with group mass, being a factor 3 lower in the highest mass group ().

Hi stacking in phase-space bins

Fig. 15 shows the Hi stacks of passive and blue galaxies in the ‘recent infall’, (ram-pressure) ‘stripping’ and ‘virialized’ regions of phase-space. These regions (in that order), roughly trace the journey of a galaxy entering a cluster (i.e. its orbital history). For these stacks we used all cluster members in the magnitude-limited sample () not associated with groups. The most significant Hi detection is that of blue galaxies in the infall region of the cluster (bottom-right panel; ). The middle and left-hand panels in the bottom row show how blue galaxies lose their Hi reservoirs as they move from the infall region of the cluster inwards. In fact, of their Hi gas is removed on first infall as they cross the RPS region of phase-space.

For passive galaxies, we only detect Hi emission in the infall region of the cluster. This is consistent with (Serra et al., 2012) that found that % of all early-type galaxies outside the Virgo cluster have Hi (with ), while in the cluster centre, only % of them show Hi emission. Moreover, the mild Hi emission we measure in passive infalling galaxies () is not far from the local Hi mass function’s , which Serra et al. (2012) measured to be in Virgo.

The small amounts of Hi gas that passive galaxies may hold before entering the cluster is also lost via ram-pressure to the point of no-detection in the ‘stripping’ or ‘virialized’ region of phase-space ( and respectively). The outcome of these stacks is summarized in panel C of Fig. 16, where it is evident that, although passive galaxies in the infall region of the cluster have Hi emission, they have 80% less Hi gas than blue ones. In other regions of phase-space passive galaxies don’t show Hi emission. This suggests that they could also experience RPS as they move in phase-space into the ‘stripping’ and ‘virialized’ regions.

3.5 Star-formation activity as a function of environment: stacking optical spectra

To understand the effect of gas stripping on the star formation activity of the cluster galaxies, in this section we complement the Hi stacking presented in Section 3.4 with stacking of the optical spectra. For homogeneity, we only used the sub-sample of galaxies observed with the MMT, shifting all the spectra to rest-frame, and normalizing the fluxes so that of the continuum between and overlap on average.

Fig. 17 shows the optical stacks for galaxies in groups, binned by galaxy colour and group mass, as in Fig. 14. Again, we did not make magnitude cuts due to the low number of galaxies. Regardless of the group mass, passive group galaxies (top) show spectra dominated by old stellar populations, with strong absorption features. Conversely, the stacks of blue group galaxies (bottom) display some differences with varying group mass. Most notably, the stacked spectra of the blue galaxies in the lower mass groups (A) show much stronger emission lines (EW(), see Table 2) and a blue continuum, than the blue galaxies in the intermediate-mass groups B or in the most massive group C.

To test the strength of the stacking results presented in the bottom panel of Fig. 17 (i.e. stacking of blue group galaxies, see also Table 2), we stacked five randomly selected blue galaxies from groups A and obtained an EW() equal or lower than that measured in groups B only of the time. Moreover, we could never reproduce the EW() of groups A from the members of group B as non of their individual EW() as as high as .

These results indicate that low-mass groups still hold a significant population of star-forming galaxies, while groups more massive than already show signs of star formation quenching.

Figure 17: Optical stacks of passive (top) and blue (bottom) galaxies inside low (Groups ; blue) intermediate (Groups ; green) and high mass (Group ; red) groups, as labelled. The number of galaxies in each stack is indicated inside the parenthesis in the top-right corner. Main spectral features are also indicated.
Figure 18: Optical stacks of passive (top) and blue (bottom) galaxies brighter than in different regions of phase-space: the ‘recent infall’ (blue curves), ‘stripping’ (green) and ‘virialized’ regions (red), as labelled. As in Fig. 17, the main absorption and emission-lines are marked, and the number of galaxies in each stack is indicated.
Optical stacking sample No. galaxies EW()
By group mass:
Groups A 6
Groups B 5
Group C 6
By location in phase-space:
Recent infalls 54
Stripping 8
Virialized 5
Table 2: This table lists the results from the optical spectra stacks performed to blue (star-forming) galaxies in different regions of the cluster, namely groups of different mass and location in phase-space, as indicated in the fist column. The second column indicates the number of galaxies used in each stack and the third the measured EW() with its bootstrapping error.

In Fig. 18 we show the stacked spectra of galaxies in the ‘recent infall’ (blue spectrum), ‘stripping’ (green) and ‘virialized’ (red) regions of phase-space. As for the Hi stacks of Fig. 15, we used all cluster members in the magnitude-limited sample () not associated with groups. When separating red passive galaxies from blue star-forming ones, we find that, the passive galaxies in the different phase-space bins have a very similar appearance, typical of old stellar populations. Star-forming galaxies display a much larger variation. The stacked spectra of star-forming galaxies in the ‘stripping’ zone shows a significantly bluer continuum and mildly stronger emission lines (EW(H)=21.5 ) than the stack of galaxies in the ‘recent infall’ region (although within errors, see Table 2). This difference is bigger for the ‘virialized’ region (EW(H)=9.6 ).

We computed the likelihood of getting a stacked spectrum with an EW(H) as low or lower than the one measured in the ‘virialized’ region by stacking at random, five blue galaxies from the ‘recent infall’ region. This exercise yielded a probability of of producing a stack with EW(H). In addition, we tested the reality of the strong blue continuum seen in the blue galaxies in the ‘stripping’ zone by looking at their individual spectra. We find that show this feature, while only 10 out of 54 () blue galaxies in the ‘recent infall’ region do.

Our results suggest that there is a brief starburst associated to the first passage through the RPS zone in phase-space, that leads to a latter quenching of their star formation. Such an enhancement in the star formation of galaxies falling into clusters has been reported in the literature (e.g. Gavazzi & Jaffe, 1985; Porter et al., 2008; Mahajan, Raychaudhury & Pimbblet, 2012).

Finally, although not shown here, we also stacked the optical spectra of star-forming galaxies (, ) inside and outside the inner part of the cluster (), split into those individually detected in Hi and those with no Hi detection. In the cluster outskirts, the stacked spectra of the Hi -detected and Hi non-detections have indistinguishable continua and emission line strengths (both have EW()). In the inner region of the cluster () the stacked spectra of the three star-forming galaxies detected in Hi still shows similarly strong emission (EW()) and a blue continuum. However, the stack of the remaining 21 star-forming cluster galaxies not detected individually in Hi shows a redder continuum and weaker lines (EW()) a factor of 2 weaker than their counterparts in the cluster outskirts. This suggests that the bulk of star-forming galaxies within are currently in the process of having their star formation quenched, indicative of a slow quenching process (see e.g. Haines et al., 2013, for a discussion), and that this is occurring only after their Hi gas reservoirs have been significantly depleted.

4 Summary and discussion

We combine HI, IR, optical, UV, and X-ray data from BUDHIES, SDSS, and LoCuSS in and around A963_1, a massive cluster at , to study the effect of (current and past) environment on the transformation of galaxies. In this section, our main findings are first summarized and then discussed.

  • We combine XMM-Newton observations and hundreds of new MMT redshifts from LoCuSS with BUDHIES spectroscopic and Hi data to better characterize A963_1. The X-ray data, together with a dynamical analysis of the cluster galaxies, reveal that A963_1 is a cluster with a large amount of substructure. We identify seven groups with masses between and at . We further use SDSS data to study the large scale structure around the cluster (beyond ), and find that at least two filamentary structures are feeding the cluster. We speculate that the groups within the cluster were accreted via such filaments.

  • When examining the properties of the Hi -detected galaxies in A963_1, we find that they have blue colour (4), sSFRs typically above , and are preferentially located outside the cluster core. Passive galaxies on the other hand, dominate the cluster core and become increasingly less frequent at larger clustercentric radii. However, the incidence of passive galaxies outside is significant.

  • Groups host at least of the passive galaxies in the outskirts of the cluster. Groups more massive than have of their galaxies quenched, and lower mass groups only . The increase of passive galaxies with group mass happens at the expense of Hi detections.

  • From a position versus velocity phase-space diagram we find that: (i) the infall region of the cluster (‘recent infall’ zone) contains a mixed population of passive, star-forming and Hi -detected galaxies; (ii) the ‘stripping’ zone, where RPS is most intense, hosts passive galaxies and a few star-forming galaxies not detected in HI; (iii) the ‘virialized’ zone, where the galaxies that have been in the cluster for over a cluster-crossing time are, is heavily populated with passive galaxies not detected in HI.

  • We use Subaru and -band images from LoCuSS to visually inspect the morphologies of the galaxies in the cluster and its sub-groups. We find many peculiar galaxies, with signatures of galaxy-galaxy interactions, and in very few cases, ‘jellyfish’ features indicative of intense RPS. Peculiar galaxies are common in and around the low-mass groups indentified inside the cluster, while less common (but still present) in the most massive group. Interestingly, peculiar galaxies are also found in the central region of the cluster. The bulk of peculiar galaxies have enhanced star formation.

  • We compute the mean Hi deficiency of blue galaxies as a function of distance from the cluster centre and find an increased deficiency (by a factor of ) in the cluster core. The radial decay of Hi deficiencies in A963_1 agrees within errors with the observations of the (nearby) Coma cluster, although close to the cluster core the cluster galaxies are a factor of 2 less HI deficient than those in Coma.

  • We stack the Hi emission in different regions of the cluster, and find that: (i) Galaxies in the cluster core have significantly lower than those in the outskirts; (ii) the Hi content of group galaxies drops significantly from to groups; (iii) when stacking in phase-space, we find that blue star-forming galaxies (likely falling into the cluster from the field), have their Hi dramatically reduced on their first passage trough the ICM (i.e. the ‘stripping’ zone in phase-space).

  • We also stack the optical spectra in different regions of the cluster. The stacked spectra of group galaxies shows that low-mass groups still host significant star formation activity, while higher mass groups are dominated by passive galaxies. Moreover, we look at the stacked spectra of galaxies in different phase-space regions and find mild evidence for a star formation enhancement in galaxies experiencing RPS as they pass though the ICM for the first time. Finally, we find signs of suppression of star formation in the cluster core when stacking star-forming Hi -poor galaxies.

Our results reveal a complex picture in which several physical mechanisms are responsible for the transformation of galaxies in clusters that are growing through the accretion of groups. We discuss the most relevant in the following.

Ram-pressure stripping: The phase-space distribution and Hi gas content of blue (star-forming galaxies) in the cluster strongly suggests that RPS by the ICM plays an important role in gas removal. In fact, Hi and optical stacking of star-forming galaxies in phase-space revealed that the galaxies’ Hi reservoirs are dramatically reduced as they first pass through the cluster, possibly causing a temporary burst of star formation. Subsequently their star formation activity decreases and eventually halts as they settle into the cluster. Our results are consistent with simulations, that show that as () galaxies experience RPS at pericentric passages in massive clusters, their star formation can be moderately enhanced (Bekki, 2014).

To fully understand the link between gas removal and star formation activity, in a future paper (Cybulski et al., in preparation), we will investigate the relative efficiency of CO versus Hi stripping in A963_1 using data from COOL BUDHIES (Cybulski et al., 2015).

Group pre-processing: Although, RPS is clearly affecting young infalling galaxies, the large fraction of passive galaxies in the infall region of the cluster suggest that RPS is not the only mechanism at play. Mass quenching (though e.g. AGN feedback or bar formation; see e.g. Gavazzi et al., 2015, and references therein) is expected to play an important role in the quenching of high-mass galaxies. For the rest, a possible quenching mechanism is ‘processing’ in group environment before falling into the cluster.

Not all galaxies arrive in the cluster from the field with largely unaffected gas reservoirs. In fact, simulations show that many of them are fed to the cluster from groups and in rare cases even from other clusters (McGee et al., 2009). Indeed, we find many groups currently being accreted into the cluster, accounting for some of the passive galaxies in the cluster outskirts. It is thus likely that at least some of the passive cluster galaxies were quenched in a previous group environment before entering the cluster (‘pre-processing’). The morphologies of the galaxies and the X-ray properties of the groups suggest that in low mass groups, gas-rich interactions among galaxies are frequent. Such interactions typically result in an enhancement of star formation that may contribute to the consumption of gas, leading to a later quenching of the group galaxies. When groups grow and approach a mass , the interactions are less frequent. However, the intensity of ram-pressure increases, contributing to the quenching of the remaining star-forming galaxies. It is possible that additional mechanisms such as tidal interactions are also at play.

Our results are consistent with the most recent multiwavelength study of galaxies in the substructures of the Virgo cluster (Boselli et al., 2014), that suggests that passive cluster galaxies have resulted from gravitational interactions in the infalling groups from which the cluster assembled. A previous BUDHIES analysis of the assembling cluster A2192_1 also supports this scenario (Jaffé et al., 2012). Simulations further show that mergers are key drivers of the morphological sequence in groups. The group galaxies in haloes are expected to frequently experience mergers and galaxy-galaxy interactions, while galaxies in groups will be affected by groups tides and RPS (e.g. strangulation).

Post-processing: In addition, we speculate that the presence of a large population of peculiar galaxies in the cluster centre arise from the passage of groups through the cluster, that “post-processes” group and cluster galaxies by enhancing RPS and inducing galaxy mergers. Simulations (Vijayaraghavan & Ricker, 2013) in fact show that the accretion of a group into a cluster can cause an enhancement of galaxy-galaxy collisions. Shock waves produced by the merger also enhance the ram-pressure intensity on group and cluster galaxies, and the increased local density leads to greater galactic tidal truncation. The combination of these phenomena can significantly transform (post-process) the cluster and group galaxies. This picture is consistent with recent WSRT observations of a merging cluster by Stroe et al. (2015), that revealed gas rich galaxies in the cluster with star formation activity triggered by the passage of the shock.

Our results highlight the entanglement of the different physical processes responsible for building up the large fraction of passive galaxies in clusters, and the need to account for the evolution of environment, imminent in a hierarchical Universe.

Acknowledgements

We thank the referee for providing constructive comments to the paper. YJ gratefully acknowledges Ricardo Demarco and FONDECYT grant no. 3130476 for their valuable support. This work was co-funded under the Marie Curie Actions of the European Commission (FP7-COFUND). CPH was funded by CONICYT Anillo project ACT-1122. MMC gratefully acknowledges support from the Dunlap Institute. The Dunlap Institute is funded through an endowment established by the David Dunlap family and the University of Toronto. We are grateful for support from a Da Vinci Professorship at the Kapteyn Institute. This work has been also supported by the Center for Galaxy Evolution Research funded by the NRF of Korea and Science Fellowship of POSCO TJ Park Foundation. XF is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-1501342. MY and RC acknowledge support from the NASA ADAP grant NNX14AF80G. GPS acknowledges support from the Royal Society and the Science and Technology Facilities Council. NO is supported by a Grant-in-Aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (26800097), the Funds for the Development of Human Resources in Science and Technology under MEXT, Japan, and Core Research for Energetic Universe in Hiroshima University (the MEXT program for promoting the enhancement of research universities, Japan). The Westerbork Synthesis Radio Telescope is operated by the ASTRON (Netherlands Institute for Radio Astronomy) with support from the Netherlands Foundation for Scientific Research (NWO). The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Participating Institutions. The Participating Institutions are The University of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, University of Pittsburgh, Princeton University, the United States Naval Observatory, and the University of Washington. This work is based [in part] on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA. Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

Appendix A Improved spectroscopic completeness

As mentioned in Section 2, we have included for this paper new redshifts in the field of Abell 963. This addition has improved our spectroscopic completeness significantly, as shown in Fig. 19, where the distribution of the targeted galaxies for spectroscopy and those with a redshift are compared as a function of -band magnitude. In particular, we are now complete (%) down to , as indicated by the lines in the bottom panel of the figure.

Figure 19: Spectroscopic completeness after inclusion of new redshifts. The open black histogram shows the R-band distribution of all the galaxies targeted for spectroscopy, while the blue histogram represents galaxies with a measured redshift. If we compare this plot with Figure 6 of Jaffé et al. (2013) it is clear that the inclusion of new redshifts (mainly from LoCuSS) has increased our spectroscopic completeness significantly (%, dashed line) for galaxies with (solid line).

Appendix B Results from the Hi stacks

This section presents Hi stacks that demonstrate the method and its reliability, and a summary of all the Hi stacking results presented in this paper.

Fig. 21 shows the Hi stacks in different bins of colour and magnitude, as labelled. The top panel includes all galaxies in each CMD bin while the bottom panel excludes direct Hi detections. It is noticeable how the 21 cm emission decreases with redder colours and fainter magnitudes. Moreover, the stacks that exclude direct Hi detections always present lower Hi fluxes, although the trend with luminosity and colour holds. This can be appreciated in Fig. 20, where the average is plotted against the three luminosity bins considered, for each colour bin, for all galaxies (solid) and excluding direct Hi detections (dashed).

Figure 20: Results of the Hi stacking in different bins of colour and magnitude, presented in Fig. 21. Solid lines correspond to the stacks that consider all galaxies in each bin, and dashed lines to those that exclude direct Hi detections. When the signal was too low, and only an upper limit on the could be determined, arrows and dotted lines were used.

Table 3 lists the resulting Hi flux and mass for all the Hi stacks presented in this paper.

Figure 21: Hi stacking in different bins of colour and -band magnitude, as labelled (see Fig. 2 as reference). In each panel, the thin grey lines correspond to individual galaxies while the thicker colour lines represent the stacked spectra. The number of galaxies that went into each stack () and number of direct Hi -detections () are indicated in the top-left corner of each panel. The top panel includes Hi -detected plus non-detected galaxies. The bottom panel shows the Hi stacks in the same bins, but excluding direct Hi -detections.
Hi stacking sample No. galaxies Flux
(see Section 3.4) (mJy km )
Colour and magnitude:
Red/passive ( or NUV undet.) and faint (19.5) 116 -0.250.31 0.19
Red/passive ( or NUV undet.) and int (R18.5) 111 0.850.23 0.170.05
Red/passive ( or NUV undet.) and bright (R) 64 0.930.42 0.25
Green (2.5NUV-R4) and faint (19.5) 41 3.620.44 0.730.09
Green (2.5NUV-R4) and int (R18.5) 72 5.980.49 1.200.10
Green (2.5NUV-R4) and bright (R) 15 12.480.58 2.510.12
Bluest (2.5NUV-R) and faint (19.5) 64 13.190.31 2.650.06
Bluest (2.5NUV-R) and int (9.5) 34 16.310.51 3.280.10
Bluest (2.5NUV-R) and bright int (18.5) 6 38.51.6 7.750.32
In/out the cluster (19.5):
Red/passive (4) and cluster core () 131 -0.040.21 0.13
Red/passive (4) and cluster outskirts () 126 1.690.34 0.340.07
Blue () and cluster core () 75 6.130.36 1.230.07
Blue () and cluster outskirts () 143 12.580.43 2.530.09
Groups by mass ()
Red/passive in Groups A 5 5.61.4 1.130.28
Red/passive in Groups B 20 -2.850.62 0.37
Red/Passive in Group C 27 -0.310.68 0.41
Blue in Groups A 11 14.690.72 2.960.14
Blue and Groups B 7 8.930.66 1.800.13
Blue and Group C 7 4.050.66 0.820.13
Phase-space regions ():
Red/Passive in ‘recent infall’ zone 81 2.370.46 0.480.09
Red/Passive in ‘stripping’ zone 28 0.080.64
Red/Passive in ‘virialized’ zone 28 -0.900.36
Blue in ‘recent infall’ zone 85 13.500.44 2.720.09
Blue in ‘stripping’ zone 8 0.940.56
Blue in ‘virialized’ zone 7 1.060.79
Table 3: This table lists the results from the Hi stacks performed in this section and in Section 3.4. Different lines correspond to different bins in colour, magnitude, regions in the cluster, groups, and phase-space zones as indicated in column 1. The second column displays the number of galaxies used in each stack, and columns 3 and 4 show the resulting Hi flux and estimated mass and errors.

Footnotes

  1. This value was measured by extending the profiles in (Martino et al., 2014) to (Haines et al. in preparation)

References

  1. Abadi M. G., Moore B., Bower R. G., 1999, MNRAS, 308, 947
  2. Abramson A., Kenney J. D. P., Crowl H. H., Chung A., van Gorkom J. H., Vollmer B., Schiminovich D., 2011, AJ, 141, 164
  3. Baldry I. K., Balogh M. L., Bower R. G., Glazebrook K., Nichol R. C., Bamford S. P., Budavari T., 2006, MNRAS, 373, 469
  4. Balogh M. L., Morris S. L., Yee H. K. C., Carlberg R. G., Ellingson E., 1999, ApJ, 527, 54
  5. Barnes J. E., Hernquist L., 1996, ApJ, 471, 115
  6. Bekki K., 2014, MNRAS, 438, 444
  7. Bekki K., Couch W. J., Shioya Y., 2002, ApJ, 577, 651
  8. Bell E. F., et al., 2005, ApJ, 625, 23
  9. Boselli A., Gavazzi G., 2006, PASP, 118, 517
  10. Boselli A. et al., 2014, A&A, 570, A69
  11. Bravo-Alfaro H., Cayatte V., van Gorkom J. H., Balkowski C., 2000, AJ, 119, 580
  12. Bravo-Alfaro H., Cayatte V., van Gorkom J. H., Balkowski C., 2001, A&A, 379, 347
  13. Brown T., Catinella B., Cortese L., Kilborn V., Haynes M. P., Giovanelli R., 2015, MNRAS, 452, 2479
  14. Byrd G., Valtonen M., 1990, ApJ, 350, 89
  15. Calzetti D. et al., 2007, ApJ, 666, 870
  16. Catinella B., Cortese L., 2015, MNRAS, 446, 3526
  17. Catinella B. et al., 2013, MNRAS, 436, 34
  18. Cayatte V., van Gorkom J. H., Balkowski C., Kotanyi C., 1990, AJ, 100, 604
  19. Chon G., Böhringer H., Smith G. P., 2012, A&A, 548, A59
  20. Chung A., van Gorkom J. H., Kenney J. D. P., Crowl H., Vollmer B., 2009, AJ, 138, 1741
  21. Chung A., van Gorkom J. H., Kenney J. D. P., Vollmer B., 2007, ApJ, 659, L115
  22. Cohn J. D., 2012, MNRAS, 419, 1017
  23. Cortese L., Catinella B., Boissier S., Boselli A., Heinis S., 2011, MNRAS, 415, 1797
  24. Crowl H. H., Kenney J. D. P., van Gorkom J. H., Vollmer B., 2005, AJ, 130, 65
  25. Cybulski R. et al., 2015, preprint (arXiv:1510.08450)
  26. Cybulski R., Yun M. S., Fazio G. G., Gutermuth R. A., 2014, MNRAS, 439, 3564
  27. Daddi E. et al., 2007, ApJ, 670, 156
  28. De Lucia G., Weinmann S., Poggianti B. M., Aragón-Salamanca A., Zaritsky D., 2012, MNRAS, 423, 1277
  29. Desai V. et al., 2007, ApJ, 660, 1151
  30. Dressler A. et al., 1997, ApJ, 490, 577
  31. Dressler A., Shectman S. A., 1988, AJ, 95, 985
  32. Duc P.-A., Bournaud F., 2008, ApJ, 673, 787
  33. Falco M., Hansen S. H., Wojtak R., Brinckmann T., Lindholmer M., Pandolfi S., 2014, MNRAS, 442, 1887
  34. Fasano G., Poggianti B. M., Couch W. J., Bettoni D., Kjærgaard P., Moles M., 2000, ApJ, 542, 673
  35. Fernández X. et al., 2013, ApJ, 770, L29
  36. Gavazzi G. et al., 2015, A&A, 580, A116
  37. Gavazzi G., Jaffe W., 1985, ApJ, 294, L89
  38. Gavazzi G. et al., 2013, A&A, 553, A90
  39. Giovanelli R., Haynes M. P., 1985, ApJ, 292, 404
  40. Gunn J. E., Gott J. R. I., 1972, ApJ, 176, 1
  41. Haines C. P., Gargiulo A., La Barbera F., Mercurio A., Merluzzi P., Busarello G., 2007, MNRAS, 381, 7
  42. Haines C. P. et al., 2013, ApJ, 775, 126
  43. Haines C. P. et al., 2009, ApJ, 704, 126
  44. Haynes M. P., Giovanelli R., 1984, AJ, 89, 758
  45. Hess K. M., Wilcots E. M., 2013, AJ, 146, 124
  46. Hudson D. S., Mittal R., Reiprich T. H., Nulsen P. E. J., Andernach H., Sarazin C. L., 2010, A&A, 513, A37
  47. Hwang H. S., Geller M. J., Diaferio A., Rines K. J., Zahid H. J., 2014, ApJ, 797, 106
  48. Jaffé Y. L., Poggianti B. M., Verheijen M. A. W., Deshev B. Z., van Gorkom J. H., 2012, ApJ, 756, L28
  49. Jaffé Y. L., Poggianti B. M., Verheijen M. A. W., Deshev B. Z., van Gorkom J. H., 2013, MNRAS, 431, 2111
  50. Jaffé Y. L., Smith R., Candlish G. N., Poggianti B. M., Sheen Y.-K., Verheijen M. A. W., 2015, MNRAS, 448, 1715
  51. Kapferer W., Sluka C., Schindler S., Ferrari C., Ziegler B., 2009, A&A, 499, 87
  52. Katayama H., Hayashida K., Takahara F., Fujita Y., 2003, ApJ, 585, 687
  53. Kenney J. D. P., van Gorkom J. H., Vollmer B., 2004, AJ, 127, 3361
  54. Kennicutt, R. C., Jr, et al., 2007, ApJ, 671, 333
  55. Leauthaud A. et al., 2010, ApJ, 709, 97
  56. Lee J., Kim S., Rey S.-C., 2015, ApJ, 807, 122
  57. Lewis I., et al., 2002, MNRAS, 334, 673
  58. McGee S. L., Balogh M. L., Bower R. G., Font A. S., McCarthy I. G., 2009, MNRAS, 400, 937
  59. Madau P., Dickinson M., 2014, ARA&A, 52, 415
  60. Mahajan S., Raychaudhury S., Pimbblet K. A., 2012, MNRAS, 427, 1252
  61. Mann A. W., Ebeling H., 2012, MNRAS, 420, 2120
  62. Martino R., Mazzotta P., Bourdin H., Smith G. P., Bartalucci I., Marrone D. P., Finoguenov A., Okabe N., 2014, MNRAS, 443, 2342
  63. Miller R. H., 1986, A&A, 167, 41
  64. Moore B., Lake G., Quinn T., Stadel J., 1999, MNRAS, 304, 465
  65. Murphy E. J. et al., 2011, ApJ, 737, 67
  66. Okabe N., Smith G. P., 2015, preprint (arXiv:1507.04493)
  67. Okabe N., Takada M., Umetsu K., Futamase T., Smith G. P., 2010, PASJ, 62, 811
  68. Peng Y., Maiolino R., Cochrane R., 2015, Nature, 521, 192
  69. Peng Y., et al., 2010, ApJ, 721, 193
  70. Poggianti B. M. et al., 2009, ApJ, 697, L137
  71. Poggianti B. M., van Gorkom J. H., 2001, in J. E. Hibbard, M. Rupen, & J. H. van Gorkom, eds, ASP Conf. Ser. Vol. 240, Gas and Galaxy Evolution. Astron. Soc. Pac., San Francisco, p. 599
  72. Porter S. C., Raychaudhury S., Pimbblet K. A., Drinkwater M. J., 2008, MNRAS, 388, 1152
  73. Quilis V., Moore B., Bower R., 2000, Science, 288, 1617
  74. Ragone-Figueroa C., Plionis M., 2007, MNRAS, 377, 1785
  75. Rieke G. H., Alonso-Herrero A., Weiner B. J., Pérez-González P. G., Blaylock M., Donley J. L., Marcillac D., 2009, ApJ, 692, 556
  76. Roediger E., 2009, Astron. Nachr., 330, 888
  77. Roediger E., Brüggen M., 2007, MNRAS, 380, 1399
  78. Salim S. et al., 2005, ApJ, 619, L39
  79. Salim S. et al., 2007, ApJS, 173, 267
  80. Sanderson A. J. R., Edge A. C., Smith G. P., 2009, MNRAS, 398, 1698
  81. Schiminovich D. et al., 2007, ApJS, 173, 315
  82. Scott T. C. et al., 2010, MNRAS, 403, 1175
  83. Scott T. C., Cortese L., Brinks E., Bravo-Alfaro H., Auld R., Minchin R., 2012, MNRAS, 419, L19
  84. Serra P. et al., 2012, MNRAS, 422, 1835
  85. Smith G. P., Kneib J.-P., Smail I., Mazzotta P., Ebeling H., Czoske O., 2005, MNRAS, 359, 417
  86. Smith R. et al., 2015, MNRAS, 454, 2502
  87. Solanes J. M., Sanchis T., Salvador-Solé E., Giovanelli R., Haynes M. P., 2002, AJ, 124, 2440
  88. Speagle J. S., Steinhardt C. L., Capak P. L., Silverman J. D., 2014, ApJS, 214, 15
  89. Stroe A., Oosterloo T., Röttgering H. J. A., Sobral D., van Weeren R., Dawson W., 2015, MNRAS, 452, 2731
  90. Tonnesen S., Bryan G. L., 2009, ApJ, 694, 789
  91. Toomre A., Toomre J., 1972, ApJ, 178, 623
  92. Verdes-Montenegro L., Yun M. S., Williams B. A., Huchtmeier W. K., Del Olmo A., Perea J., 2001, A&A, 377, 812
  93. Verheijen M. et al., 2010, Proc. of ISKAF2010 Science Meeting, A New Golden Age for Radio Astronomy.
  94. Verheijen M., van Gorkom J. H., Szomoru A., Dwarakanath K. S., Poggianti B. M., Schiminovich D., 2007, ApJ, 668, L9
  95. Vijayaraghavan R., Ricker P. M., 2013, MNRAS, 435, 2713
  96. Vollmer B., 2003, A&A, 398, 525
  97. Vollmer B., Cayatte V., Balkowski C., Duschl W. J., 2001, ApJ, 561, 708
  98. Walker I. R., Mihos J. C., Hernquist L., 1996, ApJ, 460, 121
  99. Yun M. S., Ho P. T. P., Lo K. Y., 1994, Nature, 372, 530
  100. Zabludoff A. I., Mulchaey J. S., 1998, ApJ, 496, 39
  101. Zhang Y.-Y., Finoguenov A., Böhringer H., Kneib J.-P., Smith G. P., Czoske O., Soucail G., 2007, A&A, 467, 437
  102. Zibetti S., Charlot S., Rix H.-W., 2009, MNRAS, 400, 1181
103724
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
Edit
-  
Unpublish
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel
Comments 0
Request answer
""
The feedback must be of minumum 40 characters
Add comment
Cancel
Loading ...