EDP Sciences
Free access
Issue
A&A
Volume 537, January 2012
Article Number A58
Number of page(s) 8
Section Cosmology (including clusters of galaxies)
DOI http://dx.doi.org/10.1051/0004-6361/201117973
Published online 09 January 2012

© ESO, 2012

1. Introduction

The level of star formation (SF) activity in galaxy systems is known to be suppressed relative to the field. According to the well known morphology-density relation (Dressler 1980) and the star formation rate (SFR)-density relation (Hashimoto et al. 1998; Lewis et al. 2002; Gómez et al. 2003), in the local Universe high density regions, like groups and clusters, host mostly early type galaxies characterized by a lower level of SF activity than field (mostly late-type) galaxies.

The environmental dependence of galaxy SFR may change with redshift, as galaxies in systems undergo significant evolution. In higher-redshift clusters, the fraction of blue galaxies is higher (the so-called “Butcher-Oemler” effect, Butcher & Oemler 1978, 1984; Pimbblet 2003), and so are the fractions of cluster galaxies with spectra characterized by young stellar populations (Dressler & Gunn 1983; Dressler et al. 1999; Poggianti et al. 1999) and the fraction of galaxies with late-type morphology (Dressler et al. 1997; Fasano et al. 2000; Postman et al. 2005; Smith et al. 2005). Higher-z clusters are also observed to contain a higher fraction of infrared (IR) emitting galaxies (the so called “IR Butcher-Oemler effect”, Saintonge et al. 2008; Haines et al. 2009; Temporin et al. 2009), where most of the IR emission is powered by SF.

Given this evolution and the morphology- and SFR-density relations observed locally, it is clear that a quenching of the star formation activity of galaxies in dense environments is required since z ~ 1–2, and that this quenching process must act faster in galaxy systems than in the field. When and where this quenching process takes place, is still a matter of debate. SF quenching is generally assumed to occur within the cluster environment, where processes like ram pressure (Gunn & Gott 1972), cumulative galaxy-galaxy hydrodynamic/gravitational interactions (Park & Hwang 2009), strangulation (Larson et al. 1980), and galaxy harassment (Moore et al. 1996) are particularly effective. Indeed, Tran et al. (2009) observed a higher fraction of star forming galaxies in a supergroups at z = 0.37 with respect to clusters at the same redshift. It has however also been claimed that SF quenching of cluster galaxies occur in low-mass groups or large-scale filaments prior to cluster assembly (the so-called “pre-processing” Zabludoff & Mulchaey 1998; Kodama et al. 2001; Fadda et al. 2008; Porter et al. 2008; Balogh et al. 2011). The accretion process of groups onto clusters can itself lead to star-formation quenching (Poggianti et al. 2004), perhaps caused by rapid gas consumption due to a sudden enhancement of the SF activity (Miller & Owen 2003; Coia et al. 2005; Ferrari et al. 2005). It is also not clear whether environment plays any role at all in the quenching process. According to Peng et al. (2010), mass quenching is the dominant quenching process for massive galaxies, which generally reside in massive halos like groups and clusters. Since mass quenching should occur when a galaxy reaches a limiting mass, more massive galaxies, which reside in more massive halos, should be quenched earlier than less massive ones, having reached earlier the limiting mass for quenching. This creates an environment dependent quenching of galaxy SF, which is driven by some internal, rather than external, process, such as AGN feedback. This is supported observationally by the analysis of the star-formation histories of galaxies in the Virgo cluster region (Gavazzi et al. 2002), galaxies of higher H-band luminosities being characterized by shorter timescales of SF. Additional support to this scenario comes from the analysis of chemical abundances in elliptical galaxies by Pipino & Matteucci (2004).

Another way of looking at the evolution of the SF activity in galaxy systems is to consider a global quantity such as the star formation rate per unit of halo mass, that is the sum of the SFRs of all the galaxy in a system, divided by the system total mass, Σ(SFR)/M. According to recent results (Kodama et al. 2004; Finn et al. 2004, 2005; Geach et al. 2006; Bai et al. 2009; Chung et al. 2010; Koyama et al. 2010; Hayashi et al. 2011), the evolution of the cluster Σ(SFR)/M is rapid, with a redshift dependence Σ(SFR)/M  ∝ (1 + z)α, with α ≃ 5−7. The cluster Σ(SFR)/M evolves faster up to z ~ 1 than the fraction of IR-emitting galaxies (the IR Butcher-Oemler effect) because IR-emitting galaxies are not only more numerous in higher-z clusters but also more IR luminous, given the evolution of the cluster IR luminosity function (Bai et al. 2009). In addition, there are some indications that Σ(SFR)/M anti-correlates with halo mass in galaxy systems at similar redshift (Finn et al. 2004, 2005; Homeier et al. 2005; Koyama et al. 2010).

The quantity Σ(SFR)/M thus appears to be a powerful tool for analyzing the rapid evolution of galaxies in systems of different mass. To determine Σ(SFR)/M one needs to estimate the system dynamical mass and the SFRs of its member galaxies. Robust SFR estimates can be obtained from galaxy IR luminosities, LIR, and these estimates are considered to be more reliable than the SFR based on optical emission-lines such as [OII] and Hα (Geach et al. 2006).

If one considers only the systems where the total SFR has been derived from IR observations, the Σ(SFR)/M – z relation is currently based on  ~15 galaxy clusters only (Biviano et al. 2011), all with masses1M200 ≳ 1014 M and z < 1. Masses for these systems have been derived in different ways (from X-ray or lensing data, or from kinematical analyses of the population of cluster galaxies), and at different limiting radii.

The aim of this paper is to extend the study of the Σ(SFR)/M – z relation to systems of lower masses than clusters, and to higher z, with a homogeneous estimate of the system masses and total SFRs. Specifically, we extend the mass range down to the group mass regime, and to the typical dark matter halo mass of field galaxies (1011−1012 M). We also extend the redshift range up to z ~ 1.6. We use deep far-infrared PACS data at 100 and 160 μm to accurately estimate the total SFRs of 18 systems with masses in the range 2 × 1013−3 × 1015 M, down to LIR = 1011L, i.e. the luminosity that characterizes the so-called luminous IR galaxies (LIRGs). This limit is set by the faintest luminosity observable with our data at z ~ 1.6. We use nearly complete galaxy spectroscopic samples to measure the system masses from their kinematics. To extend the Σ(SFR)/M – z relation to the field regime, i.e. to dark matter halo masses typical of field galaxies, we use published global SFR densities at different redshifts (Madau et al. 1998; Gruppioni et al. 2011; Magnelli et al. 2011).

In Sects. 2.1 and 2.2 we describe our sample and our spectroscopic and photometric data-set, in Sect. 2.3 we describe how we define system membership of the galaxies with available redshifts, and derive the system dynamical masses, in Sect. 2.4 we determine the total SFRs of the systems. The resulting Σ(SFR)/M – z relations are presented in Sect. 3 for clusters, groups (Sect. 3.1), and the field (Sect. 3.2). We summarize and discuss our results in Sect. 4.

We adopt H0 = 70 km s-1 Mpc-1, Ωm = 0.3, ΩΛ = 0.7 throughout this paper. We adopt a Salpeter IMF when calculating stellar masses and SFRs.

2. The data-set

2.1. The sample of galaxy systems

Our sample comprises 9 X-ray detected systems at 0.15 < z < 0.85, all observed in the PACS Evolutionary Probe PACS GT program (PEP, Lutz et al. 2011). The systems are listed in Table 1. In addition, we use the X-ray detected group and cluster catalog of the COSMOS field (Finoguenov et al., in prep.) to identify 27 systems at z < 0.8 with at least 10 spectroscopic members each, the minimum required for an acceptable mass estimate from the system kinematics (Girardi et al. 1993). The COSMOS field is observed as part of the PEP program. We further include in our sample 4 systems in the GOODS North and South fields at 0.7 < z < 1.6 observed as part of the PACS GOODS-Herschel Program (Elbaz et al. 2011). Of these systems, two are X-ray detected, one at z = 0.73 in GOODS-S (Gilli et al. 2003; Le Fèvre et al. 2004; Popesso et al. 2009), another one at z = 1.02 in GOODS-N (Elbaz et al. 2007), one system lies close to the Chandra CCD chip gap and it is not X-ray detected (GOODS-N system at z = 0.85 Bauer et al. 2002), and the highest redshift system (at z = 1.612) is only marginally X-ray detected in the Chandra GOODS-S map (Kurk et al. 2009). Finally, we include in our sample the Bullet cluster (Barrena et al. 2002; Markevitch et al. 2004), a well-known massive cluster undergoing a violent collision with an infalling group, which we use as a test-case for the dependence of the Σ(SFR)/M on the cluster dynamical state. PACS data for the Bullet cluster are not public yet. For this cluster we use the LIR based on Spitzer MIPS 24 μm data published by Chung et al. (2010). The considered LIRGs members lie within 1.7 Mpc from the cluster center. This is close to the r200 radius,  ~2 Mpc, estimated from the kinematics of cluster galaxies (Barrena et al. 2002).

We split our sample in groups and clusters depending on their mass estimates (see next section), groups/poor clusters with 1013 M < M200 < 3 × 1014 M and clusters with M200 ≥ 3 × 1014 M. The mass value separating groups/poor clusters from clusters approximately corresponds to the richness class 0 in the catalog of Abell et al. (1989) (using the relations between mass and velocity dispersion and between velocity dispersion and richness class for galaxy clusters; see Girardi et al. 1993; Biviano et al. 2006).

2.2. Spectral and photometric data

We take the data for the COSMOS field from the multi-wavelength and IRAC-selected catalog of Ilbert et al. (2010) complemented with public Spitzer 24 μm (Le Floc’h et al. 2009; Sanders et al. 2007) and PACS 100 and 160 μm data (Lutz et al. 2011). The association between 24 μm and PACS sources with their optical counterparts is done via a maximum likelihood method (see Lutz et al. 2011; Berta et al. 2010, for details). The photometric sources were cross-matched in coordinates with the sources for which a high-confidence spectroscopic redshift is available. For this purpose we use the public catalogs of spectroscopic redshifts complemented with other unpublished data. This catalog includes redshifts from either SDSS or the public zCOSMOS-bright data acquired using VLT/VIMOS (Lilly et al. 2007, 2009) complemented with Keck/DEIMOS (PIs: Scoville, Capak, Salvato, Sanders, Kartaltepe), Magellan/IMACS (Trump et al. 2007), and MMT (Prescott et al. 2006) spectroscopic redshifts.

Similar multi-wavelength (from UV to 24 μm) data for the GOODS-S field is provided by the MUSIC catalog (Grazian et al. 2006; Santini et al. 2009) and complemented with PACS data at 70, 100 and 160 μm and spectroscopic redshifts from the PEP catalog (see Lutz et al. 2011 for more details) and from the GMASS survey (Cimatti et al. 2008; Kurk et al. 2009).

We take the data for the GOODS-N field from the multi-wavelength catalog prepared by the PEP Team (Berta et al. 2010) including UV to MIPS 24 data, PACS 100 and 160 μm data, complemented with the spectroscopic redshift catalog of Barger et al. (2008).

We collect publicly available spectroscopic data for all clusters from the NED. For CL 0024+17 and MS 0451−03 we also collect multi-wavelength data from Moran et al. (2005, 2007). Multi-wavelength photometric data are not available for the other clusters. The references for each cluster spectroscopic sample are listed in Table 1. Each cluster sample is complemented with MIPS 24 μm data, when available, and PACS 100 and 160 μm data (Lutz et al. 2011).

The data reduction, source extraction, and a detailed explanation of the depths of the PACS maps used in this paper can be found in Lutz et al. (2011) and Elbaz et al. (2011), for the PEP and GOODS-Herschel surveys, respectively.

Table 1

Main properties of the cluster and group samples.

2.3. Membership and mass estimates of galaxy systems

To define the galaxies that are members of the systems described in Sect. 2.1, we adopt the algorithm of Mamon et al. (2010), which is based on the modeling of the mass and anisotropy profiles of cluster-sized halos extracted from a cosmological numerical simulation. This algorithm is more effective than traditional approaches (e.g. Yahil & Vidal 1977) in rejecting interlopers, while still preserving cluster members. The system membership selection depends on the location of galaxies in the system-centric distance – rest-frame velocity2, where is the system mean redshift, determined with the biweight estimator (Beers et al. 1990). diagram. The peaks of the X-ray surface brightness are adopted as centers of the X-ray detected systems. The galaxy number density maxima (estimated using an adaptive kernel technique) are adopted as centers for the other systems. The interloper rejection procedure is iterated until convergence.

thumbnail Fig. 1

Comparison of system mass estimates obtained from X-ray data with those obtained from system velocity dispersions.

Open with DEXTER

As an outcome of the procedure of Mamon et al. (2010), the system dynamical masses are obtained from biweight estimates (Beers et al. 1990) of the system velocity dispersions along the line-of-sight, σv, as in Mauduit & Mamon (2007, see Appendix A). We use the σv-based mass estimates, rather than the more traditional virial mass estimates because the σv-based mass estimates are robust against problems of spatial incompleteness, unlike the more traditional virial theorem estimator (see Biviano et al. 2006), so we do not have to worry about the possibility that in some systems of our sample the spectroscopic completeness is different in different regions.

The uncertainties in the system mass estimates are estimated from the uncertainties in the biweight estimate of σv (Beers et al. 1990) via the propagation of error analysis.

As an alternative mass estimate we consider that obtained by assuming hydrostatic equilibrium of the X-ray-emitting intra-cluster or intra-group gas. The two mass-estimates are compared in Fig. 1 for those systems which have both estimates available (see Leauthaud et al. 2010 for the COSMOS groups). The scatter is large, but there is no significant bias, on average, between the two mass estimates. There is only one system for which the two mass estimates are in disagreement by more than 3σ. The discrepancy is probably due the low number of spectroscopic members (9) available for this system. We exclude it from our sample for further analyses.

2.4. Star formation rates

We fit the spectral energy distributions (SEDs) of all system members in our sample which have multiwavelength data available from the UV to the PACS wavelengths, using a set of empirical templates of local objects (Polletta et al. 2007; Gruppioni et al. 2010) that reproduce most of the observed galaxy SEDs from the optical to the FIR wavelengths (see Rodighiero et al. 2010). We then determine the galaxy LIR by integrating the best-fit SED models from 8 to 1000 μm.

We estimate the LIR of those system members for which multiwavelength data are not available, by fitting their PACS data with SED models from Polletta et al. (2007), to extract the νLν at 100 μm rest frame, and by using the LIR −   νLν relation, as explained in Popesso et al. (in prep.).

Some of our member galaxies do not have PACS data either due to the incompleteness of the PACS photometric catalog between the 3 and 5σ levels, or because the PACS catalog are not deep enough to reach the required 1011 L limits at the redshift of the system. We determine the LIR of these galaxies from their MIPS 24 μm flux densities by using the scaling relation of Chary & Elbaz (2001), corrected as explained in Popesso et al. (in prep.) to avoid overestimations (Elbaz et al. 2011; Nordon et al. 2011). We also re-determine in the same way the LIR of the member galaxies of the Bullet cluster from the 24 μm data of Chung et al. (2010). In total, apart from the Bullet cluster, only 8 galaxies (9% of our member LIRGs) do not have their LIR estimated from PACS data: 2 in MS 0451−03, 4 in RX J0152−13 and 2 in the z ~ 1.6 GOODS-S group.

The main advantage of using PACS data in the determination of the LIR is to remove possible contamination by the emission of active galactic nuclei (AGN), since most of the rest-frame far-IR emission comes from the host galaxy (Netzer et al. 2007; Lutz et al. 2010). Thus, we do not need to exclude AGN hosts from our sample as was instead done in previous similar analyses (e.g. Geach et al. 2006; Bai et al. 2009; Biviano et al. 2011). AGN contamination is a concern, however, when estimating the SFR from 24 μm data. We do exclude the AGN hosts from the evaluation of the total SFR of the Bullet cluster, using the AGN classification of Chung et al. (2010), since PACS data are not yet publicly available for this cluster. We also cross-correlated the system members, whose LIR is based on 24 μm data only, with publicly available AGN catalogs (Johnson et al. 2003; Demarco et al. 2005; Brusa et al. 2009). None of those members turn out to be an AGN.

Finally, we estimate the galaxy SFRs from their LIR via the law of Kennicutt (1999). The LIRG LIR-limit corresponds to a limiting SFR of 17   M yr-1.

We estimate Σ(SFR), the total SFR of each system, by summing up the SFRs of the LIRGs that are members of the system. We then correct Σ(SFR) for spectroscopic incompletetness in the following way. For all systems where the PACS data reach the 1011 L limit at the system redshift, we take as a reference the PACS 100 μm catalog. We plot the 100 μm flux density vs. LIR of system member galaxies in narrow redshift bins to estimate the 100 μm flux density that corresponds to the LIRG LIR at each redshift, fLIRG(z). The spectroscopic incompleteness is defined as the fraction of PACS sources with z, among all those with flux  >fLIRG(z) within r200. The inverse of this fraction is the correction factor that we apply to the Σ(SFR) estimates to correct for incompleteness. For the systems which contain member galaxies for which the LIR-estimates are based on MIPS 24 μm flux densities, we use the same approach by using the 24 μm catalog as reference.

To test the reliability of this method, we consider a system with 100% spectroscopic completeness down to the LIR LIR level, MS 1054−03 (Bai et al. 2009). We perform 500 000 random extractions from the luminosity function of MS 1054−03 of a number of LIRGs corresponding to a given spectroscopic completeness. We repeat this test for all the spectroscopic coverage values listed in Table 1. We then estimate the uncertainty due to the incompleteness correction from the dispersion of the difference between the recovered Σ(SFR) and the real value. The uncertainty due to this correction is  ≃20% in the range of spectroscopic completeness of our sample.

The total uncertainty of the Σ(SFR) estimates is determined from the propagation of error analysis, by considering a 10% uncertainty in the LIR estimates (see Lutz et al. 2011, for further details), and the 20% uncertainty due to the completeness correction.

The spectroscopic coverage of the COSMOS field is much lower than those of the clusters and GOODS fields and spatially not uniform. Thus, to reduce the error bars of the Σ(SFR), we combine the COSMOS groups and, separately, the COSMOS clusters in four, and, respectively, two redshift bins, see Table 1. In each z-bin we then define the mean Σ(SFR), M200, and spectroscopic incompleteness, separately for the groups and the clusters contributing to the sample in that bin.

Counting the composite COSMOS groups and clusters as individual systems, our final sample comprises 9 massive clusters, 9 groups/poor clusters, and the Bullet cluster (see Table 1).

thumbnail Fig. 2

Σ(SFR)/M- redshift relation for clusters (black symbols) and groups (magenta symbols). Square symbols, triangles, and dots identify, respectively, the COSMOS composite systems, the GOODS systems, and the remaining systems. The red star identifies the Bullet cluster. The black solid line shows the best fit Σ(SFR)/M – z relation for the cluster sample, excluding the Bullet cluster. The relation fitted to the data is of the type Σ(SFR)/M  ∝ zα. The dashed line shows the relation, Σ(SFR)/M  ∝ (1 + z)α, of Bai et al. (2009) rescaled to the LIRG LIR-regime and to the r200 region. The magenta solid line shows the best fit for the sample of groups and poor clusters.

Open with DEXTER

3. Results: the Σ(SFR)/M – z relation

3.1. Clusters versus groups

The Σ(SFR)/M-redshift relation is shown in Fig. 2. Black, and magenta symbols show the Σ(SFR)/M-redshift relation for the clusters and, respectively, the groups. For both samples there is evidence for a significant Σ(SFR)/M vs. z correlation, 99% significant according to a Spearman rank correlation test. On the other hand, we do not find evidence for significant Σ(SFR)/M – M200 or M200 − z correlations within the group and cluster samples separately. Hence the observed Σ(SFR)/M – z relation within each sample must be interpreted as a genuine redshift evolution of the Σ(SFR)/M of galaxy systems.

Groups appear to be characterized by higher Σ(SFR)/M values than clusters, at all redshifts, i.e. they show a higher SF activity than massive clusters. In this sense, when the sample of groups and clusters are considered together, Σ(SFR)/M does anticorrelate with the system M200, as already found in previous studies (Finn et al. 2004, 2005; Homeier et al. 2005; Koyama et al. 2010).

The solid curves in Fig. 2 represent best-fit models to the observed Σ(SFR)/M – z relations, Σ(SFR)/M  = (66 ± 23) × z1.77±0.36 for the cluster sample and Σ(SFR)/M  = (213 ± 44) × z1.33± 0.34 for the group sample. In Fig. 2 we also show as a dashed line the Σ(SFR)/M  ∝ (1 + z)5.3 model proposed by Bai et al. (2009). We rescale their relation, which is built for Σ(SFR) computed down to a SFR limit of 2   M yr-1, to our adopted LIRG limit, by using the IR luminosity function of Bai et al. (2009). Moreover, the Σ(SFR) in Bai et al.’s relation is evaluated within a radius of r500, the radius corresponding to an overdensity equal to 500 times the critical density. To convert to our adopted radius r200, we adopt a Navarro et al. (1997) model, wich concentration c = 5 to scale for the number of galaxies inside the two radii and derive the r500/r200 ratio, and then use Fig. 7 in Bai et al. (2009) to account for the different fraction of IR-emitting galaxies at r500 and r200. The modified relation of Bai et al. (2009) appears to fail to predict the rapid evolution of cluster Σ(SFR)/M, as already suggested by Biviano et al. (2011).

3.2. Structures versus field

In Fig. 3 we compare the Σ(SFR)/M – z relation of galaxy systems with the corresponding relation for field galaxies (light blue shaded region and blue dashed line). The Σ(SFR)/M – z relation of field galaxies is obtained by dividing the observed Star Formation Rate Density (SFRD) of Magnelli et al. (2011) and the modeled SFRD of Gruppioni et al. (2011, dashed line), by the mean comoving density of the universe (Ωm × ρc where Ωm = 0.3 and ρc is the critical density of the Universe). Both SFRD have been evaluated only down to the SFR corresponding to the LIRG LIR , via the Kennicutt relation.

The field SFRD has been estimated in large comoving volumes that include galaxy systems, voids, and isolated galaxies, and is thus representative of the general field galaxy population. According to the dark halo mass function of Jenkins et al. (2001), halos of 1011−1012 M give the main contribution to the dark matter budget at all redshifts. Thus, the Σ(SFR)/M of Fig. 3 can be considered as an effective estimate of the Σ(SFR)/M of galaxy-sized dark matter halos.

Observed and modeled Σ(SFR)/M are in very good agreement within the error bars. They increase from z = 0 to z ~ 1 where they reach a plateau.

The field Σ(SFR)/M – z relation lies above both the group and the cluster relations at z ≲ 0.2. The field Σ(SFR)/M is more than an order of magnitude higher than those of Abell 2218 at z = 0.175 and of the composite COSMOS group at z = 0.1. At higher redshifts, group and field galaxy halos show comparable Σ(SFR)/M. Possibly, the group Σ(SFR)/M lies above the field relation at z ~ 1.6, consistently with a reversal of the SFR- density relation observed by Elbaz et al. (2007), Cooper et al. (2008), and Popesso et al. (2011). However, this conclusion is based on a single group Σ(SFR)/M determination. Clusters have lower Σ(SFR)/M than field galaxy halos at all redshifts up to the last measured point at z = 0.85. A blind extrapolation of the best-fit cluster Σ(SFR)/M – z relation would suggest that clusters should display a higher SFR per unit mass than field galaxies at at z ≳ 2. However, such an extrapolation is extremely uncertain with the present data and massive systems as the ones considered here should not even exist at such redshift according to ΛCDM hierarchical model predictions.

thumbnail Fig. 3

The field Σ(SFR)/M-redshift relation from Magnelli et al. (2011) (light blue shaded region) and Gruppioni et al. (2011) (dashed blue line). The shading and error bars represent 1σ confidence levels Other symbols/lines have the same meaning of the symbols as in Fig. 2.

Open with DEXTER

4. Discussion and conclusions

We find significant evidence for an evolution of the Σ(SFR)/M with z in clusters, groups, and the field. Both the cluster and the group Σ(SFR)/M increase monotonically with z and the field Σ(SFR)/M reaches a maximum at z ≳ 1.

The cluster Σ(SFR)/M evolution we find, appears to be faster, at least in the range 0 < z < 1, compared to the evolutionary relation suggested by Geach et al. (2006); Bai et al. (2007, 2009); Koyama et al. (2010). This might be due to our higher-LIR cut, if LIRGs evolve faster than less luminous IR-emitting galaxies. A faster evolution has already been suggested by Biviano et al. (2011), but only in the z ≲ 0.5 range, while in our data there is no sign of a flattening of the Σ(SFR)/M – z relation for clusters above that redshift.

We confirm earlier indications (Finn et al. 2005; Bai et al. 2007; Koyama et al. 2010) that Σ(SFR)/M is lower for systems of higher mass, and we show for the first time that this is true at all redshifts from z ~ 0 to z ~ 0.9. This is true when we compare clusters to groups and the field. On the other hand, the group Σ(SFR)/M lies below the field relation only at z ≲ 0.2. At higher redshifts, the group and the field Σ(SFR)/M become comparable. There is a hint of a reversal of the SFR activity vs. environment at z ~ 1.6, where the group Σ(SFR)/M lies above the field Σ(SFR)/M  – z relation.

Our results appear to support a scenario in which the quenching of SF occurs earlier in galaxies embedded in more massive halos, i.e. first in clusters, then in groups, and finally in the field. This would be consistent with a “halo downsizing” effect, whereby massive halos evolve more rapidly than low mass halos (Neistein et al. 2006). Although in a narrow range of halo masses, the “halo downsizing” effect has already been observed in the stellar-to-total mass ratio vs. redshift relation by Foucaud et al. (2010).

What is causing the difference between the cluster and the group/field Σ(SFR)/M  – z relation? The similarity of the field and the groups Σ(SFR)/M  – z relations suggests that a SF quenching effect is taking place mostly after galaxies enter the cluster environment, and not in groups before they merge into more massive structures, as would be predicted by the pre-processing scenario (Zabludoff & Mulchaey 1998; Kodama et al. 2001; Balogh et al. 2011). We can therefore consider the field and group Σ(SFR)/M  – z relation as indicating the normal galaxy evolution, and the cluster Σ(SFR)/M  – z relation as the accelerated evolution experienced by galaxies following their accretion to more massive structures, as individuals or in groups. However, we point out that other differences may intervene between the star forming population of the groups and the field. Indeed, in agrement with our results, Tran et al. (2009) find consistent levels of star formation in the supergroup at z = 0.37 and the field population, but they also reveal a much higher fraction of early type galaxies among the SF galaxies of the group environment with respect to the field. Similar results are obtained by Hwang & Park (2009). This could indicate a morphological transformation happening prior to the SF quenching effect.

Part of the observed evolution may come from internal, rather than external processes. In fact, most of the LIRGs in our sample are likely to be rather massive galaxies (as inferred from the almost linear, z-dependent, relation between galaxy SFRs and stellar masses of Elbaz et al. 2007), and according to Peng et al. (2010) massive galaxies evolve mostly because of an internally driven process, called “mass quenching”, caused perhaps by feedback from active galactic nuclei. But since this process is unlikely to be more efficient in quenching SF of massive galaxies in clusters than in other environments, we need to consider an additional quenching process that would operate only in the cluster environment.

Two obvious candidates for producing this additional quenching, are ram-pressure stripping (Gunn & Gott 1972) and starvation (Larson et al. 1980). Ram-pressure stripping is a fast process (Abadi et al. 1999) and is effective only in galaxy systems where the gas density is high, i.e. in clusters. However, ram-pressure stripping seems to be ruled out by the presence of LIRGs contributed by the infalling group (the “bullet”) in the Bullet cluster. The group LIRG contribution is probably the reason of the position of the Bullet cluster in the Σ(SFR)/M – z plot, intermediate between that of clusters and groups (Chung et al. 2010). If rapid ram-pressure stripping was a major star formation quenching process, the group LIRGs could have had their star formation quenched already, given that the cluster-group collision happened  ~250 Myr ago (Springel & Farrar 2007).

Starvation, caused by the removal of the hot gas halo reservoirs of galaxies, is a more likely candidate given that we see an accelerated evolution (relative to the field) both in clusters and in groups (Kawata & Mulchaey 2008; Bekki 2009). The removal of galaxy hot gas reservoirs inducing starvation can be caused by tidal galaxy-galaxy encounters or by the interaction with the intra-cluster/intra-group medium. Starvation should proceed more effectively in higher (galaxy or gas) density regions, hence it should quench SF earlier in cluster than in group galaxies, as we observe.

In order to better constrain these evolutionary scenarios it is important to extend the analysis of the Σ(SFR)/M  – z relation to even higher redshifts, where we might expect to see a reversal of the environment dependence as we approach the epoch when galaxies in dense environments underwent their first major episodes of SF. In the near future we plan to use new Herschel PACS observations of 6 poor and rich clusters at 1.4 < z < 1.9 to explore this issue.


1

The cluster virial mass, M200, is the mass contained within a sphere of radius r200, which is the radius within which the enclosed average mass density is 200 times the critical density. The two quantities are related by the expression , where G is the gravitational constant and H(z) the Hubble constant at the system redshift.

2

The galaxy rest-frame velocities with respect to the system mean velocity are obtained by the usual relation (Harrison & Noonan 1979), where is the system mean redshift.

Acknowledgments

PACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA (Germany); IFSI, OAP/AOT, OAA/CAISMI, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICYT/MCYT (Spain). We gratefully acknowledge the contributions of the entire COSMOS collaboration consisting of more than 100 scientists. More information about the COSMOS survey is available at http://www.astro.caltech.edu/~cosmos. This research has made use of NASA’s Astrophysics Data System, of NED, which is operated by JPL/Caltech, under contract with NASA, and of SDSS, which has been funded by the Sloan Foundation, NSF, the US Department of Energy, NASA, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council of England. The SDSS is managed by the participating institutions (www.sdss.org/collaboration/credits.html).

References

All Tables

Table 1

Main properties of the cluster and group samples.

All Figures

thumbnail Fig. 1

Comparison of system mass estimates obtained from X-ray data with those obtained from system velocity dispersions.

Open with DEXTER
In the text
thumbnail Fig. 2

Σ(SFR)/M- redshift relation for clusters (black symbols) and groups (magenta symbols). Square symbols, triangles, and dots identify, respectively, the COSMOS composite systems, the GOODS systems, and the remaining systems. The red star identifies the Bullet cluster. The black solid line shows the best fit Σ(SFR)/M – z relation for the cluster sample, excluding the Bullet cluster. The relation fitted to the data is of the type Σ(SFR)/M  ∝ zα. The dashed line shows the relation, Σ(SFR)/M  ∝ (1 + z)α, of Bai et al. (2009) rescaled to the LIRG LIR-regime and to the r200 region. The magenta solid line shows the best fit for the sample of groups and poor clusters.

Open with DEXTER
In the text
thumbnail Fig. 3

The field Σ(SFR)/M-redshift relation from Magnelli et al. (2011) (light blue shaded region) and Gruppioni et al. (2011) (dashed blue line). The shading and error bars represent 1σ confidence levels Other symbols/lines have the same meaning of the symbols as in Fig. 2.

Open with DEXTER
In the text